<div dir="ltr">Hi Yang,<div>This is not too bad.  I'll give explicit instructions.  First, declare a global variable, diffKTr, (along with diffKr, which is just the background, not total, vertical diffusivity, I believe) in DYNVARS.h:</div><div><br></div><div><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">      COMMON /DYNVARS_DIFFKR/</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">     &                       diffKr, diffKTr</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">      _RL  diffKr (</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures">-OLx:sNx+OLx,</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures">-OLy:sNy+OLy,Nr,nSx,nSy)</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">      _RL  diffKTr(</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures">-OLx:sNx+OLx,</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures">-OLy:sNy+OLy,Nr,nSx,nSy)</span></p></div><div><br></div><div>Then add the following lines of code to calc_3d_diffusivity.F to fill that global array with the total vertical diffusivity:</div><div><br></div><div><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">      </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">DO</span><span style="font-variant-ligatures:no-common-ligatures"> j </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">-</span><span style="font-variant-ligatures:no-common-ligatures">OLy, sNy</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">+</span><span style="font-variant-ligatures:no-common-ligatures">OLy</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">       </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">DO</span><span style="font-variant-ligatures:no-common-ligatures"> i </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">-</span><span style="font-variant-ligatures:no-common-ligatures">OLx, sNx</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">+</span><span style="font-variant-ligatures:no-common-ligatures">OLx</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">        diffKTr(i,j,k,bi,bj) </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures"> KappaRTr(i,j,k)</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">       </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">ENDDO</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">      </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">ENDDO</span></p></div><div><br></div><div>Then add the following code to diagnostics_main_init.F to add the vertical diffusivity as a diagnostic:</div><div><br></div><div><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">      diagName  </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">'KVDIFF  '</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">      diagTitle </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)"> </span><span style="font-variant-ligatures:no-common-ligatures">'Vertical diffusivity for tracers  '</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">      diagUnits </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)"> </span><span style="font-variant-ligatures:no-common-ligatures">'m^2/s           '</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">      diagCode  </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)"> </span><span style="font-variant-ligatures:no-common-ligatures">'SM      LR      '</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">        </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(213,59,211)">CALL</span><span style="font-variant-ligatures:no-common-ligatures"> DIAGNOSTICS_ADDTOLIST( diagNum,</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">     </span><span style="font-variant-ligatures:no-common-ligatures;background-color:rgb(230,230,0)">I</span><span style="font-variant-ligatures:no-common-ligatures">          diagName, diagCode, diagUnits, diagTitle, </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">0</span><span style="font-variant-ligatures:no-common-ligatures">, myThid )</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small">Then add the following lines of code to diagnostics_fill_state.F to fill that diagnostic:</span><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">        </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">IF</span><span style="font-variant-ligatures:no-common-ligatures"> ( DIAGNOSTICS_IS_ON(</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">'KVDIFF  '</span><span style="font-variant-ligatures:no-common-ligatures">,myThid) ) </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">THEN</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">         </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">DO</span><span style="font-variant-ligatures:no-common-ligatures"> bj </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures"> myByLo(myThid), myByHi(myThid)</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">          </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">DO</span><span style="font-variant-ligatures:no-common-ligatures"> bi </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures"> myBxLo(myThid), myBxHi(myThid)</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">           </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">DO</span><span style="font-variant-ligatures:no-common-ligatures"> k</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures">,Nr</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">            </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">DO</span><span style="font-variant-ligatures:no-common-ligatures"> j </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures">,sNy</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">             </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">DO</span><span style="font-variant-ligatures:no-common-ligatures"> i </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures">,sNx</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">+</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">               tmpMk(i,j,k,bi,bj) </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">=</span><span style="font-variant-ligatures:no-common-ligatures"> diffKTr(i,j,k,bi,bj)</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">             </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">ENDDO</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">            </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">ENDDO</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">           </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">ENDDO</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">          </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">ENDDO</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">         </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">ENDDO</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">         </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(213,59,211)">CALL</span><span style="font-variant-ligatures:no-common-ligatures"> DIAGNOSTICS_FILL(tmpMk,</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">'KVDIFF  '</span><span style="font-variant-ligatures:no-common-ligatures">,</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">0</span><span style="font-variant-ligatures:no-common-ligatures">,Nr,</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">0</span><span style="font-variant-ligatures:no-common-ligatures">,</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures">,</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">1</span><span style="font-variant-ligatures:no-common-ligatures">,myThid)</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">












