[MITgcm-devel] dic

[Stephanie Dutkiewicz]

Martin -

You can certainly put in the
  IF ( .NOT. usingSphericalPolarGrid ) change to insol.F

However, there is already a provision to read in the light
data instead of using insol.F - which would probably be the
best way to give the light for cubic sphere etc.
In GCHEM_OPTIONS there is a READ_PAR command - if this
is true, then light is read from from a file given as
"Filename1" in data.gchem.

I guess there should be a warning then for any other type
of grid to say: "USE READ_PAR=.TRUE. not insol.F for this
type of grid"

steph

On Wed, 15 Nov 2006, Martin.Losch at awi.de wrote:

> 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
>
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>