real.exe not working for global run

Any issues with the actual running of the WRF.

real.exe not working for global run

Postby kala » Sat May 20, 2017 8:19 pm

I'm running WRF v3.9 (I've used v3.8 before and global runs worked then, I just upgraded) and I was able to successfully run a lambert projection regional simulation through both real.exe and wrf.exe, but when I change the namelists to make it a global lat-lon projection run, I get this error in real.exe:

Code: Select all
 Domain  1: Current date being processed: 2017-05-27_00:00:00.0000, which is loop #   1 out of    2
 configflags%julyr, %julday, %gmt:        2017         147   0.00000000   
 metgrid input_wrf.F first_date_input = 2017-05-27_00:00:00
 metgrid input_wrf.F first_date_nml = 2017-05-27_00:00:00
d01 2017-05-27_00:00:00 Timing for input          2 s.
d01 2017-05-27_00:00:00          flag_soil_layers read from met_em file is  1
d01 2017-05-27_00:00:00 Turning off use of MAX WIND level data in vertical interpolation
d01 2017-05-27_00:00:00 Turning off use of TROPOPAUSE level data in vertical interpolation
Max map factor in domain 1 = *****. Scale the dt in the model accordingly.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:     564
Maybe you do not have the new map factors, try re-running geogrid
-------------------------------------------


I have thus far not been able to determine what the problem is or could be. I suspected that it had something to do with the projection or dx and dy being incorrect, but they seem to be fine. My concern is the warning about the map factor which has "*****". When searching around, other threads have floats (like 1.03, 1.09, etc) instead, but this may just be as a result of the global run and not really a problem. Has anyone else been able to do a global run with v3.9?

This is my namelist.wps:
Code: Select all
&share
 wrf_core = 'ARW',
 max_dom = 1,
 start_date = '2017-05-27_00:00:00','2006-08-16_12:00:00',
 end_date   = '2017-05-27_00:00:00','2006-08-16_12:00:00',
 interval_seconds = 21600
 io_form_geogrid = 2,
/

&geogrid
 parent_id         =   1,   1,
 parent_grid_ratio =   1,   3,
 i_parent_start    =   1,  31,
 j_parent_start    =   1,  17,
 e_we              =  463, 112,
 e_sn              =  232,  97,
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! This namelist is specific for use with the lowest resolution option for
 ! each field in the static geographic tar file.  It is mandatory to use
 ! the below settings for geog_data_res.
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !
 geog_data_res = '2deg+gtopo_10m+usgs_10m+10m+nesdis_greenfrac','2deg+gtopo_10m+usgs_10m+10m+nesdis_greenfrac',
 map_proj = 'lat-lon',
 stand_lon = 0.0
 pole_lat = 90.0
 pole_lon = 0.0
 geog_data_path = '/home/kalassak/Desktop/WRF/WPS_GEOG/'
/

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

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


and this is my namelist.input:

Code: Select all
&time_control
 run_days                            = 0,
 run_hours                           = 6,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2017, 2017, 2017,
 start_month                         = 05,   03,   03,
 start_day                           = 27,   04,   04,
 start_hour                          = 00,   00,   00,
 start_minute                        = 00,   00,   00,
 start_second                        = 00,   00,   00,
 end_year                            = 2017, 2017, 2017,
 end_month                           = 05,   03,   03,
 end_day                             = 27,   07,   07,
 end_hour                            = 06,   00,   00,
 end_minute                          = 00,   00,   00,
 end_second                          = 00,   00,   00,
 history_interval                    = 180,  60,   60,
 interval_seconds                    = 21600
 input_from_file                     = .true.,.true.,.true.,
 frames_per_outfile                  = 200,  100,  100,
 restart                             = .false.,
 restart_interval_h                  = 24,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 debug_level                         = 0
 /

 &domains
 time_step                           = 360,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 e_we                                = 463,   251,   188,
 e_sn                                = 232,   99,    182,
 e_vert                              = 30,    30,    30,
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 32,
 num_metgrid_soil_levels             = 4,
 dx                                  = 86631.8, 15000,  3333.33,
 dy                                  = 86631.8, 15000,  3333.33,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 i_parent_start                      = 1,     40,    17,
 j_parent_start                      = 1,     65,    33,
 parent_grid_ratio                   = 1,     2,     5,
 parent_time_step_ratio              = 1,     2,     5,
 feedback                            = 1,
 smooth_option                       = 00
 sfcp_to_sfcp                        = .true.
 /

 &physics
 physics_suite                       = 'tropical'
 radt                                = 30,    30,    30,
 bldt                                = 0,     0,     0,
 cudt                                = 5,     5,     5,
 icloud                              = 1,
 num_soil_layers                     = 4,
 num_land_cat                        = 24,
 sf_urban_physics                    = 0,     0,     0,
 sst_skin              = 1,
 /

 &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,     
 gwd_opt                             = 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,
 /
