[MITgcm-devel] momentum ...

[Alistair Adcroft]

P.,

I think the commenting of the IF statement arose because I was testing the
on-tile diagnostics and writing the deformation/strain fields to files (see
lower in the code). The IF should be uncommented but also an IF added to the
diagnostics lines...

A.
--
Dr Alistair Adcroft            http://www.mit.edu/~adcroft
MIT Climate Modeling Initiative        tel: (617) 253-5938
EAPS 54-1624,  77 Massachusetts Ave,  Cambridge,  MA,  USA

-----Original Message-----
From: mitgcm-devel-bounces at mitgcm.org
[mailto:mitgcm-devel-bounces at mitgcm.org] On Behalf Of Patrick Heimbach
Sent: Friday, July 16, 2004 12:40 PM
To: mitgcm-devel at mitgcm.org
Subject: [MITgcm-devel] momentum ...


Alistair, I have a few questions.

In the recent code S/R mom_fluxform has the conditions

IF (viscAstrain.NE.0. .OR. viscAtension.NE.0.) THEN

of calls to S/R MOM_CALC_TENSION, MOM_CALC_STRAIN
commented.
That IF-statement is not commented in S/R mom_vecinv.
Is there a special reason for this difference?

The other issue is that we are running a
1/4 degree setup in the Mediterranean.
Model seems to run ok if using 
vectorinvariantmomentum=.FALSE.
but crashes after about 1 month for vectorinvariantmomentum=.TRUE. Biharm.
friction in both cases is viscA4=8.E10 (typical value for this setup). Any
idea what could go wrong? Different treatment of horiz. viscosities?
Different treatment of noslip conditions?

Cheers
-Patrick


--------------------------------------------------------
Patrick Heimbach   Massachusetts Institute of Technology
FON: +1/617/253-5259                  EAPS, Room 54-1518
FAX: +1/617/253-4464             77 Massachusetts Avenue
mailto:heimbach at mit.edu               Cambridge MA 02139
http://www.mit.edu/~heimbach/                        USA

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