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

Max map factor in domain 1 = 1.01. Scale the dt in the model accordingly.

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.

dominic

New member
Hi everyone, I'm always stuck in real.exe

From the file error this is displayed:

Code:
 module_io_quilt_old.F        2931 T
 Ntasks in X            1 , ntasks in Y            1
*************************************
Configuring physics suite 'conus'

         mp_physics:      8
         cu_physics:      6
      ra_lw_physics:      4
      ra_sw_physics:      4
     bl_pbl_physics:      2
  sf_sfclay_physics:      2
 sf_surface_physics:      2
*************************************
  Domain # 1: dx =  5000.000 m
   max_ts_level must be <= number of znu half layers
   max_ts_level is reset to the number of znu half layers
REAL_EM V4.2.2 PREPROCESSOR
 *************************************
 Parent domain
 ids,ide,jds,jde            1         300           1         276
 ims,ime,jms,jme           -4         305          -4         281
 ips,ipe,jps,jpe            1         300           1         276
 *************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             766564788  bytes allocated
d01 2022-03-20_01:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.2
d01 2022-03-20_01:00:00  Input data is acceptable to use: met_em.d01.2022-03-20_01:00:00.nc
 metgrid input_wrf.F first_date_input = 2022-03-20_01:00:00
 metgrid input_wrf.F first_date_nml = 2022-03-20_01:00:00
d01 2022-03-20_01:00:00 Timing for input          0 s.
d01 2022-03-20_01:00:00          flag_soil_levels read from met_em file is  1
Max map factor in domain 1 =  1.01. Scale the dt in the model accordingly.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1208
p_top_requested < grid%p_top possible from data
-------------------------------------------
 
the metgrid files are displayed correctly and with ncview they look fine

namelist.wps

Code:
&share
 wrf_core = 'ARW',
 max_dom  = 1,
 start_date = '2022-03-20_00:00:00',
 end_date   = '2022-03-20_06:00:00',
 interval_seconds = 3600,
 io_form_geogrid  = 2,
/

&geogrid
 parent_id         =   1,
 parent_grid_ratio =   1,
 i_parent_start    =   1,
 j_parent_start    =   1,
 e_we          = 100,
 e_sn          = 99,
 dx = 3200,
 dy = 3200,
 map_proj  = 'lambert',
 ref_lat   =  44.00,
 ref_lon   =  15.00,
 truelat1  =  44.00,
 truelat2  =  44.00,
 stand_lon =  15.00,
 geog_data_res  = 'nlcd2011_30m+gtopo_30s+default',
 geog_data_path = '../WPS_GEOG',
/

&ungrib
 out_format = 'WPS',
 prefix = 'FILE',
/

&metgrid
 fg_name = 'FILE',
 io_form_metgrid = 2, 
/

namelist.input

Code:
 &time_control
 run_days                            = 0,
 run_hours                           = 6,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2022, 2000, 2000,
 start_month                         = 03,   01,   01,
 start_day                           = 20,   24,   24,
 start_hour                          = 01,   12,   12,
 end_year                            = 2022, 2000, 2000,
 end_month                           = 03,   01,   01,
 end_day                             = 20,   25,   25,
 end_hour                            = 06,   12,   12,
 interval_seconds                    = 3600
 input_from_file                     = .true.,.true.,.true.,
 history_interval                    = 180,  60,   60,
 frames_per_outfile                  = 1000, 1000, 1000,
 restart                             = .false.,
 restart_interval                    = 7200,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 /

 &domains
 time_step                           = 60,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 e_we                                = 100,    112,   94,
 e_sn                                = 99,    97,    91,
 e_vert                              = 1,    33,    33,
 p_top_requested                     = 10000,
 num_metgrid_levels                  = 1,
 num_metgrid_soil_levels             = 6,
 dx                                  = 3200,
 dy                                  = 3200,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 i_parent_start                      = 1,     31,    30,
 j_parent_start                      = 1,     17,    30,
 parent_grid_ratio                   = 1,     3,     3,
 parent_time_step_ratio              = 1,     3,     3,
 feedback                            = 1,
 smooth_option                       = 0
 /

 &physics
 physics_suite                       = 'CONUS'
 mp_physics                          = -1,    -1,    -1,
 cu_physics                          = -1,    -1,     0,
 ra_lw_physics                       = -1,    -1,    -1,
 ra_sw_physics                       = -1,    -1,    -1,
 bl_pbl_physics                      = -1,    -1,    -1,
 sf_sfclay_physics                   = -1,    -1,    -1,
 sf_surface_physics                  = -1,    -1,    -1,
 radt                                = 30,    30,    30,
 bldt                                = 0,     0,     0,
 cudt                                = 5,     5,     5,
 icloud                              = 1,
 num_land_cat                        = 21,
 sf_urban_physics                    = 0,     0,     0,
 /

 &fdda
 /

 &dynamics
 hybrid_opt                          = 2, 
 w_damping                           = 0,
 diff_opt                            = 1,      1,      1,
 km_opt                              = 4,      4,      4,
 diff_6th_opt                        = 0,      0,      0,
 diff_6th_factor                     = 0.1,   0.12,   0.12,
 base_temp                           = 290.
 damp_opt                            = 3,
 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,     
 gwd_opt                             = 1,      1,      0,
 /

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /


where can the problem be?
 
The specific error is:
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1208
p_top_requested < grid%p_top possible from data
-------------------------------------------

Can you try setting p_top_requested = 5000 (which is the default setting) and see if that works any better?
 
kwerner said:
The specific error is:
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1208
p_top_requested < grid%p_top possible from data
-------------------------------------------

Can you try setting p_top_requested = 5000 (which is the default setting) and see if that works any better?

I tried to put as value of p_top 5000 or 10000 but I always return the same error:


taskid: 0 hostname: dominic-server
module_io_quilt_old.F 2931 T
Ntasks in X 1 , ntasks in Y 1
Domain # 1: dx = 3200.000 m
REAL_EM V4.2.2 PREPROCESSOR
*************************************
Parent domain
ids,ide,jds,jde 1 234 1 237
ims,ime,jms,jme -4 239 -4 242
ips,ipe,jps,jpe 1 234 1 237
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1 , 1157918988 bytes allocated
d01 2022-03-20_01:00:00 Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.2
d01 2022-03-20_01:00:00 Input data is acceptable to use: met_em.d01.2022-03-20_01:00:00.nc
metgrid input_wrf.F first_date_input = 2022-03-20_01:00:00
metgrid input_wrf.F first_date_nml = 2022-03-20_01:00:00
d01 2022-03-20_01:00:00 Timing for input 0 s.
d01 2022-03-20_01:00:00 flag_soil_levels read from met_em file is 1
Max map factor in domain 1 = 1.00. Scale the dt in the model accordingly.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1208
p_top_requested < grid%p_top possible from data
-------------------------------------------
 
The section of code that produces this error is found in dyn_em/module_initialize_real.F.
Code:
IF ( p_top_requested .LT. grid%p_top ) THEN
               print *,'p_top_requested = ',p_top_requested
               print *,'allowable grid%p_top in data   = ',grid%p_top
               CALL wrf_error_fatal ( 'p_top_requested < grid%p_top possible from data' )
            END IF
As you can see, you should be getting printed values for p_top_requested and the allowable p_top from the data. I assume you're pasting the error message from rsl.error.0000. Can you also check in rsl.out.0000 to see if you're given those values? If you are unsure, you can attach them for me to see. Thanks!
 
Top