--- MITgcm/pkg/dic/dic_surfforcing.F	2011/04/19 21:34:32	1.25
+++ MITgcm/pkg/dic/dic_surfforcing.F	2011/05/05 21:31:58	1.26
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/dic_surfforcing.F,v 1.25 2011/04/19 21:34:32 stephd Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/dic_surfforcing.F,v 1.26 2011/05/05 21:31:58 stephd Exp $
 C $Name:  $
 
 #include "DIC_OPTIONS.h"
@@ -162,6 +162,8 @@
      &          + sca3 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
      &          + sca4 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
      &                *theta(i,j,kLev,bi,bj)
+c make sure Schmidt number isn't negative (will happen if temp>39C)
+             SchmidtNoDIC(i,j)=max(1.0 _d -2, SchmidtNoDIC(i,j))
 
 C Determine surface flux (FDIC)
 C first correct pCO2at for surface atmos pressure