[MITgcm-support] RE: a few MITgcm questions

Chris Hill cnh at mit.edu
Tue Apr 6 02:16:03 EDT 2004


Hi Ichiro,

Yes. 510*6,510,50 in a single record is what it would be for a 50 level
case. 

You can double check by looking at the size of the file in bytes from ls.
The size should be 510*6*510*50*4 for 32-bit real and 510*6*510*50*8 for
64-bit.

Chris

> -----Original Message-----
> From: Ichiro Fukumori [mailto:fukumori at jpl.nasa.gov] 
> Sent: Monday, April 05, 2004 1:54 PM
> To: cnh at MIT.EDU
> Cc: menemenlis at jpl.nasa.gov; bnc at pacific.jpl.nasa.gov
> Subject: RE: a few MITgcm questions
> 
> What is the format of the new face-based cubed sphere files?
> Is it a straight single record of a (510*6,510,50) array?
> 
> Thanks,
> Ichiro
> 
> 
> At 08:40 PM 4/3/2004 -0500, Chris Hill wrote:
> >Ichiro,
> >
> >  -r8 should not be needed, in fact g77 doesn't support it, 
> so the core 
> >code has to work ok without it. However, it may not have 
> been tested on 
> >the Altix lately. Because of g77, for the kernel code it 
> never matters, 
> >but there could be some stuff in KPP where its assumed 
> (Dimitris any thoughts?).
> >
> >  For reading files could you for now do
> >
> >   real*4           KFtheta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
> >   _RS              KFtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
> >   character*160 KFthetaFnam
> >   CALL MDSREADFIELD(KFthetaFnam,precFloat32,'RS',Nr,KFtmp, 
> 1, myThid)
> >   Kftheta = KFtmp
> >
> >   that should read a 510*6,510 file for your setup that is 
> a straight 
> >single record direct access form with 32-bit reals.
> >  Depending on
> >
> >  1 - compiler options
> >  2 - settings of -D_BYTESWAPIO
> >
> >  the input will need to be either little-endian or big-endian
> >
> >  The build options that Dimitri has for the Altix and the 
> Origin both 
> >are set for big-endian (that's the native format on MIPS but 
> not on Itanium).
> >
> >It would probably be good to have a 'R4' and 'R8' MDS arg 
> but I have to 
> >head off to Ames soon so I won't have time - I'll be in touch Monday.
> >
> >Chris
> >
> >
> >
> >
> >
> >
> > > -----Original Message-----
> > > From: Ichiro Fukumori [mailto:fukumori at jpl.nasa.gov]
> > > Sent: Saturday, April 03, 2004 6:58 PM
> > > To: cnh at MIT.EDU
> > > Cc: bnc at pacific.jpl.nasa.gov; Dimitris Menemenlis; 
> > > chenze at nas.nasa.gov
> > > Subject: a few MITgcm questions
> > >
> > > Hi Chris,
> > >
> > > Is compile option -r8 necessary?  It appears to interfere 
> with using 
> > > sgemv (a single precision BLAS matrix times vector subroutine).
> > >
> > > I have a 3D array on the cubed sphere grid that I need to read in.
> > > Could you tell me the format I should write it to and the 
> subroutine 
> > > call to read it.  Presently I have,
> > >   real*4           KFtheta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
> > >   character*160 KFthetaFnam
> > >   CALL MDSREADFIELD(KFthetaFnam,precFloat32,'RS',Nr,KFtheta,
> > > 1, myThid) Does MDSREADFIELD work with real*4 array (as 
> opposed to 
> > > _RS)?
> > >
> > > Thanks,
> > > Ichiro
> > >
> > >
> > > ===============================================================
> > > Ichiro Fukumori               | Jet Propulsion Laboratory
> > > e-mail: fukumori at jpl.nasa.gov | Mail Stop 300-323
> > > phone: +1 (818) 354-6965      | 4800 Oak Grove Drive
> > > fax: +1 (818) 393-6720        | Pasadena, California 91109, U.S.A.
> > > ===============================================================
> > >
> 
> ===============================================================
> Ichiro Fukumori               | Jet Propulsion Laboratory
> e-mail: fukumori at jpl.nasa.gov | Mail Stop 300-323
> phone: +1 (818) 354-6965      | 4800 Oak Grove Drive
> fax: +1 (818) 393-6720        | Pasadena, California 91109, U.S.A.
> ===============================================================
> 




More information about the MITgcm-support mailing list