[MITgcm-support] p /= 0 at top of atmosphere model

[Dimitris Menemenlis]

compile with:

C o Include pressure loading code
#define ATMOSPHERIC_LOADING

in CPP_OPTIONS.h

then read-in atmospheric pressure using, e.g., pkg/exf
for an example look at:
http://mitgcm.org/viewvc/MITgcm/MITgcm_contrib/llc_hires/llc_4320/input/data.exf

> On Oct 28, 2015, at 9:15 AM, Roland Young <Roland.Young at physics.ox.ac.uk> wrote:
> 
> I was wondering - is it possible (and/or straightforward) to set the pressure at the upper boundary of the atmosphere model to p /= 0? Eq. 1.54 in the documentation states p=0 at the top, but I have an unphysical temperature inversion near the top of my model, which I think I may be able to get rid of if the upper boundary pressure is not zero, so I wondered if it is possible.
> 
> Roland
> 
> 
> 
> 
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support