[MITgcm-devel] global_with_exf fails

Jean-Michel Campin jmc at ocean.mit.edu
Sat Jan 14 11:05:38 EST 2006


Hi Patrick,

The modifications you made in exf_interp.F change the 
results of the global_with_exf test (it fails at level 2).
Was it supposed to be like this (and then we need to update 
the output.txt), or is there something wrong ?

Otherwise, i think there is a potential problem
in mixing different type of argument (r8 & r4) in this call:
line 91           xg(i,j,bi,bj) = mod(xg(i,j,bi,bj),360.)
and in this one also:
line114         x_in(i) = mod(x_in(i),360.)
(xg is _RS and x_in is real*8)
and it would be safer to use the same type of argument (e.g. with
a _RS constant set to 360.  in the 1rst case and an other real*8
constant set to 360. in the 2nd case).

See you,

Jean-Michel





More information about the MITgcm-devel mailing list