[MITgcm-support] calculation of the wind stress in pkg/exf

[Wu Yang]

hi ,all
 
Does the wind stress calculate from pkg/exf:exf_bulkformulae.F ?
the Zonal/Meridional  surface wind stress is calculated through the following codes  in exf_bulkformulae.F
 
#ifdef ALLOW_ATM_WIND
           ustress(i,j,bi,bj) = 0. _d 0
           vstress(i,j,bi,bj) = 0. _d 0
#endif /* ALLOW_ATM_WIND */
           hflux  (i,j,bi,bj) = 0. _d 0
           hs(i,j,bi,bj)      = 0. _d 0
           hl(i,j,bi,bj)      = 0. _d 0
          endif
 
#else  /* ifndef ALLOW_ATM_TEMP */
#ifdef ALLOW_ATM_WIND
          umps               = sh(i,j,bi,bj)
          tmpbulk            = exf_scal_BulkCdn *
     &         ( cdrag_1/umps + cdrag_2 + cdrag_3*umps )
          ustress(i,j,bi,bj) = atmrho*tmpbulk*wspeed(i,j,bi,bj)*
     &         uwind(i,j,bi,bj)
          vstress(i,j,bi,bj) = atmrho*tmpbulk*wspeed(i,j,bi,bj)*
     &         vwind(i,j,bi,bj)
 
 
#endif
#endif /* ifndef ALLOW_ATM_TEMP */
         enddo
        enddo
       enddo
      enddo
 
 
thanks
 
Yang
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20150330/7a360fa1/attachment.htm>