WRF v4.0 w. ERA5 - "FLERCHINGER" error

Any issues with the actual running of the WRF.

WRF v4.0 w. ERA5 - "FLERCHINGER" error

Postby orndui » Mon Aug 20, 2018 7:04 am

First-time user here.

I'm trying to run WRF v4.0 using ERA5 reanalysis data as input. I've managed to successfully run WPS and real.exe but am now running into problems running wrf.exe.

I've followed this blog post and have tried to simplify my run down to no nests, and a serial build (if that's of any significance). As you can see from the WRF output log below, wrf.exe runs for 30 minutes while returning this Flerchinger warning before finally exiting with a segmentation fault:

WRF.log:
Code: Select all
*************************************
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
*************************************
WRF V4.0 MODEL
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             465766088  bytes allocated
  med_initialdata_input: calling input_input
   Input data is acceptable to use: wrfinput_d01
Timing for processing wrfinput file (stream 0) for domain        1:    0.32761 elapsed seconds
Max map factor in domain 1 =  1.00. Scale the dt in the model accordingly.
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
 LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33  CATEGORIES           2  SEASONS WATER CATEGORY =           17  SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
 LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20  CATEGORIES
 INPUT SOIL TEXTURE CLASSIFICATION = STAS
 SOIL TEXTURE CLASSIFICATION = STAS FOUND          19  CATEGORIES
ThompMP: read qr_acr_qg.dat instead of computing
ThompMP: read qr_acr_qs.dat instead of computing
ThompMP: read freezeH2O.dat instead of computing
Timing for Writing wrfout_d01_2016-12-02_00:00:00 for domain        1:    1.05070 elapsed seconds
d01 2016-12-02_00:00:00  Input data is acceptable to use: wrfbdy_d01
Timing for processing lateral boundary for domain        1:    1.11091 elapsed seconds
 Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE     74 JS      1 JE     61
WRF NUMBER OF TILES =   1
Timing for main: time 2016-12-02_00:00:10 on domain   1:    8.35096 elapsed seconds
Timing for main: time 2016-12-02_00:00:20 on domain   1:    1.22043 elapsed seconds
...
Timing for main: time 2016-12-02_00:03:40 on domain   1:    1.33022 elapsed seconds
 Flerchinger USEd in NEW version. Iterations=          10
...
 Flerchinger USEd in NEW version. Iterations=          10
Timing for main: time 2016-12-02_00:03:50 on domain   1:    1.21784 elapsed seconds
 Flerchinger USEd in NEW version. Iterations=          10
 Flerchinger USEd in NEW version. Iterations=          10
...
 Flerchinger USEd in NEW version. Iterations=          10
Timing for main: time 2016-12-02_00:30:00 on domain   1:    0.92164 elapsed seconds

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x1473b38842da in ???
#1  0x1473b3883503 in ???
#2  0x1473b2f00f1f in ???
#3  0x55b396b6e605 in ???
#4  0x55b396b8d377 in ???
#5  0x55b396bb0081 in ???
#6  0x55b396bbf8c5 in ???
#7  0x55b39673cc25 in ???
#8  0x55b3967f855a in ???
#9  0x55b396641f34 in ???
#10  0x55b396563e46 in ???
#11  0x55b395a81906 in ???
#12  0x55b395a17469 in ???
#13  0x55b395a16dbe in ???
#14  0x1473b2ee3b96 in ???
#15  0x55b395a16df9 in ???
#16  0xffffffffffffffff in ???


