[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