Derek Mallia LAIR tutorial 12 /17/ 2013

1. Weather Research and Forecasting – Stochastic Time-Inverted LagrangianTransport (WRF-STILT)model Derek Mallia LAIR tutorial 12/17/2013

2. PART I: STILT basics

3. Backward trajectory simulation using the WRF-STILT model for Hurricane Sandy

4. Basics of the STILT model • The STILT model is a 4D (x,y,z,t)Lagrangian particle dispersion model (LDPM) for atmospheric transport • http://www.stilt-model.org/pmwiki/pmwiki.php • The model was primarily developed to derive upstream influence regions for some receptor that we are trying to model for • Can also be used qualitatively to determine the origin of your air mass • Driven by meteorological models/analyses such as the ECMWF, GDAS, NAM, and WRF • However, WRF output cannot be directly used with the STILT model without converting the netcdf file format to “ARL” format which is used by the HYSPLIT model core in STILT • For optimal results, 3-hourly averaged winds are required (better mass conservation)

5. For non-WRF analyses • Go to: http://www.ready.noaa.gov/archives.php • I generally use these to get a “quick” first look • Also good for using as an outer domain • I.E if your particles leaves your outer WRF domain, you can still track them beyond that domain if you nest it with something like GDAS which is “global” • Downside of these products? • Generally, lower spatial and temporal resolution • Worse mass conservation (see Nehrkorn et at. 2010) • Not necessarily tailored to your needs GDAS domain WRF domain

6. Uinta Basin example: In the inversion… Outside of the inversion…

7. Relevant papers for details Lin, J.C., C. Gerbig, S.C. Wofsy, et al., A near-field tool for simulating the upstream influence of atmospheric observations: The Stochastic Time-Inverted Lagrangian Transport (STILT) model, J. Geophy. Res., 108(D16), 4493, doi:10.1029/2002JD003161, 2003. Nehrkorn, T., J. Eluszkiewicz, S.C. Wofsy, Lin, J.C., C. Gerbig, M. Longo, and S. Freitas. Coupled Weather Research and Forecasting--Stochastic Time-Inverted Lagrangian Transport (WRF-STILT) Model, Meteorology and Atmospheric Physics, 107, 51-64, 2010.

8. Cool model bro, where do I get this model? • First, need to download the STILT model from: • http://www.bgc-jena.mpg.de/bgc-systems/projects/stilt/pmwiki/pmwiki.php?n=FAQ.InitialSetupTasks • Also need to register • Better directions can be found here with FORTRAN hotfix: • http://home.chpc.utah.edu/~u0703457/WRF_STILT_code/STILT • Fixes random model crash… • How to run guide: • http://home.chpc.utah.edu/~u0703457/WRF_STILT_code/STILThowtoV8.pdf • Model uses name list file to run simular to that of WRF, however, can be ran using a “R” using the scripts provided in the stiltR directory • Makes the model less clunky to use, however it requires becoming slightly familiar with R (very easy to use) • Need help? Feel free to bug me (WBB 723)

9. Running STILT at the FORTRAN level • For now lets just copy the hymodelc executable (has all of the necessary bug fixes) along with my modified stiltR code and meteorological input files which can be found here: • /uufs/chpc.utah.edu/common/home/lin-group2/STILT_tutorial • Now lets get familiar with running the STILT at the FORTRAN level…

10. Two “namelist” files for STILT… • See STILThowtoV8.pdf for variable descriptions

11. Only these lines need to be modified… • See STILThowtoV8.pdf formatting of this file…

12. Type hymodelc into you command line once we have finished editing our CONTROL file… • If the CONTROL file was properly edited you code should be running! (will take a few moments)

13. Sample output file…

14. Running STILT using the stiltR interface GDAS • Lets cd into the stiltR directory and get familiar with some of the most important subroutines • Once your have your bearings lets open up automatic_trajplot_WRF.r with your favorite text editor • Lets enter R now (type R into the command line within the stiltR directory and run our tutorial code d01 d02 & 3

15. Sample output Units: PPM / μmolesm-2 s-1

16. PART II: STILT advanced settings

17. STILT model installed now what? • Need to get WRF data in ARL format and need to output specific variables that are required by the STILT model • Involves installing a converter program and editing the WRF registry + recompiling the WRF • Also had to fix some bugs as well • No worries, I supplied the ARL converter code, additional FORTRAN fixes, and easy to read instructions for WRF registry editing: http://home.chpc.utah.edu/~u0703457/WRF_STILT_code/ARL_converter_kingspeak/ • Directions for ARL converter code installation and running can be found here: http://home.chpc.utah.edu/~u0703457/WRF_STILT_code/ARL_CONVERTER • Most of the starting steps can be skipped as I did this already! • Don’t bother with having more than 100 WRF vertical levels… lots of ISSUES • ZSG code can be found here (for STEP #8) :http://home.chpc.utah.edu/~u0703457/WRF_STILT_code/Elena_codes/

18. Additional steps… • Additional all of the files found here: • http://home.chpc.utah.edu/~u0703457/WRF_STILT_code/part_drop_fix/ • This subroutines will need to be copy pasted into your merged_hysplit directory that was used to compile the original STILT model code • HELPFUL TIP, copy the merged_hysplit directory and make a new one which will be used to recompile STILT with the update code above • This is needed to fix some minor mass conservation issues and random particles dropping out for no reason (has to do with particle being dropped from memory as they leave the sub domains) • If you use stiltR, copy and paste Trajec.r to your stiltR directory which adds additional mgmin variable to prevent above issue • If not, additional line needs to be added to the SETUP.CFG file with the following MGMIN = 2000,

19. Sample output…

20. Footprints x Emissions = Footprint of emission source (mole fraction) Units: PPM / μmolesm-2 s-1 Units: micromoles / m2 s PPM /μmoles m-2 s-1 * μmoles / m2 s = PPM of stuff