Changeset 1000 for trunk/NEMO/OPA_SRC/SBC/sbcflx.F90
- Timestamp:
- 2008-05-30T11:22:27+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbcflx.F90
r888 r1000 50 50 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf ! structure of input fields (file informations, fields read) 51 51 52 REAL(wp) :: rhoa = 1.22 ! Air density kg/m3 53 REAL(wp) :: cdrag = 1.5e-3 ! drag coefficient 54 52 55 !! * Substitutions 53 56 # include "domzgr_substitute.h90" … … 96 99 INTEGER :: ierror ! return error code 97 100 REAL(wp) :: zfact ! temporary scalar 101 REAL(wp) :: ztx, zty, ztau 98 102 !! 99 103 CHARACTER(len=100) :: cn_dir ! Root directory for location of flx files … … 179 183 END DO 180 184 185 ! Estimation of wind speed as a function of wind stress 186 !CDIR NOVERRCHK 187 DO jj = 1, jpj 188 !CDIR NOVERRCHK 189 DO ji = 1, jpi 190 ztx = utau(ji,jj) * umask(ji,jj,1) 191 zty = vtau(ji,jj) * vmask(ji,jj,1) 192 ztau = SQRT( ztx * ztx + zty * zty ) 193 wndm(ji,jj) = SQRT ( ztau / ( rhoa * cdrag ) ) 194 ENDDO 195 ENDDO 196 181 197 ! control print (if less than 100 time-step asked) 182 198 IF( nitend-nit000 <= 100 .AND. lwp ) THEN
Note: See TracChangeset
for help on using the changeset viewer.