[MITgcm-devel] [MITgcm-cvs] MITgcm/pkg/seaice CVS Commit

[Jean-Michel Campin]

Martin,

The 1rst Pb should be fixed. 
And regarding AD tests, the ones that fail are:
Y Y Y Y  8>11< 2 pass  global_ocean.cs32x15.seaice_dynmix
Y Y Y Y 12> 7< 4 FAIL  lab_sea
Y Y Y Y  9> 5< 2 FAIL  lab_sea.evp
Y Y Y Y 12> 7< 4 FAIL  lab_sea.noseaicedyn
I am going to update the corresponding output_adm.*txt,
since it's likely a machine truncation problem.
And if not, we can put back the old results.

Jean-Michel

On Tue, Feb 14, 2012 at 07:19:55PM -0500, Jean-Michel Campin wrote:
> Hi Martin,
> 
> I am getting a floating point exception when trying to run
> 1D_ocean_ice_column (fwd), and also global_ocean.cs32x15.seaice (fwd)
> from the line you added:
> >      recip_HCUT = 1./HCUT
> because HCUT is zero.
> 
> Also, it looks like the changes you made change the way
> machine truncate the number so that a couple of adjoint verification
> experiment are affected and don't pass anymore (13 digits with
> gfortran on baudelaire). I am now checking this directly on baudelaire
> (was trying before on my laptop and could be due to something else).
> 
> Cheers,
> Jean-Michel
> 
> On Tue, Feb 14, 2012 at 09:52:09AM -0500, Martin Losch wrote:
> > Update of /u/gcmpack/MITgcm/pkg/seaice
> > In directory forge:/tmp/cvs-serv1373/pkg/seaice
> > 
> > Modified Files:
> > 	seaice_solve4temp.F 
> > Log Message:
> > remove one divisions by HCUT within a loop, without changing the results
> > 
> > 
> > _______________________________________________
> > MITgcm-cvs mailing list
> > MITgcm-cvs at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-cvs
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel