@@ -65,7 +65,7 @@ print_usage() {
6565 or ` conda activate unresp`
6666 or ` load your system python libraries`
6767 ^^^ these fixes can be added to .env file for bespoke Setup
68- "
68+ "
6969}
7070
7171set_viz () {
@@ -77,7 +77,7 @@ set_ffmpeg() {
7777}
7878
7979set_SO4 () {
80- setSO4=false
80+ setSO4=false
8181}
8282
8383set_model () {
@@ -98,9 +98,9 @@ while getopts 'd:n:pmsfh' flag; do
9898 m) set_model ;;
9999 f) set_ffmpeg ;;
100100 h) print_usage
101- exit 1 ;;
101+ exit 1 ;;
102102 * ) print_usage
103- exit 1 ;;
103+ exit 1 ;;
104104 esac
105105done
106106
@@ -141,18 +141,18 @@ cwd=$(pwd)
141141# ------------------- DO NOT ALTER BELOW THIS LINE------------------------#
142142# ------------------------------------------------------------------------#
143143if [ " $runmodel " = true ]; then
144- echo " ### RUNNING FORECAST SYSTEM FOR DATE " ${rundate} " ###"
144+ echo " ### RUNNING FORECAST SYSTEM FOR DATE " ${rundate} " ###"
145145fi
146146# ##TERREL###
147147if [ " $runTERREL " = true ]; then
148148 # Compile TERREL if required:
149149 cd CALPUFF_EXE
150150 if [ ! -f ./terrel_intel.exe ]; then
151- echo -n " ### COMPILING TERREL"
152- ifort -O0 -fltconsistency -w ../CALPUFF_SRC/TERREL/terrel.for -o terrel_intel.exe
153- echo " ---> FINISHED ###"
151+ echo -n " ### COMPILING TERREL"
152+ ifort -O0 -fltconsistency -w ../CALPUFF_SRC/TERREL/terrel.for -o terrel_intel.exe
153+ echo " ---> FINISHED ###"
154154 else
155- echo " ### TERREL ALREADY COMPILED ###"
155+ echo " ### TERREL ALREADY COMPILED ###"
156156 fi
157157 cd ..
158158 # Remove any old files before running:
@@ -181,11 +181,11 @@ if [ "$runCTGPROC" = true ]; then
181181 # Compile CTGPROC if required:
182182 cd CALPUFF_EXE
183183 if [ ! -f ./ctgproc_intel.exe ]; then
184- echo -n " ### COMPILING CTGPROC"
185- ifort -O0 -fltconsistency -mcmodel=medium -w ../CALPUFF_SRC/CTGPROC/ctgproc.for -o ctgproc_intel.exe
186- echo " ---> FINISHED ###"
184+ echo -n " ### COMPILING CTGPROC"
185+ ifort -O0 -fltconsistency -mcmodel=medium -w ../CALPUFF_SRC/CTGPROC/ctgproc.for -o ctgproc_intel.exe
186+ echo " ---> FINISHED ###"
187187 else
188- echo " ### CTGPROC ALREADY COMPILED ###"
188+ echo " ### CTGPROC ALREADY COMPILED ###"
189189 fi
190190 cd ..
191191 # Remove any old files before running:
@@ -214,11 +214,11 @@ if [ "$runMAKEGEO" = true ]; then
214214 # Compile MAKEGEO if required:
215215 cd CALPUFF_EXE
216216 if [ ! -f ./makegeo_intel.exe ]; then
217- echo -n " ### COMPILING MAKEGEO"
218- ifort -O0 -fltconsistency -w ../CALPUFF_SRC/MAKEGEO/makegeo.for -o makegeo_intel.exe
219- echo " ---> FINISHED ###"
217+ echo -n " ### COMPILING MAKEGEO"
218+ ifort -O0 -fltconsistency -w ../CALPUFF_SRC/MAKEGEO/makegeo.for -o makegeo_intel.exe
219+ echo " ---> FINISHED ###"
220220 else
221- echo " ### MAKEGEO ALREADY COMPILED ###"
221+ echo " ### MAKEGEO ALREADY COMPILED ###"
222222 fi
223223 cd ..
224224 # Copy data files from TERREL and CTGPROC across to the data directory
@@ -268,18 +268,18 @@ if [ "$run3DDAT" = true ]; then
268268 for i in ` seq 0 3 48` ; do
269269 hour=` printf " %02d" $i `
270270 if [ ! -f nam.t00z.afwaca${hour} .tm00.grib2 ]; then
271- echo " ### DOWNLOADING DATA FOR FORECAST HOUR " ${hour} " ###"
272- # Entire GRIB file:
273- # wget http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.${rundate}/nam.t00z.afwaca${hour}.tm00.grib2
274- # Subset of GRIB file using GRIB filter (http://nomads.ncep.noaa.gov/cgi-bin/filter_nam_crb.pl):
275- # WARNING https not http as of Jan 2019
276- curl " https://nomads.ncep.noaa.gov/cgi-bin/filter_nam_crb.pl?file=nam.t00z.afwaca" ${hour} " .tm00.grib2&" \
277- " lev_1000_mb=on&lev_100_mb=on&lev_10_mb=on&lev_150_mb=on&lev_200_mb=on&lev_20_mb=on&lev_250_mb=on&" \
278- " lev_2_mb=on&lev_300_mb=on&lev_30_mb=on&lev_400_mb=on&lev_500_mb=on&lev_50_mb=on&lev_5_mb=on&" \
279- " lev_600_mb=on&lev_700_mb=on&lev_75_mb=on&lev_7_mb=on&lev_800_mb=on&lev_850_mb=on&lev_900_mb=on&" \
280- " lev_925_mb=on&lev_950_mb=on&lev_mean_sea_level=on&var_HGT=on&var_PRMSL=on&var_RH=on&var_TMP=on&var_UGRD=on&var_VGRD=on&" \
281- " var_DZDT=on&subregion=&leftlon=272&rightlon=278&toplat=16&bottomlat=10&dir=%2Fnam." ${rundate} \
282- -o nam.t00z.afwaca${hour} .tm00.grib2
271+ echo " ### DOWNLOADING DATA FOR FORECAST HOUR " ${hour} " ###"
272+ # Entire GRIB file:
273+ # wget http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.${rundate}/nam.t00z.afwaca${hour}.tm00.grib2
274+ # Subset of GRIB file using GRIB filter (http://nomads.ncep.noaa.gov/cgi-bin/filter_nam_crb.pl):
275+ # WARNING https not http as of Jan 2019
276+ curl " https://nomads.ncep.noaa.gov/cgi-bin/filter_nam_crb.pl?file=nam.t00z.afwaca" ${hour} " .tm00.grib2&" \
277+ " lev_1000_mb=on&lev_100_mb=on&lev_10_mb=on&lev_150_mb=on&lev_200_mb=on&lev_20_mb=on&lev_250_mb=on&" \
278+ " lev_2_mb=on&lev_300_mb=on&lev_30_mb=on&lev_400_mb=on&lev_500_mb=on&lev_50_mb=on&lev_5_mb=on&" \
279+ " lev_600_mb=on&lev_700_mb=on&lev_75_mb=on&lev_7_mb=on&lev_800_mb=on&lev_850_mb=on&lev_900_mb=on&" \
280+ " lev_925_mb=on&lev_950_mb=on&lev_mean_sea_level=on&var_HGT=on&var_PRMSL=on&var_RH=on&var_TMP=on&var_UGRD=on&var_VGRD=on&" \
281+ " var_DZDT=on&subregion=&leftlon=272&rightlon=278&toplat=16&bottomlat=10&dir=%2Fnam." ${rundate} \
282+ -o nam.t00z.afwaca${hour} .tm00.grib2
283283 fi
284284 done
285285 cd ../../..
@@ -299,11 +299,11 @@ if [ "$runCALMET" = true ]; then
299299 # Compile CALMET if required:
300300 cd CALPUFF_EXE
301301 if [ ! -f ./calmet_intel.exe ]; then
302- echo -n " ### COMPILING CALMET"
303- ifort -O0 -fltconsistency -mcmodel=medium -w ../CALPUFF_SRC/CALMET/calmet.for -o calmet_intel.exe
304- echo " ---> FINISHED ###"
302+ echo -n " ### COMPILING CALMET"
303+ ifort -O0 -fltconsistency -mcmodel=medium -w ../CALPUFF_SRC/CALMET/calmet.for -o calmet_intel.exe
304+ echo " ---> FINISHED ###"
305305 else
306- echo " ### CTGPROC ALREADY COMPILED ###"
306+ echo " ### CTGPROC ALREADY COMPILED ###"
307307 fi
308308 cd ..
309309 # Remove any old data files and copy relevant new files into the data directory
@@ -321,8 +321,8 @@ if [ "$runCALMET" = true ]; then
321321 # Update input file:
322322 echo -n " ### SETTING UP CALMET INPUT FILE"
323323 sed -e " s/YYYYb/$startYear /g" -e " s/MMb/$startMonth /g" -e " s/DDb/$startDay /g" -e " s/YYYYe/$endYear /g" \
324- -e " s/MMe/$endMonth /g" -e " s/DDe/$endDay /g" -e " s/?3DDAT?/met_${rundate} .dat/g" \
325- -e " s/?NX?/$NX /g" -e " s/?NY?/$NY /g" -e " s/?DGRIDKM?/$DGRIDKM /g" ./CALPUFF_INP/calmet_template.inp > ./CALPUFF_INP/calmet.inp
324+ -e " s/MMe/$endMonth /g" -e " s/DDe/$endDay /g" -e " s/?3DDAT?/met_${rundate} .dat/g" \
325+ -e " s/?NX?/$NX /g" -e " s/?NY?/$NY /g" -e " s/?DGRIDKM?/$DGRIDKM /g" ./CALPUFF_INP/calmet_template.inp > ./CALPUFF_INP/calmet.inp
326326 echo " ---> FINISHED ###"
327327 # Run CALMET:
328328 echo " ### RUNNING CALMET"
339339if [ " $runCALPUFF " = true ]; then
340340 # Compile CALPUFF if required:
341341 if [ ! -f ./CALPUFF_EXE/calpuff_intel.exe ]; then
342- echo -n " ### COMPILING CALPUFF"
343- cd CALPUFF_SRC/CALPUFF
344- ifort -c modules.for
345- cd ../../CALPUFF_EXE
346- ifort -O0 -fltconsistency -mcmodel=medium -w ../CALPUFF_SRC/CALPUFF/calpuff.for ../CALPUFF_SRC/CALPUFF/modules.o -o calpuff_intel.exe
347- cd ..
348- echo " ---> FINISHED ###"
342+ echo -n " ### COMPILING CALPUFF"
343+ cd CALPUFF_SRC/CALPUFF
344+ ifort -c modules.for
345+ cd ../../CALPUFF_EXE
346+ ifort -O0 -fltconsistency -mcmodel=medium -w ../CALPUFF_SRC/CALPUFF/calpuff.for ../CALPUFF_SRC/CALPUFF/modules.o -o calpuff_intel.exe
347+ cd ..
348+ echo " ---> FINISHED ###"
349349 else
350- echo " ### CALPUFF ALREADY COMPILED ###"
350+ echo " ### CALPUFF ALREADY COMPILED ###"
351351 fi
352352 # Remove old and copy new CALMET data file across to the data directory
353353 echo -n " ### SETTING UP DATA DIRECTORY"
@@ -370,10 +370,10 @@ if [ "$runCALPUFF" = true ]; then
370370 echo -n " ---> NO RESTART FILE FOUND"
371371 fi
372372 sed -e " s/YYYYb/$startYear /g" -e " s/MMb/$startMonth /g" -e " s/DDb/$startDay /g" -e " s/YYYYe/$midYear /g" \
373- -e " s/MMe/$midMonth /g" -e " s/DDe/$midDay /g" -e " s/?METDAT?/calmet_${rundate} .dat/g" \
374- -e " s/?RSTARTB?/restart_$rundate .dat/g" -e " s/?RSTARTE?/restart_$middate .dat/g" \
375- -e " s/?MRES?/$mres /g" -e " s/?NX?/$NX /g" -e " s/?NY?/$NY /g" -e " s/?DGRIDKM?/$DGRIDKM /g" \
376- ./CALPUFF_INP/calpuff_template.inp > ./CALPUFF_INP/calpuff.inp
373+ -e " s/MMe/$midMonth /g" -e " s/DDe/$midDay /g" -e " s/?METDAT?/calmet_${rundate} .dat/g" \
374+ -e " s/?RSTARTB?/restart_$rundate .dat/g" -e " s/?RSTARTE?/restart_$middate .dat/g" \
375+ -e " s/?MRES?/$mres /g" -e " s/?NX?/$NX /g" -e " s/?NY?/$NY /g" -e " s/?DGRIDKM?/$DGRIDKM /g" \
376+ ./CALPUFF_INP/calpuff_template.inp > ./CALPUFF_INP/calpuff.inp
377377 echo " ---> FINISHED ###"
378378 # Run CALPUFF for first 24 hours:
379379 echo " ### RUNNING CALPUFF FOR FIRST 24 HOURS"
@@ -389,10 +389,10 @@ if [ "$runCALPUFF" = true ]; then
389389 # Set up input file for second 24hrs:
390390 echo -n " ### SETTING UP CALPUFF INPUT FILE FOR SECOND 24 HOURS"
391391 sed -e " s/YYYYb/$midYear /g" -e " s/MMb/$midMonth /g" -e " s/DDb/$midDay /g" -e " s/YYYYe/$endYear /g" \
392- -e " s/MMe/$endMonth /g" -e " s/DDe/$endDay /g" -e " s/?METDAT?/calmet_${rundate} .dat/g" \
393- -e " s/?RSTARTB?/restart_$middate .dat/g" -e " s/?RSTARTE?/restart_$enddate .dat/g" \
394- -e " s/?MRES?/1/g" -e " s/?NX?/$NX /g" -e " s/?NY?/$NY /g" -e " s/?DGRIDKM?/$DGRIDKM /g" \
395- ./CALPUFF_INP/calpuff_template.inp > ./CALPUFF_INP/calpuff.inp
392+ -e " s/MMe/$endMonth /g" -e " s/DDe/$endDay /g" -e " s/?METDAT?/calmet_${rundate} .dat/g" \
393+ -e " s/?RSTARTB?/restart_$middate .dat/g" -e " s/?RSTARTE?/restart_$enddate .dat/g" \
394+ -e " s/?MRES?/1/g" -e " s/?NX?/$NX /g" -e " s/?NY?/$NY /g" -e " s/?DGRIDKM?/$DGRIDKM /g" \
395+ ./CALPUFF_INP/calpuff_template.inp > ./CALPUFF_INP/calpuff.inp
396396 echo " ---> FINISHED ###"
397397 # Run CALPUFF for second 24 hours:
398398 echo " ### RUNNING CALPUFF FOR SECOND 24 HOURS"
@@ -419,15 +419,15 @@ if [ "$runVIS" = true ]; then
419419 mkdir ./vis/${rundate}
420420 cd Python
421421 if [ " $setSO4 " = true]; then
422- ./genmaps.py ${rundate}
422+ ./genmaps.py ${rundate}
423423 else
424- ./genmaps.py ${rundate} --SO4
424+ ./genmaps.py ${rundate} --SO4
425425 fi
426426 cd ..
427427 cd vis/${rundate}
428428 if [ ${runffmpeg} = true ]; then
429- echo " Running ffmpeg"
430- ffmpeg -f image2 -r 4 -i SO2_static_concrec0100%02d.png -vcodec mpeg4 -y -s 7680x4320 movie_${rundate} .mp4
429+ echo " Running ffmpeg"
430+ ffmpeg -f image2 -r 4 -i SO2_static_concrec0100%02d.png -vcodec mpeg4 -y -s 7680x4320 movie_${rundate} .mp4
431431 fi
432432 cd ../..
433433 echo " ---> FINISHED ###"
0 commit comments