29.01.2015 Views

CALPUFF and Postprocessors

CALPUFF and Postprocessors

CALPUFF and Postprocessors

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

integer ILANDU(nx,ny),NEARS(nx,ny)<br />

character*132 COMMENT(ncom)<br />

character*64 DATAMOD<br />

character*16 DATASET,DATAVER<br />

character*12 DATEN<br />

character*8 PMAP,DATUM<br />

character*8 CLAB1,CLAB2,CLAB3,CLAB4,CLAB5,CLAB6<br />

character*8 CLAB7,CLAB8,CLAB9,CLAB10,CLAB11,CLAB12<br />

character*4 UTMHEM<br />

logical LCALGRD<br />

CALMET.DAT File - Data Records<br />

The CALMET.DAT data records include hourly fields of winds <strong>and</strong> meteorological variables. In addition<br />

to the regular CALMET output variables, both CALGRID <strong>and</strong> <strong>CALPUFF</strong> require additional threedimensional<br />

fields of air temperature <strong>and</strong> vertical velocity. The presence of these fields in the CALMET<br />

output file is flagged by the header record logical variable, LCALGRD, having a value of TRUE.<br />

The data records contain three-dimensional gridded fields of U, V, <strong>and</strong> W wind components <strong>and</strong> air<br />

temperature, <strong>and</strong> two-dimensional fields of PGT stability class, surface friction velocity, mixing height,<br />

Monin-Obukhov length, convective velocity scale, precipitation rate (not used by CALGRID), nearsurface<br />

temperature, air density, short-wave solar radiation, relative humidity, <strong>and</strong> precipitation type<br />

codes (not used by CALGRID). A description of each variable in the data records is provided in Table F-<br />

5.<br />

Sample FORTRAN write statements for the CALMET.DAT data records are:<br />

c ---<br />

c ---<br />

c ---<br />

Write U, V, W wind components<br />

+))) Loop over vertical layers, k<br />

*<br />

* write(iunit)CLABU,NDATHR,((U(i,j,k),i=1,nx),j=1,ny)<br />

* write(iunit)CLABV,NDATHR,((V(i,j,k),i=1,nx),j=1,ny)<br />

* if(LCALGRD)write(iunit)CLABW,NDATHR((W(i,j,k+1),i=1,nx),j=1,ny)<br />

*<br />

.))) End loop over vertical layers<br />

Write 3-D temperature field<br />

if(LCALGRD.<strong>and</strong>.irtype.eq.1) then<br />

+))) Loop over vertical layers, k<br />

*<br />

* write(iunit)CLABT,NDATHR,((ZTEMP(i,j,k),i=1,nxm),j=1,nym)<br />

*<br />

.))) End loop over vertical layers<br />

endif<br />

Write 2-D meteorological fields<br />

if(irtype.eq.1) then<br />

write(iunit)CLABSC,NDATHR,IPGT<br />

write(iunit)CLABUS,NDATHR,USTAR<br />

write(iunit)CLABZI,NDATHR,ZI<br />

write(iunit)CLABL,NDATHR,EL<br />

write(iunit)CLABWS,NDATHR,WSTAR<br />

write(iunit)CLABRMM,NDATHR,RMM<br />

MAR 2006 – <strong>CALPUFF</strong><br />

F-92

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

Saved successfully!

Ooh no, something went wrong!