WRFOUT输出变量 缺少QGRAUP 已改registry

验证码yyy

WRF小白 求问 WRFOUT输出变量 有QVAPOR QCLOUD QRAIN QICE QSNOW 但是缺少QGRAUP
tate   real    -              ikjftb   moist       1         -     -    -
state   real    qv             ikjftb   moist       1         -     \
   i0rhusdf=(bdy_interp:dt)  "QVAPOR"           "Water vapor mixing ratio"      "kg kg-1"
state   real    qc             ikjftb   moist       1         -     \
   i0rhusdf=(bdy_interp:dt)  "QCLOUD"           "Cloud water mixing ratio"      "kg kg-1"
state   real    qr             ikjftb   moist       1         -     \
   i0rhusdf=(bdy_interp:dt)  "QRAIN"            "Rain water mixing ratio"       "kg kg-1"
state   real    qi             ikjftb   moist       1         -     \
   i0rhusdf=(bdy_interp:dt)  "QICE"             "Ice mixing ratio"              "kg kg-1"
state   real    qs             ikjftb   moist       1         -     \
   i0rhusdf=(bdy_interp:dt)  "QSNOW"            "Snow mixing ratio"             "kg kg-1"
state   real    qg             ikjftb   moist       1         -     rhdu
   i0rhusdf=(bdy_interp:dt)  "QGRAUP"           "Graupel mixing ratio"          "kg kg-1"
state   real    qh             ikjftb   moist       1         -     rhdu
   i0rhusdf=(bdy_interp:dt)  "QHAIL"            "Hail mixing ratio"             "kg kg-1"
state   real    -              ikjftb   dfi_moist       1         -     -    -




Registry_em_common相关代码已经改为rhdu 并且重新编译然后运行 仍然没有 QGRAUP这个变量输出变量
请问是什么问题！
以下是我的namelist.input
&time_control
run_days                            = 0,
run_hours                           = 48,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2008, 2008, 2008,
start_month                         = 07,   07,   07,
start_day                           = 16,   16,   16,
start_hour                          = 12,   12,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2008, 2008, 2008,
end_month                           = 07,   07,   07,
end_day                             = 18,   18,   18,
end_hour                            = 12,   12,   12,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 180,  60,   60,
frames_per_outfile                  = 1000, 1000, 1000,
restart                             = .false.,
restart_interval                    = 5000,
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                             = 3,
e_we                                = 79,    124,   130,
e_sn                                = 61,    100,   100,
e_vert                              = 30,    30,    30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 90000, 30000,  10000,

dy                                  = 90000, 30000,  10000,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     26,    66,
j_parent_start                      = 1,     10,    36,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 11,     11,     11,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 1,     1,     1,
radt                                = 30,    30,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                      = 1,     1,     1,
bldt                                = 0,     0,     0,
cu_physics                          = 1,     1,     0,
cudt                                = 5,     5,     5,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,     0,     0,
/

&fdda
/

&dynamics
w_damping                           = 0,
diff_opt                            = 1,      1,      1,
km_opt                              = 4,      4,      4,
diff_6th_opt                        = 0,      0,      0,
diff_6th_factor                     = 0.12,   0.12,   0.12,
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,  5000.,  5000.,

dampcoef                            = 0.2,    0.2,    0.2
khdif                               = 0,      0,      0,
kvdif                               = 0,      0,      0,
non_hydrostatic                     = .true., .true., .true.,
moist_adv_opt                       = 1,      1,      1,
scalar_adv_opt                      = 1,      1,      1,
/

&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false.,
nested                              = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

晓晓

看看用户手册，这个my_file_d01.txt的有关说明，说不定不用修改registry就能输出所有想要输出的变量~
Good luck~

验证码yyy

晓晓 发表于 2020-8-24 15:07
看看用户手册，这个my_file_d01.txt的有关说明，说不定不用修改registry就能输出所有想要输出的变量~
Good ...

你好 现在可以成功输出相关变量了 wrfout 有这个输出了 可是输出的数据都是00000请问 前辈见过这个情况么

晓晓

3层嵌套需要在namelist里有3个my_file_d0X.txt说明，这样每个domian才会有结果。。。

验证码yyy

晓晓 发表于 2020-9-15 16:34
3层嵌套需要在namelist里有3个my_file_d0X.txt说明，这样每个domian才会有结果。。。

