Changeset 1025 for trunk/NEMO/OPA_SRC/SBC/sbcflx.F90
- Timestamp:
- 2008-05-30T13:26:09+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbcflx.F90
r1000 r1025 99 99 INTEGER :: ierror ! return error code 100 100 REAL(wp) :: zfact ! temporary scalar 101 REAL(wp) :: ztx, zty, ztau 101 REAL(wp) :: ztx, zty, ztau, zcoef 102 102 !! 103 103 CHARACTER(len=100) :: cn_dir ! Root directory for location of flx files … … 183 183 END DO 184 184 185 ! Estimation of wind speed as a function of wind stress 185 ! Estimation of wind speed as a function of wind stress ( |tau|=rhoa*Cd*|U|^2 ) 186 zcoef = 0.5 / ( rhoa * cdrag ) 186 187 !CDIR NOVERRCHK 187 DO jj = 1, jpj188 DO jj = 2, jpjm1 188 189 !CDIR NOVERRCHK 189 DO ji = 1, jpi190 ztx = utau(ji,jj) * umask(ji,jj,1)191 zty = vtau(ji,jj) * vmask(ji,jj,1)190 DO ji = fs_2, fs_jpim1 ! vect. opt. 191 ztx = utau(ji-1,jj ) + utau(ji,jj) 192 zty = vtau(ji ,jj-1) + vtau(ji,jj) 192 193 ztau = SQRT( ztx * ztx + zty * zty ) 193 wndm(ji,jj) = SQRT ( ztau / ( rhoa * cdrag ) ) 194 ENDDO 195 ENDDO 194 wndm(ji,jj) = SQRT ( ztau * zcoef ) * tmask(ji,jj,1) 195 END DO 196 END DO 197 CALL lbc_lnk( wndm(:,:) , 'T', 1. ) 196 198 197 199 ! control print (if less than 100 time-step asked)
Note: See TracChangeset
for help on using the changeset viewer.