--- MITgcm/verification/hs94.128x64x5/code/external_forcing.F	2001/05/29 14:01:57	1.4
+++ MITgcm/verification/hs94.128x64x5/code/external_forcing.F	2001/06/04 20:30:49	1.5
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/verification/hs94.128x64x5/code/external_forcing.F,v 1.4 2001/05/29 14:01:57 adcroft Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/verification/hs94.128x64x5/code/external_forcing.F,v 1.5 2001/06/04 20:30:49 adcroft Exp $
 C $Name:  $
 
 #include "CPP_OPTIONS.h"
@@ -56,7 +56,8 @@
 C        termP=0.5*( rF(kLev) + min( rF(kLev+1) ,
 C    &           min(H(I,J,bi,bj),H(I,J-1,bi,bj))            ) )
          termP=0.5*( rF(kLev) + rF(kLev+1) )
-         kV=kF*MAX(0., (termP*recip_Rcol(I,J,bi,bj)-0.7)/(1.-0.7) )
+         kV=kF*MAX(0. _d 0,
+     &   (termP*recip_Rcol(I,J,bi,bj)-0.7 _d 0)/(1. _d 0-0.7 _d 0) )
          gU(i,j,kLev,bi,bj)=gU(i,j,kLev,bi,bj)
      &                      -kV*uVel(i,j,kLev,bi,bj)
         ENDIF
@@ -118,7 +119,8 @@
 C        termP=0.5*( rF(kLev) + min( rF(kLev+1) ,
 C    &           min(H(I,J,bi,bj),H(I,J-1,bi,bj))            ) )
          termP=0.5*( rF(kLev) + rF(kLev+1) )
-         kV=kF*MAX(0., (termP*recip_Rcol(I,J,bi,bj)-0.7)/(1.-0.7) )
+         kV=kF*MAX(0. _d 0,
+     &   (termP*recip_Rcol(I,J,bi,bj)-0.7 _d 0)/(1. _d 0-0.7 _d 0) )
          gV(i,j,kLev,bi,bj)=gV(i,j,kLev,bi,bj)
      &                      -kV*vVel(i,j,kLev,bi,bj)
         ENDIF
@@ -165,21 +167,22 @@
       _RL thetaLim,kT,ka,ks,term1,term2,thetaEq,termP,rSurf
 
 C--   Forcing term(s)
-      rSurf=1.E5
-      ka=1./(40.*86400.)
-      ks=1./(4. *86400.)
+      rSurf=1. _d 05
+      ka=1. _d 0/(40. _d 0*86400. _d 0)
+      ks=1. _d 0/(4. _d 0 *86400. _d 0)
       DO J=jMin,jMax
-       term1=60.*(sin(yC(1,J,bi,bj)*deg2rad)**2)
+       term1=60. _d 0*(sin(yC(1,J,bi,bj)*deg2rad)**2)
 C      termP=0.5*( rF(kLev) + min( rF(kLev+1) , H(I,J,bi,bj) ) )
        termP=0.5*( rF(kLev) + rF(kLev+1) )
-       term2=10.*log(termP/rSurf)
+       term2=10. _d 0*log(termP/rSurf)
      &          *(cos(yC(1,J,bi,bj)*deg2rad)**2)
-       thetaLim = 200. / ((termP/rSurf)**(2./7.))
+       thetaLim = 200. _d 0/((termP/rSurf)**(2. _d 0/7. _d 0))
        thetaEq=315.-term1-term2
        thetaEq=MAX(thetaLim,thetaEq)
        DO I=iMin,iMax
         kT=ka+(ks-ka)
-     &    *MAX(0., (termP*recip_Rcol(I,J,bi,bj)-0.7)/(1.-0.7) )
+     &    *MAX(0.,
+     &    (termP*recip_Rcol(I,J,bi,bj)-0.7 _d 0)/(1. _d 0-0.7 _d 0) )
      &    *COS((yC(1,J,bi,bj)*deg2rad))**4
          gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)
      &        - kT*( theta(I,J,kLev,bi,bj)-thetaEq )