Changeset 4990 for trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90
- Timestamp:
- 2014-12-15T17:42:49+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90
r4624 r4990 82 82 !! - utau, vtau i- and j-component of the wind stress 83 83 !! - taum wind stress module at T-point 84 !! - wndm 10m wind module at T-point 84 !! - wndm 10m wind module at T-point over free ocean or leads in presence of sea-ice 85 85 !! - qns, qsr non-slor and solar heat flux 86 86 !! - emp evaporation minus precipitation … … 233 233 ! Interpolate utau, vtau into the grid_V and grid_V 234 234 !------------------------------------------------- 235 235 ! Note the use of 0.5*(2-umask) in order to unmask the stress along coastlines 236 ! Note the use of MAX(tmask(i,j),tmask(i+1,j) is to mask tau over ice shelves 236 237 DO jj = 1, jpjm1 237 238 DO ji = 1, fs_jpim1 238 239 utau(ji,jj) = 0.5 * ( 2. - umask(ji,jj,1) ) * ( utau(ji,jj) * tmask(ji,jj,1) & 239 & + utau(ji+1,jj) * tmask(ji+1,jj,1) ) 240 & + utau(ji+1,jj) * tmask(ji+1,jj,1) ) & 241 & * MAX(tmask(ji,jj,1),tmask(ji+1,jj ,1)) 240 242 vtau(ji,jj) = 0.5 * ( 2. - vmask(ji,jj,1) ) * ( vtau(ji,jj) * tmask(ji,jj,1) & 241 & + vtau(ji,jj+1) * tmask(ji,jj+1,1) ) 243 & + vtau(ji,jj+1) * tmask(ji,jj+1,1) ) & 244 & * MAX(tmask(ji,jj,1),tmask(ji ,jj+1,1)) 242 245 END DO 243 246 END DO
Note: See TracChangeset
for help on using the changeset viewer.