[MITgcm-support] exf interp

[Dimitris Menemenlis]

Matt, could you send an example figure of problem?

On Oct 17, 2009, at 7:29 AM, Matthew Mazloff wrote:

> Hi Dimitris,
>
> Ah yes, I see this now.  Very nice.  Unfortunately the gridding is
> still noticeable in my plots of wind speed gradient.  I am able to
> pick out the exact locals where the forcing is prescribed -- meaning
> the interpolation is not smooth.  This is not good for plotting wind-
> stress curl. Have you noticed this on your high-res set-ups?  Do you
> think this is a problem?  Should we (can we) try a smoother interp
> method?
>
> -Matt
>
> On Oct 16, 2009, at 5:03 PM, Dimitris Menemenlis wrote:
>
>> Matt, I hardcoded bilinear interpolation for tracer fields because
>> overshoots can be problematic and bicubic for wind velocity (see
>> exf_set_uv.F) for stress fields because the second derivative
>> matters.  I would recommend to leave as is.  Other combinations can
>> (or did) cause trouble.  Dimitris
>>
>> On Oct 16, 2009, at 3:21 PM, Matthew Mazloff wrote:
>>
>>> Hello,
>>>
>>> It appears the model is equipped to do both bilinear and bicubic
>>> interpolation for exf_interp.F.  Interp_method, however, is  
>>> hardcoded
>>> to bilinear interpolation.  This does matter for my 1/16 degree set-
>>> up.  Can anyone confirm that the bicubic interpolation is ok --  
>>> and I
>>> can go ahead and use this option.
>>>
>>> -Thanks,
>>> Matt