[MITgcm-support] obcs_calc.F changes

[Martin Losch]

Ashely,

sponge layers in OBCS relax only to the same values as the boundary values (which don’t have to be tRef, please have a look at the documentation for more details, in particular <http://mitgcm.org/public/r2_manual/latest/online_documents/node236.html>, with the sections “OBCS_SPONGE” and “OBCS_PRESCRIBE_READ”)

What you want to do should be possible with the RBCS package: <http://mitgcm.org/public/r2_manual/latest/online_documents/node237.html>

maybe verifcation/exp4 can serve as an example how to use these packages, as it use both.

You can’t hack the code easily, because the initial conditions (e.g. the content of hydroThetaFile) are read into the dynamic field “theta” that is update at every timestep (so the intial conditions are lost)

Martin


> On 05 Mar 2015, at 23:19, Stroman, Ashley <acs08g at my.fsu.edu> wrote:
> 
> Hello.
> 
> I am trying to use the OBCS package to implement a sponge layer. Instead of relaxing temperature to tRef, which is the default parameter, I want to default to the temperature set up by HydrogThetaFile. Does anyone know what the command should be in obcs_calc.F instead of, for example, OBWt(j,k,bi,bj) =tRef(k)? 
> 
> Thanks
> Ashley Stroman
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support