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

[Jean-Michel Campin]

Hi Gael (and David),

Looks like we have a problem in ggl90_calc.F
(vermix testing on aces with open64 + travuv), 
I am getting a floating exception
in line 570:
      tmpVisc = MIN(tmpVisc/TKEPrandtlNumber(i,j,k),GGL90diffMax)
for i,j,k = 0,0,1
because tmpVisc is NaN 

If you can fix it quickly, or we can take a look on Monday.

Cheers,
Jean-Michel

On Fri, Aug 06, 2010 at 02:37:05PM -0400, Gael Forget wrote:
> Update of /u/gcmpack/MITgcm/pkg/ggl90
> In directory faulks.csail.mit.edu:/u/u2/gforget/MITgcm/pkg/ggl90
> 
> Modified Files:
> 	GGL90.h GGL90_TAVE.h ggl90_calc.F ggl90_calc_diff.F 
> 	ggl90_calc_visc.F ggl90_diagnostics_init.F ggl90_output.F 
> Log Message:
> Minor changes in pkg/ggl90:
>  o GGL90diffKrS was removed --> always use GGL90diffKr
>  o GGL90viscAr was removed --> replaced with GGL90viscArU, GGL90viscArV
>  o hack of mxlMaxFlag=2 --> ensure mixing between first and second level (commented out for now)
>  o change in max/min operations to ensure that smoothing is ok
>  o smoothing of GGL90viscAr was moved to ggl90_calc.F (as done for GGL90diffKr)
>  o always use diffKrNrT as background profile (i.e. never use diffKr field)
> 
> 
> _______________________________________________
> MITgcm-cvs mailing list
> MITgcm-cvs at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-cvs