/[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.11 by jmc, Mon Jul 11 19:00:34 2005 UTC revision 1.12 by jmc, Tue Feb 21 01:15:16 2006 UTC
# Line 131  C------- Line 131  C-------
131       O               albedo,       O               albedo,
132       I               myThid )       I               myThid )
133            flxSW(i,j) = flxSW(i,j)*(1. _d 0 - albedo)            flxSW(i,j) = flxSW(i,j)*(1. _d 0 - albedo)
134              siceAlb(i,j,bi,bj) = albedo
135    
136            CALL THSICE_SOLVE4TEMP(            CALL THSICE_SOLVE4TEMP(
137       I          useBulkforce, tmpflx, TFrzOce, hIce, hSnow,       I          useBulkforce, tmpflx, TFrzOce, hIce, hSnow,
# Line 150  C--    Update Sea-Ice state : Line 151  C--    Update Sea-Ice state :
151            Tice2(i,j,bi,bj)=Tice(2)            Tice2(i,j,bi,bj)=Tice(2)
152            Qice1(i,j,bi,bj)=qicen(1)            Qice1(i,j,bi,bj)=qicen(1)
153            Qice2(i,j,bi,bj)=qicen(2)            Qice2(i,j,bi,bj)=qicen(2)
           siceAlb(i,j,bi,bj) = icFrac*albedo  
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 170  C--    Update Sea-Ice state : Line 170  C--    Update Sea-Ice state :
170          CALL DIAGNOSTICS_FRACT_FILL(          CALL DIAGNOSTICS_FRACT_FILL(
171       I                   snowPrc,   iceMask,tmpFac,1,'SIsnwPrc',       I                   snowPrc,   iceMask,tmpFac,1,'SIsnwPrc',
172       I                   0,1,1,bi,bj,myThid)       I                   0,1,1,bi,bj,myThid)
173            CALL DIAGNOSTICS_FRACT_FILL(
174         I                   siceAlb,   iceMask,tmpFac,1,'SIalbedo',
175         I                   0,1,1,bi,bj,myThid)
176        ENDIF        ENDIF
177  #endif /* ALLOW_DIAGNOSTICS */  #endif /* ALLOW_DIAGNOSTICS */
178          DO j = jMin, jMax
179           DO i = iMin, iMax
180              siceAlb(i,j,bi,bj) = iceMask(i,j,bi,bj)*siceAlb(i,j,bi,bj)
181           ENDDO
182          ENDDO
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 ;
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22