Hi regcm3-users:
After converting sigma-level data to pressure levels, u and v can be displayed in grads well,but u;v can't be showed in grads. I don't know why, and sure the commands of grads I entered are right,Maybe the problem lays in the process of converting sigma-level data to pressure levels.Besides,I converted ICBCyyyymmdd to IC_Pyyyymmdd,IC_Pyyyymmdd displayed not very well.So, Any any information about converting sigma-level data to pressure levels, I'm very glad to recevive!
Any help is appreciated very very much, the settings of SIGMAtoP.f and PLEV_VAR.ctl are below:
My SIGMAtoP.f is like below:
implicit none
integer iy,jx,kx,np
parameter(iy=54,jx=92,kx=23)
parameter(np=11)
real plev(np)
data plev/1000.,925.,850.,700.,500.,400.,300.,250.,200.
& ,150.,100./
integer i,j,k,mrec,lrec
c
integer nfile
parameter (nfile=2)
character*14 inout(nfile)
data inout/'ATM.1997070100','ATM.1997100100'/ /* here two ATM files, but later just integer number(nfile) /* one data allout /'PLEV_710'/? */
data number/20,40/
character*8 allout
data allout /'PLEV_710'/
..................................
.................................
My PLEV_VAR.ctl is like below:
dset PLEV_710
title model-output on P-plane
options big_endian
undef -1.e34
pdef 92 54 lcc 33.00 89.00 46.00 27.00 30.00 60.00 89.00 45000. 45000.
xdef 263 linear 62.68 0.2027
ydef 125 linear 19.81 0.2027
zdef 11 levels 1000 925 850 700 500 400 300 250 200 150 100
tdef 124 linear 6z01jul1997 6hr
....................................
..................................
Thanks a lot!
好好学习,天天向上(good good study,day day up!)