- Timestamp:
- 2020-07-09T12:39:43+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icedyn_rhg_evp.F90
r13064 r13279 160 160 REAL(wp), DIMENSION(jpi,jpj) :: zmsk01x, zmsk01y ! dummy arrays 161 161 REAL(wp), DIMENSION(jpi,jpj) :: zmsk00x, zmsk00y ! mask for ice presence 162 REAL(wp), DIMENSION(jpi,jpj) :: zfmask , zwf! mask at F points for the ice162 REAL(wp), DIMENSION(jpi,jpj) :: zfmask ! mask at F points for the ice 163 163 164 164 REAL(wp), PARAMETER :: zepsi = 1.0e-20_wp ! tolerance parameter … … 202 202 203 203 ! Lateral boundary conditions on velocity (modify zfmask) 204 zwf(:,:) = zfmask(:,:)205 204 DO jj = 2, jpjm1 206 205 DO ji = fs_2, fs_jpim1 ! vector opt. 207 206 IF( zfmask(ji,jj) == 0._wp ) THEN 208 zfmask(ji,jj) = rn_ishlat * MIN( 1._wp , MAX( zwf(ji+1,jj), zwf(ji,jj+1), zwf(ji-1,jj), zwf(ji,jj-1) ) ) 207 zfmask(ji,jj) = rn_ishlat * MIN( 1._wp , MAX( umask(ji,jj,1), umask(ji,jj+1,1), & 208 & vmask(ji,jj,1), vmask(ji+1,jj,1) ) ) 209 209 ENDIF 210 210 END DO … … 212 212 DO jj = 2, jpjm1 213 213 IF( zfmask(1,jj) == 0._wp ) THEN 214 zfmask(1 ,jj) = rn_ishlat * MIN( 1._wp , MAX( zwf(2,jj), zwf(1,jj+1), zwf(1,jj-1) ) )214 zfmask(1 ,jj) = rn_ishlat * MIN( 1._wp , MAX( vmask(2,jj,1), umask(1,jj+1,1), umask(1,jj,1) ) ) 215 215 ENDIF 216 216 IF( zfmask(jpi,jj) == 0._wp ) THEN 217 zfmask(jpi,jj) = rn_ishlat * MIN( 1._wp , MAX( zwf(jpi,jj+1), zwf(jpim1,jj), zwf(jpi,jj-1) ) )218 217 zfmask(jpi,jj) = rn_ishlat * MIN( 1._wp , MAX( umask(jpi,jj+1,1), vmask(jpim1,jj,1), umask(jpi,jj-1,1) ) ) 218 ENDIF 219 219 END DO 220 220 DO ji = 2, jpim1 221 221 IF( zfmask(ji,1) == 0._wp ) THEN 222 zfmask(ji, 1 ) = rn_ishlat * MIN( 1._wp , MAX( zwf(ji+1,1), zwf(ji,2), zwf(ji-1,1) ) )222 zfmask(ji, 1 ) = rn_ishlat * MIN( 1._wp , MAX( vmask(ji+1,1,1), umask(ji,2,1), vmask(ji,1,1) ) ) 223 223 ENDIF 224 224 IF( zfmask(ji,jpj) == 0._wp ) THEN 225 zfmask(ji,jpj) = rn_ishlat * MIN( 1._wp , MAX( zwf(ji+1,jpj), zwf(ji-1,jpj), zwf(ji,jpjm1) ) )225 zfmask(ji,jpj) = rn_ishlat * MIN( 1._wp , MAX( vmask(ji+1,jpj,1), vmask(ji-1,jpj,1), umask(ji,jpjm1,1) ) ) 226 226 ENDIF 227 227 END DO
Note: See TracChangeset
for help on using the changeset viewer.