[Wrf-users] Restart file not produced

Yaasiin Oozeer yaasiinoozeer at gmail.com
Mon Jun 22 12:05:55 MDT 2015


Hi everyone!

I've been trying to run WRF-Chem with chem_opt = 11 (RADM2_MADE/SORGAM) on
a cluster and while the wrfout file was produced, the wrfrst file was not
produced.

The rsl files do include the line: "med_restart_out: opening
wrfrst_d01_2013-06-15_01:00:00 for writing" but unfortunately the restart
file does not get written. I've checked some other posts where the same
issue was raised. Here is the error I've been getting in the rsl file after
debugging:


<http://forum.wrfforum.com/viewtopic.php?f=6&t=9018#>d01
2013-06-15_00:15:00  output_wrf: end, fid =            5
d01 2013-06-15_00:15:00 back from outsub in open_w_dataset
d01 2013-06-15_00:15:00 calling wrf_open_for_write_commit in open_w_dataset
d01 2013-06-15_00:15:00 module_io.F: in wrf_open_for_write_commit
d01 2013-06-15_00:15:00  NetCDF error: NetCDF: One or more variable sizes
violate format constraints
d01 2013-06-15_00:15:00  NetCDF error in ext_ncd_open_for_write_commit
wrf_io.F90, line        1442
d01 2013-06-15_00:15:00 back from wrf_open_for_write_commit in
open_w_dataset
med_restart_out: opening wrfrst_d01_2013-06-15_00:15:00 for writing



I've tried the following to solve the problem:

1. Ensure that the file external/io_netcdf/wrf_io.F90 contained the
following line stat = NF_CREATE(FileName, IOR(NF_CLOBBER,NF_64BIT_OFFSET),
DH%NCID)

2. Use io_form_restart = 102 to ensure that the 'divided' restart files
would be less than 2GB

3. Use a version of netcdf that supports 64bit

Unfortunately none of the 3 things above resolved the problem. Has anyone
come across the same issue and found a working solution for it? I'd really
appreciate some help on this!



My namelist is as follows:

&time_control
 run_days                            = 1,
 run_hours                           = 0,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2013, 2013,
 start_month                         = 06, 06,
 start_day                           = 15, 15,
 start_hour                          = 00, 00,
 start_minute                        = 00, 00,
 start_second                        = 00, 00,
 end_year                            = 2013, 2013,
 end_month                           = 06, 06,
 end_day                             = 16, 16,
 end_hour                            = 00, 00,
 end_minute                          = 00, 00,
 end_second                          = 00, 00,
 interval_seconds                    = 21600,
 input_from_file                     = .true., .true.,
 history_interval                    = 60, 60,
 frames_per_outfile                  = 72, 72,
 restart                             = .False.,
 restart_interval                    = 1440,
 write_hist_at_0h_rst                = .True.,
 io_form_history                     = 2,
 io_form_restart                     = 102,
 io_form_input                       = 2,
 io_form_boundary                    = 2,
 auxinput5_inname                    = 'wrfchemi_d<domain>',
 auxinput6_inname                    = 'wrfbiochemi_d<domain>',
 auxinput7_inname                    = 'wrffirechemi_d<domain>',
 auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
 auxinput12_inname                   = 'wrf_chem_input',
 auxinput13_inname                   = 'wrfchemv_d<domain>',
 auxinput7_interval_m                = 1440, 1440,
 frames_per_auxinput7                = 1, 1,
 io_form_auxinput2                   = 2,
 io_form_auxinput5                   = 2,
 io_form_auxinput6                   = 2,
 io_form_auxinput7                   = 2,
 io_form_auxinput8                   = 2,
 io_form_auxinput12                  = 0,
 io_form_auxinput13                  = 0,
 debug_level                         = 0,
 auxinput1_inname                    = "met_em.d<domain>.<date>",
 /

 &dfi_control
/

 &domains
 time_step                           = 50,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 s_we                                = 1, 1,
 e_we                                = 250, 310,
 s_sn                                = 1, 1,
 e_sn                                = 250, 310,
 e_vert                              = 51, 51,
 eta_levels                          = 1.000, 0.998, 0.994, 0.988, 0.984,
