[MITgcm-support] Keeping track of the amount of tracer removed by sponge layers

[Martin Losch]

Hi Dan, 

I don’t think that there is a predefined diagnostic that does that for you. But implementing your own diagnostics is really not so hard. You need to define it somewhere (I would do it in rbcs_init_fixed) and then compute the diagostics and “fill” it by calling diagnostics_fill in rbcs_add_tendency. There are examples in almost every package of the model, often the diagnostics are defined in $pkg_diagnostics_init.F

There is one complication that I didn’t think off initially: rbcs_add_tendency is called per vertical level and diagnostics_fill may or may not be able to fill one vertical level at a time for a 3D diagnostic (never tried this but it may work by setting klev=k and nlev=1 for a 3D diagnostic). If it does not work, you may have to define an extra global 3D field in RBCS_FIELDS.h that you fill level by level in rbcs_add_tendency and then call diagnostics_fill with this 3D field somewhere from do_statevars_diags.F (there are examples how this is realized in other packages).

Hope that helps,

Martin


> On 09 Apr 2015, at 10:18, Dan Jones <dcjones.work at gmail.com> wrote:
> 
> Hi all,
> 
> I am using the rbcs package to implement a tracer-absorbing sponge layer in both (i.) the northern edge of my Southern Ocean model domain and (ii.) in the top 150 m. 
> 
> Does anyone know of a good way to keep track of the amount of tracer removed at each grid cell by the sponge layers?  Ideally, I'd like to use the diagnostics package for consistency with my other output, but I didn't see an obvious choice in the available diagnostics log.
> 
> Thanks,
> Dan
> 
> ---------------------------------------------------------
> 
> Dr Dan Jones
> British Antarctic Survey, NERC
> Cambridge, UK
> 
> Phone: +44 (0)1223 221505
> Fax: +44 (0)1223 362616
> Skype:  dcjones.work
> 
> ---------------------------------------------------------
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support