/[MITgcm]/MITgcm_contrib/torge/itd/code/seaice_model.F
ViewVC logotype

Diff of /MITgcm_contrib/torge/itd/code/seaice_model.F

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

revision 1.10 by torge, Tue Oct 30 21:32:54 2012 UTC revision 1.11 by torge, Wed Mar 27 18:59:53 2013 UTC
# Line 64  C     msgBuf      :: Informational/error Line 64  C     msgBuf      :: Informational/error
64  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
65        INTEGER IT        INTEGER IT
66  #endif  #endif
67  C#if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC)  #if defined(ALLOW_AUTODIFF_TAMC) || defined(SEAICE_ITD)
 #if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC) || defined(SEAICE_ITD)  
68  C>>>ToM  C>>>ToM
69        INTEGER i, j, bi, bj        INTEGER i, j, bi, bj
70  #endif  #endif
# Line 88  C--   Map thSice-variables to HEFF and A Line 87  C--   Map thSice-variables to HEFF and A
87        ENDIF        ENDIF
88  #endif /* ALLOW_THSICE */  #endif /* ALLOW_THSICE */
89    
 #ifdef SEAICE_GROWTH_LEGACY  
       IF ( .NOT.useThSice ) THEN  
 #ifdef ALLOW_AUTODIFF_TAMC  
 CADJ STORE heff  = comlev1, key=ikey_dynamics, kind=isbyte  
 CADJ STORE heffm = comlev1, key=ikey_dynamics, kind=isbyte  
 CADJ STORE area  = comlev1, key=ikey_dynamics, kind=isbyte  
 CADJ STORE hsnow = comlev1, key=ikey_dynamics, kind=isbyte  
 CADJ STORE tice  = comlev1, key=ikey_dynamics, kind=isbyte  
 #ifdef SEAICE_VARIABLE_SALINITY  
 CADJ STORE hsalt = comlev1, key=ikey_dynamics, kind=isbyte  
 #endif  
 #endif  
       DO bj=myByLo(myThid),myByHi(myThid)  
        DO bi=myBxLo(myThid),myBxHi(myThid)  
         DO j=1-OLy,sNy+OLy  
          DO i=1-OLx,sNx+OLx  
           IF ( (heff(i,j,bi,bj).EQ.0.)  
      &     .OR.(area(i,j,bi,bj).EQ.0.)  
      &     ) THEN  
            HEFF(i,j,bi,bj) = 0. _d 0  
            AREA(i,j,bi,bj) = 0. _d 0  
            HSNOW(i,j,bi,bj) = 0. _d 0  
            TICE(i,j,bi,bj) = celsius2K  
 #ifdef SEAICE_VARIABLE_SALINITY  
            HSALT(i,j,bi,bj) = 0. _d 0  
 #endif  
           ENDIF  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
       ENDIF  
 #endif  
   
90  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
91        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
92         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
93          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
94           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
 # ifdef SEAICE_GROWTH_LEGACY  
           areaNm1(i,j,bi,bj) = 0. _d 0  
           hEffNm1(i,j,bi,bj) = 0. _d 0  
 # endif  
95            uIceNm1(i,j,bi,bj) = 0. _d 0            uIceNm1(i,j,bi,bj) = 0. _d 0
96            vIceNm1(i,j,bi,bj) = 0. _d 0            vIceNm1(i,j,bi,bj) = 0. _d 0
97  # ifdef ALLOW_SITRACER  # ifdef ALLOW_SITRACER
# Line 200  C--   Apply ice velocity open boundary c Line 161  C--   Apply ice velocity open boundary c
161  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
162    
163  #ifdef ALLOW_THSICE  #ifdef ALLOW_THSICE
164        IF ( .NOT.useThSice ) THEN        IF ( useThSice ) THEN
165    #ifndef OLD_THSICE_CALL_SEQUENCE
166    #ifdef ALLOW_DEBUG
167            IF (debugMode) CALL DEBUG_CALL( 'THSICE_DO_ADVECT', myThid )
168    #endif
169            CALL THSICE_DO_ADVECT( 0, 0, myTime, myIter, myThid )
170    #endif /* OLD_THSICE_CALL_SEQUENCE */
171          ELSE
172  #endif  #endif
173  C--   Only call advection of heff, area, snow, and salt and  C--   Only call advection of heff, area, snow, and salt and
174  C--   growth for the generic 0-layer thermodynamics of seaice  C--   growth for the generic 0-layer thermodynamics of seaice
# Line 292  C     to match single category values Line 260  C     to match single category values
260         ENDDO         ENDDO
261  #endif  #endif
262  C>>>ToM  C>>>ToM
 #ifdef SEAICE_GROWTH_LEGACY  
       ELSE  
        DO bj=myByLo(myThid),myByHi(myThid)  
         DO bi=myBxLo(myThid),myBxHi(myThid)  
          DO j=1-OLy,sNy+OLy  
           DO i=1-OLx,sNx+OLx  
            areaNm1(i,j,bi,bj) = AREA(i,j,bi,bj)  
            hEffNm1(i,j,bi,bj) = HEFF(i,j,bi,bj)  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
 #endif /* SEAICE_GROWTH_LEGACY */  
263        ENDIF        ENDIF
264  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
265  CADJ STORE heffm  = comlev1, key=ikey_dynamics, kind=isbyte  CADJ STORE heffm  = comlev1, key=ikey_dynamics, kind=isbyte
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22