namelist.input
Code: Select all
 &time_control
 run_days                            = 1,
 run_hours                           = 18,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2016, 2016,
 start_month                         = 12,   12,
 start_day                           = 02,   02,
 start_hour                          = 00,   00,
 end_year                            = 2016, 2016,
 end_month                           = 12,   12,
 end_day                             = 03,   03,
 end_hour                            = 18,   18,
 interval_seconds                    = 1800,
 input_from_file                     = .true.,.true.,
 history_interval                    = 180,  60,
 frames_per_outfile                  = 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                           = 10,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 e_we                                = 74,    112,
 e_sn                                = 61,    97,
 e_vert                              = 138,    138,
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 138,
 num_metgrid_soil_levels             = 4,
 dx                                  = 30000, 10000,
 dy                                  = 30000, 10000,
 grid_id                             = 1,     2,
 parent_id                           = 1,     1,
 i_parent_start                      = 1,     31,
 j_parent_start                      = 1,     17,
 parent_grid_ratio                   = 1,     3,
 parent_time_step_ratio              = 1,     3,
 feedback                            = 1,
 smooth_option                       = 0
 /

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

 &fdda
 /

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

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /


namelist.wps:
Code: Select all
&share
 wrf_core = 'ARW',
 max_dom = 2,
 start_date = '2016-12-02_00:00:00','2016-12-02_00:00:00',
 end_date   = '2016-12-03_18:00:00','2016-12-03_18:00:00',
 interval_seconds = 151200,
 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              =  74, 112,
 e_sn              =  61,  97,
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! The default datasets used to produce the MAXSNOALB and ALBEDO12M
 ! fields have changed in WPS v4.0. These fields are now interpolated
 ! from MODIS-based datasets.
 !
 ! To match the output given by the default namelist.wps in WPS v3.9.1,
 ! the following setting for geog_data_res may be used:
 !
 ! geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default', 'maxsnowalb_ncep+albedo_ncep+default',
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !
 geog_data_res = 'default','default',
 dx = 30000,
 dy = 30000,
 map_proj = 'lambert',
 ref_lat   =  51.97,
 ref_lon   =   4.93,
 truelat1  =  40.0,
 truelat2  =  60.0,
 stand_lon =   5.0,
 geog_data_path = '/media/orndui/Bigdata/WPS_GEOG/'
/

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

&metgrid
 fg_name = 'FL','PRES',
 io_form_metgrid = 2,
/


Following a suggestion in a comment on this post I've looked at the SH2O and SMOIS variables in my wrfinput file. SMOIS looks ok but SH2O is uniformly equal to 0 over land at all levels, and 1 over sea. I don't know if this is normal? I haven't made any changes to the LSM that I know of.
I will try to compile WRF v3.9 to see if that goes any better and will post the results here for future reference.

Any tips and suggestions greatly appreciated!
orndui
 
Posts: 2
Joined: Fri Aug 17, 2018 10:41 am

Re: WRF v4.0 w. ERA5 - "FLERCHINGER" error

Postby Xun » Fri Sep 07, 2018 12:44 pm

Hi orndui,

I also had this problem before. This error is caused by the negative value in SMOIS field, which also caused artifact in the TSK and T2 in my output. In dry areas, the very small soil moisture is set to negative in ERA5. I replaced the negative value with a small positive value and then the error and the artifact were gone.

cheers!
Xun
 
Posts: 1
Joined: Fri Sep 07, 2018 12:36 pm

Re: WRF v4.0 w. ERA5 - "FLERCHINGER" error

Postby tt7fans » Fri Sep 28, 2018 3:19 am

Xun wrote:Hi orndui,

I also had this problem before. This error is caused by the negative value in SMOIS field, which also caused artifact in the TSK and T2 in my output. In dry areas, the very small soil moisture is set to negative in ERA5. I replaced the negative value with a small positive value and then the error and the artifact were gone.

cheers!


No, not caused by this.
I met the same problem when I use noah scheme, after numerous 'Flerchinger ... ' , my wrf.exe crashed. I'm still work on it, the SMOIS in wrfinput appears ok, they are not all zero nor nagative.
The SH2O looks too small but I don't know whether this is the reason.
tt7fans
 
Posts: 1
Joined: Mon Nov 13, 2017 4:36 am


Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 12 guests