NMMWRFHWRF Code Walk Through S G Gopalakrishnan Environmental

NMM-WRF/HWRF: Code Walk Through S. G. Gopalakrishnan Environmental Modeling Center NCEP/NOAA/NWS, Washington, DC Adopted from John Michalakes, NCAR

, Matthew Pyle

The Modeling Framework Real cases: Standard Initialization (WRFSI/NMMSI) Wrfv 1. /Registry. /inc. /Main. /dyn_eh or. /dyn_em or. /dyn_slt . /phys. /frame. /share. /external Post Processor: (i) Netcdf (ii) Binary

Quick walk thru’ of the real programme • http: //www. mmm. ucar. edu/wrf/WG 2/software_v 2/

http: //box. mmm. ucar. edu/wrf/WG 2/Tigers/ Welcome To the Inaugural Meeting of the WRF Software Training and Documentation Team Jan. 26 -28, 2004 NCAR, MMM Division

Is this a Jig-Saw Puzzle ? ? [/emc 2/wx 20 gk/WRFEMC/HWRFV 1/WRFV 1] % ls CVS Makefile README Registry arch clean compile configure. wrf dyn_eh dyn_em dyn_exp dyn_graps dyn_nmm dyn_slt external frame inc install_nmm. scr main phys run set_core share test tools

Directory Structure

WRF Model Directory Structure driver mediation model

#<Table> <Type> <Sym> <Dims> <Use> <Num. TLev> <Stagger> #<IO> <DNAME> <DESCRIP> <UNITS> state real t ikj dyn_nmm 1 - i 01 rh "T" "TEMPERATURE IN KELVIN" state real pd ij dyn_nmm 1 - i 01 rh “PD" " PRESSURE DEPTH IN PASCALS" state real f ij dyn_nmm 1 - irh "F" "CORIOLIS PARAMETER " state real pint ikj dyn_nmm 1 Z irh "PINT " “INTERFACE HSTATIC PRES“ -

# Exchange in the NMM core halo HALO_NMM_C dyn_nmm 8: u, v halo HALO_NMM_INTERP_DOWN dyn_nmm 120: sm, fis, t, u, v, q, pd

Boundary Arrays in the Registry • • • • state state state state real pd_bt real t_bt real q_bt real u_bt real v_bt real q 2_bt real cwm_bt ivjb ikjb ikjb ikjb dyn_nmm dyn_nmm dyn_nmm dyn_nmm 1 1 1 1 - b b b b "PDB" "PDBT" "TBT" "QBT" "UBT" "VBT" "Q 2 T" "CWMT"

Boundary Conditions • The lateral boundaries for real cases are stored as, and read from special four dimensional arrays. These arrays are specified in the Registry as follows: the state variable, say for instance U, has its boundary values stored in U_b and the boundary tencied in U_bt. The “_b” and “_bt” arrays are four dimensional (max(ide, jde), kde, spec_bdy_width, 4). • As indicated by the fourth dimensional index, the data is stored in the form of four strips, P_XSB (western), P_XEB (eastern), P_YSB (southern) and P_YEB (northern), each denoting a boundary of specified width, spec_bdy_width that is controlled by the namelist file. • While the second dimension index of these array denotes the maximum vertical dimension, the first index denotes maximum of x and y dimension. Figure below illustrates this array for a case where ide > jde • http: //www. mmm. ucar. edu/wrf/WG 2/software_2. 0/LBC_Documentation. pdf

LBC Arrays unused P_YEB jde P_XEB P_XSB A Given Domain P_YSB sp jds ids ec ide _b dy _w idt h

rconfig entries • # some examples of rconfig entries # NAMELIST 01 rconfig integer time_step namelist, namelist_01 1 # NAMELIST 02 rconfig integer start_year namelist, namelist_02 - ih "time_step“ max_domains 1993 irh "start_year" "4 DIGIT START YR" “YR" # NAMELIST 03 rconfig real dx namelist, namelist_03 max_domains 200 rconfig real dy namelist, namelist_03 max_domains 200 h h "dx" "X HORIZONTAL RES” “M” "dy" “Y HORIZONTAL RES” “M” # NAMELIST 04 rconfig integer nphs namelist, namelist_04 max_domains 12 irh "nmm input " " "

Example of derived data type # following added for nesting purpose (gopal's doing) rconfig real wbd 0 derived max_domains 0 - "wbd 0" “domain western boundary “ rconfig real sbd 0 derived max_domains 0 - "wbd 0" “domain western boundary “ ! grid configuration in mediation_integrate. F parent%wbd 0 = -(IDE -2)*parent%dx parent%sbd 0 = -((JDE-1)/2)*parent%dy

package entries # Package entries within the Registry package package rrtmscheme ra_lw_physics==1 gfdllwscheme ra_lw_physics==99 swradscheme ra_sw_physics==1 gsfcswscheme ra_sw_physics==2 gfdlswscheme ra_sw_physics==99 - - package package package sfclayscheme bl_sfclay_physics==1 myjsfcscheme bl_sfclay_physics==2 slabscheme bl_surface_physics==1 lsmscheme bl_surface_physics==2 ruclsmscheme bl_surface_physics==3 nmmlsmscheme bl_surface_physics==99 ysuscheme bl_physics==1 myjpblscheme bl_physics==2 gfsscheme bl_physics==3 mrfscheme bl_physics==99 - package package kfetascheme cu_physics==1 bmjscheme cu_physics==2 gdscheme cu_physics==3 sasscheme cu_physics==4 kfscheme cu_physics==99 - -

The solver • Quick walk thru’ of (1) WRF. F (2) module_integrate. F (3) solve_em. F http: //www. mmm. ucar. edu/wrf/WG 2/software_v 2/

How to run the wrf model [/emc 2/wx 20 gk/WRFEMC/BHWRFV 1/WRFV 2/test/nmm_real] % ls ETAMPNEW_DATA namelist. input GENPARM. TBL LANDUSE. TBL README. namelist real_nmm. scr RRTM_DATA tr 49 t 67 SOILPARM. TBL tr 49 t 85 VEGPARM. TBL tr 67 t 85 co 2_trans wrf. exe eta_micro_lookup. dat real_nmm. exe wrfinput_d 01 Wrfbdy_d 01 /nest 1/static. wrfsi. rotlat