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

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

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

revision 1.52 by heimbach, Fri Oct 10 22:56:08 2003 UTC revision 1.53 by dimitri, Thu Oct 23 07:14:49 2003 UTC
# Line 418  C--     Calculate future values on open Line 418  C--     Calculate future values on open
418          ENDIF          ENDIF
419  #endif  /* ALLOW_OBCS */  #endif  /* ALLOW_OBCS */
420    
   
421  #ifdef ALLOW_THERM_SEAICE  #ifdef ALLOW_THERM_SEAICE
422         IF (useThermSeaIce) THEN         IF (useThermSeaIce) THEN
423  #ifndef DISABLE_DEBUGMODE  #ifndef DISABLE_DEBUGMODE
# Line 430  C       including effects from ice Line 429  C       including effects from ice
429          CALL ICE_FORCING(          CALL ICE_FORCING(
430       I             bi, bj, iMin, iMax, jMin, jMax,       I             bi, bj, iMin, iMax, jMin, jMax,
431       I             myThid )       I             myThid )
        ELSE  
 #endif /* ALLOW_THERM_SEAICE */  
   
 C--     Determines forcing terms based on external fields  
 C       relaxation terms, etc.  
 #ifndef DISABLE_DEBUGMODE  
         IF ( debugLevel .GE. debLevB )  
      &    CALL DEBUG_CALL('EXTERNAL_FORCING_SURF',myThid)  
 #endif  
         CALL EXTERNAL_FORCING_SURF(  
      I             bi, bj, iMin, iMax, jMin, jMax,  
      I             myTime, myIter, myThid )  
   
 #ifdef ALLOW_THERM_SEAICE  
 C--    end of if/else block useThermSeaIce --  
432         ENDIF         ENDIF
433  #endif /* ALLOW_THERM_SEAICE */  #endif /* ALLOW_THERM_SEAICE */
434    
# Line 735  C--      Apply open boundary conditions Line 719  C--      Apply open boundary conditions
719           END IF           END IF
720  #endif   /* ALLOW_OBCS */  #endif   /* ALLOW_OBCS */
721    
 C--      Freeze water  
          IF ( allowFreezing .AND. .NOT. useSEAICE  
      &       .AND. .NOT.(useThermSeaIce.AND.k.EQ.1) ) THEN  
 #ifdef ALLOW_AUTODIFF_TAMC  
 CADJ STORE gT(:,:,k,bi,bj) = comlev1_bibj_k  
 CADJ &   , key = kkey, byte = isbyte  
 #endif /* ALLOW_AUTODIFF_TAMC */  
             CALL FREEZE( bi, bj, iMin, iMax, jMin, jMax, k, myThid )  
          ENDIF  
   
722  C--     end of thermodynamic k loop (Nr:1)  C--     end of thermodynamic k loop (Nr:1)
723          ENDDO          ENDDO
724    
Legend:
Removed from v.1.52  
changed lines
  Added in v.1.53
  ViewVC Help
Powered by ViewVC 1.1.22