[RegCNET] a question about SIGMAtoP

[yyw]

Hi,everyone:
   I used RegCM3.1 on Tibetan Plateau and got ATM.YYYYMMDDHH、RAD.YYYYMMDDHH、SRF.YYYYMMDDHH，then I convert ATM.....to PLEV_VAR by using SIGMAtoP.f.In PLEV_VAR.ctl I add a line for omega (w) after v and change pdef、xdef、ydef、tdef.But there is a problem.When I open PLEV_VAR.ctl by grads and display air temperature(t),I find the value of January is over 267 and reach 297.That is unreasonable. Does anyone encounter the problem?How to explain that?
   In SIGMAtoP.f I change the following lines :
      parameter(iy=54,jx=92,kx=23)
      parameter (nfile=12)
      data inout/'ATM.2000010100','ATM.2000020100','ATM.2000030100'
     &           ,'ATM.2000040100','ATM.2000050100','ATM.2000060100'
     &           ,'ATM.2000070100','ATM.2000080100','ATM.2000090100'
     &           ,'ATM.2000100100','ATM.2000110100','ATM.2000120100'/
      data      number/125,116,124,120,124,120,124,124,120,124,120,123/
-------------------------------------------------------------------------
PLEV_VAR.ctl is as follows:

dset ^PLEV_VAR
title model-output on P-plane
options big_endian
undef -1.e34
pdef   88   48 lcc   32.00   90.00   44.00   24.00   30.00   60.00   90.00  45000.  45000.
xdef  224 linear   65.65  0.2027
ydef  112 linear   20.20  0.2027
zdef 11 levels 1000 925 850 700 500 400 300 250 200 150 100
tdef  125 linear  6z01jan2000  6hr
vars 13
h 11 0 p-level height
t 11 0 air temperature
u 11 0 zonal wind
v 11 0 meridional wind
w 11 0 p-velocity
q 11 0 specific moisture
qc 11 0 cloud water
ps 0 99 surface pressure
slp 0 99 sea level pressure
rain 0 99 precipitation
tg 0 99 temperature of lower soil layer
swt 0 99 total soil water in mm H2O
rno 0 99 accumulated infiltration
endvars
~                              
 


　　　　　　　　yyw
　　　　　　　　ywyang at lzb.ac.cn
　　　　　　　　　　2006-10-10