129 |
CADJ STORE icflxsw(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte |
CADJ STORE icflxsw(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte |
130 |
#endif |
#endif |
131 |
|
|
132 |
|
#ifdef ALLOW_DBUG_THSICE |
133 |
DO j = jMin, jMax |
DO j = jMin, jMax |
134 |
DO i = iMin, iMax |
DO i = iMin, iMax |
135 |
IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN |
IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN |
|
#ifdef ALLOW_DBUG_THSICE |
|
136 |
IF ( dBug(i,j,bi,bj) ) THEN |
IF ( dBug(i,j,bi,bj) ) THEN |
137 |
WRITE(6,'(A,2I4,2I2)') 'ThSI_STEP_T: i,j=',i,j,bi,bj |
WRITE(6,'(A,2I4,2I2)') 'ThSI_STEP_T: i,j=',i,j,bi,bj |
138 |
WRITE(6,1010) 'ThSI_STEP_T: iceMask, hIc, hSn, Tsf =', |
WRITE(6,1010) 'ThSI_STEP_T: iceMask, hIc, hSn, Tsf =', |
139 |
& iceMask(i,j,bi,bj), iceHeight(i,j,bi,bj), |
& iceMask(i,j,bi,bj), iceHeight(i,j,bi,bj), |
140 |
& snowHeight(i,j,bi,bj), Tsrf(i,j,bi,bj) |
& snowHeight(i,j,bi,bj), Tsrf(i,j,bi,bj) |
141 |
ENDIF |
ENDIF |
142 |
|
ENDIF |
143 |
|
ENDDO |
144 |
|
ENDDO |
145 |
#endif |
#endif |
146 |
|
DO j = jMin, jMax |
147 |
|
DO i = iMin, iMax |
148 |
|
IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN |
149 |
C- surface net SW flux: |
C- surface net SW flux: |
150 |
icFlxSW(i,j,bi,bj) = icFlxSW(i,j,bi,bj) |
icFlxSW(i,j,bi,bj) = icFlxSW(i,j,bi,bj) |
151 |
& *(1. _d 0 - siceAlb(i,j,bi,bj)) |
& *(1. _d 0 - siceAlb(i,j,bi,bj)) |