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

ACM2 error :RIBX never exceeds RIC

bbbbbdst

New member
I use pxlsm and can successfully run the ACM2 scheme with three levels of nesting (36.12.04km).
But when I run 12km and 04km separately, I got this error. 36km can run successfully.
Model crashed immediately

Code:
FATAL CALLED FROM FILE:  <stdin>  LINE:     415
  RIBX never exceeds RIC, RIB(i,kte) =   0.0000000E+00 THETAV(i,1) =    300.0000     MOL=  -2.537926     TCONV =    2.587728     WST =    1.032219     KMIX =           23 UST =   0.2067498     TST =   -1.058871     U,V =   0.0000000E+00  0.0000000E+00 I,J=           1           1

I modified the timestep, but it didn't work. I changed the land type (mentioned in other posts in the forum), but it didn't work. I only use wrfbdy and wrfinput as input files, I use WRF4.3 and WPS4.4. The 04km simulation domain is located in the river plain, and there should be no very steep terrain.All other scheme can run successfully, but I really need PX lsm. I am using FNL data.

How can I fix this problem? Looking forward to your help! ! !

Code:
&time_control
 run_days                            = 28,
 run_hours                           = 0,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2020,
 start_month                         = 08,
 start_day                           = 05, 
 start_hour                          = 00, 
 end_year                            = 2020,
 end_month                           = 09,  
 end_day                             = 01,  
 end_hour                            = 00,  
 interval_seconds                    = 21600
 input_from_file                     = .true.,
 history_interval                    = 60,  60,   60,
 frames_per_outfile                  = 24,  24,   24,
 restart                             = .false.,
 restart_interval                    = 1440,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 write_hist_at_0h_rst                = .true.
 /

 &domains
 time_step                           = 8,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 e_we                                = 205,
 e_sn                                = 229,
 e_vert                              = 24,
 eta_levels                          = 1.000,0.995,0.988,0.980,0.970,0.956,     ! full sigma, bottom-up,
                                       0.938,0.916,0.893,0.868,0.839,0.808,     ! start with 1.0, end
                                       0.777,0.744,0.702,0.648,0.582,0.500,     ! with 0.0
                                       0.400,0.300,0.200,0.120,0.052,0.000,
 p_top_requested                     = 10000,
 num_metgrid_levels                  = 34,
 num_metgrid_soil_levels             = 4,
 dx                                  = 4000,
 dy                                  = 4000,
 grid_id                             = 1,   
 parent_id                           = 0,  
 i_parent_start                      = 0,    
 j_parent_start                      = 0,  
 parent_grid_ratio                   = 1,  
 parent_time_step_ratio              = 1,   
 feedback                            = 1,
 smooth_option                       = 0
 /

 &physics
 mp_physics                          = 10,
 mp_zero_out                         = 2,
 mp_zero_out_thresh                  = 1.0e-8,
 ra_lw_physics                       = 4,
 ra_sw_physics                       = 4,
 radt                                = 4,
 co2tf                               = 1, 
 sf_sfclay_physics                   = 7,
 num_soil_layers                     = 2,
 pxlsm_smois_init                    = 1,
 pxlsm_modis_veg                     = 1,
 sf_surface_physics                  = 7,
 sf_urban_physics                    = 0,
 bl_pbl_physics                      = 7,
 bldt                                = 0,
 cu_physics                          = 1,
 kfeta_trigger                       = 2
 cudt                                = 0,
 ishallow                            = 0,
 shcu_physics                        = 0,
 prec_acc_dt                         = 60,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 cu_rad_feedback                     = .true.,
 surface_input_source                = 1,
 num_land_cat                        = 20,
 num_soil_cat                        = 16,
 sst_update                          = 0,
 seaice_threshold                    = 100,
 slope_rad                           = 1,
 topo_shading                        = 1,
 shadlen                             = 25000.,
 do_radar_ref                        = 1,
 grav_settling                       = 0,
 /

 &fdda
 /

 &dynamics
 hybrid_opt                          = 2,
 w_damping                           = 1,
 diff_opt                            = 1,
 km_opt                              = 4,
 diff_6th_opt                        = 2,
 diff_6th_factor                     = 0.12,
 damp_opt                            = 3,
 base_temp                           = 290.
 zdamp                               = 5000.,
 dampcoef                            = 0.05,
 khdif                               = 0,
 kvdif                               = 0,
 non_hydrostatic                     = .true.,
 moist_adv_opt                       = 2,
 tke_adv_opt                         = 2,
 scalar_adv_opt                      = 2,
 use_theta_m                         = 1,
 /

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /
 
Last edited:
Please clarify the following questions:
(1) what WRF version did you run?
(2) What is your forcing data?
(3) How long did the model run before it crashed?
(4) Is there any difference in the options between your single-domain 4km case and the 36-12-4km nested case?
Thanks.
 
Ming Chen said:
Please clarify the following questions:
(1) what WRF version did you run?
(2) What is your forcing data?
(3) How long did the model run before it crashed?
(4) Is there any difference in the options between your single-domain 4km case and the 36-12-4km nested case?
Thanks.
1、WRF4.3 and WPS4.4
2、I don't know what forcing data is. Could you tell me how I can find it?
3、The model crashed at the very beginning, no timesteps have passed.
4、Namelist doesn't make any difference other than nested settings.
 
