Re: (no subject)

Hi folks,

This is a heads up that Dave Brown fixed this problem, and it will be
in the next release of NCL (release date Nov/early Dec). If you would
like a fix now, email me and I'll get you a test version of NCL.

--Mary

On Wed, 25 Oct 2006, Hauss Reinbold wrote:

> Hello all,
>
> I am using NCL to byteswap some fotran77-generated binary files using the
> code shown at the bottom of this email. The files are about 4.5 gigs each.
> They work fine doing one record at a time until about half-way through when
> all the values turn to zeroes. Is this because the files are too big? Is
> there a 2Gb limit on binary file size like there is for netCDF file size? I
> have developed a workaround using a fortran subroutine that reads in the
> data, but I'm curious as to why my first approach didn't work. Thanks.
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>
> begin
>
> setfileoption("bin","ReadByteOrder","BigEndian")
> setfileoption("bin","WriteByteOrder","LittleEndian")
>
> if (isfilepresent("/flame1/SPI/spidata2003-01-sw.bin")) then
> print ("removing previous file")
> system ("rm /flame1/SPI/spidata2003-01-sw.bin")
> end if
>
> irec = 0
> do nyr = 1895, 2003
> do nmo = 1, 12
>
> f = fbindirread("/flame1/SPI/spidata2003.01",irec,(/
> 1405,621/),"float")
> f = f(:,620:0)
> print ("irec "+irec+" year "+nyr+" "+nmo+"F "+f(299,299))
> fbindirwrite("/flame1/SPI/spidata2003-01-sw.bin",f)
> irec = irec + 1
> delete(f)
> end do
> end do
> end
>
> **************************************************
> Hauss Reinbold
> Assistant Research Scientist / Programmer
> Climate Ecosystem and Fire Applications
> Desert Research Institute
> Office Phone: 775-673-7386
> email: Hauss.Reinbold_at_dri.edu
> **************************************************
>
>
>
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Nov 01 2006 - 11:35:54 MST