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

Wed Oct 28 15:39:04 EDT 2015

oops, sorry, I read/replied too quickly.  please ignore my answer.  it was for ocean.

> On Oct 28, 2015, at 3:32 PM, Jean-Michel Campin <jmc at ocean.mit.edu> wrote:
> 
> 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
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support