/[MITgcm]/MITgcm/pkg/thsice/thsice_main.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_main.F

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

revision 1.37 by jmc, Thu May 2 20:11:06 2013 UTC revision 1.38 by jmc, Tue Jun 11 01:48:22 2013 UTC
# Line 102  C     (excluding overlap) and apply exch Line 102  C     (excluding overlap) and apply exch
102  CADJ STORE ocefwfx(:,:,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte  CADJ STORE ocefwfx(:,:,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte
103  CADJ STORE oceqnet(:,:,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte  CADJ STORE oceqnet(:,:,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte
104  CADJ STORE ocesflx(:,:,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte  CADJ STORE ocesflx(:,:,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte
105  # ifdef ALLOW_EXF  CADJ STORE qsw(:,:,bi,bj)     = comlev1_bibj, key=ticekey, byte=isbyte
 CADJ STORE qsw(:,:,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte  
 # endif  
 #endif  
 #ifdef ALLOW_AUTODIFF_TAMC  
106  CADJ STORE uvel (:,:,1,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte  CADJ STORE uvel (:,:,1,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte
107  CADJ STORE vvel (:,:,1,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte  CADJ STORE vvel (:,:,1,bi,bj) = comlev1_bibj, key=ticekey, byte=isbyte
108  #endif  #endif
# Line 119  CADJ STORE vvel (:,:,1,bi,bj) = comlev1_ Line 115  CADJ STORE vvel (:,:,1,bi,bj) = comlev1_
115           ENDDO           ENDDO
116          ENDDO          ENDDO
117    
118  #ifdef ALLOW_CHEAPAML  #ifndef ALLOW_AUTODIFF_TAMC
119          IF ( .NOT.useCheapAML ) THEN          IF ( .NOT.useCheapAML ) THEN
120  #endif  #endif
121             CALL THSICE_GET_OCEAN(
122          CALL THSICE_GET_OCEAN(       I                          bi, bj, myTime, myIter, myThid )
      I                         bi, bj, myTime, myIter, myThid )  
123    
124  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
125  CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj, key = ticekey
# Line 139  CADJ STORE hOceMxL(:,:,bi,bj) = comlev1_ Line 134  CADJ STORE hOceMxL(:,:,bi,bj) = comlev1_
134  CADJ STORE tOceMxL(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE tOceMxL(:,:,bi,bj) = comlev1_bibj, key = ticekey
135  CADJ STORE sOceMxL(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE sOceMxL(:,:,bi,bj) = comlev1_bibj, key = ticekey
136  CADJ STORE v2ocMxL(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE v2ocMxL(:,:,bi,bj) = comlev1_bibj, key = ticekey
137    #else
138    C-   end if not useCheapAML
139            ENDIF
140  #endif  #endif
141    
142  #ifdef OLD_THSICE_CALL_SEQUENCE  #ifdef OLD_THSICE_CALL_SEQUENCE
# Line 149  C Note: will inline this S/R once thSIce Line 147  C Note: will inline this S/R once thSIce
147       I                   bi,bj, myTime, myIter, myThid )       I                   bi,bj, myTime, myIter, myThid )
148  #endif /* OLD_THSICE_CALL_SEQUENCE */  #endif /* OLD_THSICE_CALL_SEQUENCE */
149    
150  #ifdef ALLOW_BULK_FORCE  #ifndef ALLOW_AUTODIFF_TAMC
151          IF ( useBulkforce ) THEN          IF ( useBulkforce .OR. useCheapAML ) THEN
152           CALL THSICE_GET_PRECIP(           CALL THSICE_GET_PRECIP(
153       I                  iceMask, tOceMxL,       I                  iceMask, tOceMxL,
154       O                  prcAtm(1-OLx,1-OLy,bi,bj),       O                  prcAtm(1-OLx,1-OLy,bi,bj),
# Line 159  C Note: will inline this S/R once thSIce Line 157  C Note: will inline this S/R once thSIce
157       I                  iMin,iMax,jMin,jMax, bi,bj, myThid )       I                  iMin,iMax,jMin,jMax, bi,bj, myThid )
158          ENDIF          ENDIF
159  #endif  #endif
 #ifdef ALLOW_EXF  
160          IF ( useEXF ) THEN          IF ( useEXF ) THEN
161           CALL THSICE_MAP_EXF(           CALL THSICE_MAP_EXF(
162       I                  iceMask, tOceMxL,       I                  iceMask, tOceMxL,
# Line 168  C Note: will inline this S/R once thSIce Line 165  C Note: will inline this S/R once thSIce
165       O                  icFlxSW(1-OLx,1-OLy,bi,bj),       O                  icFlxSW(1-OLx,1-OLy,bi,bj),
166       I                  iMin,iMax,jMin,jMax, bi,bj, myThid )       I                  iMin,iMax,jMin,jMax, bi,bj, myThid )
167          ENDIF          ENDIF
 #endif  
168    
169  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
170  CADJ STORE sHeating(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE sHeating(:,:,bi,bj) = comlev1_bibj, key = ticekey
171  CADJ STORE tice1(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE tice1(:,:,bi,bj) = comlev1_bibj, key = ticekey
172  CADJ STORE tice2(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE tice2(:,:,bi,bj) = comlev1_bibj, key = ticekey
173  #else  #else
174          IF ( .NOT.thSIce_skipThermo ) THEN          IF ( .NOT.( useCheapAML .OR. thSIce_skipThermo ) ) THEN
175  #endif  #endif
176            CALL THSICE_STEP_TEMP(           CALL THSICE_STEP_TEMP(
177       I                     bi, bj, iMin, iMax, jMin, jMax,       I                     bi, bj, iMin, iMax, jMin, jMax,
178       I                     myTime, myIter, myThid )       I                     myTime, myIter, myThid )
179    
 #ifdef ALLOW_CHEAPAML  
         ENDIF  
 #endif  
180  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
181  CADJ STORE empmr(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE empmr(:,:,bi,bj) = comlev1_bibj, key = ticekey
182  CADJ STORE qnet(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE qnet(:,:,bi,bj) = comlev1_bibj, key = ticekey
# Line 196  CADJ STORE Qice2(:,:,bi,bj)   = comlev1_ Line 189  CADJ STORE Qice2(:,:,bi,bj)   = comlev1_
189  CADJ STORE snowAge(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE snowAge(:,:,bi,bj) = comlev1_bibj, key = ticekey
190  CADJ STORE sHeating(:,:,bi,bj) = comlev1_bibj, key = ticekey  CADJ STORE sHeating(:,:,bi,bj) = comlev1_bibj, key = ticekey
191  #else  #else
192    C-   end if not skipThermo / useCheapAML
193          ENDIF          ENDIF
194          IF ( .NOT.thSIce_skipThermo ) THEN          IF ( .NOT.thSIce_skipThermo ) THEN
195  #endif  #endif
196            CALL THSICE_STEP_FWD(           CALL THSICE_STEP_FWD(
197       I                     bi, bj, iMin, iMax, jMin, jMax,       I                     bi, bj, iMin, iMax, jMin, jMax,
198       I                     prcAtm(1-OLx,1-OLy,bi,bj),       I                     prcAtm(1-OLx,1-OLy,bi,bj),
199       I                     snowPr, qPrcRn,       I                     snowPr, qPrcRn,
Legend:
Removed from v.1.37  
changed lines
  Added in v.1.38
  ViewVC Help
Powered by ViewVC 1.1.22