Re: POP

NCL equivalent:

   path = "./b35.012dgvm_branch07.pop.r.2315-01-01-00000"
   nrec = 254 ; NCL is zero-based
   dims = (/320,384/)
   x = fbindirread(path, nrec, dims, "double")

On 7/22/10 4:46 PM, Michael Notaro wrote:
> After some experimentation, this seemed to work.
> Still need to plot it to be sure I read it correctly.
>
>
>
> ccc f77 -byteswapio binary_pop.f
> implicit none
> integer iostat, error
> real*8 IOBUFD(320,384)
> integer recl, i, j
>
> open(unit=100, status='old',
> # file='b35.012dgvm_branch07.pop.r.2315-01-01-00000',
> # form='unformatted', access='direct', recl=8*320*384,
> # iostat=iostat)
>
> read(unit=100, rec=255) IOBUFD
>
> write(*,*) IOBUFD
>
> close(100)
>
> end
>
>
>
>
> On Jul 22, 2010, at 5:24 PM, Jaromir Jakacki wrote:
>
>> It is direct access fortran file. Please look at the header of the restart file. There is beginning record number for each varible. It easy to understand it.
>> jaromir
>>
>> On 7/22/10 11:03 PM, Michael Notaro wrote:
>>> Does anyone have either an NCL or fortran code
>>> for reading the contents of a POP ocean restart
>>> file (binary) from CCSM3?
>>>
>>> Thanks, Michael
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jul 22 17:52:18 2010