PROGRAM READSNOW
C$$$ PROGRAM  DOCUMENTATION  BLOCK
C
C PROGRAM: READSNOW         LATITUDE, LONGITUDE TO GRID COORDINATES
C   AUTHOR: GARRETT,DON      ORG: W/NMC52       DATE: 94-01-05
C
C ABSTRACT: THIS PROGRAM READS IN ONE YEAR OF UNFORMATTED WEEKLY SNOW 
C           DATA THAT IS LOCATED IN A 2X2 DEGREE GRID.  THE DATA IS READ
C           IN STARTING FROM THE SOUTH (EQUATOR) AND PROCEEDS UP TO THE
C           NORTH POLE.  NOTE THE WEEKLY DATA HAS BEEN ARCHIVED WITH ARRAY
C           DIMENSIONS OF (180,45).  WHEN THE DATA FIRST GETS READ IN IT 
C           CONVERTS THE DATA TO AN ARRAY WITH DIMENSIONS (180,45).  THIS
C           ARRAY IS THEN WRITTEN OUT TO A DIRECT ACCESS FILE, WHICH CAC
C           USES IN THE GRAPHICS PACKAGE GRADS TO DISPLAY THE DATA.   
C             
C
C PROGRAM HISTORY LOG:
C   94-01-05  GARRETT,D.
C
C   VARIABLES:
C
C     NAMES   DESCRIPTION OF VARIABLES AND TYPES
C     ------ -----------------------------------------------
C     NLAT     LATITUDE (45)
C     NLONG    LONGITUDE (180) 
C     NWEEK    # OF WEEKS IN THE PARTICULAR YEAR (52 OR 53)
C     LRECL    RECORD LENGTH (33120)
C     GRID    OUTPUT ARRAY WITH DIMENSION (180,45)
C
C
C   REMARKS: THE 2X2 DEGREE SNOW FILE WAS WRITTEN AS UNFORMATTED DATA
C             ON A SUN SPARCSTATION 1+.  IF PROBLEMS ARISE IN READING
C             THE DATA E-MAIL: wd52dg@sun1.wwb.noaa.gov
C                     OR CALL: 301(763-8227)
C
C
C$$$
C
        parameter (nlon=180,nlat=45,nweek=1)
        parameter (lrecl=33120)
c..
        dimension grid(nlon,nlat)
c..
c       open(unit=93, form='unformatted',file='/filename/year')
        open(unit=95,form='unformatted',access='direct',recl=lrecl,
     *       file='/output/filename')
c..
       do kweek=1,nweek
c
       print *,' week ',kweek
       read(93) grid
c
C  written to direct access file for GRADS use.
c
       write(95,rec=kweek) grid
c       
       enddo
c
       stop
       end