[MITgcm-support] Restoring the mean, rather than the synoptic fields

[Ryan Abernathey]

Dhruv,

There is no built-in ability to do this. It should be possible to implement
it yourself relatively easily.

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.

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.

Cheers,
Ryan



On Fri, Dec 1, 2017 at 12:28 PM, Dhruv Balwada <db194 at nyu.edu> wrote:

> Hello all,
>
> 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.
>
> Any advice or help would be appreciated.
>
> Thank you,
> Dhruv
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20171206/08f308fa/attachment.html>