<div>Thank you. It really helps me. I will read the document carefully.</div><div><div><br></div><div><br></div><div style="font-size: 12px;font-family: Arial Narrow;padding:2px 0 2px 0;">------------------ Original ------------------</div><div style="font-size: 12px;background:#efefef;padding:8px;"><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>