Changeset 7807 for branches/UKMO/r6232_HZG_WAVE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90

Timestamp:

2017-03-17T10:44:05+01:00 (7 years ago)

Author:

jcastill

Message:

Changes as in HZG wave forcing branch, but adapted to r6232

File:

• branches/UKMO/r6232_HZG_WAVE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90 (modified) (3 diffs)

branches/UKMO/r6232_HZG_WAVE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90

@@ -22 +22 @@
    USE lib_mpp         ! distribued memory computing library
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
+   USE sbcwave         ! Wave parameters
 
    IMPLICIT NONE
@@ -152 +153 @@
                zty = vtau(ji  ,jj-1) + vtau(ji,jj) 
                zmod = 0.5 * SQRT( ztx * ztx + zty * zty )
+               IF (ln_cdgw) THEN
+                  zcdrag = 1.5e-3
+                  IF (cdn_wave(ji,jj)>=0.) zcdrag = cdn_wave(ji,jj)
+                  zcoef = 1./(zrhoa*zcdrag)
+               ENDIF
                taum(ji,jj) = zmod
                wndm(ji,jj) = SQRT( zmod * zcoef )
@@ -158 +164 @@
          taum(:,:) = taum(:,:) * tmask(:,:,1) ; wndm(:,:) = wndm(:,:) * tmask(:,:,1)
          CALL lbc_lnk( taum(:,:), 'T', 1. )   ;   CALL lbc_lnk( wndm(:,:), 'T', 1. )
+
+         IF (ln_tauoc) THEN
+         utau(:,:) = utau(:,:)*tauoc_wavepar(:,:)
+         vtau(:,:) = vtau(:,:)*tauoc_wavepar(:,:)
+         taum(:,:) = taum(:,:)*tauoc_wavepar(:,:)
+         ENDIF
 
          IF( nitend-nit000 <= 100 .AND. lwp ) THEN                ! control print (if less than 100 time-step asked)