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

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

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

revision 1.7 by jmc, Thu Jul 22 22:52:59 2004 UTC revision 1.10 by jmc, Sun Jul 10 01:24:25 2005 UTC
# Line 88  C     TFrzOce   :: sea-water freezing te Line 88  C     TFrzOce   :: sea-water freezing te
88        _RL oceV2s, oceTs        _RL oceV2s, oceTs
89        _RL compact, hIce, hSnow, Tsf, Tice(nlyr), qicen(nlyr)        _RL compact, hIce, hSnow, Tsf, Tice(nlyr), qicen(nlyr)
90        _RL tmpflx(0:2), tmpdTs        _RL tmpflx(0:2), tmpdTs
91    #ifdef ALLOW_DIAGNOSTICS
92          _RL tmpFac
93    #endif
94    
95        LOGICAL dBug        LOGICAL dBug
96    
# Line 100  C-    Initialise flxAtm Line 103  C-    Initialise flxAtm
103          ENDDO          ENDDO
104         ENDDO         ENDDO
105    
106        IF ( buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN        IF ( fluidIsWater ) THEN
107         DO j = jMin, jMax         DO j = jMin, jMax
108          DO i = iMin, iMax          DO i = iMin, iMax
109  c        dBug = ( bi.EQ.3 .AND. i.EQ.15 .AND. j.EQ.11 )  c        dBug = ( bi.EQ.3 .AND. i.EQ.15 .AND. j.EQ.11 )
# Line 147  C--    Update Sea-Ice state : Line 150  C--    Update Sea-Ice state :
150            Tice2(i,j,bi,bj)=Tice(2)            Tice2(i,j,bi,bj)=Tice(2)
151            Qice1(i,j,bi,bj)=qicen(1)            Qice1(i,j,bi,bj)=qicen(1)
152            Qice2(i,j,bi,bj)=qicen(2)            Qice2(i,j,bi,bj)=qicen(2)
153  #ifdef ALLOW_TIMEAVE            siceAlb(i,j,bi,bj) = icFrac*albedo
           ice_albedo_Ave(i,j,bi,bj) = ice_albedo_Ave(i,j,bi,bj)  
      &                              + icFrac*albedo*thSIce_deltaT  
 #endif /*ALLOW_TIMEAVE*/  
154            IF ( dBug ) THEN            IF ( dBug ) THEN
155             WRITE(6,1010) 'ThSI_FWD: Tsf, Tice(1,2), frzmltMxL =',             WRITE(6,1010) 'ThSI_FWD: Tsf, Tice(1,2), frzmltMxL =',
156       &                              Tsf, Tice, frzmltMxL       &                              Tsf, Tice, frzmltMxL
# Line 164  C--    Update Sea-Ice state : Line 164  C--    Update Sea-Ice state :
164        ENDIF        ENDIF
165        dBug = .FALSE.        dBug = .FALSE.
166    
167    #ifdef ALLOW_DIAGNOSTICS
168          IF ( useDiagnostics ) THEN
169            tmpFac = 1. _d 0
170            CALL DIAGNOSTICS_FRACT_FILL(
171         I                   snowPrc,   iceMask,tmpFac,'SIsnwPrc',
172         I                   0,1,1,bi,bj,myThid)
173          ENDIF
174    #endif /* ALLOW_DIAGNOSTICS */
175    
176  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
177  C    part.2 : ice-covered fraction ;  C    part.2 : ice-covered fraction ;
178  C     change in ice/snow thickness and ice-fraction  C     change in ice/snow thickness and ice-fraction
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22