#!/bin/sh #------------------------------------------------------------- MDL="${HOME}/hysplit4" MAP="${MDL}/graphics/arlmap" MET="${HOME}/TutorialX/files" OUT="." cd $OUT if [ ! -f ASCDATA.CFG ]; then echo "-90.0 -180.0" >ASCDATA.CFG echo "1.0 1.0" >>ASCDATA.CFG echo "180 360" >>ASCDATA.CFG echo "2" >>ASCDATA.CFG echo "0.2" >>ASCDATA.CFG echo "'$MDL/bdyfiles/'" >>ASCDATA.CFG fi echo "### $0 ###" # --------------------------------------------------------------------- # --------- Create CONTROL File --------- echo "57 08 31 12 #STARTING YEAR MONTH DAY HOUR" >CONTROL echo "7 #NUMBER OF SOURCE LOCATIONS" >>CONTROL echo "37.187200 -116.067800 0.0 #SOURCE 1 LATITUDE LONGITUDE HEIGHT(m-agl)" >>CONTROL echo "37.187200 -116.067800 2666.0 #SOURCE 1 LATITUDE LONGITUDE HEIGHT(m-agl)" >>CONTROL echo "37.187200 -116.067800 5332.0 #SOURCE 1 LATITUDE LONGITUDE HEIGHT(m-agl)" >>CONTROL echo "37.187200 -116.067800 8000.0 #SOURCE 1 LATITUDE LONGITUDE HEIGHT(m-agl)" >>CONTROL echo "37.187200 -116.067800 9500.0 #SOURCE 1 LATITUDE LONGITUDE HEIGHT(m-agl)" >>CONTROL echo "37.187200 -116.067800 11000.0 #SOURCE 1 LATITUDE LONGITUDE HEIGHT(m-agl)" >>CONTROL echo "37.187200 -116.067800 12500.0 #SOURCE 1 LATITUDE LONGITUDE HEIGHT(m-agl)" >>CONTROL echo "25 #TOTAL RUN TIME (hours)" >>CONTROL echo "0 #USE MODEL VERTICAL VELOCITY" >>CONTROL echo "25000 #TOP OF MODEL DOMAIN (m-AGL)" >>CONTROL echo "1 #NUMBER OF INPUT DATA GRIDS" >>CONTROL echo "../files/" >>CONTROL echo "195708.hysplit.reanalysis.gbl" >>CONTROL echo "14 #NUMBER OF DEPOSITING POLLUTANTS" >>CONTROL echo "NGAS #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P013 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P012 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P011 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P010 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P009 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P008 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P007 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P006 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P005 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P004 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P003 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P002 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "P001 #POLLUTANT IDENTIFICATION" >>CONTROL echo "0.0 #EMISSION RATE (per hour)" >>CONTROL echo "0.0166666666666667 #HOURS OF EMISSION" >>CONTROL echo "57 08 31 12 30 #RELEASE START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "8 #NUMBER OF CONCENTRATION GRIDS" >>CONTROL echo "0.0 0.0 #CONC GRID CENTER LATITUDE LONGITUDE" >>CONTROL echo "0.03 0.03 #CONC GRID SPACING (degrees) LATITUDE LONGITUDE" >>CONTROL echo "15.0 15.0 #CONC GRID SPAN (degrees) LATITUDE LONGITUDE" >>CONTROL echo "./" >>CONTROL echo "con1hf" >>CONTROL echo "1 #NUMBER OF VERTICAL CONCENTRATION LEVELS" >>CONTROL echo "100 #HEIGHT OF EACH CONCENTRATION LEVEL (m-agl)" >>CONTROL echo "57 08 31 12 30 #SAMPLING START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "00 00 00 12 00 #SAMPLING STOP TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "0 01 00 #SAMPLING INTERVAL: TYPE HOUR MINUTE" >>CONTROL echo "0.0 0.0 #CONC GRID CENTER LATITUDE LONGITUDE" >>CONTROL echo "0.08 0.08 #CONC GRID SPACING (degrees) LATITUDE LONGITUDE" >>CONTROL echo "30.0 60.0 #CONC GRID SPAN (degrees) LATITUDE LONGITUDE" >>CONTROL echo "./" >>CONTROL echo "con1hc" >>CONTROL echo "1 #NUMBER OF VERTICAL CONCENTRATION LEVELS" >>CONTROL echo "100 #HEIGHT OF EACH CONCENTRATION LEVEL (m-agl)" >>CONTROL echo "57 08 31 12 30 #SAMPLING START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "00 00 00 00 00 #SAMPLING STOP TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "0 01 00 #SAMPLING INTERVAL: TYPE HOUR MINUTE" >>CONTROL echo "0.0 0.0 #CONC GRID CENTER LATITUDE LONGITUDE" >>CONTROL echo "0.03 0.03 #CONC GRID SPACING (degrees) LATITUDE LONGITUDE" >>CONTROL echo "15.0 15.0 #CONC GRID SPAN (degrees) LATITUDE LONGITUDE" >>CONTROL echo "./" >>CONTROL echo "dep1hf" >>CONTROL echo "1 #NUMBER OF VERTICAL CONCENTRATION LEVELS" >>CONTROL echo "0 #HEIGHT OF EACH CONCENTRATION LEVEL (m-agl)" >>CONTROL echo "57 08 31 12 30 #SAMPLING START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "00 00 00 12 00 #SAMPLING STOP TIME:YEAR MONTH DAY HOUR MINUTE " >>CONTROL echo "0 01 00 #SAMPLING INTERVAL: TYPE HOUR MINUTE" >>CONTROL echo "0.0 0.0 #CONC GRID CENTER LATITUDE LONGITUDE" >>CONTROL echo "0.08 0.08 #CONC GRID SPACING (degrees) LATITUDE LONGITUDE" >>CONTROL echo "30.0 60.0 #CONC GRID SPAN (degrees) LATITUDE LONGITUDE" >>CONTROL echo "./" >>CONTROL echo "dep1hc" >>CONTROL echo "1 #NUMBER OF VERTICAL CONCENTRATION LEVELS" >>CONTROL echo "0 #HEIGHT OF EACH CONCENTRATION LEVEL (m-agl)" >>CONTROL echo "57 08 31 12 30 #SAMPLING START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "00 00 00 00 00 #SAMPLING STOP TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "0 01 00 #SAMPLING INTERVAL: TYPE HOUR MINUTE" >>CONTROL echo "0.0 0.0 #CONC GRID CENTER LATITUDE LONGITUDE" >>CONTROL echo "0.03 0.03 #CONC GRID SPACING (degrees) LATITUDE LONGITUDE" >>CONTROL echo "15.0 15.0 #CONC GRID SPAN (degrees) LATITUDE LONGITUDE" >>CONTROL echo "./" >>CONTROL echo "conc6hf" >>CONTROL echo "2 #NUMBER OF VERTICAL CONCENTRATION LEVELS" >>CONTROL echo "0 100 #HEIGHT OF EACH CONCENTRATION LEVEL (m-agl)" >>CONTROL echo "57 08 31 12 30 #SAMPLING START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "00 00 00 12 00 #SAMPLING STOP TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "0 01 00 #SAMPLING INTERVAL: TYPE HOUR MINUTE" >>CONTROL echo "0.0 0.0 #CONC GRID CENTER LATITUDE LONGITUDE" >>CONTROL echo "0.08 0.08 #CONC GRID SPACING (degrees) LATITUDE LONGITUDE" >>CONTROL echo "30.0 60.0 #CONC GRID SPAN (degrees) LATITUDE LONGITUDE" >>CONTROL echo "./" >>CONTROL echo "conc1hc" >>CONTROL echo "2 #NUMBER OF VERTICAL CONCENTRATION LEVELS" >>CONTROL echo "0 100 #HEIGHT OF EACH CONCENTRATION LEVEL (m-agl)" >>CONTROL echo "57 08 31 12 30 #SAMPLING START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "00 00 00 00 00 #SAMPLING STOP TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "0 01 00 #SAMPLING INTERVAL: TYPE HOUR MINUTE" >>CONTROL echo "0.0 0.0 #CONC GRID CENTER LATITUDE LONGITUDE" >>CONTROL echo "0.03 0.03 #CONC GRID SPACING (degrees) LATITUDE LONGITUDE" >>CONTROL echo "15.0 15.0 #CONC GRID SPAN (degrees) LATITUDE LONGITUDE" >>CONTROL echo "./" >>CONTROL echo "dep24hf" >>CONTROL echo "1 #NUMBER OF VERTICAL CONCENTRATION LEVELS" >>CONTROL echo "0 #HEIGHT OF EACH CONCENTRATION LEVEL (m-agl)" >>CONTROL echo "57 08 31 12 30 #SAMPLING START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "00 00 00 00 00 #SAMPLING STOP TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "0 24 00 #SAMPLING INTERVAL: TYPE HOUR MINUTE" >>CONTROL echo "0.0 0.0 #CONC GRID CENTER LATITUDE LONGITUDE" >>CONTROL echo "0.08 0.08 #CONC GRID SPACING (degrees) LATITUDE LONGITUDE" >>CONTROL echo "30.0 60.0 #CONC GRID SPAN (degrees) LATITUDE LONGITUDE" >>CONTROL echo "./" >>CONTROL echo "dep24hc" >>CONTROL echo "1 #NUMBER OF VERTICAL CONCENTRATION LEVELS" >>CONTROL echo "0 #HEIGHT OF EACH CONCENTRATION LEVEL (m-agl)" >>CONTROL echo "57 08 31 12 30 #SAMPLING START TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "00 00 00 00 00 #SAMPLING STOP TIME:YEAR MONTH DAY HOUR MINUTE" >>CONTROL echo "0 24 00 #SAMPLING INTERVAL: TYPE HOUR MINUTE" >>CONTROL echo "14 #NUMBER OF DEPOSITING POLLUTANTS" >>CONTROL echo "0.0 0.0 0.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0.0 0.0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "20.0 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "45.0 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "57.5 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "70.0 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "87.5 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "112.5 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "137.5 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "162.5 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "187.5 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "225.0 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "275.0 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "350.0 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL echo "500.0 2.5 1.0 #PARTICLE:DIAMETER (um), DENSITY (g/cc), SHAPE" >>CONTROL echo "0.0 0.0 0.0 0.0 0.0 #DEP VEL (m/s), MW (g/Mole), SFC REACT. RATIO, DIFFUSIVITY RATIO, HENRY'S CONSTANT" >>CONTROL echo "0.0 0 0 #WET REMOVAL: HENRY'S (Molar/atm), IN-CLOUD (L/L), BELOW-CLOUD (1/s)" >>CONTROL echo "0.0 #RADIOACTIVE DECAY HALF-LIFE (days)" >>CONTROL echo "0.0 #POLLUTANT RESUSPENSION (1/m)" >>CONTROL #---------------------------------------------------------- echo "&SETUP " >SETUP.CFG echo "initd = 0, " >>SETUP.CFG echo "numpar = 15000, " >>SETUP.CFG echo "maxpar = 200000, " >>SETUP.CFG echo "ichem = 5, " >>SETUP.CFG echo "nbptyp = 5, " >>SETUP.CFG echo "ndump = 2, " >>SETUP.CFG echo "ncycl = 2, " >>SETUP.CFG echo "poutf = PARDUMP," >>SETUP.CFG echo "EFILE = 'EMITIMES', " >>SETUP.CFG echo "/ " >>SETUP.CFG #---------------------------------------------------------- # --------------------------------------------------------------------- # --------- Create EMITIMES emission File --------- echo "YYYY MM DD HH DURATION(hhhh) #RECORDS" >EMITIMES echo "YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w)" >>EMITIMES echo "1957 08 31 12 9999 98" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 1.24500 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.030600 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.020400 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.030600 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.020400 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.043350 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.030600 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.022950 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.015300 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.010200 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.012750 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.007650 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.005100 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 0.0 0.005100 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 2.49000 0.0 0.0 " >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.061200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.040800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.061200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.040800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.086700 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.061200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.045900 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.030600 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.020400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.025500 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.015300 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.010200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 2666.0 0.010200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 7.47000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.183600 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.122400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.183600 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.122400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.260100 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.183600 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.137700 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.091800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.061200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.076500 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.045900 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.030600 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 5332.0 0.030600 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 14.9400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.367200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.244800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.367200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.244800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.520200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.367200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.275400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.183600 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.122400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.153000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.091800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.061200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 8000.0 0.061200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 14.9400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.367200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.244800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.367200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.244800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.520200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.367200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.275400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.183600 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.122400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.153000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.091800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.061200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 9500.0 0.061200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 8.71500 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.214200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.142800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.214200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.142800 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.303450 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.214200 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.160650 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.107100 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.071400 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.089250 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.053550 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.035700 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 11000.0 0.035700 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES echo "1957 08 31 12 30 0001 37.187200 -116.067800 12500.0 0.00000 0.0 0.0" >>EMITIMES # --------------------------------------------------------------------- # --------- Prepare to run the model --------- rm -f con1hr rm -f dep1hr rm -f conc1hc rm -f conc6hf rm -f dep1hc rm -f dep1hf rm -f dep24hf rm -f dep24hc rm -f airdose.bin rm -f dose96h.bin rm -f accaird.bin rm -f tede.bin rm -f tedesum.bin rm -f tede.ps rm -f tedesum.bin rm -f airdosec.bin rm -f dose96hc.bin rm -f accairdc.bin rm -f tedec.bin rm -f tedecsum.bin rm -f tedec.ps rm -f cdosert.ps rm -f dosert.ps rm -f dose1hc.bin rm -f dose1hcsum.bin rm -f dose6h.bin rm -f dose6hsum.bin rm -f dose96h.bin rm -f dose96hc.bin rm -f dataA_coarse.txt rm -f statA_coarse.txt rm -f hysplit_smoky.txt rm -f MESSAGE rm -f message.txt rm -f PARDUMP rm -f parxplot.ps rm -f MAPTEXT.CFG ${MDL}/exec/hycs_std # --------------------------------------------------------------------- # --------- Create MAPTEXT.CFG File --------- echo " " >MAPTEXT.CFG echo " " >>MAPTEXT.CFG echo " Batch version of IND simulation" >>MAPTEXT.CFG echo " " >>MAPTEXT.CFG echo " Release: lat.: 37.1872 lon.: -116.0678 Hgt: 0 to 12500 m agl" >>MAPTEXT.CFG echo " " >>MAPTEXT.CFG echo " " >>MAPTEXT.CFG echo " Simulation: IND Fuel: U235 Particle Distribution: Glasstone" >>MAPTEXT.CFG echo " Release Quantity: 44 kT Start: 57 08 31 12 30 Duration: 00 hrs, 01 min" >>MAPTEXT.CFG echo "Model run time - 24 hrs" >>MAPTEXT.CFG echo "Meteorology: 0000Z 1 Aug 1957 - REANALYSIS" >>MAPTEXT.CFG echo "Produced by user:" >>MAPTEXT.CFG # --------------------------------------------------------------------- # --------- Create Particle Plot --------- rm -f LABELS.CFG echo 'TITLE&','### %0: cross-section ### &' >LABELS.CFG ${MDL}/exec/parxplot -iPARDUMP -k1 -m1 -z80 -j${MDL}/graphics/arlmap -n3 # --------------------------------------------------------------------- # --------- Creating 4-day TEDE fine grid maps to 12 hrs --------- ${MDL}/exec/con2rem -icon1hf -oairdose.bin -a../files/activity_IND.txt -d1 -e1 -y44 -f1 ${MDL}/exec/con2rem -idep1hf -odose96h.bin -a../files/activity_IND.txt -d1 -y44 -f1 -x96 # sum the air dose over the grid prior to adding it to the dose from deposition ${MDL}/exec/concacc -iairdose.bin -oaccaird.bin ${MDL}/exec/concadd -iaccaird.bin -bdose96h.bin -otede.bin ${MDL}/exec/concsum -itede.bin -otedesum.bin rm -f LABELS.CFG echo 'TITLE&','NOAA HYSPLIT MODEL (Fine Grid)&' >LABELS.CFG echo 'MAPID&','Total Effective Dose Equivalent (TEDE)&' >>LABELS.CFG echo 'LAYER&',' Average&' >>LABELS.CFG echo 'UNITS&','rem&' >>LABELS.CFG echo 'VOLUM&','&' >>LABELS.CFG ${MDL}/exec/concplot -itede.bin -otede.ps +n -s0 -c4 -j${MDL}/graphics/arlmap -91 -urem -k2 -v50::102051000+25::255000051+10::255102000+1::255255000+0.1::000255000 # --------------------------------------------------------------------- # --------- Creating 4-day TEDE coarse grid maps to 24 hrs --------- ${MDL}/exec/con2rem -icon1hc -oairdosec.bin -a../files/activity_IND.txt -d1 -e1 -y44 -f1 ${MDL}/exec/con2rem -idep1hc -odose96hc.bin -a../files/activity_IND.txt -d1 -y44 -f1 -x96 # sum the air dose over the grid prior to adding it to the dose from deposition ${MDL}/exec/concacc -iairdosec.bin -oaccairdc.bin ${MDL}/exec/concadd -iaccairdc.bin -bdose96hc.bin -otedec.bin ${MDL}/exec/concsum -itedec.bin -otedecsum.bin rm -f LABELS.CFG echo 'TITLE&','NOAA HYSPLIT MODEL (Coarse Grid)&' >LABELS.CFG echo 'MAPID&','Total Effective Dose Equivalent (TEDE)&' >>LABELS.CFG echo 'LAYER&',' Average&' >>LABELS.CFG echo 'UNITS&','rem&' >>LABELS.CFG echo 'VOLUM&','&' >>LABELS.CFG ${MDL}/exec/concplot -itedec.bin -otedec.ps +n -s0 -c4 -j${MDL}/graphics/arlmap -91 -urem -k2 -v50::102051000+25::255000051+10::255102000+1::255255000+0.1::000255000 # --------------------------------------------------------------------------------- # ---------Creating dose rate fine grid maps decayed 12 hours to 12 hrs ---------- ${MDL}/exec/con2rem -iconc6hf -odose6h.bin -a../files/activity_IND.txt -d0 -e1 -y44 -f1 -z12 rm -f LABELS.CFG echo 'TITLE&','NOAA HYSPLIT MODEL (Fine Grid)&' >LABELS.CFG echo 'MAPID&','Dose rate&' >>LABELS.CFG echo 'LAYER&',' Average&' >>LABELS.CFG echo 'UNITS&','rem/hr&' >>LABELS.CFG echo 'VOLUM&','&' >>LABELS.CFG ${MDL}/exec/concplot -idose6h.bin -odosert.ps +n -n1:12 -s0 -b0 -c4 -z80 -j${MDL}/graphics/arlmap -91 -urem/hr -k1 -v1::102000000+0.5::255000000+0.1::255102000+0.05::255204000+0.01::255255000+0.005::000000255 # --------------------------------------------------------------------------------- # ---------Creating dose rate coarse grid maps decayed 12 hours to 24 hrs ---------- ${MDL}/exec/con2rem -iconc1hc -odose1hc.bin -a../files/activity_IND.txt -d0 -e1 -y44 -f1 -z12 rm -f LABELS.CFG echo 'TITLE&','NOAA HYSPLIT MODEL (Coarse Grid)&' >LABELS.CFG echo 'MAPID&','Dose rate&' >>LABELS.CFG echo 'LAYER&',' Average&' >>LABELS.CFG echo 'UNITS&','rem/hr&' >>LABELS.CFG echo 'VOLUM&','&' >>LABELS.CFG ${MDL}/exec/concplot -idose1hc.bin -ocdosert.ps +n -n-2 -s0 -b0 -c4 -z80 -j${MDL}/graphics/arlmap -91 -urem/hr -k1 -v1::102000000+0.5::255000000+0.1::255102000+0.05::255204000+0.01::255255000+0.005::000000255 # ---------------------------------------------------------------------------------- # -------------- calculate coarse grid statistics with measurements ---------------- # sum the air dose over the grid prior to adding it to the dose from deposition ${MDL}/exec/concsum -idose1hc.bin -odose1hcsum.bin ${MDL}/exec/c2datem -idose1hcsum.bin -c1000.0 -ohysplit_smoky.txt -m../files/smoky_meas.txt -xi ${MDL}/exec/statmain -d../files/smoky_meas.txt -rhysplit_smoky.txt -o1 -t0 -l10.0 mv dataA.txt dataA_coarse.txt mv statA.txt statA_coarse.txt