[MITgcm-support] interpickups.m

[Baylor Fox-Kemper]

Hi Nicholas,
   It looks like rather than interpolating, it just copied to the same 
indices.  Note that interpickups tries to match x values and y values.  
That is, it reads xc, xf, etc. from the first set of pickups and from 
the pickouts, and then it interpolates from one to the other.
   Check to make sure that the x and y values cover the same range in 
pickups and pickouts.  For example, if you didn't change dxspacing and 
dyspacing when you changed resolution, the xc range wouldn't match at 
all.
   You could check quickly in matlab by making the same plots, but 
labeling the axes with xc and yc instead of just index.
   Good luck,
     -Baylor


On Sep 16, 2005, at 3:59 PM, Nicolas Wienders wrote:

>
> On Sep 16, 2005, at 2:53 PM, Baylor Fox-Kemper wrote:
>
>> Hi Nicholas,
>>   I wrote interpickups a while ago, and there have been a few changes 
>> to mnc since (sorry!).  But, I think that the basics of the program 
>> should still work.
>>   I intended that the pickout files be generated by the MITgcm.  I 
>> intended the following standard use:
>>
>> 1) Run MITgcm at coarse resolution, producing pickups.
>> 2) Run MITgcm at fine resolution for one timestep, with arbitrary 
>> initial condtions, but with the right new grid and tiling and 
>> processors, etc.  Get it to produce one set of pickups at that 
>> resolution.
>> 3) Use the pickups from 2) as pickout and use the pickups from 1) as 
>> pickin.
>>
>
>> Hi Baylor,
> Thanks a lot for the trick.
> I did try and here are the results  (field, and interpolated field).
> ....
> i wonder what i did wrong...
>
>
> <field.png><field_interp.png>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support