Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Bogus Storm STOP to 6 sec. (TC.exe)

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

L

L Alvarez

New member
Hello there,

I have followed the tutorial here->
http://www2.mmm.ucar.edu/wrf/users/docs/user_guide_V3/users_guide_chap10.htm#TCBogus

I created a bogus by modifying the namelist.input, then I executed the tc.exe and renamed the generated auxinput1_d01_<date> file output to:

Code: 
 auxinput1_d01_<date> to the name that the real.exe program expects, met_em.d01.<date>

Then I had to modify the namelist.input again by entering ->

force_use_old_data= .true.,
Click to expand...

Then I proceed to run the real.exe, which ends well but NOT wrf.exe, It always stops in the second 6 of the integration. I have painted my domain 1 with GrADS, below you can see it. I understand that the bogus assimilates it from domain 1, but something else must be activated to spread to the second, third and fourth domain? .. or the error is elsewhere.


Code: 
debug = 1000

d03 2015-08-12_00:00:06+02/03 leaving update_after_feedback_em
d03 2015-08-12_00:00:06+02/03 module_integrate: back from med_nest_feedback
d02 2015-08-12_00:00:06+02/03  DEBUG wrf_timetoa():  returning with str = [2015-08-12_00:00:06]
Timing for main: time 2015-08-12_00:00:06 on domain   2:   30.68426 elapsed seconds
d02 2015-08-12_00:00:06+02/03 in med_latbound_in
d02 2015-08-12_00:00:06+02/03 module_integrate: calling solve interface
d02 2015-08-12_00:00:06+02/03  grid spacing, dt, time_step_sound=   1800.000       6.666667               4
d02 2015-08-12_00:00:06+02/03 calling inc/HALO_EM_MOIST_OLD_E_7_inline.inc
d02 2015-08-12_00:00:06+02/03 calling inc/PERIOD_BDY_EM_MOIST_OLD_inline.inc
d02 2015-08-12_00:00:06+02/03  call rk_step_prep
d02 2015-08-12_00:00:06+02/03 calling inc/HALO_EM_A_inline.inc
d02 2015-08-12_00:00:06+02/03 calling inc/PERIOD_BDY_EM_A_inline.inc
d02 2015-08-12_00:00:06+02/03  call rk_phys_bc_dry_1
d02 2015-08-12_00:00:06+02/03  call init_zero_tendency
d02 2015-08-12_00:00:06+02/03 calling inc/HALO_EM_PHYS_A_inline.inc
d02 2015-08-12_00:00:06+02/03  call phy_prep
d02 2015-08-12_00:00:06+02/03  DEBUG wrf_timetoa():  returning with str = [2015-08-12_00:00:06]
d02 2015-08-12_00:00:06+02/03  call radiation_driver
d02 2015-08-12_00:00:06+02/03 Top of Radiation Driver
d02 2015-08-12_00:00:06+02/03 calling inc/HALO_PWP_inline.inc
d02 2015-08-12_00:00:06+02/03  call surface_driver
d02 2015-08-12_00:00:06+02/03 in MYJSFC
d02 2015-08-12_00:00:06+02/03 in NOAH DRV
d02 2015-08-12_00:00:06+02/03  call pbl_driver
d02 2015-08-12_00:00:06+02/03 in MYJPBL
d02 2015-08-12_00:00:06+02/03  call cumulus_driver
d02 2015-08-12_00:00:06+02/03 in cu_kfcup
d02 2015-08-12_00:00:06+02/03 returning from cumulus_driver
d02 2015-08-12_00:00:06+02/03  call shallow_cumulus_driver
d02 2015-08-12_00:00:06+02/03 calling inc/HALO_EM_FDDA_SFC_inline.inc
d02 2015-08-12_00:00:06+02/03  call fddagd_driver
d02 2015-08-12_00:00:06+02/03  call calculate_phy_tend
d02 2015-08-12_00:00:06+02/03  call compute_diff_metrics
d02 2015-08-12_00:00:06+02/03 calling inc/HALO_EM_TKE_C_inline.inc

Here my namelist.input ->

Code: 
&time_control
nwp_diagnostics = 1,
run_days = 0, 
run_hours = 24,
run_minutes = 0,
run_seconds = 0,
start_year = 2015,2015,2015,2015, 
start_month = 08,08,08,08, 
start_day = 12,12,12,12, 
start_hour = 00,00,00,00, 
end_year = 2015,2015,2015,2015, 
end_month = 08,08,08,08, 
end_day = 13,13,13,13, 
end_hour = 00,00,00,00, 
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,.true.,
history_interval = 180,180,60,5,
frames_per_outfile = 1000,1000,1000,1000,
restart = .false.,
restart_interval = 60,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 1000,
force_use_old_data= .true.,

