155 |
ENDDO |
ENDDO |
156 |
ENDIF |
ENDIF |
157 |
|
|
158 |
|
#ifdef SHORTWAVE_HEATING |
159 |
|
C Penetrating SW radiation |
160 |
|
DO j=jMin,jMax |
161 |
|
DO i=iMin,iMax |
162 |
|
gT(i,j,k,bi,bj) = gT(i,j,k,bi,bj) |
163 |
|
& - maskC(i,j)*Qsw(i,j,bi,bj)*swfracb(k) |
164 |
|
& *recip_Cp*recip_rhoNil*recip_dRf(k) |
165 |
|
ENDDO |
166 |
|
ENDDO |
167 |
|
#endif |
168 |
|
|
169 |
RETURN |
RETURN |
170 |
END |
END |
171 |
CStartOfInterface |
CStartOfInterface |
214 |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |
215 |
& +maskC(i,j)*( |
& +maskC(i,j)*( |
216 |
& -lambdaSaltClimRelax*(salt(i,j,kLev,bi,bj)-SSS(i,j,bi,bj)) |
& -lambdaSaltClimRelax*(salt(i,j,kLev,bi,bj)-SSS(i,j,bi,bj)) |
217 |
#ifndef USE_NATURAL_BCS |
#ifdef USE_NATURAL_BCS |
218 |
|
& +EmPmR(i,j,bi,bj)*recip_dRf(1)*salt(i,j,kLev,bi,bj) |
219 |
|
#else |
220 |
& +EmPmR(i,j,bi,bj)*recip_dRf(1)*35. |
& +EmPmR(i,j,bi,bj)*recip_dRf(1)*35. |
221 |
#endif |
#endif |
222 |
& ) |
& ) |