- Timestamp:
- 2017-03-20T17:21:42+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90
r6963 r7814 259 259 closing_net(ji,jj) = 0._wp 260 260 opning (ji,jj) = 0._wp 261 ato_i (ji,jj) = MAX( 0._wp, 1._wp - SUM( a_i(ji,jj,:) ) ) 261 262 ELSE 262 263 iterate_ridging = 1 … … 844 845 END DO 845 846 846 strength(:,:) = rn_pe_rdg * Cp * strength(:,:) / aksum(:,:) 847 strength(:,:) = rn_pe_rdg * Cp * strength(:,:) / aksum(:,:) * tmask(:,:,1) 847 848 ! where Cp = (g/2)*(rhow-rhoi)*(rhoi/rhow) and rn_pe_rdg accounts for frictional dissipation 848 849 ksmooth = 1 … … 853 854 ELSE ! kstrngth ne 1: Hibler (1979) form 854 855 ! 855 strength(:,:) = rn_pstar * vt_i(:,:) * EXP( - rn_crhg * ( 1._wp - at_i(:,:) ) ) 856 strength(:,:) = rn_pstar * vt_i(:,:) * EXP( - rn_crhg * ( 1._wp - at_i(:,:) ) ) * tmask(:,:,1) 856 857 ! 857 858 ksmooth = 1
Note: See TracChangeset
for help on using the changeset viewer.