/[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.53 by dimitri, Thu Oct 23 07:14:49 2003 UTC revision 1.54 by edhill, Fri Oct 24 05:29:35 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    
422  #ifdef ALLOW_THERM_SEAICE  #ifdef ALLOW_THERM_SEAICE
423         IF (useThermSeaIce) THEN         IF (useThermSeaIce) THEN
424  #ifndef DISABLE_DEBUGMODE  #ifndef DISABLE_DEBUGMODE
# Line 429  C       including effects from ice Line 430  C       including effects from ice
430          CALL ICE_FORCING(          CALL ICE_FORCING(
431       I             bi, bj, iMin, iMax, jMin, jMax,       I             bi, bj, iMin, iMax, jMin, jMax,
432       I             myThid )       I             myThid )
433           ELSE
434    #endif /* ALLOW_THERM_SEAICE */
435    
436    C--     Determines forcing terms based on external fields
437    C       relaxation terms, etc.
438    #ifndef DISABLE_DEBUGMODE
439            IF ( debugLevel .GE. debLevB )
440         &    CALL DEBUG_CALL('EXTERNAL_FORCING_SURF',myThid)
441    #endif
442            CALL EXTERNAL_FORCING_SURF(
443         I             bi, bj, iMin, iMax, jMin, jMax,
444         I             myTime, myIter, myThid )
445    
446    #ifdef ALLOW_THERM_SEAICE
447    C--    end of if/else block useThermSeaIce --
448         ENDIF         ENDIF
449  #endif /* ALLOW_THERM_SEAICE */  #endif /* ALLOW_THERM_SEAICE */
450    
# Line 719  C--      Apply open boundary conditions Line 735  C--      Apply open boundary conditions
735           END IF           END IF
736  #endif   /* ALLOW_OBCS */  #endif   /* ALLOW_OBCS */
737    
738    C--      Freeze water
739             IF ( allowFreezing .AND. .NOT. useSEAICE
740         &       .AND. .NOT.(useThermSeaIce.AND.k.EQ.1) ) THEN
741    #ifdef ALLOW_AUTODIFF_TAMC
742    CADJ STORE gT(:,:,k,bi,bj) = comlev1_bibj_k
743    CADJ &   , key = kkey, byte = isbyte
744    #endif /* ALLOW_AUTODIFF_TAMC */
745                CALL FREEZE( bi, bj, iMin, iMax, jMin, jMax, k, myThid )
746             ENDIF
747    
748  C--     end of thermodynamic k loop (Nr:1)  C--     end of thermodynamic k loop (Nr:1)
749          ENDDO          ENDDO
750    
Legend:
Removed from v.1.53  
changed lines
  Added in v.1.54
  ViewVC Help
Powered by ViewVC 1.1.22