[MITgcm-support] External orbiting forces!

[Angela Zalucha]

Taimaz:

If you are using the fizhi package,
subroutine ASTRO in fizhi_clockstuff.F will calculate the solar zenith 
angle as a function of obliquity and eccentricity (these constants are set 
in getcon.F).  But note that they are treated as constants, so you would 
need to modify the code if they are changing with time (I assume that's 
what you mean by precession).


I also have a subroutine to calculate solar zenith angle as a function of 
obliquity, eccentricity, and ecliptic longitude as part of my planetary 
model (not part of the public distribution).  The equation for solar 
zenith angle may be found in e.g. Peixoto and Oort, "Physics of Climate" 
textbook.


Angela




_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
Angela Zalucha, Ph.D.
A. M. Zalucha Consulting

Office located at Southwest Research Institute
1050 Walnut Street, Suite 300
Boulder, CO 80302
USA
(720) 208-7211


http://www.boulder.swri.edu/~angela/
_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-

On Fri, 14 Sep 2012, taimaz.bahadory wrote:

> Hello everybody;
> Does MITgcm support orbital parameters, such as obliquity, precession, and eccentricity? I want to run a climate model
> for, say thousands of years, in which the orbital specifications of Earth will change indeed. Is MITgcm a suitable
> package for such a case?
> Thanks
> 
>