<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">I am not sure, if this is a useful comment, but even in a non-hydrostatic simulation, the flow field is non-divergent and you coiuld diagnose wVel from the vertical integral of the horizontal divergence (in fact, the model does this in a non-hydrostatic simulation, too). In the same way, and since the operations are linear you coud diagnose dw/t from - \nabla_h du/dt, couldn’t you?<div class=""><br class=""></div><div class="">Martin<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 23. May 2022, at 15:36, 钱钰坤 <<a href="mailto:qianyk@mail3.sysu.edu.cn" class="">qianyk@mail3.sysu.edu.cn</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=GB18030" class=""><div class="">Hi,</div><div class=""><br class=""></div><div class="">I might find a workaround to get instantaneous TOTWTEND (not time-mean one)</div><div class="">by adding the following lines:</div><div class=""> <span class="pl-k" style="box-sizing: border-box; color: var(--color-prettylights-syntax-keyword); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;">CALL</span><span style="color: rgb(36, 41, 47); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;" class=""> DIAGNOSTICS_SCALE_FILL( wVel, tmpFac, </span><span class="pl-c1" style="box-sizing: border-box; color: var(--color-prettylights-syntax-constant); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;">1</span><span style="color: rgb(36, 41, 47); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;" class="">, </span><span class="pl-s" style="box-sizing: border-box; color: var(--color-prettylights-syntax-string); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;"><span class="pl-pds" style="box-sizing: border-box; color: var(--color-prettylights-syntax-string);">'</span>TOTWTEND<span class="pl-pds" style="box-sizing: border-box; color: var(--color-prettylights-syntax-string);">'</span></span><span style="color: rgb(36, 41, 47); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;" class="">,</span><span class="pl-c1" style="box-sizing: border-box; color: var(--color-prettylights-syntax-constant); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;">0</span><span style="color: rgb(36, 41, 47); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;" class="">,Nr,</span><span class="pl-c1" style="box-sizing: border-box; color: var(--color-prettylights-syntax-constant); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;">0</span><span style="color: rgb(36, 41, 47); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;" class="">,</span><span class="pl-c1" style="box-sizing: border-box; color: var(--color-prettylights-syntax-constant); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;">1</span><span style="color: rgb(36, 41, 47); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;" class="">,</span><span class="pl-c1" style="box-sizing: border-box; color: var(--color-prettylights-syntax-constant); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;">1</span><span style="color: rgb(36, 41, 47); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;" class="">,myThid )</span></div><div class="">to diagnostics_fill_state.F file, following TOTUTEND and TOTVTEND (see <a href="https://github.com/MITgcm/MITgcm/blob/20a156cdcefcb000006123cf7ddac8fc997cb603/pkg/diagnostics/diagnostics_fill_state.F#L356" class="">here</a>). And</div><div class="">also add a call to </div><div class=""> <span style="color: rgb(36, 41, 47); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;" class="">CALL </span><span style="color: rgb(36, 41, 47); font-family: ui-monospace, SFMono-Regular, "SF Mono", Menlo, Consolas, "Liberation Mono", monospace; font-size: 12px; white-space: pre;" class="">DIAGNOSTICS_ADDTOLIST()</span></div><div class="">at diagnostics_main_init.F (see <a href="https://github.com/MITgcm/MITgcm/blob/20a156cdcefcb000006123cf7ddac8fc997cb603/pkg/diagnostics/diagnostics_main_init.F#L781" class="">here</a>).</div><div class=""><br class=""></div><div class="">Then I re-compile the model and run it well. However, the tendency of w TOTWTEND</div><div class="">seems to be much larger than expected in my case:</div><div class=""><span id="cid:7239677D@602E0261.C58D8B62.png"><7239677D@602E0261.C58D8B62.png></span></div><div class=""><sign signid="0" class=""><div class=""><font class=""><br class=""></font></div><div class=""><font class="">I am not sure DIAGNOSTICS_SCALE_FILL really do the time-differencing. Is this</font></div><div class=""><font class="">the right way to get TOTWTEND? Why it looks so large in magnitude?</font></div><div class=""><font class=""><br class=""></font><div style="color:#909090;font-family:Arial Narrow;font-size:12px" class=""></div></div><div style="font-size: 14px; font-family: Verdana;" class="">
<div class=""> <font size="2" class="">Best regards <br class=""> <br class="">-----------------------------------------<br class="">Yu-Kun Qian (钱钰坤)</font><div style="font-family: Verdana;" class=""><span style="font-size: small;" class="">South China Sea Institute of Oceanology</span></div><font size="2" style="font-family: Verdana;" class="">Sun Yat-sen University<br class="">Xingang West Road, Haizhu District<br class="">Guangzhou, P.R. China<br class="">Email: <u class=""><font color="#0000ff" class=""><a href="mailto:qianyukun@scsio.ac.cn" class="">qianyukun@scsio.ac.cn</a></font></u>, <a href="mailto:qianyk@mail2.sysu.edu.cn" class="">qianyk@mail3.sysu.edu.cn</a></font></div><div class=""> <font size="2" class=""> </font> </div></div></sign></div><div class=""><includetail class=""><div class=""> </div><div style="" class=""><div style="FONT-SIZE: 12px;FONT-FAMILY: Arial Narrow;padding:2px 0 2px 0;" class="">------------------ 原始邮件 ------------------</div><div style="FONT-SIZE: 12px;background:#efefef;padding:8px;" class=""><div id="menu_sender" class=""><b class="">发件人:</b> "Matthew Mazloff"<<a href="mailto:mmazloff@ucsd.edu" class="">mmazloff@ucsd.edu</a>>;</div><div class=""><b class="">发送时间:</b> 2022年5月18日(星期三) 上午9:47</div><div class=""><b class="">收件人:</b> "MITgcm-support"<<a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a>>; <wbr class=""></div><div class=""></div><div class=""><b class="">主题:</b> Re: [MITgcm-support] vertical momentum tendency</div></div><div class=""> </div><div style="position:relative;" class=""><div id="tmpcontent_res" class=""></div>Hi Haoron<div class=""><br class=""></div><div class="">Dump snapshots at the interval you need to get dW/dt.</div><div class=""><br class=""></div><div class="">For example daily tendency would be derived from </div><div class=""> frequency(1) = -86400.0,<br class=""> timePhase(1) = 0.0,<br class=""> fields(1,1) = 'WVEL ',<br class=""> filename(1) = ‘WvelSnap’,</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Matt</div><div class=""><br class=""><div class=""><br class=""><blockquote type="cite" class=""><div class="">On May 17, 2022, at 5:47 PM, 钱钰坤 <<a href="mailto:qianyk@mail3.sysu.edu.cn" class="">qianyk@mail3.sysu.edu.cn</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class=""><font class="">Hi Haoran,</font></div><div class=""><br class=""></div><div class="">Thank you for your help. But I don't really know how these two diagnostics</div><div class="">could help in deriving the tendency of w (unit of acceleration).</div><div class=""><br class=""></div><div class="">My understanding is that these two diagnostics are useful in calculating the</div><div class="">eddy momentum flux when the output is time-average <> one. For example,</div><div class=""><w'u'> = <wu> - <w><u></div><div class="">If we output monthly mean wu, w, and u, we can get eddy flux <w'u'> for this</div><div class="">month.</div><div class=""><div class=""><font class=""><br class=""></font><div style="color:#909090;font-family:Arial Narrow;font-size:12px" class=""></div></div><div style="font-size: 14px; font-family: Verdana;" class="">
<div class=""> <font size="2" class="">Best regards <br class=""> <br class="">-----------------------------------------<br class="">Yu-Kun Qian (钱钰坤)</font><div style="font-family: Verdana;" class=""><span style="font-size: small;" class="">South China Sea Institute of Oceanology</span></div><font size="2" style="font-family: Verdana;" class="">Sun Yat-sen University<br class="">Xingang West Road, Haizhu District<br class="">Guangzhou, P.R. China<br class="">Email: <u class=""><font color="#0000ff" class=""><a href="mailto:qianyukun@scsio.ac.cn" class="">qianyukun@scsio.ac.cn</a></font></u>, <a href="mailto:qianyk@mail2.sysu.edu.cn" class="">qianyk@mail3.sysu.edu.cn</a></font></div><div class=""> <font size="2" class=""> </font> </div>
</div></div><div class=""> </div><div class=""><div class=""> </div><div class=""> </div><div style="" class=""><div style="FONT-SIZE: 12px;FONT-FAMILY: Arial Narrow;padding:2px 0 2px 0;" class="">------------------ 原始邮件 ------------------</div><div style="FONT-SIZE: 12px;background:#efefef;padding:8px;" class=""><div id="menu_sender" class=""><b class="">发件人:</b> "Haoran Xu"<<a href="mailto:xu1jian2wei3@163.com" class="">xu1jian2wei3@163.com</a>>;</div><div class=""><b class="">发送时间:</b> 2022年5月17日(星期二) 晚上7:40</div><div class=""><b class="">收件人:</b> "MITgcm-support"<<a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a>>; <wbr class=""></div><div class=""></div><div class=""><b class="">主题:</b> Re: [MITgcm-support] vertical momentum tendency</div></div><div class=""> </div><div style="position:relative;" class=""><div id="tmpcontent_res" class=""></div><div style="line-height: 1.7; font-size: 14px; font-family: Arial;" class=""><div style="margin: 0;" class=""><span style="font-family: Arial;" class="">Hi，</span></div><div style="margin: 0;" class=""><span style="font-family: Arial;" class=""> If I <span style="color: rgb(102, 102, 102); font-family: tahoma, "Microsoft YaHei", 微软雅黑, 宋体, "Malgun Gothic", sans-serif;" class="">understand correctly,these two diagnostics could help.</span></span></div><div style="margin: 0;" class=""><br class=""></div><div style="margin: 0;" class=""><br class=""></div><div style="margin: 0;" class=""><pre style="box-sizing: border-box; font-family: SFMono-Regular, Menlo, Monaco, Consolas, "Liberation Mono", "Courier New", Courier, monospace; font-size: 12px; white-space: pre; line-height: 1.4; margin-top: 0px; margin-bottom: 0px; padding: 12px; overflow: auto; color: rgb(64, 64, 64);" class="">WU_VEL |WU LR|m.m/s^2 |Vertical Transport of Zonal Momentum
WV_VEL |WV LR|m.m/s^2 |Vertical Transport of Meridional Momentum</pre></div><div style="margin: 0;" class=""><br class=""></div><div style="margin: 0;" class=""> Haoran Xu.</div><div style="position:relative;zoom:1" class=""></div><div style="margin: 0;" class=""><br class=""></div><p id="isReplyContentTip" class="">At 2022-05-17 18:30:48, "钱钰坤" <<a href="mailto:qianyk@mail3.sysu.edu.cn" class="">qianyk@mail3.sysu.edu.cn</a>> wrote:</p><blockquote id="isReplyContent" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid" class=""><div class=""><font class="">Hi all,</font></div><div class=""><font class=""><br class=""></font></div><div class=""><font class="">Just a quick question. Why there is no tendency of vertical momentum,</font></div><div class=""><font class="">like TOTWTEND, in the diagnostics even the non-hydrostatic option is</font></div><div class=""><font class="">turned on?</font></div><div class=""><font class=""><br class=""></font></div><div class=""><font class="">To close the energy budget, I guess we need this diagnostic. Is there any</font></div><div class=""><font class="">workaround to obtain it if there is no TOTWTEND?</font></div><div class=""><br class=""></div><div class=""><div class=""><div style="color:#909090;font-family:Arial Narrow;font-size:12px" class=""></div></div><div style="font-size: 14px; font-family: Verdana;" class=""><div class=""> <font size="2" class="">Best regards.<br class=""> <br class=""></font><span style="font-family: "lucida Grande", Verdana, "Microsoft YaHei";" class="">---------------------------</span><font size="2" class=""><br class="">Yu-Kun Qian (钱钰坤)</font></div><div class=""><font size="2" class="">South China Sea Institute of Oceanology<br class=""></font></div><font size="2" class="">Sun Yat-sen University <br class="">Xingang West Road, Haizhu District <br class="">Guangzhou, P.R. China<br class="">Email: <a href="mailto:qianyk@mail2.sysu.edu.cn" class="">qianyk@mail3.sysu.edu.cn</a> </font>
<div class=""> <font size="2" class=""> </font> </div></div></div><div class=""> </div><div class=""></div></blockquote></div><br class=""><br class=""><span title="neteasefooter" class=""><div class=""> <br class="webkit-block-placeholder"></div></span></div></div></div>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class=""><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=-35OiAkTchMrZOngvJPOeA&r=HSHS34ROdfL7f2oVxuKB761hrMvWN1RWA5k3SHGts_s&m=3dmk9mE5wbsK8pGhM9ik2Mywehebv53YMmvBG-OJeJ1YILefVCNpFtkCisx_V2rL&s=-KY9LHIxlYVAlDaSW3axV0lhAZkr0gMbYyogCKxqBQs&e=" class="">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=-35OiAkTchMrZOngvJPOeA&r=HSHS34ROdfL7f2oVxuKB761hrMvWN1RWA5k3SHGts_s&m=3dmk9mE5wbsK8pGhM9ik2Mywehebv53YMmvBG-OJeJ1YILefVCNpFtkCisx_V2rL&s=-KY9LHIxlYVAlDaSW3axV0lhAZkr0gMbYyogCKxqBQs&e=</a> <br class=""></div></blockquote></div><br class=""></div>
</div></div><!--<![endif]--></includetail></div>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br class=""></div></blockquote></div><br class=""></div></body></html>