GPCP 1DD data.

Hi Prince:

You need read data following closely the Fortran code from the website as
below:

   do iday = 1,mday(ii) ; day cycles here
        do nn = 0,nlat-1
           index = nn+(iday-1)*nlat+nstart
           fin(count,:,nn) = fbindirread(fils(ii),index,(/nlon/),"float")
        end do
  end do

Notice nstart = 1, as the first recode is description of the GPCP 1DD data.

I am reading the data for a specific year at once. "count" is looping over
the day, "fils(ii)" are all the monthly data files.

For a month, you can change my code as you need.

Yi Wang, PhD
Pacific Northwest National Lab
PO BOX 999, MSIN K9-24
Richland, Washington 99352
USA
Phone: (509) 375-2386
Fax: (509) 375-6448

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Nov 26 2007 - 12:59:00 MST