[MITgcm-support] Vertical Biharmonic Diffusivity & a Passing Solid Object.

[Martin Losch]

Hi Jason,

1) vertical viscousity is only harmonic/laplacian. Vertical viscous fluxes are computed in mom_common/mom_[u,v]_rviscflux.F
This routine is called from either mom_fluxform/mom_fluxform.F or mom_vecinv/mom_vecinv.F
You might want to check out both routines to see how the horizontal biharmonic fluxes are computed.
2) Not quite clear to me what you want to do. Initial conditions (for velocity) can be specified in "data&PARM05". In general there is external_forcing.F that contains subroutines external_forcing_[u,v]. You can modify these to add extra tendencies to gU and gV.

Martin

On Sep 21, 2011, at 8:48 PM, Flanagan, Jason FORNATL, EI wrote:

> Hi all,
>  
> We are currently interested in modeling Thermohaline staircases and the effects of a passing solid object on vertical transport (temperature/salinity fluxes etc). To this end we would greatly appreciate advice on the following;
>  
> 1)      We have had some success in implementing the horizontal operator for biharmonic diffusivity (‘diffK4’ in our input/data file). Is there an equivalent vertical operator available in MITgcm? If so, how can we access it? If not, any advice on how to implement this operator in MITgcm would be most welcome.
> 2)      Where in the source code would be the best place to introduce a disruption to the velocity components that would be ellipsoidal in nature and non-stationary?
>  
> Many Thanks,
> Jason.
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support