[MITgcm-support] clever pickup too clever for me

[Dimitris Menemenlis]

Matt, if it's helpful, in
http://ecco2.jpl.nasa.gov/data1/cube/matlab/mk_pickup78.m
there are short descriptions of some of the more recent and common  
pickup formats
and code to go from older to newer formats:

% pickup before cube76 plus cube77 had 403 horizontal levels
%   1:50  uVel
%  51:100 guNm1
% 101:150 vVel
% 151:200 gvNm1
% 201:250 theta
% 251:300 gtNm1
% 301:350 salt
% 351:400 gsNm1
%     401 etaN
%     402 dEtaHdt
%     403 etaH
% pickup for cube76 and for cube78 and later have 303 levels
%   1:50  uVel
%  51:100 vVel
% 101:150 theta
% 151:200 salt
% 201:250 guNm1
% 251:300 gvNm1
%     301 etaN
%     302 dEtaHdt
%     303 etaH
pickup75to76mapping=[1:50 101:150 201:250 301:350 51:100 151:200 401  
402 403];

% pickup_seaice before cube66 plus cube67 to cube69 had 15 levels
%    1  TICES
%    2  YNEG
%    3  HSNOW
%  4:6  UICE
%  7:9  VICE
% 10:12 HEFF
% 13:15 AREA
% pickup_seaice for cube66 plus cube70 to cube75 had 21 levels
%  1:7  TICES
%    8  YNEG
%    9  HSNOW
% 10:12 UICE
% 13:15 VICE
% 16:18 HEFF
% 19:21 AREA
% pickup_seaice for cube76 to cube78 have 22 horizontal levels
%  1:7  TICES
%    8  YNEG
%    9  HSNOW
% 10:12 UICE
% 13:15 VICE
% 16:18 HEFF
% 19:21 AREA
%    22 HSALT
% pickup_seaice for cube78 and later have 13 horizontal levels
%  1:7  TICES
%    8  AREA
%    9  HEFF
%    10 HSNOW
%    11 HSALT
%    12 UICE
%    13 VICE
pickup65to76mapping=[1 1 1 1 1 1 1 2:15 0];
pickup75to76mapping=[1:21 0];






Dimitris Menemenlis
DMenemenlis at gmail.com

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