[Wrf-users] ndown not going

BasitAli Khan BasitAli.Khan at kaust.edu.sa
Mon Oct 8 04:38:28 MDT 2012


Hi guys,

I am trying to do a one-way nested run with ndown using wrfchemV3.3.1.  I followed the procedure given in http://www.mmm.ucar.edu/wrf/OnLineTutorial/CASES/NestRuns/ndown1-2.htm. I am facing problems on step 4 (Run ndown.exe). When i run ndown.exe it crashes in the very beginning but i don't get any error message.  I ran ndown.exe with debug_level =100 and the output of rsl.out.0000 is as follows:

 NDOWN_EM V3.3.1 PREPROCESSOR
   ndown_em: calling alloc_and_configure_domain coarse
  *************************************
  Parent domain
  ids,ide,jds,jde  1 661 1 400
  ims,ime,jms,jme  -4 90 -4 57
  ips,ipe,jps,jpe  1 83 1 50
  *************************************
 DYNAMICS OPTION: Eulerian Mass Coordinate
    alloc_space_field: domain  1 ,  520047484  bytes allocated
    setup_timekeeping:  set xtime to  0.0000000000E+00
    setup_timekeeping:  set julian to  139.0000000
   setup_timekeeping:  returning...
   ndown_em: calling model_to_grid_config_rec
   ndown_em: calling set_scalar_indices_from_config
   ndown_em: calling init_wrfio
   Entering ext_gr1_ioinit
   Entering ext_gr2_ioinit

After this line "Entering ext_gr2_joinit" ndown.exe crashes. I tried many things, like turned off the chemistry altogether, but nothing seems to work. I will be grateful if  someone could could  help in spotting where i am going wrong and how to fix it. Copied below is the namelist.input  that i used for ndown.exe.


&time_control

 run_days                            = 3,

 run_hours                           = 00,

 run_minutes                         = 0,

 run_seconds                         = 0,

 start_year                          = 2006, 2006, 2006,

 start_month                         = 05,   05,   05,

 start_day                           = 20,   20,   24,

 start_hour                          = 00,   00,   12,

 start_minute                        = 00,   00,   00,

 start_second                        = 00,   00,   00,

 end_year                            = 2006, 2006, 2000,

 end_month                           = 05,   05,   01,

 end_day                             = 23,   23,   25,

 end_hour                            = 00,   00,   12,

 end_minute                          = 00,   00,   00,

 end_second                          = 00,   00,   00,

 interval_seconds                    = 3600,

 input_from_file                     = .true.,.true.,.false.,  ! v.imp param for nested doms; should be .true.for each active dom.

 history_interval                    = 60,  60,   60,

 frames_per_outfile                  = 1000, 1000, 1000,

 restart                             = .false.,

 restart_interval                    = 4320,

 io_form_history                     = 2,

 io_form_restart                     = 2,

 io_form_input                       = 2,

 io_form_boundary                    = 2,

!!! auxinput5_inname                    = 'wrfchemi_d<domain>_<date>',

!!! auxinput7_inname                    = 'wrffirechemi_d<domain>_<date>',

!!! auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>_<date>',

! auxinput12_inname                  = 'wrf_chem_input',

! auxinput13_inname                  = 'wrfchemv_d<domain>',

!!! auxinput5_interval_m                = 60, 60, 31680, 60,

!!! auxinput7_interval_m                = 60, 60, 31680, 60,

!!! auxinput8_interval_m                = 60, 60, 31680, 60,

! auxinput13_interval_m              = 60, 60, 31680, 60,

 io_form_auxinput2                   = 2,

!!! io_form_auxinput5                   = 2,

!!! io_form_auxinput7                   = 2,                  ! 0, for the final wrfrun

!!! io_form_auxinput8                   = 2,                  ! 0, for the final wrfrun

!!! frames_per_auxinput5                = 1,1,

!!! frames_per_auxinput7                = 1,1,

!!! frames_per_auxinput8                = 1,1,

 debug_level                         = 0,

 auxinput1_inname                    = "met_em.d<domain>.<date>",

