[MITgcm-devel] RE: Higher order dissip on the cube
Alistair Adcroft
adcroft at MIT.EDU
Sun Feb 8 13:40:46 EST 2004
I recently "upgraded" the way that global_sum() is used in mon_stats_r*.F -
it got past the verification experiments so not sure what's going on? Try
down grading mon_stats_*.F routines - maybe I messed up something that
involves MPI?
cvs update -D "5 days ago" pkg/monitor/mon_stats*.F
A.
--
Dr Alistair Adcroft http://www.mit.edu/~adcroft
MIT Climate Modeling Initiative tel: (617) 253-5938
EAPS 54-1523, 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 Dimitris Menemenlis
Sent: Sunday, February 08, 2004 1:27 PM
To: Alistair Adcroft
Cc: mitgcm
Subject: RE: [MITgcm-devel] RE: Higher order dissip on the cube
Alistair, some good news with "non-slippery" coasts.
The hi-res cube integration has gotten past the 1996
Equatorial Pacific sea-mount blow-up with viscA4Grid=1e-3
Anyone knows what happened to the monitor routines?
When I test-report global_ocean.cs32x15, I get a
bunch of NANs for mean and del2's and zeros for sd fields:
(PID.TID 0000.0001) %MON dynstat_eta_mean =
NAN
(PID.TID 0000.0001) %MON dynstat_eta_sd =
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_del2 =
NAN
(PID.TID 0000.0001) %MON dynstat_uvel_max =
1.8517049457442E-01
(PID.TID 0000.0001) %MON dynstat_uvel_min =
-2.8835551845614E-01
(PID.TID 0000.0001) %MON dynstat_uvel_mean =
NAN
(PID.TID 0000.0001) %MON dynstat_uvel_sd =
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_del2 =
NAN
(PID.TID 0000.0001) %MON dynstat_vvel_max =
1.7299533857874E-01
(PID.TID 0000.0001) %MON dynstat_vvel_min =
-2.2520988227009E-01
(PID.TID 0000.0001) %MON dynstat_vvel_mean =
NAN
(PID.TID 0000.0001) %MON dynstat_vvel_sd =
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_del2 =
NAN
(PID.TID 0000.0001) %MON dynstat_wvel_max =
1.0345009107981E-04
(PID.TID 0000.0001) %MON dynstat_wvel_min =
-2.0499843500452E-04
(PID.TID 0000.0001) %MON dynstat_wvel_mean =
NAN
(PID.TID 0000.0001) %MON dynstat_wvel_sd =
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_del2 =
NAN
(PID.TID 0000.0001) %MON dynstat_theta_max =
3.0924013283713E+01
(PID.TID 0000.0001) %MON dynstat_theta_min =
-1.9000000000000E+00
(PID.TID 0000.0001) %MON dynstat_theta_mean =
NAN
(PID.TID 0000.0001) %MON dynstat_theta_sd =
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_theta_del2 =
NAN
(PID.TID 0000.0001) %MON dynstat_salt_max =
4.1446610064801E+01
(PID.TID 0000.0001) %MON dynstat_salt_min =
1.8388095152689E+01
(PID.TID 0000.0001) %MON dynstat_salt_mean =
NAN
(PID.TID 0000.0001) %MON dynstat_salt_sd =
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_salt_del2 =
NAN
--
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Insitute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656; fax: 818-393-6720
_______________________________________________
MITgcm-devel mailing list
MITgcm-devel at mitgcm.org http://dev.mitgcm.org/mailman/listinfo/mitgcm-devel
More information about the MITgcm-devel
mailing list