!auxhist2_outname=’afwa_d<domain>_<date>’
!io_form_auxhist2 = 5,
!auxhist2_interval = 21600,21600,60,5,
!frames_per_auxhist2 = 1,1,1,1,
/

&tc  
insert_bogus_storm = .true.  
remove_storm = .false.  
latc_loc = 28.557  
lonc_loc = -16.484  
vmax_meters_per_second = 45  
rmax = 90000.0  
vmax_ratio = 0.9  
/

&afwa
afwa_diag_opt = 1,1,1,1,  
afwa_ptype_opt = 1,1,1,1,
afwa_vil_opt = 1,1,1,1,  
afwa_radar_opt = 1,1,1,1,
afwa_severe_opt = 1,1,1,1,
afwa_vil_opt = 1,1,1,1,
afwa_icing_opt = 1,1,1,1,  
afwa_vis_opt = 1,1,1,1,    
afwa_cloud_opt = 1,1,1,1,  
afwa_therm_opt = 1,1,1,1,  
afwa_turb_opt = 1,1,1,1,   
afwa_buoy_opt  = 1,1,1,1,   
afwa_ptype_ccn_tmp = 264.15, 
afwa_ptype_tot_melt = 50,    
/


&domains
use_adaptive_time_step = .false.,
step_to_output_time = .true.,
target_cfl = 1.2,1.2,1.2,1.2,
max_step_increase_pct = 5,51,51,51,
starting_time_step = -1,
max_time_step = -1,
min_time_step = -1,
adaptation_domain = 4,
time_step = 20,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 4,
e_we = 250,250,250,250,  
e_sn = 250,250,250,190,
!e_vert = 35,35,35,35,
e_vert = 45,45,45,45,
p_top_requested = 5000,
!sfcp_to_sfcp = .true.,
auto_levels_opt = 2,
max_dz = 1000,
dzbot = 50,
dzstretch_s = 1.5,
dzstretch_u = 1.3,
!eta_levels = 1.000, 0.998, 0.994, 0.989, 0.983,
!             0.890, 0.850, 0.823, 0.800, 0.750,
!             0.700, 0.650, 0.600, 0.550, 0.500,
!             0.450, 0.400, 0.350, 0.300, 0.250,
!             0.150, 0.100, 0.050, 0.000,
!eta_levels                           = 1.000, 0.990, 0.978, 0.964, 0.946,
!                                       0.922, 0.894, 0.860, 0.817, 0.766,
!                                       0.707, 0.644, 0.576, 0.507, 0.444,
!                                       0.380, 0.324, 0.273, 0.228, 0.188,
!                                       0.152, 0.121, 0.093, 0.069, 0.048,
!                                       0.029, 0.014, 0.000,
num_metgrid_levels = 18,
num_metgrid_soil_levels = 2,
dx = 5400,1800,600,200,  
dy = 5400,1800,600,200, 
grid_id = 1,2,3,4,
parent_id = 1,1,2,3,
i_parent_start = 1,90,75,100,
j_parent_start = 1,80,90,150,
parent_grid_ratio = 1,3,3,3,
parent_time_step_ratio = 1,3,3,3,
feedback = 1,
smooth_option = 0,
!numtiles = 2, 
!nproc_x = 2, 
!nproc_y = 8, 
/

&physics
physics_suite = 'CONUS'
mp_physics = 7,7,7,7,
cu_physics = 10,10,10,10,
!ra_lw_physics = 1,1,1,1,
!ra_sw_physics = 1,1,1,1,
bl_pbl_physics = 2,2,2,2,
!sf_sfclay_physics = 1,1,1,1,
!sf_surface_physics = 2,2,2,2,
radt = 45,45,45,45,
gsfcgce_hail = 1,
gsfcgce_2ice = 1,
topo_wind = 1,
ysu_topdown_pblmix = 1,
bldt = 0, 0, 0, 0,
kfeta_trigger = 2,
cu_rad_feedback = true,
sst_update = 0,
surface_input_source = 3,
do_radar_ref = 1,

lightning_option = 1, 1, 1, 1,
lightning_start_seconds = 600,600,600,600,
lightning_dt = 60, 20, 20, 20,
flashrate_factor =1.0,0.95,0.9,0.8,
cellcount_method =1,1,1,1,
cldtop_adjustment = 2, 2, 2, 2,
iccg_method = 4, 4, 4, 4,
iccg_prescribed_num = 3.,
iccg_prescribed_den = 1.,
/

 &fdda
 /

