[MITgcm-devel] balancing ctrl-obcs-flow

[Matthew Mazloff]

Hello,

I updated the ALLOW_OBCS_CONTROL_MODES code and am testing it now.

I have yet to get to revising the I/O yet, but its on my to-do list.

-Matt

ps> Martin, I like the new ctrl_getobcs?.f code -- much cleaner!



On Mar 16, 2011, at 6:25 AM, Jean-Michel Campin wrote:

> 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
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel