(no subject)

From: Hauss Reinbold <Hauss.Reinbold_at_nyahnyahspammersnyahnyah>
Date: Wed, 25 Oct 2006 15:31:36 -0700

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"



    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,(/
         f = f(:,620:0)
         print ("irec "+irec+" year "+nyr+" "+nmo+"F "+f(299,299))
         irec = irec + 1
       end do
    end do

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
Received on Wed Oct 25 2006 - 16:31:36 MDT

This archive was generated by hypermail 2.2.0 : Wed Nov 01 2006 - 11:37:05 MST