Run Forecasting with GFS Initial Data

Introduction to the forum

Run Forecasting with GFS Initial Data

Postby SussanaDiaz » Fri Jul 27, 2018 12:45 pm

Hi, everyone!
I'm beginner in WRF real data case, I have been running some simulations with GFS data and everthing gone OK, but now I need to run a real 72 hr forecasting.

My GFS data was download from:
https://nomads.ncdc.noaa.gov/data/gfs4/201807/20180723/?C=M;O=A
and
http://rda.ucar.edu/data/ds083.2/

Initial date:2018-07-23_00:00 and end date=2018-07-26_00:00 UTC
so the namelist.wps is:

Code: Select all
&share
 wrf_core = 'ARW',
 max_dom = 2,
 start_date = '2018-07-23_00:00:00', '2018-07-23_00:00:00',
 end_date   = '2018-07-26_00:00:00', '2018-07-26_00:00:00',
 interval_seconds = 21600,
 io_form_geogrid = 2,
 opt_output_from_geogrid_path = '/home/wrf/WRF/Build_WRF/Domain_Wizard/FORECAST_TUTO/',
 debug_level = 0,
/

&geogrid
 parent_id         = 1,1,
 parent_grid_ratio = 1,3,
 i_parent_start    = 1,36,
 j_parent_start    = 1,23,
 e_we          = 100,91,
 e_sn          = 82,70,
 geog_data_res = '10m','10m',
 dx = 59700,
 dy = 59700,
 map_proj =  'lambert',
 ref_lat   = 24.584,
 ref_lon   = -100.541,
 truelat1  = 24.584,
 truelat2  = 24.584,
 stand_lon = -100.541,
 geog_data_path = '/home/wrf/WRF/Build_WRF/WPS_GEOG',
 opt_geogrid_tbl_path = '/home/wrf/WRF/Build_WRF/Domain_Wizard/FORECAST_TUTO/',
 ref_x = 50.0,
 ref_y = 41.0,
/

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

&metgrid
 fg_name = 'FILE',
 io_form_metgrid = 2,
 opt_output_from_metgrid_path = '/home/wrf/WRF/Build_WRF/Domain_Wizard/FORECAST_TUTO/',
 opt_metgrid_tbl_path = '/home/wrf/WRF/Build_WRF/Domain_Wizard/FORECAST_TUTO/',
/

&mod_levs
 press_pa = 201300 , 200100 , 100000 ,
             95000 ,  90000 ,
             85000 ,  80000 ,
             75000 ,  70000 ,
             65000 ,  60000 ,
             55000 ,  50000 ,
             45000 ,  40000 ,
             35000 ,  30000 ,
             25000 ,  20000 ,
             15000 ,  10000 ,
              5000 ,   1000
 /


&domain_wizard
 grib_data_path = '/home/wrf/WRF/Build_WRF/DATA/FIRST_FORECAST',
 grib_vtable = 'Vtable.GFS',
 dwiz_name    =FORECAST_TUTO
 dwiz_desc    =
 dwiz_user_rect_x1 =284
 dwiz_user_rect_y1 =248
 dwiz_user_rect_x2 =611
 dwiz_user_rect_y2 =489
 dwiz_show_political =true
 dwiz_center_over_gmt =true
 dwiz_latlon_space_in_deg =10
 dwiz_latlon_linecolor =-8355712
 dwiz_map_scale_pct =12.5
 dwiz_map_vert_scrollbar_pos =0
 dwiz_map_horiz_scrollbar_pos =0
 dwiz_gridpt_dist_km =59.7
 dwiz_mpi_command =null
 dwiz_tcvitals =null
 dwiz_bigmap =Y
/


