/[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.9 by torge, Fri Oct 26 15:07:52 2012 UTC revision 1.10 by torge, Tue Oct 30 21:32:54 2012 UTC
# Line 63  C     msgBuf      :: Informational/error Line 63  C     msgBuf      :: Informational/error
63  #endif  #endif
64  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
65        INTEGER IT        INTEGER IT
 CToM a constant divergence rate for an otherwise purely thermodynamic 1-D test case (not a permanent addition)  
       _RL divergence  
66  #endif  #endif
67  C#if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC)  C#if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC)
68  #if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC) || defined(SEAICE_ITD)  #if defined(SEAICE_GROWTH_LEGACY) || defined(ALLOW_AUTODIFF_TAMC) || defined(SEAICE_ITD)
# Line 313  CADJ STORE heffm  = comlev1, key=ikey_dy Line 311  CADJ STORE heffm  = comlev1, key=ikey_dy
311  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
312    
313  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
 c      create some open water  
        divergence = 1.0/(5*365.*86400./SEAICE_deltaTtherm)  
        DO bj=myByLo(myThid),myByHi(myThid)  
         DO bi=myBxLo(myThid),myBxHi(myThid)  
          DO IT=1,nITD  
           DO J=1-OLy,sNy+OLy  
            DO I=1-OLx,sNx+OLx  
             AREAITD(I,J,IT,bi,bj)=MAX(0.0,AREAITD(I,J,IT,bi,bj)  
      &                           -divergence)  
             HEFFITD(I,J,IT,bi,bj)=HEFFITD(I,J,IT,bi,bj)  
      &                           -divergence*HEFFITD(I,J,IT,bi,bj)  
            ENDDO  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
314  #ifdef SEAICE_DEBUG  #ifdef SEAICE_DEBUG
315  C     ToM: generate some test output  C     ToM: generate some test output
316         WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F8.4)'         WRITE(HlimitMsgFormat,'(A,I2,A)') '(A,',nITD,'F8.4)'
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22