<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<div style="font-family: '微软雅黑'; font-size: 10.5pt; color: #000000">
<p class="MsoNormal"><span lang="EN-US">Hello everyone,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> I’m running an ocean model with a spatial varying vertical diffusivity field, by setting “<span style="font-family: Georgia;"><i>diffKrFile=diffKr.bin</i></span>”. Meanwhile, The KPP package is disabled and the GMRedi
(code is attached below) package is enabled. <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> I wonder if the vertical diffusivity field would change with the model run, or it just keeps being its preset value. And how to output the vertical diffusivity field?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> I tried to look into the diffusive flux (<span style="font-family: Georgia;">DFrE_TH, DFrI_TH, DFrE_SLT
</span>and <span style="font-family: Georgia;">DFrI_SLT</span>) produced by diagnostics package, but they have no pattern of the preset diffusivity. Is there any way to calculate the diffusivity by diffusive flux (As the unit of diffusive flux is
</span>℃·<span lang="EN-US">m^3</span>·<span lang="EN-US">s, it seems the diffusivity cannot be calculated by k = F/grad_z(theta))? Besides, I noticed the DFrE_SLT and DFrI_SLT are all zeros. Does that mean the salinity field has no diffusion? If so, how to
set the salinity diffusivity?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US">Many thanks for any advice!<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> </span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US">Jianing Li<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> </span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US">data.gmredi:<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"># GM+Redi package parameters:<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> </span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US">#-from MOM :<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"># GM_background_K: G & Mc.W diffusion coefficient<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"># GM_maxSlope : max slope of isopycnals<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"># GM_Scrit : transition for scaling diffusion coefficient<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"># GM_Sd : half width scaling for diffusion coefficient<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"># GM_taper_scheme: slope clipping or one of the tapering schemes<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"># GM_Kmin_horiz : horizontal diffusion minimum value
<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> </span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US">#-Option parameters (needs to "define" options in GMREDI_OPTIONS.h")<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"># GM_isopycK : isopycnal diffusion coefficient (default=GM_background_K)<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"># GM_AdvForm : turn on GM Advective form (default=Skew flux form)<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> </span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> &GM_PARM01
<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> GM_background_K = 1.e+3,<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> GM_taper_scheme = 'gkw91',<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> GM_maxSlope = 1.e-2,<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> GM_Kmin_horiz = 0.,<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> GM_Scrit = 4.e-3,<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> GM_Sd = 1.e-3,<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:10.5pt;mso-char-indent-count:1.0"><span lang="EN-US"> &<o:p></o:p></span></p>
</div>
</body>
</html>