arwpost生成的ctl中z层数与interp_levels的对应关系

一颗梅子在南京

grads画图的时候遇到一个问题，ARWpost生成的高度层为29层，但是namelist.arwpost里面的inter_levels好像没有设置这么多层，不LL知道这两者的对应关系是什么样的：
CTL文件如下：
title  OUTPUT FROM WRF V3.5 MODEL
pdef  111  96 lcc  30.464  118.766   56.000   48.500  60.00000  30.00000  115.00000  10000.000  10000.000
xdef  291 linear  112.43001   0.04504504
ydef  219 linear   25.40103   0.04504504
zdef   29 linear 1 1  
tdef   15 linear 00Z28JUN2016     180MN      

namelist.arwpost文件如下：红色部分应该就是inter_levels的设置吧，数了一下应该是19层
&datetime
start_date = '2013-10-05_12:00:00',
end_date   = '2013-10-07_06:00:00',
interval_seconds = 3600,
tacc = 0,
debug_level = 0,
/
&io
input_root_name = '../WRFV3/test/em_real/wrfout_d02_2013-10-*'
output_root_name = './wrfout.d02'
plot = 'all_list'
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'
mercator_defs = .true.
/
split_output = .true.
frames_per_outfile = 2

plot = 'all'
plot = 'list'
plot = 'all_list'
! Below is a list of all available diagnostics
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'

&interp
interp_method = 0,
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
/
extrapolate = .true.
interp_method = 0,     ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,

求大神指导一下，谢谢！

孤独卡齐莫多

看说明,意思是说要设置了interp_method=1才会用你自己设置的层数.

男紫汗

楼上正解。另外你重新生成后能不能画一下高度场，看看等高线的情况

一颗梅子在南京

孤独卡齐莫多 发表于 2016-7-28 15:04
看说明,意思是说要设置了interp_method=1才会用你自己设置的层数.

我有试过改成1后生成的ctl还是29层

一颗梅子在南京

孤独卡齐莫多 发表于 2016-7-28 15:04
看说明,意思是说要设置了interp_method=1才会用你自己设置的层数.

extrapolate = .true.
interp_method = 0,     ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,
这段应该是只是起到一个解释的作用吧！
我再问一个问题：当interp_method设为0的时候，对应的interp_levels是不是对应的设置为：,interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0,11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,
而当interp_method设为1时，对应的interp_levels是不是对应的设置为：interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,

孤独卡齐莫多

在那遥远的地方 发表于 2016-7-28 15:39
感觉这两个应该只是高度的表示方式不一样吧，实际用的时候要注释掉一行的，有一种

我记得我以前是设为1的,你可以试试,我也记不太清楚

一颗梅子在南京

在那遥远的地方 发表于 2016-7-28 17:21
是的，设置成1是自定义

恩，我搞定了，设置成1以后就可以了