<div dir="ltr">Dhruv,<div><br></div><div>There is no built-in ability to do this. It should be possible to implement it yourself relatively easily.</div><div><br></div><div>The general approach is to modify the function EXTERNAL_FORCING_T in model/src/external_forcing.F. Here you can add custom fortran code to apply an arbitrary, problem-specific forcing.</div><div><br></div><div>The one big challenge is how to calculate T_bar, the zonal mean temperature. This is not defined locally to a tile, so it will involve some sort of MPI call to sum the temperature field globally along the i-axis. I don't know how to do this. Maybe Jean Michel can give a hint.</div><div><br></div><div>Cheers,</div><div>Ryan</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Dec 1, 2017 at 12:28 PM, Dhruv Balwada <span dir="ltr"><<a href="mailto:db194@nyu.edu" target="_blank">db194@nyu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello all,<br>
<br>
I was wondering if there is a standard package, or some modified code that someone has written which would all me to do the following task. With RBCS package I can restore the stratification (temp and salt), to some prescribed profiles. However the way this works is by relaxing the synoptic fields, such as 1/(tau)*(T(x,y,z,t) - T_ref). Instead I want to not restore the synoptic fields, but instead some zonal mean profiles, such as 1/(tau)*(T_bar(y,z,t) - T_ref). The goal is not overly damp the large eddies, but instead just nudge the mean state.<br>
<br>
Any advice or help would be appreciated.<br>
<br>
Thank you,<br>
Dhruv<br>
<br>
<br>
<br>
______________________________<wbr>_________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
</blockquote></div><br></div>