New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 1025 for trunk/NEMO/OPA_SRC/SBC/sbcflx.F90 – NEMO

Ignore:
Timestamp:
2008-05-30T13:26:09+02:00 (16 years ago)
Author:
cetlod
Message:

adding wind speed module variable, see ticket 172

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/SBC/sbcflx.F90

    r1000 r1025  
    9999      INTEGER  ::   ierror       ! return error code 
    100100      REAL(wp) ::   zfact        ! temporary scalar 
    101       REAL(wp) ::   ztx, zty, ztau 
     101      REAL(wp) ::   ztx, zty, ztau, zcoef 
    102102      !! 
    103103      CHARACTER(len=100) ::  cn_dir                               ! Root directory for location of flx files 
     
    183183      END DO 
    184184 
    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 )  
    186187!CDIR NOVERRCHK 
    187       DO jj = 1, jpj 
     188      DO jj = 2, jpjm1 
    188189!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) 
     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)  
    192193            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. ) 
    196198 
    197199      ! control print (if less than 100 time-step asked) 
Note: See TracChangeset for help on using the changeset viewer.