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

Jean-Michel Campin jmc at ocean.mit.edu
Tue Feb 14 19:19:55 EST 2012


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



More information about the MITgcm-devel mailing list