[gradsusr] Help new to grads

[Huddleston, John]

Bob,

Use a ftp client and go to ftp.cpc.ncep.noaa.gov<ftp://ftp.cpc.ncep.noaa.gov/> and then into the wd52dg/snow directory.

There is a wkly_89x89 directory. Download wk2010.

There is a progs directory and the readsnow.f file. Download readsnow.f file. Uncomment the open(unit=93 line and edit it to read the file wk2010 in your local path. Edit the open(unit=95 line to save the output file of a name of your choosing. Compile the application using 'gfortran -o readsnow readsnow.f'. Run it './readsnow'

There is a grads_snw directory and a control file. Download plt_snw.ctl and edit the top line to be the same name as your output file from unit-95 above.

grads -l
open plt_snw.ctl
d snow

If any of this makes no sense you need to take some basic IT and programming courses.

John Huddleston, PhD, PE

________________________________
From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] On Behalf Of bob duker [dukerbob at gmail.com]
Sent: Tuesday, May 04, 2010 11:08 PM
To: gradsusr at gradsusr.org
Subject: [gradsusr] Help new to grads


Hello All,
I'm doing snow Research and looking for a program to covert the climate prediction center (CPC) Northern Hemosphere raw weekly snow data (in 1's and 0's) into a useable number. I don't have much coding experience or computer science background, so i'm not sure how to covert the data myself, or what program to use. Also I don't have any fortran compiler. The CPC does covert the data into a nice table format, but they don't keep it up todate. Their website is http://www.cpc.ncep.noaa.gov/data/snow/. On their website they have a link to grads.

The program they give for reading weekly 89*89 snow data is

PROGRAM AREA
C*********************************************************************
C
C THIS PROGRAM SIMPLY READS IN THE SNOW DATA AND
C WRITES IT OUT FOR A GIVEN YEAR.
C
C VARIABLES:
C
C    ISNDAT - 89 X 89 SNOW DATA
C
C FILES:
C
C    A YEAR OF WEEKLY SNOW DATA IS STORED IN A FILE WITH
C    THE NAMING CONVENTION:
C
C wk19xx   xx - is the year (ex. wk1973)
C
C*********************************************************************

      DIMENSION ISNDAT(89,89)
      CHARACTER*1 ICHAR1,ICHAR2
      OPEN(11,FILE='INFILE',
     *   FORM='FORMATTED',ACCESS='SEQUENTIAL')
      OPEN(12,FILE='OUTFILE',
     *   FORM='FORMATTED',ACCESS='SEQUENTIAL')
       DO 199 LP = 1,52
C*********************************************************************
C
C READS IN THE WEEKLY SNOW DATA.
C
C*********************************************************************
        DO 200 I = 1,89
        READ(11,100,END=10) IYEARS,IWKS,IROW,ICHAR1,(ISNDAT(I,J),J=1,45)
  100    FORMAT (I4,2I2,A1,45I1)
        READ(11,101,END=10)IYEARS,IWKS,IROW,ICHAR2,(ISNDAT(I,J),J=46,89)
  101    FORMAT (I4,2I2,A1,44I1)
C*********************************************************************
C
C WRITE OUT THE WEEKLY SNOW DATA.
C
C*********************************************************************
         WRITE(12,100) IYEARS,IWKS,IROW,ICHAR1,(ISNDAT(I,J),J=1,45)
         WRITE(12,101) IYEARS,IWKS,IROW,ICHAR2,(ISNDAT(I,J),J=46,89)
  200   CONTINUE
READ(11,102) ICONST
  102    FORMAT(I4)
WRITE(12,102) ICONST
C
  199   CONTINUE
   10   CONTINUE
999     STOP
         END

I'm using windows 7 64bit operating system. My question is how do i use grads to view the snow data or is there a fortran compiler the above program will work on that i should use instead?
Thanks,
Bob R Duke
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100505/5049913d/attachment-0003.html