你好 请问是要看WRFV3.1中文手册不分么
我从那篇文章上看到
在垂直插值气象要素场的过程中，real程序要求在每个相同的层次上和其它气象要素一样都有一个3-D的气压场。calc_ecmwf_p.exe应用程 序可以被用于利用ECMWF的sigma坐标数据集来生成这样一个气压场。考虑了地表气压场（地面气压场的记录）及一些列系数A和 B，calc_ecmwf_p.exe可以在ECMWF的sigma坐标中的k层的（i，j）点计算出气压，即Pijk = Ak Bk*Psfcij。这些系数可sigma高度-系数的对应表中找到。这个文件被写在普通的txt文件，ecmwf_coeffs中。例如，16个sigma层的emcwf_coeffs文件会包含

这一部分 有个txt.... 其它地方我没看到又说namelist my_file_d0X的相关说明啊 是我照的章节不对还是文章不对呢

xinshou6

他给你说了 你看新版除非你用的旧版那可能还没有 为啥出不来 我也不懂 是不是物理机制不好 前面有个chem的改个东西就行了  还有为啥是rh没有i呢 是不需要吗----搜下好像grib也有他，改i01试试？Registry.EM_COMMON里都有应该能直接抄吧 phy里没看到他应该不涉及什么机制把

诸葛皮蛋

晓晓 发表于 2020-8-24 15:07
看看用户手册，这个my_file_d01.txt的有关说明，说不定不用修改registry就能输出所有想要输出的变量~
Good ...

你好，我用WRF4.0，QVAPOR QCLOUD QRAIN QICE QSNOW QGRAUP这几个变量都没输出，我按照说明手册添加了my_file_d01.txt，但仍然没输出，还报错了，请问这个方法难到也要重新编译才可以吗？

晓晓

诸葛皮蛋 发表于 2023-4-28 10:31
你好，我用WRF4.0，QVAPOR QCLOUD QRAIN QICE QSNOW QGRAUP这几个变量都没输出，我按照说明手册添加了my_ ...

不可能吧！！！这几个变量都是wrfout里默认输出的变量，
float QVAPOR(Time, bottom_top, south_north, west_east) ;
                QVAPOR:FieldType = 104 ;
                QVAPOR:MemoryOrder = "XYZ" ;
                QVAPOR:description = "Water vapor mixing ratio" ;
                QVAPOR:units = "kg kg-1" ;
                QVAPOR:stagger = "" ;
                QVAPOR:coordinates = "XLONG XLAT XTIME" ;
        float QCLOUD(Time, bottom_top, south_north, west_east) ;
                QCLOUD:FieldType = 104 ;
                QCLOUD:MemoryOrder = "XYZ" ;
                QCLOUD:description = "Cloud water mixing ratio" ;
                QCLOUD:units = "kg kg-1" ;
                QCLOUD:stagger = "" ;
                QCLOUD:coordinates = "XLONG XLAT XTIME" ;
        float QRAIN(Time, bottom_top, south_north, west_east) ;
                QRAIN:FieldType = 104 ;
                QRAIN:MemoryOrder = "XYZ" ;
                QRAIN:description = "Rain water mixing ratio" ;
                QRAIN:units = "kg kg-1" ;
                QRAIN:stagger = "" ;
                QRAIN:coordinates = "XLONG XLAT XTIME" ;
        float QICE(Time, bottom_top, south_north, west_east) ;
                QICE:FieldType = 104 ;
                QICE:MemoryOrder = "XYZ" ;
                QICE:description = "Ice mixing ratio" ;
                QICE:units = "kg kg-1" ;
                QICE:stagger = "" ;
                QICE:coordinates = "XLONG XLAT XTIME" ;
        float QSNOW(Time, bottom_top, south_north, west_east) ;
                QSNOW:FieldType = 104 ;
                QSNOW:MemoryOrder = "XYZ" ;
                QSNOW:description = "Snow mixing ratio" ;
                QSNOW:units = "kg kg-1" ;
                QSNOW:stagger = "" ;
                QSNOW:coordinates = "XLONG XLAT XTIME" ;
        float QGRAUP(Time, bottom_top, south_north, west_east) ;
                QGRAUP:FieldType = 104 ;
                QGRAUP:MemoryOrder = "XYZ" ;
                QGRAUP:description = "Graupel mixing ratio" ;
                QGRAUP:units = "kg kg-1" ;
                QGRAUP:stagger = "" ;
                QGRAUP:coordinates = "XLONG XLAT XTIME" ;
   。。。。

// global attributes:
                :TITLE = " OUTPUT FROM WRF V4.0 MODEL" ;

按照手册来操作肯定可以输出想要的数据，自己再确认一下吧

诸葛皮蛋

QAQ,我再去找找原因，我说呢，我看注册文件里也是默认输出的量呀，可能是我没读取出来，我再试试。谢谢你的答复！

nightfalls

请问楼主解决了吗？我也遇到了同样的问题，注册表里SWUPTCLN的输出已经是rhdu了，但是wrfout就没有。