[MITgcm-support] Re: mdsio slice routines

Tue May 27 15:02:39 EDT 2008

Hi Matt,
Thanks for your reply. The error message of the optimization run (only
with southern or western boundary) is
--------------------------------------------------------------------------------------------------------
forrtl: severe (67): input statement requires too much data, unit 20,
file /tmp/mitgcm_so_suneet/JHU_ctrl _MIT_CE_000.opt0000
--------------------------------------------------------------------------------------------------------

The interesting points to notice are that:

(i) when I don't use mdsio slice routines and copy obcs control vector
files from different nodes manually on the master node, my
optimization run (optim.x) works fine for all the boundaries without
any error, whereas, on using those mdsio routines I get abovesaid
error message. Now, the problem here is that I can not manually copy
obcs control vector files every time while running assimilation
experiment in which model is run several times in a loop.

(ii) On using mdsio slice routines for northern and eastern boundary,
my optimization run (optim.x) again works fine but not with southern
and western boundary.

Therefore, personally, I think there is some bug in the mdsio routines
because the present problem arises as soon as I try to use them for
southern/western boundaries?

Though I tried hard but I am not able to figure out the likely bug or
inconsistency either in linesearch or mdsio slice routines.

Please help me sort out this problem.
Regards,
Suneet

On Tue, May 27, 2008 at 2:01 PM, Matthew Mazloff <mmazloff at mit.edu> wrote:
> Hi Suneet,
>
> I have never checked.  I only use northern boundary.   I do not think there
> is a problem with your model set-up because you said the model runs.  Its
> the linesearch, optim.x, that doesn't work, correct?  So there is likely a
> bug or inconsistency in there somewhere...what error does the optimization
> give you?
>
> -Matt
>
>
> On May 27, 2008, at 1:36 PM, Suneet Dwivedi wrote:
>
>> Hi Matt,
>>
>> May I request you to please let me know whether your mdsio slice
>> routines work for you (for your model setup) for southern or western
>> boundary? I am not sure whether my current problem with optimization
>> is associated with mdsio slice routines or there is any other problem
>> with my model setup.
>>
>> Suneet
>>
>> On Fri, May 23, 2008 at 4:06 PM, Suneet Dwivedi <suneetdwivedi at gmail.com>
>> wrote:
>>>
>>> Hi Matt,
>>> Just to let you know that when I don't use those mdsio slice routines
>>> and copy obcs control vector files from different nodes manually on
>>> the master node, my optimization run (optim.x) works fine for all the
>>> boundaries. But, this is not a practical solution for running
>>> assimiation experiment. So, I am sure that there is some discrepancy
>>> between call routines for southern and western boundaries.
>>> Cheers,
>>> Suneet
>>>
>>> On Fri, May 23, 2008 at 2:28 PM, Suneet Dwivedi <suneetdwivedi at gmail.com>
>>> wrote:
>>>>
>>>> yes everytime when i make any modification in code (including change
>>>> in ECCO_CPPOPTIONS.h) , I always recompile optim.x.
>>>> Thanks,
>>>> Suneet
>>>>
>>>> On Fri, May 23, 2008 at 2:08 PM, Matthew Mazloff <mmazloff at mit.edu>
>>>> wrote:
>>>>>
>>>>> just checking, but did you recompile optim.x when you added these
>>>>> controls
>>>>> (i.e. changed ECCO_CPPOPTIONS.h)?
>>>>>
>>>>>
>>>>> On May 23, 2008, at 1:49 PM, Suneet Dwivedi wrote:
>>>>>
>>>>>> Hi Matt,
>>>>>> Thanks for your reply. My do not 'work' means that the
>>>>>> assimilation/optimization run fails when I try to allow
>>>>>> southern/western open boundary control. So the model runs absolutely
>>>>>> fine (as before) but optimization (optim.x) fails with south/west OB
>>>>>> control.
>>>>>> Cheers,
>>>>>> Suneet
>>>>>>
>>>>>> On Fri, May 23, 2008 at 12:16 PM, Matthew Mazloff <mmazloff at mit.edu>
>>>>>> wrote:
>>>>>>>
>>>>>>> Hi Suneet,
>>>>>>>
>>>>>>> Northern and southern boundaries should both use same I/O
>>>>>>> routines.....and
>>>>>>> eastern and western should both use same I/O routines.  If the
>>>>>>> northern
>>>>>>> boundary works but the southern doesnt than than there must be a
>>>>>>> discrepancy
>>>>>>> between call routines....that should not be the case....I had a quick
>>>>>>> look
>>>>>>> for this discrepancy but couldnt find it....not sure whats going on.
>>>>>>>  What
>>>>>>> do you mean the routines do not 'work' for the south and west
>>>>>>> boundaries?
>>>>>>>
>>>>>>> -Matt
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> On May 22, 2008, at 5:18 PM, Suneet Dwivedi wrote:
>>>>>>>
>>>>>>>> Hi Matt,
>>>>>>>> Thanks for your reply. But, what is the difference between northern
>>>>>>>> and southern boundary and similarly between eastern and western
>>>>>>>> boundary?
>>>>>>>> As per your advice, I modified MDSREADFIELDYZ and MDSWRITEFIELDYZ.
>>>>>>>> Now
>>>>>>>> the routines work for eastern boundary in addition to the northern
>>>>>>>> boundary but still not for south and west.
>>>>>>>> Waiting for your response,
>>>>>>>> Cheers,
>>>>>>>> Suneet
>>>>>>>>
>>>>>>>> On Thu, May 22, 2008 at 3:04 PM, Matthew Mazloff <mmazloff at mit.edu>
>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>> Hi Suneet,
>>>>>>>>>
>>>>>>>>> The northern (southern) boundary is xz.  I think you just need to
>>>>>>>>> apply
>>>>>>>>> these changes to the eastern (western) boundary = yz.  so for
>>>>>>>>> example
>>>>>>>>> just
>>>>>>>>> change (appropriately) MDSREADFIELDYZ to look like MDSREADFIELDXZ.
>>>>>>>>>
>>>>>>>>> -Matt
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On May 22, 2008, at 2:39 PM, Suneet Dwivedi wrote:
>>>>>>>>>
>>>>>>>>>> Hi Patrick and Matt,
>>>>>>>>>> I wish to modify mdsio slice routines (mdsio_gl_slice.F,
>>>>>>>>>> mdsio_slice.F, mdsio_slice_loc.F, mdsio_readvector.F,
>>>>>>>>>> mdsio_writevector.F) sent to me by Matt so as to make them work
>>>>>>>>>> for
>>>>>>>>>> other boundaries in addition to the northern boundary. Matt told
>>>>>>>>>> me
>>>>>>>>>> that these routines are written only for the northern boundary. I
>>>>>>>>>> am
>>>>>>>>>> not sure which variables in the routines denote the northern
>>>>>>>>>> boundary.
>>>>>>>>>> This information is necessary for me to extend the routines for
>>>>>>>>>> any
>>>>>>>>>> other boundary. May I request you to please tell me the same.
>>>>>>>>>> Cheers,
>>>>>>>>>> Suneet
>>>>>>>>>
>>>>>>>>>
>>>>>>>
>>>>>>> _______________________________________________
>>>>>>> MITgcm-support mailing list
>>>>>>> MITgcm-support at mitgcm.org
>>>>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>>>>>
>>>>>> _______________________________________________
>>>>>> MITgcm-support mailing list
>>>>>> MITgcm-support at mitgcm.org
>>>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>>>
>>>>> _______________________________________________
>>>>> MITgcm-support mailing list
>>>>> MITgcm-support at mitgcm.org
>>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>>>
>>>>
>>>
>
>