Homework on CALPUFF

[ Home | Homework | ]


Homework Assignments


Date 2024 Details


11 Mar

Stull and/or Tim covers Panoply installation and use.  

If time, Stull will get a head-start on CALPUFF technical info (Users Guide v5, section 2.1 ) 



13 Mar

Tim Chui leads us in the installation of Calpuff, Calmet, Calpost, CalWRF

Bring your laptop to class for CALPUFF model set up.
.
18 Mar in class

 Tim will finish guiding us on CALPUFF installation.

Skim the CALPUFF-version 6 User Instructions.
In particular, see
Table 1-1 for CALMET major features. (p1-21, or p 32 of the pdf file).
Table 1-2 for CALPUFF major features.  (p 1-26, or p 37 of the pdf file).
Table 1-3 for required meteorol. data. (p 1-29, or p40 of the pdf file)
Table 1-5 summary of input data. (p1-34 , or p45 of the pdf file) 
Read section 7.5 on CALWRF Preprocessor  (starting on p 7-56, or p242 of the pdf file).

Stull will continue discussion of the Users Guide v5 section 2.1 on Puff Averaging and Slugs. 

Be prepared to discuss in class, and compare with AERMOD and HYSPLIT.
What version of the system can read WRF files as meteorological input?
What language is CALPUFF written in?
What operating systems does it run on?

20 Mar

Schedule of presentations, based on sections from the CALPUFF User Guide-v5: see Piazza to sign up for sections to present.
Please create a presentation (roughly 10 - 15 minutes) that does more than just regurgitate the eqs from the Users Manual.  Instead, explain how the eqs work or show how they were derived if it helps our understanding.  Also create opportunities for the class to participate and debate, such as we did with the puff reflections.  You do NOT need to create PowerPoint presentations - - you can use the whiteboard and/or resources from the internet.

Note to ALL: please have your presentation ready before your scheduled date, in case others before you finish early. 

.

Calpuff Run Assignment

=====
Part A.
Create a folder on your own local laptop desktop called yourname-cal. (e.g., stull-cal).
Inside that folder, create subdirectories: calpuff, calmet, calpost, calwrf.
These folders will hold the files that you download.

In part B, you will download some files from Optimum to your local computer, where you can use friendlier tools to view large files.

In part C, you will answer questions about those files that you downloaded.

=====
Part B.
Download files from Optimum to your local laptop, to view them more easily. Note, if you are off campus, you might need to run your VPN first.

1) CALPUFF:
Open a terminal window on your own laptop (do not log into optimum) and use cd to navigate into your own folder that you just created called calpuff.

First, use secure copy to put onto your own local laptop a copy of the:
- input control file INP
- output LST file from the first CALPUFF run that we made
Examples of code to do this will be provided via Piazza.

2) CALMET:
In the terminal window on your own laptop (not on optimum), use cd to navigate into your own folder that you had created called calmet.

Use secure copy to put onto your own local laptop a copy of the
- input control file INP
- output LST file
Examples of code to do this will be provided via Piazza.

3) CALPOST
In the terminal window on your own laptop (not on optimum), use cd to navigate into your own folder that you had created called calpost.

Use secure copy to put onto your own local laptop a copy of the
- input control file INP
- output LST file
Examples of code to do this will be provided via Piazza.

4) CALWRF
In the terminal window on your own laptop (not on optimum), use cd to navigate into your own folder that you had created called calwrf.

Use secure copy to put onto your own local laptop a copy of the
- output WRF.m2d file
Examples of code to do this will be provided via Piazza.

=====
Part C.

1. Each individual student turn submit (via Canvas if possible) the 3 .LST files you downloaded from Part B.

2. Work in groups (at each table in our classroom) to turn in one submission for the whole group, being sure to indicate which students are in your group. Explain what types of info are shown in each of the files, starting at the lines indicated here:
a) calpuff.lst: start at line 2046. (just after the last !END! command)
b) calmet.lst start at line 927 (just after the last !END! command)
c) calpost.lst start at line 833 (just after the last !END! command)
d) WRF.m2d (all lines)
In your description, you can give explanations relative to line numbers.
Example: Lines 1-832 are an echo of the input files.
Lines xxx - yyy show the rain rate.
Etc.
Work together, and submit one document from each group, via email to me. Make sure that everyone in your group know what types of info are in all of the parts of all the .LST and .m2d files

Hints:
Compare the .LST output files with the .INP input files to help you understand what you are looking at.
The input files contain mostly explanation for human users, explaining what the options are.
The portion of input actually used by the programs is contained between the ! ! signs. Look for them.
Take advantage of CALPUFF User Guide v6, and other resources that are available.

-end-


.  This assignment assumes that you have run all the CALPUFF demos on the Optimum computer as taught in class, and that you already created all the input and output files for CALPUFF, CALMET, CALPOST, and CALWRF. Note, if you used different file names than we specified in the demo, then you will likely need to modify the scripts that we provide for this assignment.

25 Mar

Finish presentations, based on sections from the CALPUFF User Guide-v5:

Note to ALL: please have your presentation ready before your scheduled date, in case others before you finish early. 

25 Mar

HW on Slugs (2024), due today (but OK to turn in late).
The code below is written in "R", but you can write your program in any language.
You might be able to utilize snippets from your Gaussian plume homework answers, to enable
you to solve this homework more quickly.

# Slugs
# R. Stull, 17 Feb 2024

# HW Exercise Statement =======
# Theme: Slugs (Gaussian packets of pollutants stretched in the along-wind direction)
#
# Use eq 2-14 (page 2-7 of Calpuff User Guide v5 technical details, written by
# Scire, Strimaitis & Yamartino, 2000, by Earth Tech, Inc.)
# https://www.src.com/calpuff/download/CALPUFF_UsersGuide.pdf
#
# Given:
# Meteorology. Steady state (assume constant in time & space).
U = 8.          # (m/s) mean wind, positive toward the east
tL = 60.         # (s) Lagrangian time scale in both vertical and crosswind directions
sigma_v = 1.5      # (m/s) crosswind turbulence standard deviation
sigma_w = 0.5      # (m/s) vertical turbulence standard deviation
h = 1000.         # (m) boundary layer depth

# Crosswind spread sigma_y: is given by Stull (Pract.Meteor.) eq. (19.13a),
# Assume sigma_x = sigma_y.
# Vertical spread sigma_z: is given by Stull (Pract.Meteor.) eq. (19.13b),

# Pollutant source info.
q = 200.       # (g/s) emission rate of pollutant
He = 50.       # (m) pollutant source emission height after any plume rise
t1minutes = 20.      # Slug initial point (1) was emitted 20 minutes ago.
t2minutes = 10.      # Slug final point (2) was emitted 10 minutes ago.

# domain, assuming origin (x, y, z) = (0,0,0) is a base of smoke stack
xmax = 15000.         # x downwind domain size for x = 0 to xmax (m)
ymax = 1000.         # y crosswind domain size for y = -ymax to +ymax (m)
# All receptors are at the surface, z = 0.

# spatial resolution for calculations of concentration
delx = 200.         # x increment (m)
dely = 50.         # y increment (m)

EXERCISES:
# 1) Plot a graph of the "causality" function F vs. downwind distance x.
# Hints:
# Use Stull's clarification (see lecture notes) of CALPUFF v5 User Guide eq (2-15).
# Compute the horiz. distances x1 and x2 of the initial and final slug points from the source
# In the causality eq, for any location x, the distance da2 = x - x2
# In the causality eq, for any location x, the distance da1 = x - x1

# 2) What is the value of the vertical term "g", for an x location corresponding to the middle of the slug.
# Hints:
# See eq (2-2) of the CALPUFF User Guide v5 for the definition of g .
# This eq already assume reflection from the ground, and assume receptors at the ground z=0.
# For this exercise, assume NO reflection from the top of the boundary layer,
# Hence, do NOT sum over all n (where n = number of reflections).
# Instead, assume n = 0 (i.e., no reflection from top of boundary layer).

# 3) Draw a contour plot of instantaneous concentration C at the surface,
# for a grid of receptors that span the region near the slug where
# concentration is significantly different from zero.
# Hints: In eq (2-14) in the CALPUFF User Guide v5, dc is just crosswind distance
# from the slug centerline. Namely: dc = y .
# Don't forget that sigma_y varies with downwind distance x
# But you can let g be constant, with the value you found in part (2).
# Use a contour interval of 20 µg/m3.

# 4) (a) Discuss how the one slug you plotted would be combined with earlier and later slugs to give concentration along the whole smoke plume.
(b) What are the advantages and disadvantages of using slugs instead of puffs?
(c) The contours you plotted in part (3) were "instantaneous" concentration values. What would you need to do to create time-averaged concentrations, such as are used in the ambient air quality standards.?

 

Use any computer language.
Let me know if I forgot to give you a needed variable. 


 --> switch to using the CMAQ homework schedule now. 


Other Info





[ Home | Homework | ]

http://www.eos.ubc.ca/courses/atsc507/ADM/
Last updated 23 Mar 2024
Copyright © 2016, 2018, 2021, 2024 by Roland Stull.