/[MITgcm]/MITgcm/model/src/swfrac.F
ViewVC logotype

Diff of /MITgcm/model/src/swfrac.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.14 by jmc, Tue Sep 4 16:49:44 2007 UTC revision 1.15 by jmc, Tue Aug 4 21:51:53 2009 UTC
# Line 69  c     INTEGER count0, count1 Line 69  c     INTEGER count0, count1
69  c     INTEGER jerl(12)  c     INTEGER jerl(12)
70  c     DATA jerl / 2 , 2 , 2 , 3 , 3 , 3 , 4 , 4 , 4 , 4 , 3 , 2 /  c     DATA jerl / 2 , 2 , 2 , 3 , 3 , 3 , 4 , 4 , 4 , 4 , 3 , 2 /
71  #endif /* ALLOW_CAL */  #endif /* ALLOW_CAL */
 C  
 C     Jerlov water type :  I       IA      IB      II      III  
 C                jwtype    1       2       3       4       5  
 C  
       DATA rfac         /  0.58 ,  0.62 ,  0.67 ,  0.77 ,  0.78 /  
       DATA a1           /  0.35 ,  0.6  ,  1.0  ,  1.5  ,  1.4  /  
       DATA a2           / 23.0  , 20.0  , 17.0  , 14.0  ,  7.9  /  
72  CEOP  CEOP
73    
74    C     Jerlov water type :
75    C                  I          IA         IB         II         III
76    C     jwtype :     1          2          3          4          5
77          DATA rfac / 0.58 _d 0, 0.62 _d 0, 0.67 _d 0, 0.77 _d 0, 0.78 _d 0/
78          DATA a1   / 0.35 _d 0, 0.6  _d 0, 1.0  _d 0, 1.5  _d 0, 1.4  _d 0/
79          DATA a2   / 23.0 _d 0, 20.0 _d 0, 17.0 _d 0, 14.0 _d 0, 7.9  _d 0/
80    
81  #ifdef ALLOW_CAL  #ifdef ALLOW_CAL
82  ceh3 this should have an IF ( useCALENDAR ) THEN  ceh3 this should have an IF ( useCALENDAR ) THEN
83  CML(  CML(
# Line 98  CML) Line 98  CML)
98  #endif /* ALLOW_CAL */  #endif /* ALLOW_CAL */
99    
100        DO i = 1,imax        DO i = 1,imax
101           facz = fact*swdk(i)          facz = fact*swdk(i)
102           IF (facz .LT. (-200.)) THEN          IF ( facz .LT. -200. _d 0 ) THEN
103              swdk(i) = 0.            swdk(i) = 0. _d 0
104           ELSE          ELSE
105              swdk(i) =       rfac(jwtype)  * exp(facz/a1(jwtype))            swdk(i) =       rfac(jwtype)  * exp( facz/a1(jwtype) )
106       $                + (1.-rfac(jwtype)) * exp(facz/a2(jwtype))       &       + (1. _d 0 - rfac(jwtype)) * exp( facz/a2(jwtype) )
107           ENDIF          ENDIF
108        ENDDO        ENDDO
109    
110        RETURN        RETURN
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22