[MITgcm-devel] balancing ctrl-obcs-flow

[Martin Losch]

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