Re: GPCP 1dd binary read

Hi Prince,

Greetings from Korea (how r things?).

There is the possibility that the first record is empty. GPCP data may have
data gaps. You may have to catch those. I do the following (my data is monthly
mean data though).

;--------------%< code

setfileoption("bin","ReadByteOrder","BigEndian")

iyrec=0
do iyr = yr1,yr2
  infile="gpcp_v2_psg."+iyr
  irec=1
  do imon=1,nmon
    do j = 1,nlat
      prec(iyrec,j-1,:) = fbindirread(infile, irec, nlon, "float")
      irec=irec+1
    end do
      if all(ismissing(prec(iyrec,:,:)))
       print("All Missing Values during month "+imon)
      break
      end if
    iyrec=iyrec+1
  end do
end do
;--------------%< code

saji
...

* Prince K. XAVIER <xavier_at_lmd.ens.fr> [2007-11-26 17:59:42 +0100]:

> Hi all,
>
> I am trying to read the GPCP 1DD data at
> http://www1.ncdc.noaa.gov/pub/data/gpcp/1dd/data/
> using the options:
>
> setfileoption("bin","ReadByteOrder","BigEndian")
> prec = fbindirread("gpcp_1dd_p1d.199707",1,(/180,360/),"float")
>
> But I find get the first record empty. I think I am doing something wrong
> with the record length. If someone has read this data using NCL, could you
> help to set the correct file options?
>
> thanks in advance
> Prince
>
> --
> Prince K. XAVIER
> Laboratoire de Meteorologie Dynamique
> Ecole Normale Superieure
> 24 Rue Lhomond 75231 Paris cedex 05 FRANCE
> Tel: +33 (0)1 44 32 22 32
> Fax: +33 (0)1 43 36 83 92
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
Saji N. Hameed
APEC Climate Center          				+82 51 668 7470
National Pension Corporation Busan Building 12F         
Yeonsan 2-dong, Yeonje-gu, BUSAN 611705			saji_at_apcc21.net
KOREA
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk