[MITgcm-devel] balancing ctrl-obcs-flow

[Jean-Michel Campin]

Hi Martin and Matt,

There are some simple I/O S/R for reading vector (not the one used 
for tape): 
in pkg/mdsio: MDS_READVEC_LOC & MDS_WRITEVEC_LOC 
in pkg/rw: READ_GLVEC_RL,RS & WRITE_GLVEC_RL,RS 
 which call the previous one (but only support limited choice)
And I think they could be used to read in Nr*Nr*Nr array (not a tiled array)

Cheers,
Jean-Michel

On Wed, Mar 16, 2011 at 10:48:29AM +0100, Martin Losch wrote:
> Hi Matt,
> 
> I am not really qualified to answer your question, but I think that all I/O is somewhat taylored to the horizontal grid dimensions and you'll have organize the I/O of an array(Nr,Nr,Nr) yourself.
> You could try to copy some of the mdsio stuff about opening the file etc, but that actual read statement will be very different from the mdsio-way.
> 
> Martin
> 
> PS. try to avoid single (lonely, as in "not pairs") apostrophes ("'") in comments (e.g. as in "don't" -> "do not" ore "dont"), they break some scripts (for the html pages, as far as I remember).
> 
> 
> On Mar 15, 2011, at 6:09 PM, Matthew Mazloff wrote:
> 
> > Hi Martin,
> > 
> > Sounds good -- I can update my code by the end of the week.
> > 
> > I have one question though.  As the code is written the user must provide a mode file in ctrl_init_obcs_variables.F  The explanation is in the code.....   Can you recommend a way to read this NRxNRxNR file in.  It is currently coded:
> > 
> >         open(117, file='baro_invmodes.bin', access='direct',
> >      &     recl=2*WORDLENGTH*Nr*Nr, status='old')
> >         do j = 1,Nr
> >            read(117,rec=j) ((modesv(k,i,j), k=1,Nr), i=1,Nr)
> >         end do
> > 
> > but I am sure there is a much better way to do this, e.g. using the mdsio routines
> > 
> > Thanks
> > Matt
> > 
> > 
> > http://mitgcm.org/viewvc/MITgcm/MITgcm_contrib/SOSE/BoxAdj/code_ad/ctrl_init_obcs_variables.F
> > 
> > 
> > 
> > On Mar 15, 2011, at 9:45 AM, Martin Losch wrote:
> > 
> >> Matt,
> >> 
> >> this is now checked-in:
> >> 
> >>> Matt, when you find the time, please update your BoxAdj experiment (but you might want to wait until I have sorted out the "xx_fldperiod" stuff).
> >> 
> >> so you can update your BoxAdj experiment, if you like.
> >> 
> >> Martin
> >> 
> >> 
> >> _______________________________________________
> >> MITgcm-devel mailing list
> >> MITgcm-devel at mitgcm.org
> >> http://mitgcm.org/mailman/listinfo/mitgcm-devel
> > 
> > _______________________________________________
> > MITgcm-devel mailing list
> > MITgcm-devel at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-devel
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel