[MITgcm-devel] RE: Higher order dissip on the cube
Chris Hill
cnh at mit.edu
Mon Feb 9 09:17:17 EST 2004
All,
Here is the latest from the road. I created a _rx form for the exch2_
changes and tried to undo as many of my debugging steps. Unfortunately I
didn't think to keep a separate directory of CVS checkout so I'm not sure I
rembered everything I had done. Nevertheless here is what it looks like:
1. exch2_uv needs the wider loop mods that were in the exch2_uv_xyz_rl.F I
posted to JM and devel yesterday.
2. to exactly reproduce the exch1_ behavior with wider loops it also needs
to have the null region zeroing that is in that code too.
3. with 1. and 2. I can get the exact same numbers with exch1_ or exch2_
provided I turn off non-linear free surface stuff.
4. when I turn on non-lin free surface everything breaks!
The slightly weird thing I haven't resolved is
o with exch2_uv_xy_rl.F I seem to get a different way of crashing with
non-lin free surf than with exch1 even though
I have written out a full 2d field after the two exchanges and compared
them in octave and they appear the same.
JM/Alistair can you get Andy to double check 3. and 4. and then help with
figuring out whether it's a problem in calc_rstar or a problem in
exch2_uv_xy
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 2:01 PM
> To: MITgcm-devel at mitgcm.org
> Subject: RE: [MITgcm-devel] RE: Higher order dissip on the cube
>
> 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
>
> _______________________________________________
> 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