[MITgcm-devel] CALL IMPLDIFF

[Jean-Michel Campin]

Hi Patrick,

I would like to change the argument "tracerId" in the 4 IMPLDIFF calls
that are in dynamics.F:
for e.g., -1 with gU ; -2 with gV ; and keep 0 for vVelD & uVelD
( or 0,-1,-2,-3 with gU,gV,vVelD,uVelD )

The purpose is to add a diagnostic of vertical viscous flux
for U & V.

Is this going to cause problems for the adjoint (if
this tracerId is used in storage directive ?)

Thanks,

Jean-Michel