WRF Derecho case. (Experiment 4 at end)

Size: px

Start display at page:

Download "WRF Derecho case. (Experiment 4 at end)"

Melvin Page
5 years ago
Views:

1 WRF Derecho case (Experiment 4 at end) 1

2 Set up WRF environment mkdir DERECHO! cd DERECHO! cp /home/c115-test/derecho/make_all_links.csh.! cp /home/c115-test/derecho/namelist.*.! cp /home/c115-test/derecho/control_file.*.! cp /home/c115-test/derecho/plot_derecho_horiz.gs.! cp /home/c115-test/derecho/vtable.!!!!!!./make_all_links.csh!! 2

3 Outline Set up domain: geogrid.exe! Unpack inibalizabon data: ungrib.exe! Prepare inibalizabon data: metgrid.exe! IniBalize WRF model: real.exe! Run WRF model: wrf.exe! 3

4 Two namelist files namelist.wps Used for geogrid, ungrib, metgrid namelist.input Used for real, wrf 4

5 Tour of namelist.wps Example for a 2- domain run, inibalized with North American Regional Reanalysis (NARR) data 5

6 &share secbon &share! wrf_core = 'ARW',! max_dom = 2,! start_date = ' _00:00:00', ' _00:00:00,! end_date = ' _12:00:00', ' _12:00:00,! interval_seconds = 10800! io_form_geogrid = 2,! /! Don t worry about greyed text 2 domains will be defined Start date 00Z 28 June 2012 End date 12Z 30 June 2012 IniBalizaBon data interval 3 h (10800 sec) 6

7 &geogrid secbon &geogrid! parent_id = 1, 1,! parent_grid_ratio = 1, 3,! i_parent_start = 1, 11,! j_parent_start = 1, 10,! e_we = 48, 76,! e_sn = 35, 43,! geog_data_res dx = 54000,! dy = 54000,! map_proj = 'lambert',! ref_lat = 41.,! ref_lon = -86.0,! truelat1 = 30.0,! truelat2 = 60.0,! stand_lon = -86.0,! = 'modis_30s+10m','modis_30s+2m',! geog_data_path = '/home/c115-test/geogv36'! /! Parent domain 48x35 at 54 km resolubon Nest is 76x43 At grid_rabo = 3, that is 18 km spacing Nest lower le] corner at (11,10) relabve to parent Lambert projecbon, with central labtude 41 N, central longitude 86 W If ref_lon stand_lon, domain is rotated relabve to north Topo databases (do not change) 7

8 Other secbons &ungrib! out_format = 'WPS',! prefix = 'FILE',! /!! &metgrid! constants_name = '/home/fovell/fixed_narr',! fg_name = 'FILE'! io_form_metgrid = 2,! /! constants_name needed only for NARR inibalizabon data The rest rarely needs to be changed 8

9 Prepare domain: geogrid namelist.wps already edited 9

10 geogrid steps Run geogrid.exe! Creates files geo_em.d01.nc, geo_em.d02.nc Run plotgrids.exe! Creates simple domain picture gmeta Run idt gmeta! Displays domain 10

11 Telescoping domain D01 D02 (11,10) 11

12 Extract data: ungrib 12

13 ungrib step IniBalize with NARR data, stored at /home/c115- test/derecho/narrdata link_grib.csh /home/c115-test/derecho/ NARRDATA/merged*! Result: a large number of files called GRIBFILE* ungrib.exe! Result: a large number of files called FILE* 13

14 Prepare data: metgrid 14

15 metgrid step metgrid.exe! Result: a sequence of files called met_em.d01* and met_em.d02* 15

16 Tour of namelist.input Used for real.exe and wrf.exe! 16

