[MITgcm-support] Bug in exf_interp routine?

[Jean-Michel Campin]

Hi Nora,

Thanks for pointing to this.
I think you are right regarding this line of code (South pole treatment) and 
it should be changed to minic what is there for the N.pole.
Will fix this later today.

The S.pole treatment has probably not be carefully checked since S.pole is 
generally far from any ocean, whereas I remember doing some tests for
near N.pole interpolation.

Cheers,
Jean-Michel

On Mon, Sep 19, 2016 at 11:26:01AM +0000, Nora Loose wrote:
> Hi all,
> 
> 
> I have a question regarding the exf_interp routine, probably easiest to answer for the code developers.
> 
> Is it possible that there is a small bug in line 151 of exf_interp.F (where 2 rows are added at the southern end of the input grid and care is taken of their latitude assignments)? I am using version c65u, and the exf_interp routine has not been updated since then (as far as I can see).
> 
> The code reads as
> 
> |y_in(0)| > yPole and |y_in(-1)| > yPole,
> 
> where I think it should be
> 
> |y_in(0)| < yPole and |y_in(-1)| > yPole.
> 
> I might be wrong??? but when applying the same procedure at the north pole (line 169), something equivalent to the suggested version is done.
> 
> 
> Thanks for your help,
> 
> Nora
> 
> ---------------------------------------------
> Nora Loose, PhD research fellow
> Department of Earth Science
> Bjerknes Centre for Climate Research
> 
> University of Bergen
> Allegaten 70, 5007 Bergen, Norway
> 
> Room 104
> Tel +47 55 58 32 54
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support