ndown from d03 to d04 - SMOIS not carried over from wrfout?

Discuss the nesting capability within the model itself and any problems you might have run into.

ndown from d03 to d04 - SMOIS not carried over from wrfout?

Postby Wyrmux » Fri Mar 24, 2017 12:37 pm

Hi!

I'm attempting to nest a fine LES domain within a coarser setup with three domains in a valley. I've already ran the three coarser domains and obtained wrfout_d03 files at a rate of 30 min as input for ndown (I've renamed them to d01, as per the instructions from the ndown manual).

However, I come across a 'problem'. I'm not sure if it is actually a 'problem' or the actual way in which ndown works... When I run ndown and obtain my wrfinput_d01 and wrfbdy_01 files for the fourth nested LES domain, I see that the soil moisture fields (contained in wrfinput_d01) are the same as those from the reanalysis data that was used to drive the previous coarser run:

http://imgur.com/a/TZd1o
http://imgur.com/a/EQTCL

NB: The reason the moisture is much lower than the one from the reanalysis product (NARR in this case) is because NARR gave values much larger than observations, so I divided the moisture by a factor of 3 (just within the valley) before the coarser run to reduce some of the cold bias I was observing beforehand because of inappropriate H/LE partitioning.

On the other hand, when I compare e.g. the 10-m wind speed from the first minute of the LES run and the corresponding time from the wrfout_d01 file (not shown here), the agreement is there, so the initial conditions were carried down as expected (from wrfout_d01 files and not the reanalysis).

My question would be: Is SMOIS supposed to be carried down from the reanalysis and not from the wrfout_d01 files? It makes more sense to do it from the latter, since the model has spun up (for 21 h) and is now somewhat closer to reality. As it is right now, it seems like I need to divide the SMOIS fields by 3 AGAIN just to get reasonable values, which doesn't make any sense.

As for performing ndown, I also noticed that currently two schools of thought exist. Specifically, one which modifies the met_em files:

viewtopic.php?f=10&t=1729#p19757
http://www2.mmm.ucar.edu/wrf/OnLineTuto ... /ndown.htm

... and one which does not:

http://www2.mmm.ucar.edu/wrf/users/docs ... way_Nested
http://www2.mmm.ucar.edu/wrf/OnLineTuto ... /ndown.htm

I tried both, and I still get the too-moist SMOIS field carried down from NARR. Additionally, the former approach yielded distortions in the topography of d04, probably similar to this post:

viewtopic.php?f=10&t=3550

My namelist files:

namelist.wps
Code: Select all
&share
 wrf_core = 'ARW',
 max_dom = 4,
 start_date = '2006-04-29_00:00:00', '2006-04-29_00:00:00', '2006-04-29_00:00:00', '2006-04-29_21:00:00',
 end_date   = '2006-04-30_00:00:00', '2006-04-30_00:00:00', '2006-04-30_00:00:00', '2006-04-30_00:00:00',
 interval_seconds = 10800
 io_form_geogrid = 2,
/

&geogrid
 parent_id         =   1,  1,  2,  3,
 parent_grid_ratio =   1,  3   3,  5,
 i_parent_start    =   1,  51, 50, 132,
 j_parent_start    =   1,  34, 52, 107,
 e_we              =   150,  178,  220, 131,
 e_sn              =   120,  175,  223, 116,
 geog_data_res     = 'nlcd2006_30s+gtopo_2m','nlcd2006_30s+gtopo_30s','nlcd2006_30s+gtopo_30s','nlcd2011_1s+gtopo_30s',
 dx = 12150,
 dy = 12150,
 map_proj = 'mercator',
 ref_lat   =  36.50,
 ref_lon   =  -119.00,
 truelat1  =  0.0,
 truelat2  =  0.0,
 stand_lon =  -119.00,
 geog_data_path = '/static/geog_data/geog/'
/

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

&metgrid
 fg_name = 'FILE',
 constants_name = './NARR_FF:1979-11-08_00',
 io_form_metgrid = 2,
/


namelist.input.coarse
Code: Select all
 &time_control
 run_days                            = 0,
 run_hours                           = 3,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2006, 2006, 2006,  2006,
 start_month                         = 04,   04,  04,   04,
 start_day                           = 29,   29,  29,   29,
 start_hour                          = 00,   00,  00,   21,
 start_minute                        = 00,   00,   00,   00,
 start_second                        = 00,   00,   00,   00,
 end_year                            = 2006, 2006, 2006,    2006,
 end_month                           = 04,   04,  04,   04,
 end_day                             = 30,   30,  30,  30,
 end_hour                            = 00,   00,  00,  00,
 end_minute                          = 00,   00,   00,  00,
 end_second                          = 00,   00,   00,  00,
 interval_seconds                    = 10800
 input_from_file                     = .true.,.true.,.true.,.true.,
 fine_input_stream                   = 0,0,0,0,
 history_interval                    = 180,  180, 30, 30,
 frames_per_outfile                  = 1000, 1000, 2,  2,
 restart                             = .false.,
 restart_interval                    = 5000,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 debug_level                         = 0
 /

 &domains
 time_step                           = 27,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 4,
 s_we                  = 1,     1,    1, 1,
 e_we                                = 150, 178, 220, 131,
 s_sn                  = 1,     1,    1, 1,
 e_sn                                = 120, 175, 223, 116,
 e_vert                              = 65,    65,   65,  65,
 p_top_requested                     = 10000,
 num_metgrid_levels                  = 30,
 num_metgrid_soil_levels             = 4,
 dx                                  = 12150, 4050, 1350, 270,
 dy                                  = 12150, 4050, 1350, 270,
 grid_id                             = 1,     2,     3,   4,
 parent_id                           = 1,     1,     2,   3,
 i_parent_start                      = 1,  51,  50, 132,
 j_parent_start                      = 1,  34,  52, 107,
 parent_grid_ratio                   = 1,     3,   3,  5,
 parent_time_step_ratio              = 1,     3,   3,  15,
 feedback                            = 0,
 smooth_cg_topo                      = .true.,
 smooth_option                       = 0,
 numtiles                            = 2,
 /

 &physics
 mp_physics                          = 0,     0,     0,  0,
 ra_lw_physics                       = 4,     4,     4,  4,
 ra_sw_physics                       = 4,     4,     4,  4,
 radt                                = 9,    9,    9,  9,
 sf_sfclay_physics                   = 1,     1,     1,  1,
 sf_surface_physics                  = 2,     2,     2,  2,
 bl_pbl_physics                      = 1,     1,     1,  0,
 bldt                                = 0,     0,     0,  0,
 cu_physics                          = 1,     0,     0,  0,
 cudt                                = 5,     0,     0,  0,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_land_cat                        = 40,
 num_soil_layers                     = 4,
 sf_urban_physics                    = 0,     0,     0,  0,
 slope_rad                           = 1,     1,     1,  1,
 topo_shading                        = 1,     1,     1,  1,
 shadlen                             = 25000,
 /

 &fdda
 /

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

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /



namelist.input.ndown

Code: Select all
  &time_control
 run_days                            = 0,
 run_hours                           = 3,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2006, 2006,
 start_month                         = 04,   04,
 start_day                           = 29,   29,
 start_hour                          = 21,   21,
 start_minute                        = 00,   00,
 start_second                        = 00,   00,
 end_year                            = 2006, 2006,
 end_month                           = 04,   04,
 end_day                             = 30,   30,
 end_hour                            = 00,   00,
 end_minute                          = 00,   00,
 end_second                          = 00,   00,
 interval_seconds                    = 1800,
 input_from_file                     = .true.,.true.,
 fine_input_stream                   = 0,0,
 history_interval                    = 180,  1,
 frames_per_outfile                  = 1000, 1,
 restart                             = .false.,
 restart_interval                    = 5000,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 io_form_auxinput2                   = 2
 debug_level                         = 0
 /

 &domains
 time_step                           = 3,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 2,
 s_we                  = 1,    1,
 e_we                                = 220,  131,
 s_sn                  = 1,    1,
 e_sn                                = 223,  116,
 e_vert                              = 65,   65,
 p_top_requested                     = 10000,
 num_metgrid_levels                  = 30,
 num_metgrid_soil_levels             = 4,
 dx                                  = 1350, 270,
 dy                                  = 1350, 270,
 grid_id                             = 1,     2,
 parent_id                           = 1,     1,
 i_parent_start                      = 1,  132,
 j_parent_start                      = 1,  107,
 parent_grid_ratio                   = 1,  5,
 parent_time_step_ratio              = 1,  15,
 feedback                            = 0,
 smooth_cg_topo                      = .true.,
 smooth_option                       = 0,
 numtiles                            = 2,
 /

 &physics
 mp_physics                          = 0,     0,
 ra_lw_physics                       = 4,     4,
 ra_sw_physics                       = 4,     4,
 radt                                = 9,     9,
 sf_sfclay_physics                   = 1,     1,
 sf_surface_physics                  = 2,     2,
 bl_pbl_physics                      = 1,     0,
 bldt                                = 0,     0,
 cu_physics                          = 0,     0,
 cudt                                = 0,     0,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_land_cat                        = 40,
 num_soil_layers                     = 4,
 sf_urban_physics                    = 0,     0,     
 slope_rad                           = 1,     1,
 topo_shading                        = 1,     1,
 shadlen                             = 25000,
 /

 &fdda
 /

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

 &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 = 0,
 nio_groups = 1,
 /




namelist.input.les (This one is used when finally running the LES domain on its own, as well as when running real.exe when I change the met_em files, of course with history_interval set to 10800).

Code: Select all
 &time_control
 start_year                          = 2006,
 start_month                         = 04,
 start_day                           = 29,
 start_hour                          = 21,
 start_minute                        = 00,
 start_second                        = 00,
 end_year                            = 2006,
 end_month                           = 04,
 end_day                             = 30,
 end_hour                            = 00,
 end_minute                          = 00,
 end_second                          = 00,
 interval_seconds                    = 1800,
 input_from_file                     = .true.,
 history_interval_m                  = 0,
 history_interval_s                  = 10,
 frames_per_outfile                  = 1,
 restart                             = .false.,
 restart_interval                    = 5000,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 io_form_auxinput2                   = 2
 debug_level                         = 0
 /

 &domains
 time_step                           = 0,
 time_step_fract_num                 = 2,
 time_step_fract_den                 = 10,
 max_dom                             = 1,
 s_we                  = 1,
 e_we                                = 131,
 s_sn                  = 1,
 e_sn                                = 116,
 e_vert                              = 65,
 p_top_requested                     = 10000,
 num_metgrid_levels                  = 30,
 num_metgrid_soil_levels             = 4,
 dx                                  = 270,
 dy                                  = 270,
 grid_id                             = 1,
 parent_id                           = 0,
 i_parent_start                      = 132,
 j_parent_start                      = 107,
 parent_grid_ratio                   = 1,
 parent_time_step_ratio              = 1,
 feedback                            = 0,
 smooth_cg_topo                      = .true.,
 smooth_option                       = 0,
 numtiles                            = 2,
 /

 &physics
 mp_physics                          = 0,
 ra_lw_physics                       = 4,
 ra_sw_physics                       = 4,
 radt                                = 9,
 sf_sfclay_physics                   = 1,
 sf_surface_physics                  = 2,
 bl_pbl_physics                      = 0,
 bldt                                = 0,
 cu_physics                          = 0,
 cudt                                = 0,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_land_cat                        = 40,
 num_soil_layers                     = 4,
 sf_urban_physics                    = 0,     
 slope_rad                           = 1,
 topo_shading                        = 1,
 shadlen                             = 25000,
 /

 &fdda
 /

 &dynamics
 w_damping                           = 0,
 diff_opt                            = 2,
 km_opt                              = 2,
 diff_6th_opt                        = 2,
 diff_6th_factor                     = 0.12,
 base_temp                           = 290.
 damp_opt                            = 0,
 zdamp                               = 5000.,
 dampcoef                            = 0.2,
 khdif                               = 0,
 kvdif                               = 0,
 non_hydrostatic                     = .true.,
 moist_adv_opt                       = 1,   
 scalar_adv_opt                      = 1,     
 sfs_opt                             = 2,
 mix_isotropic                       = 1,
 mix_upper_bound                     = 0.1,
 /

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /

Wyrmux
 
Posts: 1
Joined: Tue Dec 03, 2013 10:34 am

Return to Nesting

Who is online

Users browsing this forum: No registered users and 3 guests