[MITgcm-devel] RE: viscosity diagnostics

[Alistair Adcroft]

Andrea,

Think I found the problem - there isn't one!  I changed the freq from 10 to
2 and
The numbers are so small that they look like zeros in the netcdf file - set
the freq to 3 and I get 1.E-25 in a straight line that looks exactly like an
index shift. I think it's OK.

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 Alistair Adcroft
Sent: Tuesday, July 27, 2004 10:06 AM
To: MITgcm-devel at mitgcm.org
Cc: 'Andrea Molod'
Subject: [MITgcm-devel] RE: viscosity diagnostics


Andrea,

The shifting did not happen at one point in time. I have a version from
shortly after the first implementation of the viscA4 diagnostic that works
fine. I'll start from there and update incrementally until it breaks...

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: Andrea Molod [mailto:molod at ocean.mit.edu] 
Sent: Monday, July 26, 2004 6:47 PM
To: Alistair Adcroft
Subject: viscosity diagnostics


hi alistair,
have been playing with this diagnostics stuff a bit today.
i did see your 'shifting' behavior in your diagnostics, and started to try
to see what's going on. i made some temporary arrays in mom_vi_hdissip, of
all sorts of shapes and sizes to use the other 'flag's in fill_diagnostics,
and got the same behavior no matter what i did. then i tested the bibjflg 3
flag in the held-suarez with uVel, and it worked just fine. then i filled
the visc array with fake numbers, ones that should show the shift (i think),
like setting every i,j to the value of j, and it worked just fine. is it
possible that the visca4_d array HAS that shifted pattern in it? please say
yes........

if not, i will get busy again....

andrea

------------------------------------------------------------------
Andrea M. Molod, PhD.
Department of Earth, Atmospheric & Planetary Sciences 
MIT
Tel: (617) 253-5458, Email: amolod at mit.edu 
WWW:http://paoc.mit.edu/paoc/people/person.asp?position=Postdoc&who=molod
----------------------------------------------------------------


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