[MITgcm-devel] RE: Higher order dissip on the cube

Alistair Adcroft adcroft at MIT.EDU
Sun Feb 8 16:00:44 EST 2004


C., Can you try down-grading the mon_stats*.F to Thursday - I modified them
on Friday.

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 Chris Hill
Sent: Sunday, February 08, 2004 2:42 PM
To: MITgcm-devel at mitgcm.org
Subject: RE: [MITgcm-devel] RE: Higher order dissip on the cube


I get the monitor NAN's too on a checkout from yesterday and running
global_ocean.cs single process. I had thought it might be something I had
added locally so I was ignoring them, but maybe it eas in the main branch I
checked out. Note - the model runs ok despite the monitor NAN's e.g.
cg2d_init residuals are fine.

chris

> -----Original Message-----
> From: mitgcm-devel-bounces at mitgcm.org 
> [mailto:mitgcm-devel-bounces at mitgcm.org] On Behalf Of Alistair Adcroft
> Sent: Sunday, February 08, 2004 1:41 PM
> To: MITgcm-devel at mitgcm.org
> Subject: RE: [MITgcm-devel] RE: Higher order dissip on the cube
> 
> 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
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-devel
> 

_______________________________________________
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