[MITgcm-support] R/W pickup

[Martin Losch]

Hi Chen

you can just use rdmds.m (and rdmeta.m) in utils/matlab
p = rdmds(‘pickup’,yourIter);
will give you a 3D field (nx,ny,nn)
then you’ll have to determine, which of the nn levels belong to which field (see pickup*.meta for a list of files in pickup*.data)
after modifying p, you can write it like this:
fid=fopen(‘pickup.withNewName.data’,’w’,’b’);
fwrite(fid,p,’real*8’);
fclose(fid);

Martin


> On 06 Jan 2016, at 14:59, Chen Cheng <cc2713206 at 163.com> wrote:
> 
> Dear all,
>  
> I am trying to set up a regional coupled ice shelf-ocean-sea ice experiment in the Prydz Bay with MITgcm. However, the model has blown up after nearly 3-years integration. Right now, I want to make some modifications on the pickup files to continue running. If possible, could you provide the relevant executable source codes encapsulated in m-files for Matlab to read or write the pickup files.
>  
> Thanks, and happy new year.
>  
> Chen Cheng
> 
> 
>  
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support