[MITgcm-devel] dic

[Martin.Losch at awi.de]

Hi Steph,
you are right, once one decides to read the radiation, there shouldn't be any problem. However, it's nice to have radiation 
from astronomical parameters. I will check in a fix to insol.F like this
#ifndef READ_PAR
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
#endif /* READ_PAR */

and also add something like this in insol.F (it should really be in a routine such as dic_check.F but we don't have that, 
yet).
#ifndef READ_PAR
C     latitute in radians, backed out from coriolis parameter
C     (makes latitude independent of grid)
      IF ( usingCurvilinearGrid ) THEN
       STOP 'S/R INSOL: use #define READ_PAR with usingCurvilinearGrid'
      ENDIF
#endif /* READ_PAR */

OK?

Martin


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



----- Original Message -----
From: Stephanie Dutkiewicz <stephd at ocean.mit.edu>
Date: Wednesday, November 15, 2006 6:24 pm
Subject: Re: [MITgcm-devel] dic

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