[MITgcm-support] EXF period with interpolation

[Jean-Michel Campin]

Hi Renske,

This "STOP" in exf_set_uv.F has been added not too long ago (June 2014)
so I suspect it is necessary (does catch a real problem). 

The issue has to do with the fact that the vector version (EXF_INIT_UV) that 
would do the initial interpolation of a 2 compoinent vector field is missing (not yet
written), whereas, for a scalar field, we have an initialisation routine (EXF_INIT_GEN) 
in addition to the in-time-step loop routine EXF_SET_GEN (both in same source file: 
exf_set_gen.F).

There might be a simple way to go arround this but would need to check carefully.

In the mean time, if you are not in this case: usingCurvilinearGrid .OR. rotateGrid
you could try to set:
      uvInterp_stress = .FALSE.
      uvInterp_wind   = .FALSE.
      uvInterp_climstr= .FALSE.
by commenting out lines 517 & 521 in exf_init_fixed.F:
c # ifdef EXF_USE_OLD_INTERP_POLE
....
c # endif

Cheers,
Jean-Michel

On Mon, Oct 12, 2015 at 02:29:24PM +0000, Renske Gelderloos wrote:
> Hi all,
> 
> 
> I'm trying to use the EXF package with on-the-fly interpolation and (initially) non-time-varying forcing. If I have understood this correctly, I set repeatPeriod to 0 and that should work. However, there is a safe guard built into exf_set_uv.F, that if USE_EXF_INTERPOLATION is set and uvecperiod .LE. 0 the model stops. Now, in my case uvecperiod is 0 because I have time-independent forcing. Did I misinterpret something or should this safe guard be .LT. and not .LE.?
> 
> 
> Renske

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support