[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