/[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.8 by jmc, Tue Oct 19 02:42:04 2004 UTC revision 1.9 by jmc, Mon Jan 31 19:37:06 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 tmpFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
93          LOGICAL  DIAGNOSTICS_IS_ON
94          EXTERNAL DIAGNOSTICS_IS_ON      
95    #endif
96    
97        LOGICAL dBug        LOGICAL dBug
98    
# Line 147  C--    Update Sea-Ice state : Line 152  C--    Update Sea-Ice state :
152            Tice2(i,j,bi,bj)=Tice(2)            Tice2(i,j,bi,bj)=Tice(2)
153            Qice1(i,j,bi,bj)=qicen(1)            Qice1(i,j,bi,bj)=qicen(1)
154            Qice2(i,j,bi,bj)=qicen(2)            Qice2(i,j,bi,bj)=qicen(2)
155  #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*/  
156            IF ( dBug ) THEN            IF ( dBug ) THEN
157             WRITE(6,1010) 'ThSI_FWD: Tsf, Tice(1,2), frzmltMxL =',             WRITE(6,1010) 'ThSI_FWD: Tsf, Tice(1,2), frzmltMxL =',
158       &                              Tsf, Tice, frzmltMxL       &                              Tsf, Tice, frzmltMxL
# Line 164  C--    Update Sea-Ice state : Line 166  C--    Update Sea-Ice state :
166        ENDIF        ENDIF
167        dBug = .FALSE.        dBug = .FALSE.
168    
169    #ifdef ALLOW_DIAGNOSTICS
170          IF ( useDiagnostics ) THEN
171    
172           IF ( DIAGNOSTICS_IS_ON('SIsnwPrc',myThid) ) THEN
173            DO j=1,sNy
174             DO i=1,sNx
175              tmpFld(i,j) = iceMask(i,j,bi,bj)*snowPrc(i,j,bi,bj)
176             ENDDO
177            ENDDO
178            CALL DIAGNOSTICS_FILL(tmpFld,'SIsnwPrc',0,1,2,bi,bj,myThid)
179           ENDIF
180    
181          ENDIF
182    #endif /* ALLOW_DIAGNOSTICS */
183    
184  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
185  C    part.2 : ice-covered fraction ;  C    part.2 : ice-covered fraction ;
186  C     change in ice/snow thickness and ice-fraction  C     change in ice/snow thickness and ice-fraction
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22