17 &Bme_control secbon (edited) &time_control! run_days = 0,! run_hours run_seconds = 24,! = 0,! run_minutes start_year = 0,! = 2012, 2012,! We will run 24 hours Star2ng 12Z 29 June start_month = 06, 06,! start_day = 29, 29,! start_hour = 12, 12,! Ending 12Z 30 June end_year = 2012, 2012,! end_month = 06, 06,! end_day = 30, 30,! end_hour = 12, 12,! interval_seconds = 10800! From inibalizabon s interval input_from_file =.true.,.true.,! history_interval = 60, 60,! frames_per_outfile = 1, 1,! restart =.false.,! restart_interval = 99999,! /! Save history data every 60 min for both domains Each history Bme wriken to separate file Not a restart run 17

18 &domains secbon (edited) &domains! time_step = 120,! max_dom = 2,! e_we = 48, 76,! e_sn = 35, 43,! e_vert = 51, 51,! num_metgrid_levels = 30! num_metgrid_soil_levels = 4! dx = 54000, 18000,! dy = 54000, 18000,! grid_id = 1, 2,! parent_id = 0, 1,! i_parent_start = 1, 11,! j_parent_start = 1, 10,! parent_grid_ratio = 1, 3,! parent_time_step_ratio = 1, 3,! feedback = 1,! p_top_requested = 10000! /! 120 second Bme step (for outer domain; inner domain Bme step is 1/3 of this) Domain dimensions must agree with namelist.wps!!! 51 verbcal levels will be selected by real.exe (you can override) 30 metgrid levels and 4 soil levels (determined by inibalizabon data source) Grid spacing must agree, too!!! Nest posibon also must agree!!! Parent to nest grid spacing and Bme step rabos both 3 Feedback = 1 means nest influences parent Top of model pressure Pa = 100 mb (dictated by data source) 18

19 &physics secbon (edited) &physics! NUM_LAND_CAT = 20,! no_mp_heating = 0,! mp_physics = 3, 3,! ra_lw_physics = 4, 4,! ra_sw_physics = 4, 4,! radt = 10, 10,! 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,! num_soil_layers = 5,! sf_urban_physics = 0,! mp_zero_out = 0,! /! MODIS landuse database used WSM3 microphysics (mp = 3) RRTMG LW and SW radiabon (opbon 4) RadiaBon called every 10 min (radt) Surface layer (sfclay) and PBL scheme usually go together. This uses YSU PBL Noah land surface model (sf_surface = 2) Call PBL scheme every Bme step (bldt = 0) Kain- Fritsch cumulus (cu_physics = 1) used in both domains, called every 5 min (cudt) 19

20 &dynamics secbon &dynamics! w_damping = 0,! diff_opt = 1, 1,! km_opt = 4, 4,! diff_6th_opt = 1,! diff_6th_factor = 0.12,! base_temp = 290.! damp_opt = 0,! zdamp = 5000., 5000.,! dampcoef = 0.01, 0.01,! khdif = 0, 0,! kvdif = 0, 0,! non_hydrostatic =.true.,.true.,! moist_adv_opt = 2, 2,! scalar_adv_opt = 2, 2,! /! For real data runs, this secbon will rarely be modified. Same for &bdy_control (not shown) 20

21 IniBalize and run WRF namelist.input already edited 21

22 Run real.exe and wrf.exe! Run real.exe! Creates files wrrdy_d01, wrfinput_d01, wrfinput_d02 Run wrf.exe! ulimit s unlimited! nohup wrf.exe > wrf.out &! To monitor progress of model run tail -f wrf.out! Creates wrfout_d01*, wrfout_d02* files 22

23 Postprocessing WRF output wrf2grads 23

24 wrf2grads Creates the outputs examined in the sea- breeze experiment Edit control_file.z file to include list of wrfout* files (for only one domain) WRF verbcal coordinate is terrain- following This configurabon interpolates WRF output to a nice set of constant height levels, based on lowest height in domain being analyzed w2g control_file.z exper_name! Creates exper_name.ctl, exper_name.dat 24

25 29-30 June 2012 derecho 25

