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 />

Sample FORTRAN write st<strong>at</strong>ements for the CALMET.DAT header records are:<br />

c ---<br />

Header record #1 - File Declar<strong>at</strong>ion -- 24 words<br />

write(iomet) DATASET,DATAVER,DATAMOD<br />

c ---<br />

Header record #2 - Number <strong>of</strong> comment lines -- 1 word<br />

write(iomet) NCOM<br />

c ---<br />

Header record #3 to NCOM+2 (Comment record section) -- 33 words each<br />

write(iomet) COMMENT<br />

c ---<br />

Header record #NCOM+3 - run control parameters -- 39 words<br />

write(iomet) IBYR,IBMO,IBDY,IBHR,IBSEC,<br />

2 IEYR,IEMO,IEDY,IEHR,IESEC,AXTZ,IRLG,IRTYPE,<br />

3 NX,NY,NZ,DGRID,XORIGR,YORIGR, IWFCOD,NSSTA,<br />

4 NUSTA,NPSTA,NOWSTA,NLU,IWAT1,IWAT2,LCALGRD<br />

5 PMAP,DATUM,DATEN,FEAST,FNORTH,UTMHEM,IUTMZN,<br />

6 RNLAT0,RELON0,XLAT1,XLAT2<br />

c ---<br />

Header record #NCOM+4 - cell face heights (NZ + 1 words)<br />

write(iomet)CLAB1,IDUM,IDUM,IDUM,IDUM,ZFACEM<br />

c ---<br />

c ---<br />

Header records #NCOM+5 & 6 - x, y coordin<strong>at</strong>es <strong>of</strong> surface st<strong>at</strong>ions<br />

(NSSTA words each record)<br />

if(nssta.ge.1)then<br />

write(iomet)CLAB2,IDUM,IDUM,IDUM,IDUM,XSSTA<br />

write(iomet)CLAB3,IDUM,IDUM,IDUM,IDUM,YSSTA<br />

endif<br />

c ---<br />

c ---<br />

Header records #NCOM+7 & 8 - x, y coordin<strong>at</strong>es <strong>of</strong> upper air st<strong>at</strong>ions<br />

(NUSTA words each record)<br />

if(nusta.ge.1)then<br />

write(iomet)CLAB4,IDUM,IDUM,IDUM,IDUM,XUSTA<br />

write(iomet)CLAB5,IDUM,IDUM,IDUM,IDUM,YUSTA<br />

endif<br />

c ---<br />

c ---<br />

Header records #NCOM+9 & 10 - x, y coordin<strong>at</strong>es <strong>of</strong> precipit<strong>at</strong>ion st<strong>at</strong>ions<br />

(NPSTA words each record)<br />

if(npsta.ge.1)then<br />

write(iomet)CLAB6,IDUM,IDUM,IDUM,IDUM,XPSTA<br />

write(iomet)CLAB7,IDUM,IDUM,IDUM,IDUM,YPSTA<br />

endif<br />

c ---<br />

Header record #NCOM+11 - surface roughness lengths (NX * NY words)<br />

write(iomet)CLAB8,IDUM,IDUM,IDUM,IDUM,Z0<br />

c ---<br />

Header record #NCOM+12 - land use c<strong>at</strong>egories (NX * NY words)<br />

write(iomet)CLAB9,IDUM,IDUM,IDUM,IDUM,ILANDU<br />

c ---<br />

Header record #NCOM+13 - elev<strong>at</strong>ions (NX * NY words)<br />

write(iomet)CLAB10,IDUM,IDUM,IDUM,IDUM,ELEV<br />

c ---<br />

Header record #NCOM+14 - leaf area index (NX * NY words)<br />

write(iomet)CLAB11,IDUM,IDUM,IDUM,IDUM,XLAI<br />

c ---<br />

c ---<br />

Header record #NCOM+15 - nearest surface st<strong>at</strong>ion no. to each<br />

grid point (NX * NY words)<br />

9-122

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

Saved successfully!

Ooh no, something went wrong!