No rain in output

Topics related to current and future physics in the WRF as well as any problems you may have.

Re: No rain in output

Postby olver.hernandez » Mon Nov 26, 2012 7:26 am

Hi,
I have the same problem as no rain for Colombia in a case where was raining over a lot of areas over the domain.
I'm just a beginner with WRF and I'm trying to run it with the default configuration for three nested domains, using the GFS outputs as boundary and initialization (GFS time: 2012/11/15 00:00).

This is my namelist.input
Code: Select all
 &time_control
 run_days                            = 1,
 run_hours                           = 0,   
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2012,     2012,   2012,
 start_month                         = 11,       11,       11,
 start_day                           = 15,       15,       15,
 start_hour                          = 12,       12,       12,
 start_minute                        = 00,       00,       00,
 start_second                        = 00,       00,       00,
 tstart                              = 00,       
 end_year                            = 2012,     2012,     2012,
 end_month                           = 11,       11,        11,
 end_day                             = 16,       16,        16,
 end_hour                            = 12,       12,       12,
 end_minute                          = 00,       00,       00,
 end_second                          = 00,       00,       00,
 interval_seconds                    = 10800,
 history_interval                    = 180,       60,      15,
 frames_per_outfile                  = 1,        1,       1,
 restart                             = .false.,
 restart_interval                    = 5400,
 reset_simulation_start              = F,
 io_form_input                       = 2,
 io_form_history                     = 2,
 io_form_restart                     = 2,
 io_form_boundary                    = 2,
 io_form_auxinput1                   = 2,
 debug_level                         = 1,
 /

 &domains
 time_step                           = 60,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 3,
 e_we                                = 95,       41,     68,
 e_sn                                = 190,      118,    136,
 e_vert                              = 35,       35,     35,
 num_metgrid_levels                  = 27,
 dx                                  = 0.16190330633, 0.053967769,   0.053967769,   
 dy                                  = 0.16190330633, 0.053967769,   0.053967769,   
 p_top_requested                     = 5000.,
 ptsgm                               = 42000.,
 grid_id                             = 1,        2,       3,
 parent_id                           = 0,        1,       2,
 i_parent_start                      =  1,      40,      8,
 j_parent_start                      =  1,      78,     55,
 parent_grid_ratio                   = 1,        3,      3,
 parent_time_step_ratio              = 1,        3,      3,
 eta_levels                          =          1.0000, 0.9974, 0.9940, 0.9900, 0.9854, 0.9796, 0.9723, 0.9635,
                                                0.9528, 0.9401, 0.9252, 0.9079, 0.8882, 0.8659, 0.8410, 0.8133,
                                                0.7828, 0.7494, 0.7133, 0.6742, 0.6323, 0.5878, 0.5406, 0.4915,
                                                0.4409, 0.3895, 0.3379, 0.2871, 0.2378, 0.1907, 0.1465, 0.1056,
                                                0.0682, 0.0332, 0.0000,

/

 &physics
 mp_physics                          = 5,        5,      5,
 ra_lw_physics                       = 99,       99,     99,
 ra_sw_physics                       = 99,       99,     99,
 nrads                               = 20,       60,     180,
 nradl                               = 20,       60,     180,
 co2tf                               = 1,
 sf_sfclay_physics                   = 2,        2,      2,
 sf_surface_physics                  = 2,        2,      2,
 bl_pbl_physics                      = 2,        2,      2,
 nphs                                = 2,        4,      6,
 cu_physics                          = 2,        2,      2,
 ncnvc                               = 2,        4,      6,
 tprec                               = 3,        3,      3,
 theat                               = 6,        6,      6,
 tclod                               = 6,        6,      6,
 trdsw                               = 6,        6,      6,
 trdlw                               = 6,        6,      6,
 tsrfc                               = 6,        6,      6,
 pcpflg                              = .false.,  .false.,    .false.,
 num_soil_layers                     = 4,
 mp_zero_out                         = 0
 gwd_opt                             = 0
 /

 &dynamics
 coac                                = 1.6,
 codamp                              = 6.4,
 slophc                              = 0.0064,
 euler_adv                           = .true.,
 idtadt                              = 2,
 idtadc                              = 1
 /

 &bdy_control
 spec_bdy_width                      = 1,
 specified                           = .true.,
 nested                              = .false.
 /

 &fdda
 /

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1
 /


Tha output files donesn't contains any precipitacion output. You can see that through the next command:
Code: Select all
ncdump -h wrfout_d01_2012-11-15_21\:00\:00 |grep RAIN