0.981, 0.976, 0.971, 0.964, 0.958, 0.947, 0.943, 0.932, 0.925, 0.903,
0.879, 0.853, 0.824, 0.792,                   0.758, 0.737, 0.723, 0.685,
0.647, 0.608, 0.568, 0.527, 0.487, 0.447, 0.408, 0.370, 0.333, 0.297,
0.263, 0.231, 0.200, 0.171, 0.144,                   0.119, 0.096, 0.075,
0.062, 0.056, 0.047, 0.039, 0.031, 0.024, 0.018, 0.011, 0.005, 0.000,
 num_metgrid_levels                  = 61,
 num_metgrid_soil_levels             = 4,
 dx                                  = 9000, 3000,
 dy                                  = 9000, 3000,
 grid_id                             = 1, 2,
 parent_id                           = 1, 1,
 i_parent_start               = 1, 76
 j_parent_start               = 1, 78
 parent_grid_ratio                   = 1, 3,
 parent_time_step_ratio              = 1, 3,
 p_top_requested                     = 5000,
 numtiles                 = 16,
 feedback                            = 1,
 smooth_option                       = 0
 zap_close_levels                    = 50
 interp_type                         = 1
 t_extrap_type                       = 2
 force_sfc_in_vinterp                = 0
 use_levels_below_ground             = .true.
 use_surface                         = .true.
 lagrange_order                      = 1
 /
 sfcp_to_sfcp                        = .true.

 &physics
 mp_physics                          = 10, 10,
 progn                               = 1, 1,
 ra_lw_physics                       = 4, 4,
 ra_sw_physics                       = 4, 4,
 radt                                = 27, 27,
 sf_sfclay_physics                   = 2, 2,
 sf_surface_physics                  = 2, 2,
 bl_pbl_physics                      = 2, 2,
 bldt                                = 0, 0,
 cu_physics                          = 5, 0,
 cugd_avedx                       = 1,
 cu_diag                             = 1, 0,
 cudt                                = 0, 0, 0,
 ishallow                            = 1,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 num_land_cat                        = 20,
 sf_urban_physics                    = 0,
 mp_zero_out                         = 2,
 mp_zero_out_thresh                  = 1.e-12
 maxiens                             = 1,
 maxens                              = 3,
 maxens2                             = 3,
 maxens3                             = 16,
 ensdim                              = 144,
 cu_rad_feedback                     = .true., .false.,
 /


 &fdda
 grid_fdda                    =  1, 0,
 gfdda_inname                 = "wrffdda_d<domain>",
 gfdda_end_h                  =  24,24,
 gfdda_interval_m             = 360, 360,
 fgdt                         =   0, 0,
 if_no_pbl_nudging_uv         =   0, 0,
 if_no_pbl_nudging_t          =   0, 0,
 if_no_pbl_nudging_q          =   0, 0,
 if_zfac_uv                   =   0, 0,
  k_zfac_uv                   =  18, 18,
 if_zfac_t                    =   0, 0,
  k_zfac_t                    =  18, 18,
 if_zfac_q                    =   0, 0,
  k_zfac_q                    =  18, 18,
 guv                          = 0.0003, 0.0003,
 gt                           = 0.0003, 0.0003,
 gq                           = 0.0003, 0.0003,
 if_ramping                   = 1,
 dtramp_min                   = 60.0,
 io_form_gfdda                = 2,
 /

 &dynamics
 rk_ord                              = 3,
 w_damping                           = 1,
 diff_opt                            = 1,
 km_opt                              = 4,
 diff_6th_opt                        = 0,
 diff_6th_factor                     = 0.12,
 base_temp                           = 290.
 damp_opt                            = 0,
 zdamp                               = 5000., 5000.,
 dampcoef                            = 0.01, 0.01,
 khdif                               = 0, 0,
 kvdif                               = 0, 0,
 non_hydrostatic                     = .True., .True.,
 moist_adv_opt                       = 2, 2,
 scalar_adv_opt                      = 2, 2,
 chem_adv_opt                        = 2, 2,
 tke_adv_opt                         = 2, 2,
 time_step_sound                     = 4, 4,
 h_mom_adv_order                     = 5, 5,
 v_mom_adv_order                     = 3, 3,
 h_sca_adv_order                     = 5, 5,
 v_sca_adv_order                     = 3, 3,
 /

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 2,
 nio_groups = 1,
 /

  &chem
 kemit                               = 1,
 chem_opt                            = 11, 11,
 bioemdt                             = 0, 0,
 photdt                              = 30, 30,
 chemdt                              = 2.5, 2.5,
 io_style_emissions                  = 1,
 emiss_inpt_opt                      = 1, 1,
 emiss_opt                           = 5, 5,
 emiss_opt_vol                       = 0, 0,
 emiss_ash_hgt                       = 20000.,
 chem_in_opt                         = 1, 0,
 phot_opt                            = 1, 1,
 gas_drydep_opt                      = 1, 1,
 aer_drydep_opt                      = 1, 1,
 bio_emiss_opt                       = 0, 0,
 dust_opt                            = 1,
 dmsemis_opt                         = 1,
 seas_opt                            = 2,
 gas_bc_opt                          = 1, 1,
 gas_ic_opt                          = 1, 1,
 aer_bc_opt                          = 1, 1,
 aer_ic_opt                          = 1, 1,
 gaschem_onoff                       = 1, 1,
 aerchem_onoff                       = 1, 1,
 wetscav_onoff                       = 1, 1,
 cldchem_onoff                       = 1, 1,
 vertmix_onoff                       = 1, 1,
 chem_conv_tr                        = 1, 0,
 biomass_burn_opt                    = 1, 1,
 plumerisefire_frq                   = 120, 120,
 have_bcs_chem                       = .True., .false.,
 aer_ra_feedback                     = 1, 1,
 aer_op_opt                          = 1, 1,
 opt_pars_out                        = 0,
 diagnostic_chem                     = 0,
 /

 &tc
 /

 &noah_mp
 /

 &scm
 /

 &fire
 /

  &diags
 /



Thanks everyone!

Yaasiin
The University of Nottingham
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20150623/8e903fce/attachment-0001.html 


More information about the Wrf-users mailing list