29.01.2015 Views

Table of Contents - The Atmospheric Studies Group at TRC

Table of Contents - The Atmospheric Studies Group at TRC

Table of Contents - The Atmospheric Studies Group at TRC

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Section 9: CALPUFF<br />

CONC.DAT File - D<strong>at</strong>a Records<br />

<strong>The</strong> CONC.DAT d<strong>at</strong>a records consist <strong>of</strong> "NSPOUT+2" sets <strong>of</strong> records for each hour <strong>of</strong> the CALPUFF run<br />

and for each source for each hour <strong>of</strong> the CALPUFF run if the source contribution option is selected<br />

(NSPOUT is the number <strong>of</strong> output species in the CALPUFF run). <strong>The</strong> first record <strong>of</strong> each set contains<br />

the d<strong>at</strong>e and time period, and the second contains the identity <strong>of</strong> the source <strong>of</strong> the d<strong>at</strong>a in the records<br />

which follow it. <strong>The</strong> next "NSPOUT" record groups contain the predicted concentr<strong>at</strong>ions in g/m 3 , for<br />

each species flagged for output in the control file. If the source contribution option is active, there will be<br />

one output record set for each source, followed by an additional record set for the total due to all sources.<br />

<strong>The</strong> source name used for the total due to all sources is “TOTAL”, the type is “0”, and the number is “1”.<br />

<strong>Table</strong> 9-53 describes the variables.<br />

Sample FORTRAN read st<strong>at</strong>ements for one d<strong>at</strong>a record set (in uncompressed form<strong>at</strong>) are:<br />

READ(iunit)nyrb,njulb,nhrb,nsecb,nyre,njule,nhre,nsece<br />

READ(iunit) istype,isnum,sname,sxkm,sykm<br />

┌─<br />

│<br />

│<br />

│<br />

│<br />

│<br />

│<br />

│<br />

│<br />

│<br />

│<br />

└─<br />

LOOP OVER OUTPUT SPECIES<br />

GRIDDED RECEPTOR CONCENTRATIONS<br />

IF(LSGRID)READ(iunit)CSPECG,CONCG<br />

DISCRETE RECEPTOR CONCENTRATIONS<br />

IF(NDREC.GT.0)READ(iunit)CSPECD,CONCD<br />

COMPLEX TERRAIN RECEPTOR CONCENTRATIONS<br />

IF(NCTREC.GT.0)READ(iunit)CSPECCT,CONCCT<br />

END LOOP OVER OUTPUT SPECIES<br />

where the following declar<strong>at</strong>ions apply:<br />

and<br />

Character*15 CSPECG,CPSECD,CSPECCT<br />

Character*16 SNAME<br />

Real CONCG(nxg,nyg),CONCD(NDREC),CONCCT(NCTREC)<br />

nxg = IESAMP - IBSAMP+1<br />

nyg = JESAMP - JBSAMP+1<br />

CALPUFF contains an option to compress the d<strong>at</strong>a by replacing strings <strong>of</strong> zeroes with a coded repetition<br />

factor. <strong>The</strong> factor is a neg<strong>at</strong>ive number whose absolute value indic<strong>at</strong>es the number <strong>of</strong> consecutive zeroes<br />

th<strong>at</strong> have been replaced by the repetition factor. This method is especially useful in reducing the size <strong>of</strong><br />

the output file when large segments <strong>of</strong> the receptor arrays lie upwind <strong>of</strong> the puffs during an hour, thereby<br />

producing long strings <strong>of</strong> zeroes in the output arrays. For example, the following record with d<strong>at</strong>a for 20<br />

receptors requires 20 unpacked words:<br />

0.0, 0.0, 0.0, 0.0, 0.0, 1.2, 3.5, 0.0, 0.0, 0.0,<br />

0.0, 0.0, 0.0, 0.7, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0<br />

<strong>The</strong>se d<strong>at</strong>a in packed form would be represented in six words:<br />

9-253

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!