And the namelist.intput:
Code: Select all
&time_control           
run_days                 = 0,
run_hours                = 84,
run_minutes              = 0,
run_seconds              = 0,
start_year               = 2018,     2018,
start_month              = 7,        7,
start_day                = 23,       23,
start_hour               = 12,       12,
start_minute             = 00,       00,
start_second             = 00,       00,
end_year                 = 2018,     2018,
end_month                = 7,        7,
end_day                  = 26,       26,
end_hour                 = 12,       12,
end_minute               = 00,       00,
end_second               = 00,       00,
interval_seconds         = 21600,
input_from_file          = .true.,   .true.,
history_interval         = 180,       60,
frames_per_outfile       = 1000,     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                 
eta_levels               = 1.000, 0.994, 0.987, 0.979, 0.97,
                0.96, 0.949, 0.937, 0.924, 0.909,
                0.892, 0.873, 0.851, 0.826, 0.798,
                0.768, 0.736, 0.702, 0.666, 0.629,
                0.5915, 0.5536, 0.5153, 0.4773, 0.44,
                0.404, 0.3695, 0.3375, 0.3085, 0.2845,
                0.2645, 0.2465, 0.2305, 0.2165, 0.2035,
                0.1915, 0.1792, 0.1667, 0.1539, 0.1407,
                0.1272, 0.1134, 0.0995, 0.0855, 0.0713,
                0.0571, 0.0429, 0.0287, 0.0145, 0.000,
time_step                = 358,
time_step_fract_num      = 0,
time_step_fract_den      = 1,
max_dom                  = 2,
e_we                     = 100,       91,
e_sn                     = 82,       70,
e_vert                   = 50,       50,
p_top_requested          = 6000,
num_metgrid_levels       = 50,
num_metgrid_soil_levels  = 4,
dx                       = 59700,    19900,
dy                       = 59700,    19900,
grid_id                  = 1,        2,
parent_id                = 1,        1,
i_parent_start           = 1,       36,
j_parent_start           = 1,       23,
parent_grid_ratio        = 1,        3,
parent_time_step_ratio   = 1,        3,
feedback                 = 1,
smooth_option            = 0,
/

&physics                 
mp_physics               = 3,        3,
ra_lw_physics            = 1,        1,
ra_sw_physics            = 1,        1,
radt                     = 30,       30,
sf_sfclay_physics        = 1,        1,
sf_surface_physics       = 2,        2,
bl_pbl_physics           = 1,        1,
bldt                     = 0,        0,
cu_physics               = 1,        1,
cudt                     = 5,        5,
isfflx                   = 1,
ifsnow                   = 0,
icloud                   = 1,
surface_input_source     = 1,
num_soil_layers          = 4,
sf_urban_physics         = 0,        0,
maxiens                  = 1,
maxens                   = 3,
maxens2                  = 3,
maxens3                  = 16,
ensdim                   = 144,
/

&fdda                   
/

&dynamics               
w_damping                = 0,
diff_opt                 = 1,
km_opt                   = 4,
diff_6th_opt             = 0,        0,
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,
/

&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,
/


The problem is when I run Ungrib.exe:

Code: Select all
2018-07-27 11:34:11.323 ---  *** StarStart_date =  2018-07-23_00:2018-07-27 11:34:11.324 --- Start_date =  20output format is WPS
Path to intermediate files is ./

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

Inventory for date = 2018-07-23 00:00:00

PRES   TT       UU       VV       RH       HGT      PSFC     PMSL     SM000010 SM010040 SM040100 SM100200 SM010200 ST000010 ST010040 ST040100 ST100200 ST010200 SEAICE   LANDSEA  LANDN    SOILHGT  SKINTEMP SNOW     SNOWH    UMAXW    VMAXW    PMAXW    PMAXWNN  TMAXW    HGTMAXW  UTROP    VTROP    PTROP    PTROPNN  TTROP    HGTTROP
-------------------------------------------------------------------------------
2013.0  O        O        O        O        O        O        X        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O     
2001.0  X        X        X        X        O        X        O        X        X        X        X        O        X        X        X        X        O        X        X        X        X        X        X        O        X        X        X        O        X        X        X        X        X        O        X        X     
1000.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 975.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 950.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 925.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 900.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 850.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 800.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 750.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 700.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 650.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 600.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 550.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 500.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 450.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 400.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 350.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 300.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 250.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 200.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 150.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                                                                                                                                                                                                                                                                                             
  20.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
  10.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
   7.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
   5.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
   3.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
   2.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
   1.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
