- Timestamp:
- 2017-09-26T15:24:17+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_da.F90
r8534 r8563 135 135 136 136 ! Contribution to salt flux 137 sfx_lam_1d(ji) = sfx_lam_1d(ji) + rhoic * h t_i_1d(ji) * zda * sm_i_1d(ji) * r1_rdtice137 sfx_lam_1d(ji) = sfx_lam_1d(ji) + rhoic * h_i_1d(ji) * zda * sm_i_1d(ji) * r1_rdtice 138 138 139 139 ! Contribution to heat flux into the ocean [W.m-2], (<0) 140 hfx_thd_1d(ji) = hfx_thd_1d(ji) - zda * r1_rdtice * ( h t_i_1d(ji) * r1_nlay_i * SUM( e_i_1d(ji,1:nlay_i) ) &141 + h t_s_1d(ji) * r1_nlay_s * SUM( e_s_1d(ji,1:nlay_s) ) )140 hfx_thd_1d(ji) = hfx_thd_1d(ji) - zda * r1_rdtice * ( h_i_1d(ji) * r1_nlay_i * SUM( e_i_1d(ji,1:nlay_i) ) & 141 + h_s_1d(ji) * r1_nlay_s * SUM( e_s_1d(ji,1:nlay_s) ) ) 142 142 143 143 ! Contribution to mass flux 144 wfx_lam_1d(ji) = wfx_lam_1d(ji) + zda * r1_rdtice * ( rhoic * h t_i_1d(ji) + rhosn * ht_s_1d(ji) )144 wfx_lam_1d(ji) = wfx_lam_1d(ji) + zda * r1_rdtice * ( rhoic * h_i_1d(ji) + rhosn * h_s_1d(ji) ) 145 145 146 146 ! new concentration 147 147 a_i_1d(ji) = a_i_1d(ji) - zda 148 148 149 ! ensure that h t_i = 0 where a_i = 0149 ! ensure that h_i = 0 where a_i = 0 150 150 IF( a_i_1d(ji) == 0._wp ) THEN 151 h t_i_1d(ji) = 0._wp152 h t_s_1d(ji) = 0._wp151 h_i_1d(ji) = 0._wp 152 h_s_1d(ji) = 0._wp 153 153 ENDIF 154 154 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.