小白转wrf出了.dat和.ctl 用grads画图降水量全为负值

蚊子

按照网上找的资料转wrf。出了.dat和.ctl 用grads画图降水量全为负值
namelist.wps如下

&share

wrf_core = 'ARW',

max_dom = 2,

start_date = '2011-05-01_00:00:00','2011-05-01_00:00:00',

end_date   = '2011-05-01_12:00:00','2011-05-01_12:00:00',

interval_seconds = 21600

io_form_geogrid = 2,

/

&geogrid

parent_id         =   1,   1,

parent_grid_ratio =   1,   3,

i_parent_start    =   1,  30,

j_parent_start    =   1,  19,

s_we              =   1,   1,

e_we              =  74, 112,

e_sn              =  59, 100,

geog_data_res     = '10m','2m',

dx = 30000,

dy = 30000,

map_proj = 'lambert',

ref_lat   =  31.5,

ref_lon   =  103,

truelat1  =  30.0,

truelat2  =  60.0,

stand_lon =  103,

geog_data_path = '/home/geog/geog_complete/geog'

/

&ungrib

out_format = 'WPS',

prefix = 'FILE',

/

&metgrid

fg_name = 'FILE'

io_form_metgrid = 2,

/

&mod_levs

press_pa= 201300 , 200100 , 100000 , 95000 , 90000 , 85000 , 80000 , 75000 , 70000 , 65000 , 60000 , 55000 , 50000 , 45000 , 40000 , 35000 , 30000 , 25000 , 20000 , 15000 , 10000 , 5000 , 1000

/

namelist.input如下

&time_control

run_days                            = 0,

run_hours                           = 12,

run_minutes                         = 0,

run_seconds                         = 0,

start_year                          = 2011, 2011,

start_month                         = 05,   05,   

start_day                           = 01,   01,   

start_hour                          = 00,   00,   

start_minute                        = 00,   00,   

start_second                        = 00,   00,   

end_year                            = 2011, 2011,

end_month                           = 05,   05,   

end_day                             = 01,   01,   

end_hour                            = 12,   12,   

end_minute                          = 00,   00,   

end_second                          = 00,   00,   

interval_seconds                    = 21600

input_from_file                     = .true.,.true.,

history_interval                    = 180,  60,   

frames_per_outfile                  = 1000, 1000,

restart                             = .false.,

restart_interval                    = 500,

io_form_history                     = 2

io_form_restart                     = 2

io_form_input                       = 2

io_form_boundary                    = 2

debug_level                         = 0

/

&domains

time_step                           = 180,

time_step_fract_num                 = 0,

time_step_fract_den                 = 1,

max_dom                             = 2,

e_we                                = 74,    112,   

e_sn                                = 59,    100,   

e_vert                              = 30,    30,   

p_top_requested                     = 5000,

num_metgrid_levels                  = 27,

num_metgrid_soil_levels             = 4,

dx                                  = 30000, 10000,  

dy                                  = 30000, 10000,  

grid_id                             = 1,     2,     

parent_id                           = 1,     1,     

i_parent_start                      = 1,     30,   

j_parent_start                      = 1,     19,   

parent_grid_ratio                   = 1,     3,     

parent_time_step_ratio              = 1,     3,     

feedback                            = 1,

smooth_option                       = 0

/

&physics

mp_physics                          = 2,     2,   

ra_lw_physics                       = 1,     1,     

ra_sw_physics                       = 1,     1,     

radt                                = 10,    10,   

sf_sfclay_physics                   = 1,     1,     

sf_surface_physics                  = 2,     2,   

bl_pbl_physics                      = 1,     1,     

bldt                                = 0,     0,     

cu_physics                          = 2,     2,     

cudt                                = 0,     0,   

isfflx                              = 1,

ifsnow                              = 0,

icloud                              = 1,

surface_input_source                = 1,

num_soil_layers                     = 4,

sf_urban_physics                    = 0,     0,     

maxiens                             = 1

maxens                              = 3

maxens2                             = 3

maxens3                             = 16

ensdim                              = 144

/

&fdda

/

&dynamics

w_damping                           = 0,

diff_opt                            = 1,      1,      

km_opt                              = 4,      4,      

diff_6th_opt                        = 0,      0,      

diff_6th_factor                     = 0.12,   0.12,   

base_temp                           = 290.

damp_opt                            = 0,

zdamp                               = 5000.,  5000.,  

dampcoef                            = 0.2,    0.2,   

khdif                               = 0,      0,      

kvdif                               = 0,      0,      

non_hydrostatic                     = .true., .true.,

moist_adv_opt                       = 1,      1,           

scalar_adv_opt                      = 1,      1,         

/

&bdy_control

spec_bdy_width                      = 5,

spec_zone                           = 1,

relax_zone                          = 4,

specified                           = .true., .false.,

nested                              = .false., .true.,

/

&grib2

/

&namelist_quilt

nio_tasks_per_group = 0,

nio_groups = 1,

/

namelist.ARWpost如下

&datetime

start_date = '2011-05-01_00:00:00',

end_date   = '2011-05-01_12:00:00',

interval_seconds = 3600,

tacc = 0,

debug_level = 0,

/

&io

input_root_name = '/home/mode/WRFV3/run/wrfout_d01_2011-05-01_00:00:00'

output_root_name = '/home/mode/data/test'

plot = 'all_list'

fields = 'height,pressure,tk,tc'

output_type = 'grads'

mercator_defs = .true.

/

split_output = .true.

frames_per_outfile = 1

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 = 1,

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,

麻烦前辈们看看是不是wrf编译时出现问题，刚学不太懂，降水量为什么会为负？

男紫汗

这是缺测值吧。。

蚊子

嗯嗯我到会去看，发现wrfout  do1和do2 0时到8时21分好像没出全。被该是12小时的。met是12个小时。应该是namelist.input错了。预处理时用的二重嵌套，望指点。