<div>Hi Yang,</div><div> </div><div>So far as I know, GM parameterization is an adiabatic scheme, which means that this parameterization do not change water properties (e.g., density).  You can read their paper in 1995 which showed that GM scheme is effectively an advection rather than diffuion.  So neither GM90 nor GM_VisbK (Visbeck et al. 1997) is diapycnal.</div><div> </div><div>In MITgcm, diagnosing diapycnal diffusion is not straightforward as it is a z-coordinate model.  But you can use the resorting method by Griffies et al. (2000) to diagnose the 'effective' diapycnal diffusivity.  Effective means those processes act to homogenize water density, such as the explicit diapycnal diffusivity you set in the model (diffKr), spurious numerical diffusivity associated with advection scheme, those caused by other parameterizations (e.g., KPP, tidal mixing) etc.</div><div> </div><div>You can refer to Griffies's paper for details.  Hope this helps.</div><div> </div><div>References:</div><div>Griffies, S. M., R. C. Pacanowski, and R. W. Hallberg, 2000: Spurious Diapycnal Mixing Associated with Advection in a z-Coordinate Ocean Model. Mon. Wea. Rev., 128, 538-564.</div><div>Visbeck, M., J. Marshall, and T. Haine, 1997: Specification of Eddy Transfer Coefficients in Coarse-Resolution Ocean Circulation Models. J. Phys. Oceanogr., 27, 381-402.</div><sign signid="0"><div> </div><div> </div><div><br>------------------</div><div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px;"><div> <font size="2">Best regards <br> <br>Yu-Kun Qian (钱钰坤) <br>Center for Monsoon and Environment Research</font> </div> <font size="2">Department of Atmospheric Sciences<br>School of Environmental Science and Engineering<font size="2"> <br> </font>Sun Yat-sen University <br>No. 135 Xingang West Road, Haizhu District <br>Guangzhou, 510275, P.R. China <br>Tel; 020-84115227 <br>Email: <a href="mailto:qianyk@mail2.sysu.edu.cn">qianyk@mail3.sysu.edu.cn</a> </font> 
<div> <font size="2"> </font> </div></div></div></sign><div> </div><div> </div><div><includetail><div> </div><div> </div><div style="color: rgb(0, 0, 0);"><div style="padding: 2px 0px; font-family: Arial Narrow; font-size: 12px;">------------------ 原始邮件 ------------------</div><div style="background: rgb(239, 239, 239); padding: 8px; font-size: 12px;"><div id="menu_sender"><b>发件人:</b> "yzb"<853245241@qq.com>;</div><div><b>发送时间:</b> 2018年7月16日(星期一) 晚上8:23</div><div><b>收件人:</b> "mitgcm-support"<mitgcm-support@mitgcm.org>; <wbr></div><div></div><div><b>主题:</b> Re: [MITgcm-support] How to diagnose the diapycnal diffusivity?</div></div><div> </div><div style="position: relative;"><div>Thank you. It really helps me. I will read the document carefully.</div><div><div><br></div><div><br></div><div style="padding: 2px 0px; font-family: Arial Narrow; font-size: 12px;">------------------ Original ------------------</div><div style="background: rgb(239, 239, 239); padding: 8px; font-size: 12px;"><div><b>From: </b> "Martin Losch"<Martin.Losch@awi.de>;</div><div><b>Date: </b> Mon, Jul 16, 2018 08:20 PM</div><div><b>To: </b> "MITgcm Support"<mitgcm-support@mitgcm.org>;<wbr></div><div></div><div><b>Subject: </b> Re: [MITgcm-support] How to diagnose the diapycnal diffusivity?</div></div><div><br></div>As far as I know, the parameters \kappa_{\rho} and \kappa_{GM} are constant parameters that you specify in data.gmredi (from GMREDI.h:<br>C     GM_isopycK       :: Isopycnal diffusivity [m^2/s] (Redi-tensor)<br>C     GM_background_K  :: Thickness diffusivity [m^2/s] (GM bolus transport)<br>). You  can have variable \kappa_{GM} (documentation: <https://mitgcm.readthedocs.io/en/latest/phys_pkgs/gmredi.html#variable> ), but you need to turn it on explicitly with the appropriate CPP- and runtime-flags, and then GM_VisbK is the approriate diagnostics for this particular version of variable \kappa_{GM}.<br>In recent code, there are all sorts of other options that I am not familiar with, but they maybe worth exploring.<br><br>Martin<br><br><br>> On 16. Jul 2018, at 13:57, yzb <853245241@qq.com> wrote:<br>> <br>> <br>> Thanks for you reply. Sorry for the confusion of my question. This is something I should have written more about the diapycnal diffusivity. The diffusivity I want to calculate is the diapycnal mixing ( ,unit is m^2/s), maybe can be regarded as eddy coefficient. I am not sure the  which is in the Gmredi pkg can be represented for it. There is a diagnostics option in Gmredi pkg which is called 'GM_VisbK' and this may be relate to my question.<br>> <br>> ------------------ Original ------------------<br>> From:  "Martin Losch"<Martin.Losch@awi.de>;<br>> Date:  Mon, Jul 16, 2018 07:35 PM<br>> To:  "MITgcm Support"<mitgcm-support@mitgcm.org>;<br>> Subject:  Re: [MITgcm-support] How to diagnose the diapycnal diffusivity?<br>> <br>> Yang,<br>> <br>> when you run your model with a default debugLevel (=1), then you’ll a file called “available_diagnostics.log”, that will show you all diagnostics for your particular configuration. If you have turned on GM (useGMredi=.TRUE. in data.pkg), then you the GM-diagnostics should be part of that list. I am not sure if there is a diagnostic for diapycnal diffusivity (do you mean diffusive flux?) …<br>> <br>> Martin<br>> <br>> > On 16. Jul 2018, at 11:12, yzb <853245241@qq.com> wrote:<br>> > <br>> > Hello,<br>> > <br>> >      I want to diagnose the diapycnal diffusivity directly from model. But I don't know the diagnostic name of it. Could anyone tell me the diagnostic name? Meanwhile, I found the package of Gmredi may contain this diagnosis. If so, what's the diagnostic name in that package? <br>> > <br>> >     Thanks,<br>> >     Yang<br>> > _______________________________________________<br>> > MITgcm-support mailing list<br>> > MITgcm-support@mitgcm.org<br>> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br>> <br>> _______________________________________________<br>> MITgcm-support mailing list<br>> MITgcm-support@mitgcm.org<br>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br>> _______________________________________________<br>> MITgcm-support mailing list<br>> MITgcm-support@mitgcm.org<br>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br><br>_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br></div></div></div><!--<![endif]--></includetail></div>