[MITgcm-support] clever pickup too clever for me

Matthew Mazloff mmazloff at MIT.EDU
Thu Sep 4 16:34:31 EDT 2008


Ahhh, JMC, i committed the ultimate support line sin of not first  
googling http://forge.csail.mit.edu/pipermail/mitgcm-support/ !!!   
sorry about that
well thanks for pointing it out :)

and now i have the pickup stuff sorted out...the initial fields look  
good (including initial seaice fields).

But the model is still blowing up...

Any idea...this is so wierd?  I'll dig a little deeper

-matt

ps> thanks for the link Dimitris

On Sep 4, 2008, at 6:15 AM, Jean-Michel Campin wrote:

> Hi Matt,
>
> Please have a look at :
> http://forge.csail.mit.edu/pipermail/mitgcm-support/2007-October/ 
> 005080.html
> Also, if you did not keep the meta files,
> there are few verification experiment with "minimal" pickup*.meta
> file that allow to restart from an old pickup.
> (and this is also explained in the email mentionned above).
>
> Jean-Michel
>
> On Wed, Sep 03, 2008 at 10:35:10PM -0700, Matthew Mazloff wrote:
>> Hi Patrick,
>>
>> Actually it seems that if you provide no meta file at all, which  
>> is what
>> I did, the model will try to read the pickup in a "new" way and  
>> crash.  I
>> think you need the old meta file and also need
>> pickupStrictlyMatch=.FALSE.,
>>
>> I think that with an an old meta file nbFields = 0 and is read in  
>> as how
>> use to be:
>>
>>           CALL READ_REC_3D_RL( fn, fp, Nr, uVel,  1, myIter,myThid )
>>           CALL READ_REC_3D_RL( fn, fp, Nr, guNm1, 2, myIter,myThid )
>>           CALL READ_REC_3D_RL( fn, fp, Nr, vVel,  3, myIter,myThid )
>>           CALL READ_REC_3D_RL( fn, fp, Nr, gvNm1, 4, myIter,myThid )
>>           CALL READ_REC_3D_RL( fn, fp, Nr, theta, 5, myIter,myThid )
>>           CALL READ_REC_3D_RL( fn, fp, Nr, gtNm1, 6, myIter,myThid )
>>           CALL READ_REC_3D_RL( fn, fp, Nr, salt,  7, myIter,myThid )
>>           CALL READ_REC_3D_RL( fn, fp, Nr, gsNm1, 8, myIter,myThid )
>>           CALL READ_REC_3D_RL( fn,fp,1, etaN, nj*Nr+1,  
>> myIter,myThid )
>>
>> if no meta file is presented, as I did, nbfields = -1, and reads in
>> pickup file in the new way (which goes uVel, vVel, theta, salt,  
>> guNm1,
>> gvNm1, gtNm1, gsNm1, etaN), which isn't consistent, and this  
>> causes the
>> model to crash.
>>
>> So I'll make up some meta files and throw them in there...(or I could
>> hardcode nbFields to 0)....and let you know if it works
>>
>> -Matt
>>
>> ps> of course, i don't know what's gonna happen with seaice,  
>> should be
>> fun :)
>>
>>
>>
>>
>>
>>
>>
>>
>> On Sep 3, 2008, at 9:56 PM, Patrick Heimbach wrote:
>>
>>> Hi Matt,
>>>
>>> what has changed is that you now
>>> also provide the .meta file, and its content is checked for
>>> consistency
>>> (might be important for seaice where it's not clear which options/
>>> grids etc you're using).
>>> Previously, the .meta files weren't evaluated.
>>> If you have an old pickup and didn't keep the .meta file,
>>> then just add in file data under PARM03 following option (as you  
>>> did)
>>>   pickupStrictlyMatch=.FALSE.,
>>> In that case the model warns you that the pickup might not be
>>> consistent,
>>> but continues running (assuming you know what you're doing).
>>> So ideally, keep the .meta files for your new runs and provide them
>>> when re-starting from pickup.
>>>
>>> -p.
>>>
>>>
>>>
>>> On Sep 3, 2008, at 11:47 PM, Matthew Mazloff wrote:
>>>
>>>> Hello JMC,
>>>>
>>>> I recently updated. It has been quite awhile since last update.
>>>> When trying to run the model it immediately blew up.  I strongly
>>>> suspect that it is the  "clever pickup" as the model blows up in 2
>>>> time steps.  I think fields are being improperly initialized.
>>>>
>>>> I don't have the pickup meta files so I use
>>>>  pickupStrictlyMatch=.FALSE.,
>>>> and I get the following warnings
>>>> (PID.TID 0001.0001) WARNING >> MDS_READ_META: file: pickup.
>>>> 0000000001.meta , pickup.0000000001.001.001.meta
>>>> (PID.TID 0001.0001) WARNING >> MDS_READ_META: Files DO not exist
>>>> (PID.TID 0001.0001) WARNING >> READ_PICKUP:  no field-list found
>>>> (PID.TID 0001.0001) WARNING >>  try to read pickup as currently
>>>> written
>>>> (PID.TID 0001.0001) WARNING >> MDS_READ_META: file: pickup_seaice.
>>>> 0000000001.meta , pickup_seaice.0000000001.001.001.meta
>>>> (PID.TID 0001.0001) WARNING >> MDS_READ_META: Files DO not exist
>>>> (PID.TID 0001.0001) WARNING >> SEAICE_READ_PICKUP:  no field-list
>>>> found
>>>> (PID.TID 0001.0001) WARNING >>  try to read pickup as currently
>>>> written
>>>>
>>>> So what has changed...I don't mind restructuring pickup.data or
>>>> making a meta file...but can you please let me know the best way to
>>>> do this.
>>>>
>>>> Or is it possibly not the clever pickup?
>>>>
>>>> -Thanks
>>>> Matt
>>>>
>>>>
>>>>
>>>> _______________________________________________
>>>> MITgcm-support mailing list
>>>> MITgcm-support at mitgcm.org
>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>
>>> ---
>>> Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
>>> MIT | EAPS 54-1518 | 77 Massachusetts Ave | Cambridge MA 02139 USA
>>> FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach
>>>
>>>
>>> _______________________________________________
>>> 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




More information about the MITgcm-support mailing list