This is an old (1 month) thread ...

I find the latest postproc.f useful for compressing large model
output on CDF files, so I have started debugging it
(initialization of some xxx3d_p arrays, and declaration of some
variable types). I'm now using the plv option correctly on ATM fields.

Unfortunately, I have also checked that using the daily average of postproc.f 
is only very-similar-to, but not identical-to, using the combination of

  DAY_month.f
  SIGMAtoP

suggested by Bi in another recent post on this subject.

 I have not touched the RAD/CHEM processing and the usgs option yet
 (no check on usgs so it is considered always .true. even if is .false. in user.in).
It is a starting point but it works.

I'm not posting the code since I'm still fixing it. However
I'll be glad to know if anybody else is working on this.

Sandro



The version on the web doesn't have plv option tested. Both usgs and plv 
options are set to false for now in  postproc.in (and left for later use) - 
This means at the moment, sigma to pressure should not be used. I left this 
option incomplete because of some issues in RegCM output e.g I wanted PS 
(surface pressure) variable in both RAD and CHEM files, which is needed for 
pressure level interpolation - since PS is not there in RAD/CHEM files, I found 
it least useful either to make postprocessing of RAD and CHEM dependent on the 
presenece of ATM files (for reading PS form there) or to let this option 
working for ATM only. Moreover, there are very easy few lines scripts in 
ferret/ncl which could do  this sigma to pressure interpolation. 

However, if you think that its still useful to have at least ATM files working 
with  plv option, I can test the code for that.


Moet 


Quoting kdkorte at ncsu.edu:

> All,
> 
> I am having problems when changing from sigma levels to pressure levels
> using the plv option within postproc.in
> 
> I try to view various variables at several different levels and times and
> get no value errors.  Previously I used postproc with plv set to false and
> had no problems.
> 
> Any help would be greatly appreciated.
> 
> Thanks,
> 
> Kurt Korte
> kdkorte at ncsu.edu
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
> 


-- 
====================================================
Sandro Calmanti                             ENEA-ACS
Tel:  +39 0630484072    Via Anguillarese, 301, 00060 
Fax:  +39 0630484264       S. Maria di Galeria (ROMA)
http://clima.casaccia.enea.it/staff/calmanti