</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">        </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(205,121,35)">ENDIF</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><br></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small">Finally, add the following lines to input/data.diagnostics (change the frequency or output file number if you need to):</span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">  frequency(18) = 2635200.0,</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">   fields(1:2,18) = </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">'KVDIFF  '</span><span style="font-variant-ligatures:no-common-ligatures">,</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(195,55,32)">'GM_Kwz  '</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">

</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">   filename(18) = </span><span style="font-variant-ligatures:no-common-ligatures">'diags/kappa_3d_set1'</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">,</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)"><br></span></p><p style="margin:0px;font-stretch:normal;line-height:normal"><span style="font-variant-ligatures:no-common-ligatures">In the output files, diags/kappa_3d_set1*, you just need to subtract the second field (GM_Kwz) from the first (KVDIFF) to get the diapycnal diffusivity.  That's my understanding.</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small">-David</span></span></p></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jul 16, 2018 at 7:39 AM, 钱钰坤 <span dir="ltr"><<a href="mailto:qianyk@mail3.sysu.edu.cn" target="_blank">qianyk@mail3.sysu.edu.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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><u></u><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" target="_blank">qianyk@mail3.sysu.edu.cn</a> </font> 
<div> <font size="2"> </font> </div></div></div><u></u><div> </div><div> </div><div><u></u><div> </div><div> </div><div style="color:rgb(0,0,0)"><div style="padding:2px 0px;font-family:Arial Narrow;font-size:12px">------------------ 原始邮件 ------<wbr>------------</div><div style="background:rgb(239,239,239);padding:8px;font-size:12px"><div id="m_-4829972403956366898menu_sender"><b>发件人:</b> "yzb"<<a href="mailto:853245241@qq.com" target="_blank">853245241@qq.com</a>>;</div><div><b>发送时间:</b> 2018年7月16日(星期一) 晚上8:23</div><div><b>收件人:</b> "mitgcm-support"<<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-<wbr>support@mitgcm.org</a>>; </div><div></div><div><b>主题:</b> Re: [MITgcm-support] How to diagnose the diapycnal diffusivity?</div></div><div> </div><div><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"<<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>>;</div><div><b>Date: </b> Mon, Jul 16, 2018 08:20 PM</div><div><div class="h5"><div><b>To: </b> "MITgcm Support"<<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@<wbr>mitgcm.org</a>>;</div><div></div><div><b>Subject: </b> Re: [MITgcm-support] How to diagnose the diapycnal diffusivity?</div></div></div></div><div><div class="h5"><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: <<a href="https://mitgcm.readthedocs.io/en/latest/phys_pkgs/gmredi.html#variable" target="_blank">https://mitgcm.readthedocs.<wbr>io/en/latest/phys_pkgs/gmredi.<wbr>html#variable</a>> ), 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 <<a href="mailto:853245241@qq.com" target="_blank">853245241@qq.com</a>> 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"<<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>>;<br>> Date:  Mon, Jul 16, 2018 07:35 PM<br>> To:  "MITgcm Support"<<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@<wbr>mitgcm.org</a>>;<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 <<a href="mailto:853245241@qq.com" target="_blank">853245241@qq.com</a>> 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>> > ______________________________<wbr>_________________<br>> > MITgcm-support mailing list<br>> > <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>> > <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>> <br>> ______________________________<wbr>_________________<br>> MITgcm-support mailing list<br>> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>> ______________________________<wbr>_________________<br>> MITgcm-support mailing list<br>> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br><br>______________________________<wbr>_________________<br>MITgcm-support mailing list<br><a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br></div></div></div></div></div><u></u></div><br>______________________________<wbr>_________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
<br></blockquote></div><br></div>