<div>Hi,</div><div><br></div><div>I might find a workaround to get instantaneous TOTWTEND (not time-mean one)</div><div>by adding the following lines:</div><div> <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;"> 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;">, </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;">,</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;">,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;">,</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;">,</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;">,myThid )</span></div><div>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">here</a>). And</div><div>also add a call to </div><div> <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;">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;">DIAGNOSTICS_ADDTOLIST()</span></div><div>at diagnostics_main_init.F (see <a href="https://github.com/MITgcm/MITgcm/blob/20a156cdcefcb000006123cf7ddac8fc997cb603/pkg/diagnostics/diagnostics_main_init.F#L781">here</a>).</div><div><br></div><div>Then I re-compile the model and run it well. However, the tendency of w TOTWTEND</div><div>seems to be much larger than expected in my case:</div><div><img src="cid:7239677D@602E0261.C58D8B62.png"></div><div><sign signid="0"><div><font><br></font></div><div><font>I am not sure DIAGNOSTICS_SCALE_FILL really do the time-differencing. Is this</font></div><div><font>the right way to get TOTWTEND? Why it looks so large in magnitude?</font></div><div><font><br></font><div style="color:#909090;font-family:Arial Narrow;font-size:12px"></div></div><div style="font-size:14px;font-family:Verdana;color:#000;">
<div> <font size="2">Best regards <br> <br>-----------------------------------------<br>Yu-Kun Qian (钱钰坤)</font><div style="font-family: Verdana;"><span style="font-size: small;">South China Sea Institute of Oceanology</span></div><font size="2" style="font-family: Verdana;">Sun Yat-sen University<br>Xingang West Road, Haizhu District<br>Guangzhou, P.R. China<br>Email: <u><font color="#0000ff">qianyukun@scsio.ac.cn</font></u>, <a href="mailto:qianyk@mail2.sysu.edu.cn">qianyk@mail3.sysu.edu.cn</a></font></div><div> <font size="2"> </font> </div></div></sign></div><div><includetail><div> </div><div style="font:Verdana normal 14px;color:#000;"><div style="FONT-SIZE: 12px;FONT-FAMILY: Arial Narrow;padding:2px 0 2px 0;">------------------ 原始邮件 ------------------</div><div style="FONT-SIZE: 12px;background:#efefef;padding:8px;"><div id="menu_sender"><b>发件人:</b> "Matthew Mazloff"<mmazloff@ucsd.edu>;</div><div><b>发送时间:</b> 2022年5月18日(星期三) 上午9:47</div><div><b>收件人:</b> "MITgcm-support"<mitgcm-support@mitgcm.org>; <wbr></div><div></div><div><b>主题:</b> Re: [MITgcm-support] vertical momentum tendency</div></div><div> </div><div style="position:relative;"><div id="tmpcontent_res"></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><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="">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= <br class=""></div></blockquote></div><br class=""></div>
</div></div><!--<![endif]--></includetail></div>