26 Some useful websites Storm PredicBon Center storm reports (can access past storm data) hkp:// Daily weather maps since 2003 hkp:// index.html Episodes page at NCAR (displays past radar mosaic and satellite imagery) hkp://locust.mmm.ucar.edu 26

27 hkp:// Storm PredicBon Center storm reports 27

28 Almost 1200 total wind reports 28

29 Highest wind gust report SW FORT WAYNE ALLEN IN MEASURED 91 MPH GUST AT FORT WAYNE INTERNATIONAL AIRPORT (IWX) 91 mph (41 m/s) wind gust 1854Z at 5 miles SW of Fort Wayne, IN (labtude N, longitude W) 1195 total reports, 37 exceeding 65 kt (75 mph, 33.4 m/s) 29

30 hkp:// 29 June

31 Radar composite 1624Z 29 June 2012 hkp://locust.mmm.ucar.edu 31

32 Radar composite 1923Z 29 June 2012 hkp://locust.mmm.ucar.edu 32

33 Radar composite 2154Z 29 June 2012 hkp://locust.mmm.ucar.edu 33

34 Control WRF simulabon 34

35 Control simulabon setup IniBalized 12Z 29 June 2012, run 24 h IniBalized with NARR reanalysis data 2 domains: 54 and 18 km resolubon WSM3 microphysics, Kain- Fritsch cumulus, Noah land surface model, YSU PBL 35

36 Using plot_derecho_horiz.gs for t = 14 H 36

37 Common WRF opbons for low resolubon (54 & 18 km) simulabons Microphysics (mp_physics) mp_physics = 3 WSM3 mp_physics = 5 WSM5 Cumulus schemes (cu_physics) cu_physics = 1 (Kain- Fritsch) cu_physics = 6 (Tiedtke) PBL schemes (bl_pbl_physics) bl_pbl_physics = 1 (YSU, with sf_sfclay_physics = 1) bl_pbl_physics = 2 (MYJ, with sf_sfclay_physics = 2) bl_pbl_physics = 5 (MYNN2, with sf_sfclay_physics = 1 or 2) bl_pbl_physics = 7 (ACM2, with sf_sfclay_physics = 1 or 7) LSMs (sf_surface_physics) sf_surface_physics = 1 (TD scheme) sf_surface_physics = 2 (Noah scheme) sf_surface_physics = 3 (RUC scheme) sf_surface_physics = 7 (PX scheme) 37

38 Experiment 4 38

39 QuesBon #1 Make the control run simulabon and verify you can create an image using plot_derecho_horiz.gs at t = 14 like that shown earlier QualitaBvely compare the storm to radar composite imagery from the Episodes web page. Comment on the accuracy of the Bme at which the simulated storm arrives in NE Indiana. Too early, too late, or on Bme? If late or early, by how many hours? 39

40 QuesBon #2 NoBce the high sea- level pressure ( slvl ) area of approximately 1014 mb at this Bme. WRF does not permit us to easily separate buoyancy from dynamic pressure. However, do you think buoyancy pressure is contribubng to this H pressure? Do you think dynamic pressure is contribubng? Why or why not? JusBfy your answer. Note you do not have to compute or esbmate the magnitude of buoyancy or dynamic pressure. Only the sign. 40

41 QuesBon #3 You can start the model simulabon at earlier Bmes, as early as 00Z on 28 June Does giving the model more Bme to spin up improve the forecast? Provide evidence for your conclusion. Note #1: Whatever the start Bme, run the model unbl 12Z on June 30. This means you have to alter start_day and start_hour, as well as run_hours! Note #2: If you do not edit control_file.z, then (say) t=14 will be the same Bme in each GrADS output file, independent of the actual model start Bme 41

WRF Pre-Processing (WPS)

WRF Pre-Processing (WPS) NCAR Earth System Laboratory National Center for Atmospheric Research NCAR is Sponsored by NSF and this work is partially supported by the Willis Research Network and the Research Partnership to Secure