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

29.01.2015 Views

Section 9: CALPUFF Table 9-54 (Concluded) Unformatted DFLX.DAT file - Header Header Record NCOM+15 - Volume Source Names (Variable Emissions File) (Included only if NVL2 > 0) No. Variable Type Description 1 ISTYPE integer Source type (8) 2 CNAMVL2 C*16 array Source names C*16 = Character*16 DFLX.DAT File - Data Records The DFLX.DAT data records consist of a set of "NDFOUT+1" records for each hour of the CALPUFF runs (NDFOUT is the number of species flagged as being stored in the output file). The first record of each set contains the date and hour of the data in the records which follow it. The next "NDFOUT" records contain predicted one-hour averaged dry deposition fluxes in g/m 2 /s for each relevant species (see Table 9-55). Sample FORTRAN read statements for the data records (in uncompressed format) are: READ(iunit)nyrb,njulb,nhrb,nsecb,nyre,njule,nhre,nsece READ(iunit) istype,isnum,sname,sxkm,sykm ┌─ │ │ │ │ │ │ │ └─ LOOP OVER DRY DEPOSITED SPECIES STORED ON DISK GRIDDED RECEPTOR DRY FLUXES IF(LSGRID)READ(iunit)CDFG,DFLXG DISCRETE RECEPTOR DRY FLUXES IF(NDREC.GT.0)READ(iunit)CDFD,DFLXD END LOOP OVER DRY DEPOSITED SPECIES STORED ON DISK where the following declarations apply: Character*15 CDFG,CDFD Real DFLXG(nxg,nyg),DFLXD(NDREC) and 9-274

Section 9: CALPUFF nxg = IESAMP - IBSAMP+1 nyg = JESAMP - JBSAMP+1 CALPUFF contains an option to compress the data by replacing strings of zeroes with a coded repetition factor. The factor is a negative number whose absolute value indicates the number of consecutive zeroes that have been replaced by the repetition factor. This method is especially useful in reducing the size of the output file when large segments of the receptor arrays lie upwind of the puffs during an hour, thereby producing long strings of zeroes in the output arrays. For example, the following record with data for 20 receptors requires 20 unpacked words: 0.0, 0.0, 0.0, 0.0, 0.0, 1.2, 3.5, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.7, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 These data in packed form would be represented in six words: -5., 1.2, 3.5, -6., 0.7, -6. Table 9-55: Unformatted DFLX.DAT File - Data Records (Records 1 and 2 of each set) No. Variable Type Description (Record 1) 1 NYRB integer Year (4 digits) at start of averaging period 2 NJULB integer Julian day at start of averaging period 3 NHRB integer Hour (00-23 LST) at start of averaging period 4 NSECB integer Second (0000-3600 LST) at start of averaging period 5 NYRE integer Year (4 digits) at end of averaging period 6 NJULE integer Julian day at end of averaging period 7 NHRE integer Hour (00-23 LST) at end of averaging period 8 NSECE integer Second (0000-3600 LST) at start of averaging period No. Variable Type a Description (Record 2) 1 ISTYPE integer Source type 2 ISNUM integer Source number of this type 3 SNAME C*16 Source name 4 SXKM real Source X-coordinate (km) in the modeling map projection 5 SYKM real Source Y-coordinate (km) in the modeling map projection 9-275

Section 9: CALPUFF<br />

<strong>Table</strong> 9-54 (Concluded)<br />

Unform<strong>at</strong>ted DFLX.DAT file - Header<br />

Header Record NCOM+15 - Volume Source Names (Variable Emissions File)<br />

(Included only if NVL2 > 0)<br />

No. Variable Type Description<br />

1 ISTYPE integer Source type (8)<br />

2 CNAMVL2 C*16 array Source names<br />

C*16 = Character*16<br />

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

<strong>The</strong> DFLX.DAT d<strong>at</strong>a records consist <strong>of</strong> a set <strong>of</strong> "NDFOUT+1" records for each hour <strong>of</strong> the CALPUFF<br />

runs (NDFOUT is the number <strong>of</strong> species flagged as being stored in the output file). <strong>The</strong> first record <strong>of</strong><br />

each set contains the d<strong>at</strong>e and hour <strong>of</strong> the d<strong>at</strong>a in the records which follow it. <strong>The</strong> next "NDFOUT"<br />

records contain predicted one-hour averaged dry deposition fluxes in g/m 2 /s for each relevant species (see<br />

<strong>Table</strong> 9-55).<br />

Sample FORTRAN read st<strong>at</strong>ements for the d<strong>at</strong>a records (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 />

LOOP OVER DRY DEPOSITED SPECIES STORED ON DISK<br />

GRIDDED RECEPTOR DRY FLUXES<br />

IF(LSGRID)READ(iunit)CDFG,DFLXG<br />

DISCRETE RECEPTOR DRY FLUXES<br />

IF(NDREC.GT.0)READ(iunit)CDFD,DFLXD<br />

END LOOP OVER DRY DEPOSITED SPECIES STORED ON DISK<br />

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

Character*15 CDFG,CDFD<br />

Real DFLXG(nxg,nyg),DFLXD(NDREC)<br />

and<br />

9-274

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

Saved successfully!

Ooh no, something went wrong!