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！）