kala
 
Posts: 5
Joined: Wed Mar 08, 2017 3:18 am

Re: real.exe not working for global run

Postby kala » Sat May 20, 2017 10:59 pm

Turns out the problem is in module_initialize_real.F. The code for assigning map factors was messed up (there were two if statements in the same block testing for the same thing, when they should be testing for different things). Basically, it was always testing whether the map projection was polar, and there was no test for the lat-lon projection at all. I copied the section of code from v3.8.1 into my 3.9 installation and recompiled, and voila, it worked.

I'm just wondering if anyone else has run into the problem in the past month, or do very few people run global sims?
kala
 
Posts: 5
Joined: Wed Mar 08, 2017 3:18 am

Re: real.exe not working for global run

Postby cyclone6 » Sun Jul 30, 2017 5:04 am

Hi,

I'm curious to know what you are doing global simulations for? I'm thinking about running some operationally, because I'm a meteorologist who just really enjoys playing with this.

Can I ask where exactly those lines of code are?

ted
cyclone6
 
Posts: 15
Joined: Tue Mar 06, 2012 2:16 am

Re: real.exe not working for global run

Postby kala » Wed Aug 02, 2017 11:19 pm

I am setting up WRF to run global simulations of an Earth-like planet that has different properties from Earth. I am a meteorology student and I'm interested in tropical cyclones so that was mostly what I was focusing on in my simulations. I'm mostly messing around with it for fun, too.

I don't exactly remember what line of code I fixed, and now I realize my solution explanation was subpar, but I'll try to explain what I think it was.

Around line 472 in v3.9's module_initialize_real.F, the ELSE IF statement is probably the same thing as the one above it, which it shouldn't be, so change it so that it's similar to v3.8's (or even the file here https://github.com/yyr/wrf/blob/master/ ... ize_real.F) like so:

Code: Select all
      ELSE IF ( ( .NOT. config_flags%map_proj .EQ. PROJ_CASSINI ) .AND. ( flag_mf_xy .EQ. 1 ) ) THEN
         IF ( grid%msfvx(its,jts) .EQ. 0 ) THEN
            CALL wrf_error_fatal ( 'Maybe you do not have the new map factors, try re-running geogrid' )
         END IF
kala
 
Posts: 5
Joined: Wed Mar 08, 2017 3:18 am

Re: real.exe not working for global run

Postby cyclone6 » Wed Aug 09, 2017 11:57 am

Thanks for your response! I just got back from SE Asia, and WRF did a nice job dealing with a few tropical cyclones that moved into the domain during my daily runs. Those runs were not global - just a couple of large domains (36-km mother and a 12-km nest).

Looking forward to seeing how WRF does over 3 or 4 days with just a single initialization hour. No nudging allowed there - its just WRF doing the forecasting...

Fingers crossed.
cyclone6
 
Posts: 15
Joined: Tue Mar 06, 2012 2:16 am


Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 3 guests