[MITgcm-support] change hFac in mom_vecinv package

[Alistair Adcroft]

All sounds fine - I'll let you know how it changes all my tests that do have
VI with topography. I'm running MITgcm in shallow water mode at the moment
examining energy and enstrophy conservation (using Sadourny coriolis) and
these are all sensitive to how the hFac' are done - including the upwinding
of fixed topography - which, as we discussed before, breaks these
properties.

A.
--
Dr Alistair Adcroft              email: Alistair.Adcroft at noaa.gov
NOAA/GFDL, Princeton University               Tel: (609) 452-6538
Forrestal Campus, U.S. Route 1, P.O. Box 308, Princeton, NJ 08542


-----Original Message-----
From: mitgcm-support-bounces at mitgcm.org
[mailto:mitgcm-support-bounces at mitgcm.org] On Behalf Of Jean-Michel Campin
Sent: Saturday, August 02, 2003 2:08 PM
To: mitgcm-support at mitgcm.org
Subject: [MITgcm-support] change hFac in mom_vecinv package



Hi,

I will check-in few changes in the mom_vecinv package:
In order to use it with partial cell, I will replace few
S/R with the modified source from aim_5l_cs/code.
a) coriolis: no hFac
b) flux of relative vorticity: I keep the original form
(that was not working very well with partial cell)
within "IF ( use_original_hFac ) THEN / ENDIF" block, with
   PARAMETER ( use_original_hFac=.FALSE. )
and use the aim_5l_cs/code version as default. I don't have time now to test
other options (as, e.g, a formulation with no hFac).

I had to say that the new version is doing well on different tests (exp4,
enstrophy and energy conservation with z* in a SSH adjustment problem with
rotation or solid body rotation background flow) and has been used for the
atmosphere (both AIM & HS) with topography for more than a year now without
any problem.

These changes do not affect any verification/exp/ results because none that
uses VI has topography.

There might be some issue with TAF & the adjoint (sorry, Patrick).

Jean-Michel

_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org
http://dev.mitgcm.org/mailman/listinfo/mitgcm-support