/[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.10 by edhill, Thu Oct 9 04:19:18 2003 UTC revision 1.11 by heimbach, Thu Nov 6 22:01:43 2003 UTC
# Line 60  C     max number of different water type Line 60  C     max number of different water type
60        _RL facz        _RL facz
61        _RL rfac(nwtype),a1(nwtype),a2(nwtype)        _RL rfac(nwtype),a1(nwtype),a2(nwtype)
62        INTEGER i        INTEGER i
63  #ifdef ALLOW_CALENDAR  #ifdef ALLOW_CAL
64        _RL     fac        _RL     fac
65        logical first, changed        logical first, changed
66        integer count0, count1        integer count0, count1
67        integer myiter        integer myiter
68        integer jerl(12)        integer jerl(12)
69        data jerl / 2 , 2 , 2 , 3 , 3 , 3 , 4 , 4 , 4 , 4 , 3 , 2 /        data jerl / 2 , 2 , 2 , 3 , 3 , 3 , 4 , 4 , 4 , 4 , 3 , 2 /
70  #endif /* ALLOW_CALENDAR */  #endif /* ALLOW_CAL */
71  C  C
72  C     Jerlov water type :  I       IA      IB      II      III  C     Jerlov water type :  I       IA      IB      II      III
73  C                jwtype    1       2       3       4       5  C                jwtype    1       2       3       4       5
# Line 77  C Line 77  C
77        DATA a2           / 23.0  , 20.0  , 17.0  , 14.0  ,  7.9  /        DATA a2           / 23.0  , 20.0  , 17.0  , 14.0  ,  7.9  /
78  CEOP  CEOP
79    
80  #ifdef ALLOW_CALENDAR  #ifdef ALLOW_CAL
81  ceh3 this should have an IF ( useCALENDAR ) THEN  ceh3 this should have an IF ( useCALENDAR ) THEN
82        myiter=0        myiter=0
83        call  cal_GetMonthsRec(        call  cal_GetMonthsRec(
84       O     fac, first, changed, count0, count1,       O     fac, first, changed, count0, count1,
85       I     mytime, myiter, mythid )       I     mytime, myiter, mythid )
86        jwtype=jerl(count0)        jwtype=jerl(count0)
87  #else /* ALLOW_CALENDAR undef */  #else /* ALLOW_CAL undef */
88        jwtype=2        jwtype=2
89  #endif /* ALLOW_CALENDAR */  #endif /* ALLOW_CAL */
90    
91        DO i = 1,imax        DO i = 1,imax
92           facz = fact*swdk(i)           facz = fact*swdk(i)
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22