<div dir="ltr"><div><div><div>I have tried to compute it myself, and it works!  Since thie error message is showing up in the STDERR files, I thought there might be a problem somewhere else...<br><br></div>And I think you're right, my vertical diffusivity should be constant everywhere. Do I understand that if a diagnostics is all zeros, the " WARNING -   has not been filled (ndiag=  0 )" is written in STDERRs?<br></div><br>--<br></div>Pascal<br></div><br><div class="gmail_quote"><div dir="ltr">Le mar. 10 juil. 2018 à 11:42, Bertrand Louis Rene Delorme <<a href="mailto:bdelorme@stanford.edu" target="_blank">bdelorme@stanford.edu</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div lang="FR" link="#0563C1" vlink="#954F72">
<div class="m_7312205381034102362m_-7010446484631517999WordSection1">
<p class="MsoNormal"><span lang="EN-US">Salut Pascal,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">I am not sure about the
</span><span lang="EN-US">“momVort3” diagnostic: have you tried to compute this term yourself with UVEL and VVEL to see what solution you get? For the “DFrI_TH” diagnostic, maybe your vertical diffusivity is just constant
 everywhere? <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Cheers,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Bertrand<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">De : </span></b><span style="font-size:12.0pt;color:black">MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" target="_blank">mitgcm-support-bounces@mitgcm.org</a>> au nom de Pascal Bourgault <<a href="mailto:pascal.bourgault@gmail.com" target="_blank">pascal.bourgault@gmail.com</a>><br>
<b>Répondre à : </b>"<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>" <<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>><br>
<b>Date : </b>mardi 10 juillet 2018 à 8:07 AM<br>
<b>À : </b>"<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>" <<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>><br>
<b>Objet : </b>[MITgcm-support] Diagnostics have not been filled<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Hi, <u></u><u></u></p>
</div>
<p class="MsoNormal">In my data.diagnostics, I ask the model to output the "momVort3" and the "DFrI_TH" diagnostic, and many others. I am using mpi and the following packages are compiled in my mitgcmuv : gfp, flt, mnc, mdsio. Also, momVort3 and DFrI_TH both
 appear in the available_diagnostics.log file.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">However, the output values for those are all zeros in the final files. For example, in my dynDiag*.nc files, where I save UVEL, VVEL, WVEL and momVort3, the three velocities look good, but momVort3 is only zeros.<u></u><u></u></p>
</div>
<p class="MsoNormal">In the STDERR.* of the model I have this<br>
"""<br>
(PID.TID 0000.0001) - WARNING - from DIAGNOSTICS_OUT at iter=     17286<br>
(PID.TID 0000.0001) - WARNING -   diag.#   117 : DFrI_TH  (#   7 ) in outp.Stream: tempDiag<br>
(PID.TID 0000.0001) - WARNING -   has not been filled (ndiag=  0 )<br>
(PID.TID 0000.0001) WARNING DIAGNOSTICS_OUT  => write ZEROS instead<br>
(PID.TID 0000.0001) - WARNING - from DIAGNOSTICS_OUT at iter=     17286<br>
(PID.TID 0000.0001) - WARNING -   diag.#   174 : momVort3 (#   4 ) in outp.Stream: dynDiag<br>
(PID.TID 0000.0001) - WARNING -   has not been filled (ndiag=  0 )<br>
 """<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">repeating numberous times.<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Are there some options in some *_OPTIONS.h files I need to turn on to save those fields? I don't really understand the error...<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">I can provide any other config info needed, I just don't know what is useful in this case...<u></u><u></u></p>
</div>
<p class="MsoNormal">Thanks!<br>
<br>
--<u></u><u></u></p>
</div>
<p class="MsoNormal">Pascal<u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>

_______________________________________________<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" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>