calpuf.inp

Changes made to calpuf.inp

• METDAT – Change to data/calmet.dat
• RSTARTE – Uncomment and set to restart_so2_e.dat. This is the name given to the generated restart file at simulation end.
• OZDAT – Comment out
• DEBUG – Uncomment and set to debug.dat (we want an output debug file)
• BALDAT – Uncomment and set to massbal.dat (we want an output mass balance file)
• LCFILES – Set to T
• IBYR, IBMO, IBDY, IBHR – Set to required simulation start date/- time (e.g. 2017, 12, 04, 0)
• IEYR, IEMO, IEDY, IEHR – Set to required simulation end date/- time (e.g. 2017, 12, 06, 0)
• ABTZ – Set to UTC+0000 (we want times in UTC)
• NSPEC – Set to 2 (we have SO2 and SO4)
• NSE – Set to 1 (only one species actually emitted; SO4 is created in the atmosphere)
• MRESTART – Set to 2 for main run (i.e. write a restart file at the end of the run)
• MTRANS – Set to 0 (don’t model transitional plume rise)
• MTIP – Set to 0 (don’t model stack downwash)
• MBDW – (not needed?) Set to 2 (use PRIME building downwash method)
• MSHEAR – Set to 1 (vertical wind shear modelled)
• MCHEM – Set to 1 (use MESOPUFF II cheical scheme)
• MWET – Set to 0 (don’t model wet removal)
• MPARTL – Set to 0 (don’t model partial plume penetration)
• MPARTLBA – Set to 0 (don’t model partial plume penetration)
• CSPEC – Delete all but SO2 and SO4 from both the list and the table
• IUTMZN – Set to 16
• DATUM – Set to WGS-84
• NX – Set to 91 (1km res) or 901 (100m res)
• NY – Set to 55 (1km res) or 541 (100m res)
• DGRIDKM – Set to 1.0 (1km res) or 0.1 (100m res)
• NZ – Set to 24
• ZFACE – Set to "0,20,60,140,260,540,1000,2000,3000,4000,5000,6000,7000, 8000,9000,10000,11000,12000,13000,14000,15000,16000,17000,18000,19000"
• XORIGKM – Set to 525.09417
• YORIGKM – Set to 1294.23926
• IECOMP – Set to 91 (1km res) or 901 (100m res) (X index of UR corner)
• JECOMP – Set to 55 (1km res) or 541 (100m res) (J index of UR corner)
• LSAMP – Set to F (no gridded receptors)
• IDRY – Set to 0 (no dry fluxes output)
• IWET – Set to 0 (no wet fluxes output)
• IVIS – Set to 0 (no RH output)
• LCOMPRS – Set to F (no compression to output file)
• IPFTRAK – Set to 0 (no puff-tracking output)
• INRISE – Set to 1 (create plume properties output file rise.dat)
• ICFRQ – Set to 24 (concentration print interval in timesteps)
• IPRTU – Set to 1 (use g/m3 for printing)
• SPECIES LIST – Delete all but SO2 and SO4, and set both to "1, 1, 0, 0, 0, 0, 0" (just print and save concentrations)
• LDEBUG – Set to T (debug output wanted) NOTE: Set this to F to dramatically reduce .lst file size
• NPFDEB – Set to 1044 (track this many puffs in debug)
• NN2 – Set to 24 (stop debug output after this number of timesteps)
• XHILL2M – (not needed?) Set to 0 (factor to convert horizontal dimensions to meters if MHILL is 1)
• ZHILL2M – (not needed?) Set to 0 (factor to convert vertical dimensions to meters if MHILL is 1)
• INPUT GROUP 7 – (not needed?) Delete all but SO2 line
• INPUT GROUP 8 – (not needed?) Delete all but SO4 line
• INPUT GROUP 10 – (not needed?) Delete all but SO2 line, and change liquid precip to 2.61E-05 and frozen precip to 3.0E-05
• MOZ – Set to 0 (use monthly background ozone value)
• BCKO3 – Set all values to 30 (monthly background ozone)
• BCKNH3 – Set all values to 1.0 (monthly background ammonia)
• XMXLEN – Set to 0.5 (maximum length of a slug)
• NCOUNT – Set to 5 (number of iterations when computing transport wind)
• SZCAP_M – Set to -1.0 (disable)
• WSCALM – Set to 0.01 (minimum wind speed for non-calm conditions)
• XMINZI – Set to 20.0 (minimum mixing height)
• SL2PF – Set to 5.0 (slug-to-puff transition criterion factor)
• SYSPLITH – Set to 0.8 (min sig-y before split)
• SHSPLITH – Set to 1.0 (min puff elongation rate before split)
• RSAMPBC – Set to 15.0 (search radius for sampling nearest BC puff)
• MDEPBC – Set to 0 (no adjustment for depletion)
• NPT1 – Set to 0 (no point sources to follow)
• NSPT1 – Set to 0 (no variable emissions factors)
• Subgroup 13b – Delete all sources
• Subgroup 13d – Delete all sources
• NAR1 – Set to 1 (one area source to follow)
• Subgroup 14b – Add following lines (source parameters):

! SRCNAM = MASAYA2 !
! X = 0, 520, 50, 0.004 !
!END!

• Subgroup 14c – Add following lines (source coordinates):

1 ! SRCNAM = MASAYA2 !
1 ! XVERT = 590.495, 590.740, 590.702, 590.435!
1 ! YVERT = 1324.876, 1324.831, 1325.058, 1325.078!
!END!

• NREC – Set to 18000 (number of non-gridded receptor points to follow)
• NRGRP – Set to 0 (no groups of receptors)
• Subgroup 20b – Delete all groups
• Subgroup 20c – Paste in 18000 non-gridded receptor points from Sara’s input file.