[MITgcm-support] KPP

Thu Mar 1 03:31:20 EST 2007

Hi Riema,
this one fall into the category RTFM:
In order to enable packages, you have to edit packages.conf and NOT  
CPP_OPTIONS.h
KPP is  special in that you have to set #define SHORTWAVE_HEATING in  
CPP_OPTIONS.h in addtion. Check kpp_readparms.F for the namelist of  
parameters you can set at runtime. the Richardson (buoyancy divided  
vertical velocity shear) number is computed in the KPP package. I  
don't see, why you want to change that. I am guessing that you mean  
Ricr (the "critical" Richardson number), which is part of the  
namelist defined in kpp_readparms.F

Hope that helps,
Martin

On 1 Mar 2007, at 05:24, Riema Rachmayani wrote:

> 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
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support