[MITgcm-devel] ptracer conservation

[Matthew Mazloff]

Hi Martin

Unfortunately this is causing major recomputations. I was able to fix one with an extra ptracer store in thermodynamics.F. Analogous to lines 140 to 142 but for ptracer, and starting on line 155.

But I still have a lot of the nature:
TAF RECOMPUTATION LOOP WARNING DOLOOP_STMT ad_input_code.f:476756 in the_main_loop
 extensive recomputations are required.
 The responsible variables are : totsurfcorptr

I tried putting a store of “totsurfcorptr” in ptracers_ad_check_lev?_dir.h but that didn’t help

I have to run but can work on this more tonight

-Matt


> On Aug 22, 2016, at 7:43 AM, Martin Losch <Martin.Losch at awi.de> wrote:
> 
> 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
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel