[MITgcm-devel] dic

[Martin.Losch at awi.de]

Steph,

Hezi and I are trying to use pkg/dic for the Gulf of Eilat. For a long time we couldn't get it to work because of grid scale 
noise. However, now we found the problem: it's the insolation (insol.F), which is hardwired to latitude=YC, but because 
we are using a cartesian grid YC is in meters and sin/cos(YC*pi/180)  produces garbage. Do you want that fixed in the 
repository. I could do something like this:
c latitude in radians
          lat=YC(1,j,1,bj)/180.d0*3.1416 ! your original code
C     latitute in radians, backed out from coriolis parameter
C     (makes latitude independent of grid)
          IF ( .NOT. usingSphericalPolarGrid )
     &         lat = asin( fCori(1,j,1,bj)/(2. _d 0 *omega) )*deg2rad

What do you think? Because you do this insol thing only for varying J, it will still not work for the cubes sphere or rotated 
grids, but at least for cartesian grids with the y-axis in the north south direction. Or do it right and go for a 2d insolation 
field.

Martin


Martin Losch
Alfred Wegener Institute 
Postfach 120161, 27515 Bremerhaven, Germany; 
Tel./Fax: ++49(0471)4831-1872/1797