- Timestamp:
- 2015-03-03T11:15:02+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r5079 r5116 184 184 ENDIF 185 185 186 ! -- Energy needed to bring ocean surface layer until its freezing (<0, J.m-2) --- !186 ! --- Energy needed to bring ocean surface layer until its freezing (<0, J.m-2) --- ! 187 187 zqfr = tmask(ji,jj,1) * rau0 * rcp * fse3t_m(ji,jj) * ( t_bo(ji,jj) - ( sst_m(ji,jj) + rt0 ) ) 188 188 189 ! --- Energy from the turbulent oceanic heat flux (W/m2) --- ! 190 zfric_u = MAX( SQRT( ust2s(ji,jj) ), zfric_umin ) 191 fhtur(ji,jj) = MAX( 0._wp, rswitch * rau0 * rcp * zch * zfric_u * ( ( sst_m(ji,jj) + rt0 ) - t_bo(ji,jj) ) ) ! W.m-2 192 fhtur(ji,jj) = rswitch * MIN( fhtur(ji,jj), - zqfr * r1_rdtice / MAX( at_i(ji,jj), epsi10 ) ) 193 ! upper bound for fhtur: the heat retrieved from the ocean must be smaller than the heat necessary to reach 194 ! the freezing point, so that we do not have SST < T_freeze 195 ! This implies: - ( fhtur(ji,jj) * at_i(ji,jj) * rtdice ) - zqfr >= 0 196 189 197 !-- Energy Budget of the leads (J.m-2). Must be < 0 to form ice 190 qlead(ji,jj) = MIN( 0._wp , zqld - zqfr )198 qlead(ji,jj) = MIN( 0._wp , zqld - ( fhtur(ji,jj) * at_i(ji,jj) * rdt_ice ) - zqfr ) 191 199 192 200 ! If there is ice and leads are warming, then transfer energy from the lead budget and use it for bottom melting … … 198 206 ENDIF 199 207 ! 200 !-- Energy from the turbulent oceanic heat flux --- !201 !clem zfric_u = MAX ( MIN( SQRT( ust2s(ji,jj) ) , zfric_umax ) , zfric_umin )202 zfric_u = MAX( SQRT( ust2s(ji,jj) ), zfric_umin )203 fhtur(ji,jj) = MAX( 0._wp, rswitch * rau0 * rcp * zch * zfric_u * ( ( sst_m(ji,jj) + rt0 ) - t_bo(ji,jj) ) ) ! W.m-2204 ! upper bound for fhtur: we do not want SST to drop below Tfreeze.205 ! So we say that the heat retrieved from the ocean (fhtur+fhld) must be < to the heat necessary to reach Tfreeze (zqfr)206 ! This is not a clean budget, so that should be corrected at some point207 fhtur(ji,jj) = rswitch * MIN( fhtur(ji,jj), - fhld(ji,jj) - zqfr * r1_rdtice / MAX( at_i(ji,jj), epsi10 ) )208 209 208 ! ----------------------------------------- 210 209 ! Net heat flux on top of ice-ocean [W.m-2]
Note: See TracChangeset
for help on using the changeset viewer.