[MITgcm-devel] dic

[Stephanie Dutkiewicz]

Martin -
You won't need the #ifndef READ_PAR

as insol.F is only called #ifndef READ_PAR
(see bio_export.F).

But otherwise your fixes/checks sound great.

steph

On Thu, 16 Nov 2006, Martin.Losch at awi.de wrote:

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