[MITgcm-support] How to add a new equation of state in mitgcm

[赖燕红]

Hi Martin,

Thanks very much for your detailed instructions on how to add a new EOS or a parameter in MITgcm. I have tried by myself in mitgcm and your instructions really help a lot. Thanks again!

Yanhong


> -----原始邮件-----
> 发件人: mitgcm-support-request at mitgcm.org
> 发送时间: 2021-12-18 01:00:05 (星期六)
> 收件人: mitgcm-support at mitgcm.org
> 抄送: 
> 主题: MITgcm-support Digest, Vol 222, Issue 10
> 
> Send MITgcm-support mailing list submissions to
> 	mitgcm-support at mitgcm.org
> 
> To subscribe or unsubscribe via the World Wide Web, visit
> 	http://mailman.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. How to add a new equation of state in mitgcm (???)
>    2. Re: How to add a new equation of state in mitgcm (Martin Losch)
> 
> 
> ----------------------------------------------------------------------
> 
> Message: 1
> Date: Fri, 17 Dec 2021 17:57:35 +0800 (GMT+08:00)
> From: ??? <yhlai at pku.edu.cn>
> To: mitgcm-support at mitgcm.org
> Subject: [MITgcm-support] How to add a new equation of state in mitgcm
> Message-ID: <75210593.80e02.17dc7d35583.Coremail.yhlai at pku.edu.cn>
> Content-Type: text/plain; charset="utf-8"
> 
> Hi, 
> 
> I want to add a new equation of state in MITgcm. Particularly, this equation of state is assumed to be a piecewise function, that means the density has a different dependence on pressure and temperature under different temperatures. I found some files in /models/ that might relate to the settings of the equation of state, for example, /models/src/ini_eos.F & Find_rho.F & find_alpha.F, /models/inc/EOS.h & PARMS.h. But I'm still very confused on how to realize a piecewise equation of state that depends on temperature. Does anyone know or ever try how to realize this? 
> 
> 
> 
> 
> Also, if I want to add a new parameter in /input/data PARM01, where do I need to define this new parameter? Or do I need to modify other files at the same time?
> 
> 
> 
> 
> Thanks in advance!
> 
> 
> 
> 
> Yanhong 
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20211217/1bb72e23/attachment-0001.html>
> 
> ------------------------------
> 
> Message: 2
> Date: Fri, 17 Dec 2021 15:48:59 +0100
> From: Martin Losch <Martin.Losch at awi.de>
> To: MITgcm Support <mitgcm-support at mitgcm.org>
> Subject: Re: [MITgcm-support] How to add a new equation of state in
> 	mitgcm
> Message-ID: <27040E38-2884-42AE-B299-58C1EA1BFD28 at awi.de>
> Content-Type: text/plain; charset="utf-8"
> 
> Yanhong,
> 
> if you want to add your own eos, then the right place to start is ?model/src/find_rho.F?, where density (actually rho-rhoConst) is computed for level ?k?. The routine computes pressure ?locPres? as a fucntion of depth (level k). and with that you should have all ingredients to compute your density.
> For some mixed layer models, you need alpha and beta (thermal expansion and haline contraction coefficients), which are computed in find_alpha.F and find_beta.F
> 
> If you want to add a new parameter, you need to define it in model/inc/PARAMS.h (plenty of example), don?t forget to put it into the appropriate common block, set the default in model/src/set_defaults.F and add it to the namelist in model/src/ini_parms.F
> For your own safety it makes sense to print its value in model/src/config_summary.F
> 
> Martin
> 
> > On 17. Dec 2021, at 10:57, ??? <yhlai at pku.edu.cn> wrote:
> > 
> > Hi, 
> > 
> > I want to add a new equation of state in MITgcm. Particularly, this equation of state is assumed to be a piecewise function, that means the density has a different dependence on pressure and temperature under different temperatures. I found some files in /models/ that might relate to the settings of the equation of state, for example, /models/src/ini_eos.F & Find_rho.F & find_alpha.F, /models/inc/EOS.h & PARMS.h. But I'm still very confused on how to realize a piecewise equation of state that depends on temperature. Does anyone know or ever try how to realize this? 
> > 
> > 
> > 
> > Also, if I want to add a new parameter in /input/data PARM01, where do I need to define this new parameter? Or do I need to modify other files at the same time?
> > 
> > 
> > 
> > Thanks in advance!
> > 
> > 
> > 
> > Yanhong 
> > 
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> 
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: smime.p7s
> Type: application/pkcs7-signature
> Size: 1665 bytes
> Desc: not available
> URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20211217/e655901b/attachment-0001.p7s>
> 
> ------------------------------
> 
> Subject: Digest Footer
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> 
> 
> ------------------------------
> 
> End of MITgcm-support Digest, Vol 222, Issue 10
> ***********************************************