[MITgcm-devel] ptracer conservation

[Bruno Deremble]

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
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ptracers_readparms.F
Type: application/octet-stream
Size: 9439 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20160727/1f134c63/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ptracers_calc_wsurf_tr.F
Type: application/octet-stream
Size: 2633 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20160727/1f134c63/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ptracers_apply_forcing.F
Type: application/octet-stream
Size: 3378 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20160727/1f134c63/attachment-0005.obj>