[Wrf-users] Wrf-users Digest, Vol 139, Issue 35

Thorsten Beisiegel tbeisiegel at unal.edu.co
Thu Mar 31 07:28:43 MDT 2016


Dear Yaoping,

as I see from your rsl. file, instability occurs from vertical velocity. I
suggest: Put the upper level damping flag damp_opt=1 and dampcoef to
something between 0.01 and 0.1. For me it works fine with 0.01 even in
steep mountains. If this still doesn't work you can try additionally to put
w_damping=1.

I hope this will help to get your run stable.

with best regards,

Thorsten

2016-03-30 21:12 GMT-05:00 <wrf-users-request at ucar.edu>:

> Send Wrf-users mailing list submissions to
>         wrf-users at ucar.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         http://mailman.ucar.edu/mailman/listinfo/wrf-users
> or, via email, send a message with subject or body 'help' to
>         wrf-users-request at ucar.edu
>
> You can reach the person managing the list at
>         wrf-users-owner at ucar.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Wrf-users digest..."
>
>
> Today's Topics:
>
>    1. vertical velocity as output variable in ARWpost
>       (Jose Augusto Paix?o Veiga)
>    2. Exceeded cfl=2 at the corners of the domain (Wang, Yaoping)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 30 Mar 2016 15:49:32 -0400
> From: Jose Augusto Paix?o Veiga <veiga.uea at gmail.com>
> Subject: [Wrf-users] vertical velocity as output variable in ARWpost
> To: "wrf-users at ucar.edu" <wrf-users at ucar.edu>
> Message-ID:
>         <CAK=
> jcYOTLUJyLSVJtqGPuboyT4ByGiwkEHYhf9dBkwdSa_4how at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Dear users,
>
> someone knows how to include omega vertical velocity (dp/dt) as an output
> of ARWpost?
>
> I found the following list of available output diagnostic variables in WRF.
> But I see no vertical velocity included in.
>
>
> *cape* - 3d cape
> *cin* - 3d cin
> *mcape* - maximum cape
> *mcin* - maximum cin
>
> *clfr - *low/middle and high cloud fraction
> *dbz* - 3d reflectivity
> *max_dbz* - maximum reflectivity
>
> *geopt* - geopotential
> *height* - model height in km
>
> *lcl* - lifting condensation level
> *lfc* - level of free convection
> *pressure* - full model pressure in hPa
> *rh* - relative humidity
> *rh2* - 2m relative humidity
> *theta* - potential temperature
> *tc* - temperature in degrees C
> *tk* - temperature in degrees K
> *td* - dew point temperature in degrees C
> *td2* - 2m dew point temperature in degrees C
>
> *slp* - sea level pressure
>
> *umet* and *vmet* - winds rotated to earth coordinates
> *u10m* and *v10m* - 10m winds rotated to earth coordinates
> *wdir* - wind direction
> *wspd* - wind speed coordinates
> *wd10* - 10m wind direction
> *ws10* - 10m wind speed
>
> Any help on this will be very welcome .
>
> Thanks in advance?
>
>
>
> Jos? Augusto P. Veiga
> CREA n.: 24161
>
> ======================================
> Universidade do Estado do Amazonas
> Departamento de Meteorologia
> Escola Superior de Tecnologia (EST)
> ======================================
> Av. Darcy Vargas, 1200, Manaus-AM Brasil
> Work phone:    (92) 3878 4317, Ramal 4317
> Skype: veiga_j.a.p.
> ======================================
> CV: http://lattes.cnpq.br/4027612512091565
> URL: https://scientificmet.wordpress.com
> ======================================
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/wrf-users/attachments/20160330/cd04e447/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Thu, 31 Mar 2016 01:04:52 +0000
> From: "Wang, Yaoping" <wang.3866 at buckeyemail.osu.edu>
> Subject: [Wrf-users] Exceeded cfl=2 at the corners of the domain
> To: "wrfhelp at ucar.edu" <wrfhelp at ucar.edu>
> Cc: "wrf-users at ucar.edu" <wrf-users at ucar.edu>
> Message-ID:
>         <
> BY2PR01MB057625398F4EA0C39A69A13B2990 at BY2PR01MB057.prod.exchangelabs.com>
>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hello,
>
>
> I am new to WRF and am setting up a run over the North American domain at
> 30km resolution. But I keep getting "exceeded cfl=2". After testing various
> time steps it seems that the error keep occurring at the corners of my
> domain. I also checked the wind speeds and divergence in the input
> meteorological data, but did not find any clear outlier models. Could
> anyone give some suggestions as to what can be done to solve the problem?
>
>
> The trials and errors:
>
> time step = 180s
>
> d01 1967-01-01_00:00:00         5537 points exceeded cfl=2 in domain d01
> at time 1967-01-01_00:00:00 hours
> d01 1967-01-01_00:00:00  MAX AT i,j,k:          283          42
>  2 vert_cfl,w,d(eta)=  2.4396536E+07            NaN  9.9999905E-03
> d01 1967-01-01_00:00:00         1906 points exceeded cfl=2 in domain d01
> at time 1967-01-01_00:00:00 hours
> d01 1967-01-01_00:00:00  MAX AT i,j,k:           21         190
> 20 vert_cfl,w,d(eta)=   1102.228      -6444.851      3.5938144E-02
> Timing for main: time 1967-01-01_00:03:00 on domain   1:   15.48090
> elapsed seconds
>
>
> time step = 90s
>
> d01 1967-01-01_00:00:00         6080 points exceeded cfl=2 in domain d01
> at time 1967-01-01_00:00:00 hours
> d01 1967-01-01_00:00:00  MAX AT i,j,k:          283          46
> 29 vert_cfl,w,d(eta)=   6955.001     -1.8047905E+15  9.2214663E-03
> d01 1967-01-01_00:00:00         1969 points exceeded cfl=2 in domain d01
> at time 1967-01-01_00:00:00 hours
> d01 1967-01-01_00:00:00  MAX AT i,j,k:          278          56
> 21 vert_cfl,w,d(eta)=   2879.670       56544.66      3.2001421E-02
> Timing for main: time 1967-01-01_00:01:30 on domain   1:   13.09157
> elapsed seconds
>
>
> time step = 30s, namelist.input &domains  smooth_cg_topo = .true.
>
> d01 1967-01-01_00:00:00         5537 points exceeded cfl=2 in domain d01
> at time 1967-01-01_00:00:00 hours
> d01 1967-01-01_00:00:00  MAX AT i,j,k:          283          42
>  2 vert_cfl,w,d(eta)=  2.4396536E+07            NaN  9.9999905E-03
> d01 1967-01-01_00:00:00         1906 points exceeded cfl=2 in domain d01
> at time 1967-01-01_00:00:00 hours
> d01 1967-01-01_00:00:00  MAX AT i,j,k:           21         190
> 20 vert_cfl,w,d(eta)=   1102.228      -6444.851      3.5938144E-02
>
>
> time step = 10s, namelist.input &domains  smooth_cg_topo = .true.
>
> d01 1967-01-01_00:00:00         2548 points exceeded cfl=2 in domain d01
> at time 1967-01-01_00:00:00 hours
> d01 1967-01-01_00:00:00  MAX AT i,j,k:          285         179
> 23 vert_cfl,w,d(eta)=  8.4951393E+14  2.3136813E+23  2.4293438E-02
> d01 1967-01-01_00:00:00          987 points exceeded cfl=2 in domain d01
> at time 1967-01-01_00:00:00 hours
> d01 1967-01-01_00:00:00  MAX AT i,j,k:          285         196
> 14 vert_cfl,w,d(eta)=   24310.53       372558.9      6.8435431E-02
> Timing for main: time 1967-01-01_00:00:10 on domain   1:   15.84628
> elapsed seconds
>
>
> The namelist.input file:
>
>  &time_control
>  run_days                            = 3,
>  run_hours                           = 0,
>  run_minutes                         = 0,
>  run_seconds                         = 0,
>  start_year                          = 1967,
>  start_month                         = 01,
>  start_day                           = 01,
>  start_hour                          = 00,
>  start_minute                        = 00,
>  start_second                        = 00,
>  end_year                            = 1967,
>  end_month                           = 01,
>  end_day                             = 04,
>  end_hour                            = 00,
>  end_minute                          = 00,
>  end_second                          = 00,
>  interval_seconds                    = 21600
>  input_from_file                     = .true.
>  history_interval                    = 60,
>  frames_per_outfile                  = 1000,
>  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                           = 10,
>  time_step_fract_num                 = 0,
>  time_step_fract_den                 = 1,
>  max_dom                             = 1,
>  s_we                                = 1,
>  s_sn                                = 1,
>  s_vert                              = 1,
>  e_we                                = 290,
>  e_sn                                = 200,
>  e_vert                              = 30,
>  p_top_requested                     = 5000,
>  num_metgrid_levels                  = 27,
>  num_metgrid_soil_levels             = 4,
>  dx                                  = 30000,
>  dy                                  = 30000,
>  grid_id                             = 1,
>  parent_id                           = 1,
>  i_parent_start                      = 1,
>  j_parent_start                      = 1,
>  smooth_option                       = 0
>  smooth_cg_topo = .true.
>  /
>
>  &physics
>  mp_physics                          = 3,
>  ra_lw_physics                       = 1,
>  ra_sw_physics                       = 1,
>  radt                                = 30,
>  sf_sfclay_physics                   = 1,
>  sf_surface_physics                  = 2,
>  bl_pbl_physics                      = 1,
>  bldt                                = 0,
>  cu_physics                          = 1,
>  cudt                                = 5,
>  isfflx                              = 1,
>  ifsnow                              = 1,
>  icloud                              = 1,
>  surface_input_source                = 1,
>  num_soil_layers                     = 4,
>  sf_urban_physics                    = 0,
>  /
>
>  &fdda
>  /
>
>  &dynamics
>  w_damping                           = 0,
>  diff_opt                            = 1,
>  km_opt                              = 4,
>  diff_6th_opt                        = 0,
>  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,
>  /
>
>  &bdy_control
>  spec_bdy_width                      = 5,
>  spec_zone                           = 1,
>  relax_zone                          = 5,
>  specified                           = .true.,
>  nested                              = .false.,
>  /
>
>  &grib2
>  /
>
>  &namelist_quilt
>  nio_tasks_per_group = 0,
>  nio_groups = 1,
>  /
>
>
> Thank you very much,
>
> Yaoping Wang
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/wrf-users/attachments/20160331/e7303f9c/attachment.html
>
> ------------------------------
>
> _______________________________________________
> Wrf-users mailing list
> Wrf-users at ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/wrf-users
>
>
> End of Wrf-users Digest, Vol 139, Issue 35
> ******************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20160331/52df642d/attachment-0001.html 


More information about the Wrf-users mailing list