69 |
_RL HCORR (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL HCORR (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
70 |
C SEAICE_SALT contains m of ice melted (<0) or created (>0) |
C SEAICE_SALT contains m of ice melted (<0) or created (>0) |
71 |
_RL SEAICE_SALT(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL SEAICE_SALT(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
72 |
C actual ice thickness |
C actual ice thickness with upper and lower limit |
73 |
|
_RL hIceLoc (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
74 |
|
|
75 |
_RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
76 |
C actual snow thickness |
C actual snow thickness |
77 |
_RL hSnwLoc(1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL hSnwLoc(1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
163 |
DO I=1,sNx |
DO I=1,sNx |
164 |
cph need to adjoint-store AREA again before using it in further init. |
cph need to adjoint-store AREA again before using it in further init. |
165 |
cph (all these initialisations involving AREA are nasty "non-linear") |
cph (all these initialisations involving AREA are nasty "non-linear") |
166 |
|
C COMPUTE ACTUAL ICE THICKNESS AND PUT MINIMUM/MAXIMUM |
167 |
|
C ON ICE THICKNESS FOR BUDGET COMPUTATION |
168 |
HICE(I,J) = HEFF(I,J,2,bi,bj)/areaLoc(I,J) |
HICE(I,J) = HEFF(I,J,2,bi,bj)/areaLoc(I,J) |
169 |
|
HICE(I,J) = MAX(HICE(I,J),0.05 _d +00) |
170 |
|
HICE(I,J) = MIN(HICE(I,J),9.0 _d +00) |
171 |
hSnwLoc(I,J) = HSNOW(I,J,bi,bj)/areaLoc(I,J) |
hSnwLoc(I,J) = HSNOW(I,J,bi,bj)/areaLoc(I,J) |
172 |
ENDDO |
ENDDO |
173 |
ENDDO |
ENDDO |