/



 &dfi_control

/


 &domains

 time_step                           = 30,

 time_step_fract_num                 = 0,

 time_step_fract_den                 = 1,

 max_dom                             = 2,

 s_we                                =   1,    1,     1,

 e_we                                = 661,    601,

 s_sn                                =   1,    1,     1,

 e_sn                                = 400,    601,

 e_vert                              = 51,    51,    51,

 eta_levels                          = 1, 0.9965, 0.993, 0.988, 0.983, 0.9765,

                                       0.97, 0.962, 0.954, 0.944, 0.934,

                                       0.9215, 0.909, 0.8945, 0.88, 0.8647463,

                                       0.8494927, 0.8342391, 0.8189855, 0.8037319, 0.7884783,

                                       0.7579709, 0.7020623, 0.6493446, 0.5996661, 0.5528811,

                                       0.5088496, 0.4674375, 0.4285163, 0.3919626, 0.357658,

                                       0.3254893, 0.295348, 0.2671299, 0.2407355, 0.2160695,

                                       0.1930407, 0.1715619, 0.1515494, 0.1329236, 0.1156082,

                                       0.09953026, 0.08462016, 0.07081139, 0.05804047, 0.04624676,

                                       0.0353724, 0.02536213, 0.0161632, 0.007725209, 0

 num_metgrid_levels                  = 38,

! num_metgrid_soil_levels            = 4,

 dx                                  = 7500, 1500,  3333.33,

 dy                                  = 7500, 1500,  3333.33,

 grid_id                             = 1,     2,     3,

 parent_id                           = 0,     1,     2,

 i_parent_start                      = 1,     175,    30,

 j_parent_start                      = 1,     220,    30,

 parent_grid_ratio                   = 1,     5,     3,

 parent_time_step_ratio              = 1,     5,     3,

 p_top_requested                     = 5000,

 feedback                            = 0,

 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.

!

! OPTIONS FOR ADAPTIVE TIME STEP

 use_adaptive_time_step              = .true.,

 step_to_output_time                 = .true.,

 target_cfl                          = 1.2, 1.2,1.2,

 max_step_increase_pct               = 5, 51,51,

 starting_time_step                  = -1, -1,-1,

 max_time_step                       = -1, -1,-1,

 min_time_step                       = -1, -1,-1,

 adaptation_domain                   = 1,

/


 &physics

 mp_physics                          = 2,     2,     2,  ! wetscav_onoff =1 works only with option 2

 progn                               = 1, ! prognostic cloud droplet number

 ra_lw_physics                       = 4,     4,     1,

 ra_sw_physics                       = 4,     4,     1,

 radt                                = 5,     5,      5,   ! 1 min/km of dx (10 for 10km grid); same value for all nests

 sf_sfclay_physics                   = 2,     2,     1,    ! surface layer scheme

 sf_surface_physics                  = 3,     3,     1,    ! Land surface scheme

 bl_pbl_physics                      = 2,     2,     1,    ! boundary layer scheme

 bldt                                = 0,     0,     0,    ! mins b/w bl phys if =0 call every time step

 cu_physics                          = 5,     0,     0,    ! cumulus option; 3 & 5 works with cu_diag and cu_rad_feedback

 cu_diag                             = 1,     0,     0,    ! time average cumulus clouds

 cudt                                = 0,     0,     0,    ! mins b/w cu_phys calls, =0 call every time step

 ! cugd_avedx                         = 1,     1,          ! NEVER TURN ON THIS OPTION. Num of neighbouring grid boxes over which subsidence will spread fm conv param.

 ishallow                            = 0,

 isfflx                              = 1,

 ifsnow                              = 1,

 icloud                              = 1,

 surface_input_source                = 1,

 num_soil_layers                     = 6,

 sf_urban_physics                    = 0,     0,     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.,            ! feedback fm param convection to radiation and photolysis schemes

