[RegCNET] something about SIGMAtoP of the latest r egcm3

Tue Aug 15 12:01:57 CEST 2006

Hi,

Shekhar wrote an email to RegCNET about this problem back in July. There
is an additional variable in the atmospheric output for the most recent
RegCM release - omega. This is accounted for in the SIGMAtoP.f program,
but not in the sample PLEV.CTL file in Commons/tools.

Therefore, you need to add a line for omega (w) after v in PLEV_VAR.ctl.
Also, don't forget to change the vars to 13 instead of 12. So, your
PLEV_VAR.ctl should look something like this (with your own PDEF, XDEF,
and YDEF):

set ^PLEV_VAR
title model-output on P-plane
options big_endian
undef -1.e34
pdef 49 32 lcc 45.39 13.48 24.50 16.00 30.00 60.00 13.48 60000. 60000.
xdef 175 linear -9.60 0.2703
ydef 78 linear 34.70 0.2703
zdef 11 levels 1000 925 850 700 500 400 300 250 200 150 100
tdef 60 linear 06z01Jul1994 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 omega
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

Let me know if this doesn't solve your problem!

cheers,
sara

jwzhou66 wrote:
> Hello,everyone!
> I have a question about the SIGMAtoP of the latest regcm3(2006).
> After the simulating by the regcm3,I do SIGMAtoP in output
> directory,and copy pdef,xdef,ydef,tdef from ATM.********.ctl file,but
> while displaying the PLEV_VAR results in
> GrADS,the result is not correct and illogical.I also display the
> viriable in ATM.********.ctl in GrADS,it seems reasonable,the
> following is what I find:
> ga-> open PLEV_VAR.ctl
> Scanning description file: PLEV_VAR.ctl
> Data file PLEV_VAR is open as file 1
> LON set to 33.9 172.834
> LAT set to -3.68 61.7326
> LEV set to 1000 1000
> Time values set: 1998:5:25:0 1998:5:25:0
> Notice: Implied interpolation for file PLEV_VAR.ctl
> Interpolation will be performed on any data displayed from this file
> ga-> set t 1 19
> Time values set: 1998:5:25:0 1998:5:29:12
> ga-> d rain
> Notice: Automatic Grid Interpolation Taking Place
> Constant field. Value = 0
> Contouring: 0.002 to 0.02 interval 0.002
> Contouring: -0.006 to 0 interval 0.002
> Too many contour levels -- adjusting cint
> Contouring: -24 to 32 interval 2
> Contouring: -26 to 28 interval 2
> Contouring: 266 to 306 interval 2
> Too many contour levels -- adjusting cint
> Contouring: 1260 to 1560 interval 20
> Contouring: 0 to 0 interval 20
> Contouring: 20 to 0 interval 20
> Contouring: -0 to 0 interval 20
> Contouring: -20 to 40 interval 20
> Contouring: -0 to 40 interval 20
> Contouring: 220 to 220 interval 20
> Contouring: 11560 to 12560 interval 20
> Contouring: 280 to 300 interval 20
> Contouring: 0 to 0 interval 20
> Contouring: 20 to 0 interval 20
> Contouring: -0 to 0 interval 20
> Contouring: -20 to 20 interval 20
> Press enter to continue:::::
> ga-> open ATM.1998052500.ctl
> Scanning description file: ATM.1998052500.ctl
> Data file ATM.1998052500 is open as file 1
> LON set to 33.9 172.834
> LAT set to -3.68 61.7326
> LEV set to 1008.43 1008.43
> Time values set: 1998:5:25:0 1998:5:25:0
> Notice: Implied interpolation for file ATM.1998052500.ctl
> Interpolation will be performed on any data displayed from this file
> ga-> set t 1 19
> Time values set: 1998:5:25:0 1998:5:29:12
> ga-> d tpr
> Notice: Automatic Grid Interpolation Taking Place
> Constant field. Value = 0
> Contouring: -0 to 350 interval 50
> Contouring: -0 to 1950 interval 50
> Contouring: -0 to 1100 interval 50
> Contouring: -0 to 1100 interval 50
> Contouring: -0 to 900 interval 50
> Contouring: -0 to 950 interval 50
> Contouring: -0 to 1300 interval 50
> Contouring: -0 to 850 interval 50
> Contouring: -0 to 800 interval 50
> Contouring: -0 to 1200 interval 50
> Contouring: -0 to 1550 interval 50
> Contouring: -0 to 1350 interval 50
> Contouring: -0 to 1050 interval 50
> Contouring: -0 to 1400 interval 50
> Contouring: -0 to 1550 interval 50
> Contouring: -0 to 1900 interval 50
> Contouring: -0 to 900 interval 50
> Contouring: -0 to 850 interval 50
> Press enter to continue::::: ga->
> and I also do the SIGMAtoP for ICBC in Input directory,its results is
> right.I ever used the regcm3 of 2005and find no problem in SIGMAtoP of
> output,so I wondered!
> Please tell me !
> At last,can anyone give me the source of regcm3 in 2005,thank you !
> Hope for you reply!
>
>
>
>
>
>
>
> 张 靓 影 TCL 敢 梦 敢 想 广 告剧 本 征 集
> 为 张 靓 影 编 写 广告 剧 本 ，赢 取 24 台 笔 记 本 和 3 万 元 现 金
> ！
> <http://adclient.163.com/event.ng/Type=click&FlightID=65619&AdID=67427&TargetID=2095&Values=31,43,51,60,72,78,85,92,100,110,312,332,499,589,700,702,733,734&Redirect=http://tcl.163.com>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
>   

-- 

Dr. Sara A. Rauscher
Earth System Physics Section 
The Abdus Salam International Centre for Theoretical Physics
Strada Costiera 11 
34014 Trieste ITALY
ph: +39 040 2240 225 fax: +39 040 2240 449
email: srausche at ictp.it