&dynamics
hybrid_opt = 2,
rk_ord = 3,
w_damping = 1,
diff_opt = 1,1,1,1,
km_opt = 4,4,4,4,
diff_6th_opt = 2,2,2,2,
diff_6th_factor = 0.12,0.12,0.12,0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 10000.,10000,10000,10000,
dampcoef = 0.2,0.2,0.2,0.2,
khdif = 0,0,0,0,
kvdif = 0,0,0,0,
smdiv = 0.1,
emdiv = 0.01,
epssm = 0.1,0.3,0.5,0.8,
non_hydrostatic = .true.,true.,true.,true.,
gwd_opt = 1,
h_mom_adv_order = 5,5,5,5,
h_sca_adv_order = 5,5,5,5,
v_mom_adv_order = 3,5,5,5,
v_sca_adv_order = 3,3,3,3,
moist_adv_opt = 1,1,1,1,
scalar_adv_opt = 1,1,1,1,
chem_adv_opt = 0
tke_adv_opt = 1
/

 &bdy_control
 spec_bdy_width                      = 5,
 specified                           = .true.
 /

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /
 

Attachments

  • bogus_scheme.png
    bogus_scheme.png
    49.1 KB · Views: 1,705
Hi,
Unfortunately we are aware of a problem with the TC Bogus program beginning in V3.7. The problem causes a change to the soil information, making it corrupt. There is an option to get around it - you can take all the necessary fields from the new auxinput file and overwrite them in the original met_em file so that they maintain soil information. You would need to use some other program, such as NCL to do this, though.

A couple of other things to note: I notice that you're using 4 domains. The Users Guide mentions that users with multiple domains should horizontally-interpolate the generated meteorological fields to the fine-grid domains. Did you do that? It also mentions that the tc.exe program should be built serially. Was your program built serially?
 
Hi,

At first I thought I could not stand any physics, then lower the time step, and always with the same fault.
This is usually what the model says when it blew-up -> "Error_in_finding_100_hPa_up"
I will download a 3.6.1 version and I will do a couple of tests.

http://www2.mmm.ucar.edu/wrf/users/docs/user_guide_V3/users_guide_chap10.htm#TCBogus

The Users Guide mentions that users with multiple domains should horizontally-interpolate the generated meteorological fields to the fine-grid domains.
Click to expand...

I didn't do this. How it is performed?.. thank you.


START PROCESSING DATA
Processing time --- 2015-08-12_00:00:00
Found the right date - continue
Processing time --- 2015-08-12_00:05:00
Found the right date - continue
Troubles finding level 100 above ground.
Problems first occur at ( 2, 2)
Surface pressure = NaN hPa.
Error_in_finding_100_hPa_up
Click to expand...
 
Hi kwerner,

I have already downloaded and compiled WRF 3.6.1 model version...

The Users Guide mentions that users with multiple domains should horizontally-interpolate the generated meteorological fields to the fine-grid domains.
Click to expand...

How it is performed?.. ->

In a @wrfhelp response with date; March 21, 2009 ->

If you have more than one domain, just set fine_input_stream = 2 for the nests. This will tell the model to use the high resolution static data for the nests, but interpolate the met data from domain 1.
http://forum.wrfforum.com/viewtopic.php?f=12&t=929
Click to expand...

Would I have to do this in namelist.input? ->

Code: 
&time_control
input_from_file = .true.,.true.,.true.,.true.,
fine_input_stream   = 0,2,2,2,

Would I also have to do this? ->

Code: 
&bdy_control
spec_bdy_width = 5,
specified  = .T.,.F.,.F.,F.,
nested  = .F.,.T.,.T.,.T.,

Thank you..

L Alvarez
 
Hi,
I apologize for the delay. We are preparing for an overseas tutorial for next week and have gotten a bit behind with forum posts. I am not very familiar with the tc bogus process, but the response you received from wrfhelp (likely someone other than me) seems reasonable. I don't think you would need to make any changes to the &bdy_control section, though. Perhaps just adding 'fine_input_stream = 0, 2, 2, 2' to the &time_control section would be enough.
 
kwerner said:
Hi,
Unfortunately we are aware of a problem with the TC Bogus program beginning in V3.7. The problem causes a change to the soil information, making it corrupt. There is an option to get around it - you can take all the necessary fields from the new auxinput file and overwrite them in the original met_em file so that they maintain soil information. You would need to use some other program, such as NCL to do this, though.

A couple of other things to note: I notice that you're using 4 domains. The Users Guide mentions that users with multiple domains should horizontally-interpolate the generated meteorological fields to the fine-grid domains. Did you do that? It also mentions that the tc.exe program should be built serially. Was your program built serially?
Click to expand...


Hi, may I know how exactly do you do this one? Which NCL commands because unfortunately, I cannot also run my ./tc.exe. I tried to overwrite my aux_met_em to met_em but it stops after a while. I do not have problem in simulating without bogussing.
 
You must log in or register to reply here.
Top