飑线理想试验

irish

请问有人做过wrf自带飑线理想试验吗？默认是二维的，我把y方向格点设多一点 变成三维的，就会没有对流生成，降水都是零，请问有人知道如何修改吗？是源代码里面热泡的问题吗？还有请问热泡的初始位温是在哪个位置设的呀？

ljchen1989

理想实验初始条件在程序里面hard-wired的
看看生成初始条件的子程序dyn_em/module_initialize_squall2d_x.F哈
热泡也在里面

irish

ljchen1989 发表于 2018-9-13 09:25
理想实验初始条件在程序里面hard-wired的
看看生成初始条件的子程序dyn_em/module_initialize_squall2d_x. ...

DO J = jts, min(jde-1,jte)
    yrad = config_flags%dy*float(j-nyc)/4000.
!     yrad = 0.
    DO I = its, min(ide-1,ite)
      xrad = config_flags%dx*float(i-nxc)/4000.
!     xrad = 0.
      DO K = 1, kte-1

!  put in preturbation theta (bubble) and recalc density.  note,
!  the mass in the column is not changing, so when theta changes,
!  we recompute density and geopotential

        zrad = 0.5*(grid%ph_1(i,k,j)+grid%ph_1(i,k+1,j)  &
                   +grid%phb(i,k,j)+grid%phb(i,k+1,j))/g
        zrad = (zrad-1500.)/1500.
        RAD=SQRT(xrad*xrad+yrad*yrad+zrad*zrad)
        IF(RAD <= 1.) THEN
           grid%t_1(i,k,j)=grid%t_1(i,k,j)+delt*COS(.5*PI*RAD)**2
           grid%t_2(i,k,j)=grid%t_1(i,k,j)
           qvf = 1. + rvovrd*moist(i,k,j,P_QV)
           grid%alt(i,k,j) = (r_d/p1000mb)*(grid%t_1(i,k,j)+t0)*qvf* &
                        (((grid%p(i,k,j)+grid%pb(i,k,j))/p1000mb)**cvpm)
           grid%al(i,k,j) = grid%alt(i,k,j) - grid%alb(i,k,j)
        ENDIF
      ENDDO
这个是程序里关于热泡那一块，没有找到关于初始位温扰动的值是在哪块设定的呢？？还有热启动是加热泡在里面，请问知道冷启动是怎么回事儿？谢谢~

ljchen1989

grid%t_1(i,k,j)=grid%t_1(i,k,j)+delt*COS(.5*PI*RAD)**2
这一行就是啊，开头声明变量的时候有一个参数，好像是delt，生命热泡振幅的
热启动是循环同化里面的，和这个没关系，理想实验都是冷启动

irish

ljchen1989 发表于 2018-9-14 19:26
grid%t_1(i,k,j)=grid%t_1(i,k,j)+delt*COS(.5*PI*RAD)**2
这一行就是啊，开头声明变量的时候有一个参数， ...

嗯嗯 谢谢~还有一个问题请教一下，请问理想试验有类似真实个例那种rsl.out吗？因为我在namelist里面设的模拟时长6小时，10分钟输出一次，可是出来的wrfout里面数据只有12个时次，也不知道是哪里出了问题。。。。求解答~~万分感谢！！

ljchen1989

要不你把namelist或者rsl贴出来看看？有可能崩了
而且你要看一眼你生成的ic是不是满足你需求，你改程序不一定改的对，需要验证

irish

ljchen1989 发表于 2018-9-14 22:51
要不你把namelist或者rsl贴出来看看？有可能崩了
而且你要看一眼你生成的ic是不是满足你需求，你改程序不 ...

&time_control
run_days                            = 0,
run_hours                           = 6,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 0001,
start_month                         = 01,
start_day                           = 01,
start_hour                          = 00,
start_minute                        = 00,
start_second                        = 00,
end_year                            = 0001,
end_month                           = 01,
end_day                             = 01,
end_hour                            = 06,
end_minute                          = 00,
end_second                          = 00,
history_interval                    = 30,
frames_per_outfile                  = 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                           = 9,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
s_we                                = 1,
e_we                                = 402,
s_sn                                = 1,
e_sn                                = 120,
s_vert                              = 1,
e_vert                              = 81,
dx                                  = 500,
dy                                  = 500,
ztop                                = 20000.,
/

&physics
mp_physics                          = 2,
ra_lw_physics                       = 0,
ra_sw_physics                       = 0,
radt                                = 0,
sf_sfclay_physics                   = 0,
sf_surface_physics                  = 0,
bl_pbl_physics                      = 0,
bldt                                = 0,
cu_physics                          = 0,
cudt                                = 0,
num_soil_layers                     = 5,
/

&fdda
/

&dynamics
rk_ord                              = 3,
diff_opt                            = 2,      2,      2,
km_opt                              = 2,      2,      2,
damp_opt                            = 2,
dampcoef                            = .003,
zdamp                               = 5000.,
khdif                               = 300,
kvdif                               = 1,
smdiv                               = 0.1,
emdiv                               = 0.01,
epssm                               = 0.1,
time_step_sound                     = 6,
h_mom_adv_order                     = 5,
v_mom_adv_order                     = 3,
h_sca_adv_order                     = 5,
v_sca_adv_order                     = 3,
mix_full_fields                     = .true.,
non_hydrostatic                     = .true.,
moist_adv_opt                       = 1,      1,      1,     
scalar_adv_opt                      = 1,      1,      1,     
chem_adv_opt                        = 1,      1,      1,     
tke_adv_opt                         = 1,      1,      1,     
/

&bdy_control
periodic_x                          = .false.,
symmetric_xs                        = .false.,
symmetric_xe                        = .false.,
open_xs                             = .true.,
open_xe                             = .true.,
periodic_y                          = .true.,
symmetric_ys                        = .false.,
symmetric_ye                        = .false.,
open_ys                             = .false.,
open_ye                             = .false.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
我只是把他原来模拟时长1小时改成了6小时，把interval从60改成5000，然后格点数增多了一点，出来的时次就是不对，有时候是12，有时候是1，很奇怪。interval是60的时候会有wrfrst文件生成，改成5000后没有了，但也没有对流了。。。很奇怪。。求解答~~万分感谢

ljchen1989

history_interval是30，半小时输出一次，算上IC总共输出13次？
cu_physics=0，把积云对流关了，只保留mp，只有浅对流，rainnc？

ljchen1989

而且你模式稳定么？
dx=0.5km，dt=9s
dt/dx=18
官方建议小于6么，不是？

irish

ljchen1989 发表于 2018-9-16 12:33
而且你模式稳定么？
dx=0.5km，dt=9s
dt/dx=18

谢谢您了，问了别人，说是我用的模块不对，我编译的是em_squall2d_x 说是这个模拟不了三维的，我准备换em_quarter_ss试一下。就是探空要换一下，换成飑线的