[MITgcm-devel] ptracer conservation

[Martin Losch]

Hi Bruno,

I just added your code (with a slight modification in 
ptracers_readparms.F). Please test, if you find this OK.

The verification results are not changed, but we may want to change one 
of the experiements to use this flag (actually: linFSConserveTr). Does 
global_ocean.90x40x15.dwnslp make sense? There, one ptracer is used to 
mimic salinity and we have linear free surface.

Martin

On 07/27/2016 09:27 PM, Bruno Deremble wrote:
>
> Hi,
>
> I was looking with a student at an offline configuration for tracer
> advection. The (u,v) velocity field comes form another model (Hycom) and
> has a divergent component. With an implicit free surface, it seems that
> the tracer field is not conserved.
>
> For some reason the correction term needed to conserve the ptracer with
> a free surface is only partly implemented (cf. flag PTRACERS_calcSurfCor
> in thermodynamics.F)
>
> I finished the implementation in a similar way than what is currently
> done for temperature and salinity.
>
> the modified files are (cf. attached) :
> ptracers_calc_wsurf_tr.F
> ptracers_apply_forcing.F
> ptracers_readparms.F
>
> and the associated flag is PTRACERS_calcSurfCor (set by default to false)
>
> these are minor modifications but I think it could be helpful to add it
> to the main mitgcm.
>
> thanks
> bruno
>
>
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>

-- 
Martin Losch
Alfred Wegener Institute for Polar and Marine Research
Postfach 120161, 27515 Bremerhaven, Germany;
Tel./Fax: ++49(0471)4831-1872/1797