[MITgcm-support] KPP

Wed Feb 28 23:24:11 EST 2007

hi all,

now i'd like to ask about Richardson number in \pkg\kpp\kpp_readparms.F, it is enable when KPP is enable too?? so i wrote #define ALLOW_KPP in CPP_OPTIONS.h, but i have an error message :
In ./CPP_OPTIONS.h there is an illegal line: #define ALLOW_KPP, what is that??i'd like to change value of Ri in kpp_readparms.F,what i have to do then??

thank you,
regards
rima

----- Original Message ----
From: "mitgcm-support-request at mitgcm.org" <mitgcm-support-request at mitgcm.org>
To: mitgcm-support at mitgcm.org
Sent: Saturday, February 24, 2007 12:00:10 AM
Subject: MITgcm-support Digest, Vol 44, Issue 17

Send MITgcm-support mailing list submissions to
    mitgcm-support at mitgcm.org

To subscribe or unsubscribe via the World Wide Web, visit
    http://mitgcm.org/mailman/listinfo/mitgcm-support
or, via email, send a message with subject or body 'help' to
    mitgcm-support-request at mitgcm.org

You can reach the person managing the list at
    mitgcm-support-owner at mitgcm.org

When replying, please edit your Subject line so it is more specific
than "Re: Contents of MITgcm-support digest..."

Today's Topics:

   1. nonhydrostatic (Riema Rachmayani)
   2. Re: nonhydrostatic (Martin Losch)

----------------------------------------------------------------------

Message: 1
Date: Fri, 23 Feb 2007 00:34:57 -0800 (PST)
From: Riema Rachmayani <imoth_22 at yahoo.com>
Subject: [MITgcm-support] nonhydrostatic
To: mitgcm-support at mitgcm.org
Message-ID: <20070223083457.65898.qmail at web53802.mail.yahoo.com>
Content-Type: text/plain; charset="iso-8859-1"

hi, all...

when we define our model for nonhydrostatic...we just describe :
#define ALLOW_NONHYDROSTATIC in code/CPP_OPTIONS.h and activated nonHydrostatic=.TRUE. in input/data and automatically MITgcm will compute the experiment by it self  specially for vertical velocity right??? just want to make it clear and my curiousity about it...

is theta' in rho=rho0(1-alfatheta*theta') ...eq (3.12) in chapter 3 is define as Tref in input/data???
Tref is reference potential temperature (theta0) or fluctuation potential temperature (theta')?? or total potential temperature (theta=theta0+theta') ?? or we use eq 3.13 to find Tref?? 

thank you,
regards
-rima-

Send instant messages to your online friends http://uk.messenger.yahoo.com 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://forge.csail.mit.edu/pipermail/mitgcm-support/attachments/20070223/ccfadb81/attachment-0001.html

------------------------------

Message: 2
Date: Fri, 23 Feb 2007 13:52:46 +0100
From: Martin Losch <Martin.Losch at awi.de>
Subject: Re: [MITgcm-support] nonhydrostatic
To: mitgcm-support at mitgcm.org
Message-ID: <91141BF8-406B-4F1E-AA1D-7D1A5B99D369 at awi.de>
Content-Type: text/plain; charset=US-ASCII; delsp=yes; format=flowed

Riema,

to turn on non-hydrostatic code you do exactly what you described.  
However, you can also set additional parameters for the cg3d solver  
(in analogy to cg2d), for example cg3dMaxIters, and  
cg3dTargetResidual, otherwise you'll use some defaults for those  
parameters (150 and 1e-7).

if (and only if) eosType = 'LINEAR', then
rho = rhoNil*(1-talpha*(theta-tRef) + sbeta*(salt-sRef))
default for rhoNil=999.8, but you can set that in in "data",  
talpha=2e-4,sbeta=7.4e-4 by default and you can set those parameters  
in "data" as well. tRef(k) and sRef(k) on the other hand, you have to  
set in "data". They are meant to be a reference temperature/salinity  
for the model level k. How you choose them is up to you.
You can also use other eosTypes, where the density is a nonllinear  
function of theta and salt and pressure (JMD95Z, JMD95P, MDJWF),  
based on observations.

Martin
On 23 Feb 2007, at 09:34, Riema Rachmayani wrote:

> hi, all...
>
> when we define our model for nonhydrostatic...we just describe :
> #define ALLOW_NONHYDROSTATIC in code/CPP_OPTIONS.h and activated  
> nonHydrostatic=.TRUE. in input/data and automatically MITgcm will  
> compute the experiment by it self  specially for vertical velocity  
> right??? just want to make it clear and my curiousity about it...
>
> is theta' in rho=rho0(1-alfatheta*theta') ...eq (3.12) in chapter 3  
> is define as Tref in input/data???
> Tref is reference potential temperature (theta0) or fluctuation  
> potential temperature (theta')?? or total potential temperature  
> (theta=theta0+theta') ?? or we use eq 3.13 to find Tref??
>
> thank you,
> regards
> -rima-
>
>
>
> Send instant messages to your online friends http:// 
> uk.messenger.yahoo.com
> _______________________________________________
> 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

End of MITgcm-support Digest, Vol 44, Issue 17
**********************************************

Send instant messages to your online friends http://uk.messenger.yahoo.com 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://forge.csail.mit.edu/pipermail/mitgcm-support/attachments/20070226/69d10d47/attachment-0001.html

------------------------------

_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org
http://mitgcm.org/mailman/listinfo/mitgcm-support

End of MITgcm-support Digest, Vol 44, Issue 18
**********************************************

Send instant messages to your online friends http://uk.messenger.yahoo.com 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20070228/f2e9ce9d/attachment.htm>