! chem_diag     = 0,                  ! chemical tendeny diagnostic; does'nt work with dust-only setup

 /


 &fdda

 /


 &dynamics

 rk_ord                              = 3,

 w_damping                           = 1,

 diff_opt                            = 1,

 km_opt                              = 4,

 diff_6th_opt                        = 0,

 diff_6th_factor                     = 0.12, 0.12

 base_temp                           = 290.

 damp_opt                            = 0,

 zdamp                               = 5000.,  5000.,  5000.,

 dampcoef                            = 0.01,   0.01,   0.01

 khdif                               = 0,      0,      0,

 kvdif                               = 0,      0,      0,

 non_hydrostatic                     = .true., .true., .true.,

 moist_adv_opt                       = 2,      2,      0,

 scalar_adv_opt                      = 2,      2,      0,

!* chem_adv_opt                        = 2,      2,      0,

 tke_adv_opt                         = 2,      2,      0,

 time_step_sound                     = 4,      4,      4,

 h_mom_adv_order                     = 5,      5,      5,

 v_mom_adv_order                     = 3,      3,      3,

 h_sca_adv_order                     = 5,      5,      5,

 v_sca_adv_order                     = 3,      3,      3,

 /


 &bdy_control

 spec_bdy_width                      = 5,

 spec_zone                           = 1,

 relax_zone                          = 4,

 specified                           = .true., .false.,.false.,

 nested                              = .false.,.false.,.false.,

 /


 &grib2

 /


 &namelist_quilt

 nio_tasks_per_group = 0,

 nio_groups = 1,

 /


! * &chem

! * kemit                               = 1,

! * chem_opt                            = 0,

! bioemdt                            = 0,          0,    ! update time intv bioemissions,mins

! * phot_opt                            = 0,          1,    ! photolysis

! * photdt                              = 0,                ! update interval of photolysis in mins

! * chemdt                              = 5,                ! time step used by chemistry in mins

! * io_style_emissions                  = 2,

! * emiss_inpt_opt                      = 1,

! * emiss_opt                           = 5,          0,

! emiss_opt_vol                      = 0,          0,

! emiss_ash_hgt                      = 20000.,

! * chem_in_opt                         = 0,          0,

! * gas_drydep_opt                      = 1,          1,

! * aer_drydep_opt                      = 1,          1,

! * bio_emiss_opt                       = 0,          0,

! * dust_opt                            = 1,

! * dmsemis_opt                         = 1,

! * seas_opt                            = 1,

! * gas_bc_opt                          = 0,          0,

! * gas_ic_opt                          = 0,          0,

! * aer_bc_opt                          = 0,          0,

! * aer_ic_opt                          = 0,          0,

! * gaschem_onoff                       = 1,          0,

! * aerchem_onoff                       = 1,          0,

! * wetscav_onoff                       = 0,          0,    ! only works with chem_opt 9 -12

! * cldchem_onoff                       = 0,          0,    ! only works with chem_opt 9 -12

! * vertmix_onoff                       = 1,          0,

! * chem_conv_tr                        = 1,          0,

! * biomass_burn_opt                    = 1,          0,

! * plumerisefire_frq                   = 120,         0,

! * have_bcs_chem                       = .false., .false., .false.,

! * aer_ra_feedback                     = 1,                ! feedback fm aero to rad schemes; for direct/indirect effect

! * aer_op_opt                          = 1,                ! Aero optical properties

! * opt_pars_out                        = 1,                ! output optical properties

! * diagnostic_chem                     = 0,

! * /

----


Best regards,

Basit A. Khan, Ph.D.
Postdoctoral Research Fellow
Division of Physical Sciences & Engineering
Office# 3204, Level 3, Building 1,
King Abdullah University of Science & Technology
4700 King Abdullah Blvd, Box 2753, Thuwal 23955 –6900,
Kingdom of Saudi Arabia.

Office: +966(0)2 808 0276,  Mobile: +966(0)5 0860 3617
E-mail: basitali.khan at kaust.edu.sa<mailto:basitali.khan at kaust.edu.sa>
Skype name: basit.a.khan

________________________________
This message and its contents including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20121008/9bf81364/attachment-0001.html 


More information about the Wrf-users mailing list