Ming Chen said:
Please clarify the following questions:
(1) what WRF version did you run?
(2) What is your forcing data?
(3) How long did the model run before it crashed?
(4) Is there any difference in the options between your single-domain 4km case and the 36-12-4km nested case?
Thanks.
I change nested case namelist to single-domian case directly. wrfinut and wrfbdy used in single-case are producted by nested namelist (real.exe).
This is my nested case namelist:
Code:
 &time_control
 run_days                            = 28,
 run_hours                           = 0,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2020, 2020, 2020,
 start_month                         = 08,   08,   08,
 start_day                           = 05,  05,   05,
 start_hour                          = 00,   00,   00,
 end_year                            = 2020, 2020, 2020,
 end_month                           = 09,   09,   09,
 end_day                             = 01,   01,   01,
 end_hour                            = 00,   00,   00,
 interval_seconds                    = 21600
 input_from_file                     = .true.,.true.,.true.,
 history_interval                    = 60,  60,   60,
 frames_per_outfile                  = 24,  24,   24,
 restart                             = .false.,
 restart_interval                    = 1440,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 write_hist_at_0h_rst                = .true.
 /

 &domains
 time_step                           = 120,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 3,
 e_we                                = 186,   148,   205,
 e_sn                                = 149,   241,   229,
 e_vert                              = 24,    24,    24,
 eta_levels                          = 1.000,0.995,0.988,0.980,0.970,0.956,     ! full sigma, bottom-up,
                                       0.938,0.916,0.893,0.868,0.839,0.808,     ! start with 1.0, end
                                       0.777,0.744,0.702,0.648,0.582,0.500,     ! with 0.0
                                       0.400,0.300,0.200,0.120,0.052,0.000,
 p_top_requested                     = 10000,
 num_metgrid_levels                  = 34,
 num_metgrid_soil_levels             = 4,
 dx                                  = 36000, 12000, 4000,
 dy                                  = 36000, 12000, 4000,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 i_parent_start                      = 1,     84,    65,
 j_parent_start                      = 1,     32,    69,
 parent_grid_ratio                   = 1,     3,     3,
 parent_time_step_ratio              = 1,     3,     3,
 feedback                            = 1,
 smooth_option                       = 0
 /

 &physics
 mp_physics                          = 10,   10,   10,
 cu_physics                          = 1,     1,     0,
 ra_lw_physics                       = 4,     4,     4,
 ra_sw_physics                       = 4,     4,     4,
 bl_pbl_physics                      = 7,     7,     7,
 sf_sfclay_physics                   = 7,     7,     7,
 sf_surface_physics                  = 7,     7,     7,
 radt                                = 20,    20,    20,
 bldt                                = 0,     0,     0,
 cudt                                = 0,     0,     0,
 icloud                              = 1,
 num_land_cat                        = 20,
 sf_urban_physics                    = 0,     0,     0,
 /

 &fdda
 / 

 &dynamics
 hybrid_opt                          = 2, 
 w_damping                           = 1,
 diff_opt                            = 1,      1,      1,
 km_opt                              = 4,      4,      4,
 diff_6th_opt                        = 2,      2,      2,
 diff_6th_factor                     = 0.12,   0.12,   0.12,
 base_temp                           = 290.
 damp_opt                            = 3,
 zdamp                               = 5000.,  5000.,  5000.,
 dampcoef                            = 0.05,    0.05,    0.05, 
 khdif                               = 0,      0,      0,
 kvdif                               = 0,      0,      0,
 non_hydrostatic                     = .true., .true., .true.,
 moist_adv_opt                       = 2,      2,      2,     
 scalar_adv_opt                      = 2,      2,      2,     
 use_theta_m                         = 1,
 gwd_opt                             = 1,
 /

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /
 
If the model crashed immediately, it often indicates that the input data is wrong. By forcing data, I mean the data you ungrib, which are used to produce initial and boundary conditions for WRF. Please clarify what did you use.
 
Ming Chen said:
If the model crashed immediately, it often indicates that the input data is wrong. By forcing data, I mean the data you ungrib, which are used to produce initial and boundary conditions for WRF. Please clarify what did you use.
I use NCEP FNL data (GRIB2 6 HOURLY FILES 2007.12.06 to CURRENT )
 
If the model crashed immediately, it often indicates that the input data is wrong. By forcing data, I mean the data you ungrib, which are used to produce initial and boundary conditions for WRF. Please clarify what did you use.
Not sure have you got my reply?~ Could the FNL data be the cause of the error?
 
I am not sure whether FNL data could be the reason for the failed case. But very often if the case failed immediately, it may indicate the input data is wrong.
Another issue is that, you run a triply-nested case with 36-12-4km resolution. Please turn on cumulus scheme for the first two domains and turn it off for the 3rd domain.
 
I am not sure whether FNL data could be the reason for the failed case. But very often if the case failed immediately, it may indicate the input data is wrong.
Another issue is that, you run a triply-nested case with 36-12-4km resolution. Please turn on cumulus scheme for the first two domains and turn it off for the 3rd domain.
I set cu_physics = 1, 1, 0, . Do I need to change anything?
 
Sorry I had the wrong impression that cu_physics = 0, 0, 0. ow that you set it to cu_physics = 1, 1, 0, I believe this is the correct setting.
You may also try to run the case with single domain and set radt = 36. If it works, then run 2-domain nested case , triply- nested case...
Such tests can tell other the data is fine, whether the large difference between the 4km resolution and the 1-degreee FNL is an issue..
 
Top