Changeset 9121
- Timestamp:
- 2017-12-18T21:28:47+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90
r9119 r9121 236 236 DO jk = 1, nlay_s 237 237 WHERE( v_s(:,:,:) > epsi20 ) !--- icy area 238 t_s(:,:,jk,:) = MAX( -100._wp , & 239 & MIN( r1_cpic * ( -r1_rhosn * ( e_s(:,:,jk,:) / v_s(:,:,:) * zlay_s ) + lfus ) , & 240 & 0._wp ) ) & 241 & + rt0 238 t_s(:,:,jk,:) = rt0 + MAX( -100._wp , & 239 & MIN( r1_cpic * ( -r1_rhosn * ( e_s(:,:,jk,:) / v_s(:,:,:) * zlay_s ) + lfus ) , 0._wp ) ) 242 240 ELSEWHERE !--- no ice 243 241 t_s(:,:,jk,:) = rt0
Note: See TracChangeset
for help on using the changeset viewer.