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

[Jean-Michel Campin]

Hi Roland,

I have the impression that there is some confusion here:
#define ATMOSPHERIC_LOADING (as well as pkg/exf) is specific to an ocean set-up 
and does not do much in an atmospheric set-up (which I though was the 
 context of your question).

And regarding how to set-up to top atmospheric pressure to tTop > 0:
 since pTop = Ro_SeaLevel + SUM( delR )_k=1,Nr
you can adjust Ro_SeaLevel or the level pressure-thickness delR to 
achieve pTop > 0.

It should work just fine with pure pressure-coord., but there are some parts
of the p*-coord code (I am thinking of geopotential calculation)
with select_rStar > 0 that assumes pTop = 0 
and would not work is pTop is larger or same magnitude as delR(Nr).

Cheers,
Jean-Michel

On Wed, Oct 28, 2015 at 11:10:39AM -0400, Dimitris Menemenlis wrote:
> 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
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support