[MITgcm-devel] RE: Higher order dissip on the cube
Alistair Adcroft
adcroft at MIT.EDU
Mon Feb 9 09:14:38 EST 2004
I'm confused. To clarify, we have two outstanding issues:
1) NaN's appearing in the monitor output
2) An inconsistency between exch1 and ech2 when using higher order
dissipation
I'm not sure that they are related, though Chris's last email to Andy
suggests otherwise. Am I missing something?
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: Monday, February 09, 2004 9:06 AM
To: MITgcm-devel at mitgcm.org
Subject: RE: [MITgcm-devel] RE: Higher order dissip on the cube
Andy,
Things have moved on - see next message.
Chris
> -----Original Message-----
> From: mitgcm-devel-bounces at mitgcm.org
> [mailto:mitgcm-devel-bounces at mitgcm.org] On Behalf Of Andrew Eichmann
> Sent: Sunday, February 08, 2004 3:46 PM
> To: MITgcm-devel at mitgcm.org
> Subject: Re: [MITgcm-devel] RE: Higher order dissip on the cube
>
> Do you still want the exchanges checked like you mentioned on Friday?
>
> Andy
>
>
> * Chris Hill <cnh at mit.edu> wrote on Sun, Feb 08, 2004 at 02:41:31PM
> -0500
> > 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
> _______________________________________________
> 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