<div dir="ltr"><div><div></div>So I did some tests and setting vectorInvariantMomentum and implicitDiffusion to TRUE did activate the computation and filling of momVort3 and DFrI_TH. But it changes the equations used, so computing momVort3 afterwards, or commenting the line in mom_fluxform.F, makes more sense here.<br><br>Merci à tous!<br><br>--<br></div>Pascal<br></div><br><div class="gmail_quote"><div dir="ltr">Le mer. 11 juil. 2018 à 03:32, Bruno Deremble <<a href="mailto:bruno.deremble@ens.fr">bruno.deremble@ens.fr</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I have had the same issue with the vorticity. If you are using the flux<br>
form, the vorticity is only computed if you use variable viscosity<br>
<br>
a quick fix is to comment<br>
IF ( useVariableVisc ) THEN<br>
in mom_fluxform.F<br>
<br>
bruno<br>
<br>
On Tuesday, Jul 10 2018, 17:41:47, Bertrand Louis Rene Delorme <<a href="mailto:bdelorme@stanford.edu" target="_blank">bdelorme@stanford.edu</a>> wrote:<br>
> Salut Pascal,<br>
><br>
> I am not sure about the “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?<br>
><br>
> Cheers,<br>
> Bertrand<br>
><br>
><br>
> De : 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>
> Répondre à : "<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>
> Date : mardi 10 juillet 2018 à 8:07 AM<br>
> À : "<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>
> Objet : [MITgcm-support] Diagnostics have not been filled<br>
><br>
> Hi,<br>
> 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.<br>
><br>
> 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.<br>
> 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>
> """<br>
> repeating numberous times.<br>
> 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...<br>
> I can provide any other config info needed, I just don't know what is useful in this case...<br>
> Thanks!<br>
_______________________________________________<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>