is:
Code: Select all
   float F_RAIN(Time, bottom_top, south_north, west_east) ;
      F_RAIN:FieldType = 104 ;
      F_RAIN:MemoryOrder = "XYZ" ;
      F_RAIN:description = "Rain fraction of liquid part of CWM" ;
      F_RAIN:units = "" ;
      F_RAIN:stagger = "" ;
   float AVRAIN(Time) ;
      AVRAIN:FieldType = 104 ;
      AVRAIN:MemoryOrder = "0  " ;
      AVRAIN:description = "  of times gridscale precip called in NHEAT steps" ;
      AVRAIN:units = "" ;
      AVRAIN:stagger = "" ;
   float QRAIN(Time, bottom_top, south_north, west_east) ;
      QRAIN:FieldType = 104 ;
      QRAIN:MemoryOrder = "XYZ" ;
      QRAIN:description = "Rain water mixing ratio" ;
      QRAIN:units = "kg kg-1" ;
      QRAIN:stagger = "" ;
      CUPPT:description = "ACCUMULATED CONVECTIVE RAIN SINCE LAST CALL TO THE RADIATION" ;

Then not precipitacion outputs :(

Any suggestions?
olver.hernandez
 
Posts: 2
Joined: Mon Nov 26, 2012 6:54 am

Re: No rain in output

Postby alainaketh » Wed Oct 05, 2016 12:38 pm

Hi everyone,

I just solved this problem. I am using ECMWF multimodel seasonal forecast data for a simulation of the Katrina hurricane in 2005.

I'll skip the details of this dataset and go straight to the solution: my input data didn't have the relative humidity fields required by the Vtable, but did have specific humidity that is not included in the Vtable. Just like jimmyc pointed out, this is vital information to produce rainfall since (unless you're working with the WRF3D or 4D-VAR) you're not actually entering rainfall information into the model.

So I modified my Vtable (which by the way is Vtable.ECMWF) from this:

Code: Select all
GRIB | Level| Level| Level| metgrid  |  metgrid | metgrid                                  |
Code | Code |   1  |   2  | Name     |  Units   | Description                              |
-----+------+------+------+----------+----------+------------------------------------------+
 129 | 100  |   *  |      | GEOPT    | m2 s-2   |                                          |
     | 100  |   *  |      | HGT      | m        | Height                                   |
 130 | 100  |   *  |      | TT       | K        | Temperature                              |
 131 | 100  |   *  |      | UU       | m s-1    | U                                        |
 132 | 100  |   *  |      | VV       | m s-1    | V                                        |
 157 | 100  |   *  |      | RH       | %        | Relative Humidity                        |
 165 |  1   |   0  |      | UU       | m s-1    | U                                        | At 10 m
 166 |  1   |   0  |      | VV       | m s-1    | V                                        | At 10 m
 167 |  1   |   0  |      | TT       | K        | Temperature                              | At 2 m
 168 |  1   |   0  |      | DEWPT    | K        |                                          | At 2 m
     |  1   |   0  |      | RH       | %        | Relative Humidity                        | At 2 m
 172 |  1   |   0  |      | LANDSEA  | 0/1 Flag | Land/Sea flag                            |
 129 |  1   |   0  |      | SOILGEO  | m2 s-2   |                                          |
     |  1   |   0  |      | SOILHGT  | m        | Terrain field of source analysis         |
 134 |  1   |   0  |      | PSFC     | Pa       | Surface Pressure                         |
 151 |  1   |   0  |      | PMSL     | Pa       | Sea-level Pressure                       |
 235 |  1   |   0  |      | SKINTEMP | K        | Sea-Surface Temperature                  |
  31 |  1   |   0  |      | SEAICE   | 0/1 Flag | Sea-Ice-Flag                             |
  34 |  1   |   0  |      | SST      | K        | Sea-Surface Temperature                  |
 141 |  1   |   0  |      | SNOW_EC  | m        |                                          |
     |  1   |   0  |      | SNOW     | kg m-2   |Water Equivalent of Accumulated Snow Depth|
 139 | 112  |   0  |   7  | ST000007 | K        | T of 0-7 cm ground layer                 |
 170 | 112  |   7  |  28  | ST007028 | K        | T of 7-28 cm ground layer                |
 183 | 112  |  28  | 100  | ST028100 | K        | T of 28-100 cm ground layer              |
 236 | 112  | 100  | 255  | ST100255 | K        | T of 100-255 cm ground layer             |
  39 | 112  |   0  |   7  | SM000007 | fraction | Soil moisture of 0-7 cm ground layer     |
  40 | 112  |   7  |  28  | SM007028 | fraction | Soil moisture of 7-28 cm ground layer    |
  41 | 112  |  28  | 100  | SM028100 | fraction | Soil moisture of 28-100 cm ground layer  |
  42 | 112  | 100  | 255  | SM100255 | fraction | Soil moisture of 100-255 cm ground layer |
-----+------+------+------+----------+----------+------------------------------------------+
#
#  Grib codes are from Table 128
#  http://www.ecmwf.int/products/data/technical/GRIB_tables/table_128.html

# snow depth is converted to the proper units in rrpr.F


and simply added the line

Code: Select all
 133 | 100  |   *  |      | SPECHUMD | kg kg-1  | Specific Humidity            |


between the CODE 132 and 157 lines. I took that additional line from this website: https://wrfandme.wordpress.com/2016/02/05/converting-ecmwf-to-wrf-input/

To check that the model is reading humidity, check your log.ungrib file after running ungrib. After adding that line, my log file now includes the SPECHUMD column and looks like this:

Code: Select all
 *** Starting program ungrib.exe ***
Start_date =  2015-08-01_00:00:00 ,      End_date = 2015-08-06_00:00:00
output format is WPS
Path to intermediate files is ./

###############################################################################

Inventory for date = 2015-08-01 00:00:00

PRES   GEOPT    HGT      TT       UU       VV       SPECHUMD RH       DEWPT    LANDSEA  SOILGEO  SOILHGT  PSFC     PMSL     SKINTEMP SEAICE   SST      SNOW_DEN SNOW_EC  SNOW     SNOWH    ST000007 ST007028 ST028100 ST100255 SM000007 SM007028 SM028100 SM100255
-------------------------------------------------------------------------------
1000.0  X                 X        X        X        X                                                                                                                                                                                                           
 925.0  X                 X        X        X        X                                                                                                                                                                                                           
 850.0  X                 X        X        X        X                                                                                                                                                                                                           
 700.0  X                 X        X        X        X                                                                                                                                                                                                           
 500.0  X                 X        X        X        X                                                                                                                                                                                                           
 400.0  X                 X        X        X        X                                                                                                                                                                                                           
 300.0  X                 X        X        X        X                                                                                                                                                                                                           
 200.0  X                 X        X        X        X                                                                                                                                                                                                           
 100.0  X                 X        X        X        X                                                                                                                                                                                                           
  70.0  X                 X        X        X        X                                                                                                                                                                                                           
  50.0  X                 X        X        X        X                                                                                                                                                                                                           
  30.0  X                 X        X        X        X                                                                                                                                                                                                           
  10.0  X                 X        X        X        X                                                                                                                                                                                                           
-------------------------------------------------------------------------------

###############################################################################

Inventory for date = 2015-08-02 00:00:00

PRES   GEOPT    HGT      TT       UU       VV       SPECHUMD RH       DEWPT    LANDSEA  SOILGEO  SOILHGT  PSFC     PMSL     SKINTEMP SEAICE   SST      SNOW_DEN SNOW_EC  SNOW     SNOWH    ST000007 ST007028 ST028100 ST100255 SM000007 SM007028 SM028100 SM100255
-------------------------------------------------------------------------------
1000.0  X                 X        X        X        X                                                                                                                                                                                                           
 925.0  X                 X        X        X        X                                                                                                                                                                                                           
 850.0  X                 X        X        X        X                                                                                                                                                                                                           
 700.0  X                 X        X        X        X                                                                                                                                                                                                           
 500.0  X                 X        X        X        X                                                                                                                                                                                                           
 400.0  X                 X        X        X        X                                                                                                                                                                                                           
 300.0  X                 X        X        X        X                                                                                                                                                                                                           
 200.0  X                 X        X        X        X                                                                                                                                                                                                           
 100.0  X                 X        X        X        X                                                                                                                                                                                                           
  70.0  X                 X        X        X        X                                                                                                                                                                                                           
  50.0  X                 X        X        X        X                                                                                                                                                                                                           
  30.0  X                 X        X        X        X                                                                                                                                                                                                           
  10.0  X                 X        X        X        X                                                                                                                                                                                                           
-------------------------------------------------------------------------------


For all my time steps.

So yeah, check that you have Geopotential height, U and V wind speed and relative and/or specific humidity.

Hope this helps,

Alaina.
alainaketh
 
Posts: 26
Joined: Fri Jun 26, 2015 12:30 pm

Previous

Return to Physics

Who is online

Users browsing this forum: No registered users and 2 guests

cron