-------------------------------------------------------------------------------

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

Inventory for date = 2018-07-23 00:00:00

PRES   TT       UU       VV       RH       HGT      PSFC     PMSL     SM000010 SM010040 SM040100 SM100200 SM010200 ST000010 ST010040 ST040100 ST100200 ST010200 SEAICE   LANDSEA  LANDN    SOILHGT  SKINTEMP SNOW     SNOWH    UMAXW    VMAXW    PMAXW    PMAXWNN  TMAXW    HGTMAXW  UTROP    VTROP    PTROP    PTROPNN  TTROP    HGTTROP
-------------------------------------------------------------------------------
2013.0  O        O        O        O        O        O        X        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O        O     
2001.0  X        X        X        X        O        X        O        X        X        X        X        O        X        X        X        X        O        X        X        X        X        X        X        X        X        X        X        X        X        X        X        X        X        X        X        X     
1000.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 975.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 950.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 925.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 900.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 850.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 800.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 750.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 700.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 650.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 600.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 550.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 500.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 450.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 400.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 350.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 300.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 250.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 200.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
 150.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                                                                                                                                                                                                                                                                                             
  20.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
  10.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
Subroutine DATINT: Interpolating 3-d files to fill in any missing data...
Looking for data at time 2018-07-23_00
Found file:      FILE:2018-07-23_00
Looking for data at time 2018-07-23_06
ERROR: Data not found: 2018-07-23_06:00:00.0000
 Name of source model =>NCEP GFS Analysis       GRID   4
   7.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
   5.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
   3.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
   2.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
   1.0  X        X        X        X        X                                                                                                                                                                                                                                                                                             
-------------------------------------------------------------------------------
   -         -         -         -
2018-07-27 11:34:14.253 ---     15000         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.254 ---     10000         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.254 ---      7000         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.255 ---      5000         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.255 ---      3000         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.255 ---      2000         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.256 ---      1000         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.256 ---       700         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.257 ---       500         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.257 ---       300         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.257 ---       200         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.258 ---       100         X         X         X         X         X         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -         -
2018-07-27 11:34:14.258 --- -------------------------------------------------
2018-07-27 11:34:14.628 --- Subroutine DATINT: Interpolating 3-d files to fill in any missing data...
2018-07-27 11:34:14.629 --- Looking for data at time 2018-07-23_00
2018-07-27 11:34:14.629 --- Found file:      FILE:2018-07-23_00
2018-07-27 11:34:14.629 --- Looking for data at time 2018-07-23_06
2018-07-27 11:34:14.629 --- ERROR: Data not found: 2018-07-23_06:00:00.0000/


I don't understand how run a forecasting just a simulation, please help me!
Regards,

Ps. Sorry for my poor english!
como dijo Celia "my english is not very good looking".
SussanaDiaz
 
Posts: 1
Joined: Fri Jul 27, 2018 12:05 pm

Re: Run Forecasting with GFS Initial Data

Postby dcvz » Mon Jul 30, 2018 8:21 pm

A forecast and simulation are essentially the same thing, except the forecast uses lateral boundary conditions from a forecast model (in your case a GFS forecast).

The first thing you need to check is to make sure you downloaded the correct files from nomads (the .grb2 files, not the .inv files). Then you need to make sure the files are named correctly (the 000-h file contains the 000-h forecast and the 006-h file contains the 006-h forecast). Then, make sure that GRIBFILE.AAA points to the 000-h file and GRIBFILE.AAB points to the 006-h file, etc.).

I suggest setting debug_level = 100 in your namelist.wps to give you more printout.
The first couple of lines in your ungrib printout looks messed up. It was probably due to cut and paste, but we need to see what ungrib thinks is the start and end time.

The only problems I see in your namelists, is the end_date for domain 2 should be the same as the start_date. There is no need for metgrid to process forecasts on domain 2. Also, your namelist.input has a start hour of 12 and not 0.
dcvz
 
Posts: 172
Joined: Tue Apr 15, 2008 12:02 am


Return to Welcome to the WRF User's Forum

Who is online

Users browsing this forum: No registered users and 2 guests