Changeset 7807 for branches/UKMO/r6232_HZG_WAVE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.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/sbcblk_core.F90 (modified) (10 diffs)

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

@@ -17 +17 @@
    !!            3.7  !  2014-06  (L. Brodeau) simplification and optimization of CORE bulk
    !!----------------------------------------------------------------------
-
+   !!            3.6     2016-06  (V.Alari) Added wave dependent drag
+   !coefficient and ocean side stress
    !!----------------------------------------------------------------------
    !!   sbc_blk_core    : bulk formulation as ocean surface boundary condition (forced mode, CORE bulk formulea)
@@ -41 +42 @@
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
    USE prtctl          ! Print control
-   USE sbcwave, ONLY   :  cdn_wave ! wave module
+!   USE sbcwave, ONLY   :  cdn_wave ! wave module (V.Alari 27.06.2016) We get
+!   info from sbc_oce
    USE sbc_ice         ! Surface boundary condition: ice fields
    USE lib_fortran     ! to use key_nosignedzero
@@ -183 +185 @@
          !
          lhftau = ln_taudif                        ! do we use HF tau information?
+
+
+
          jfld = jpfld - COUNT( (/.NOT. lhftau/) )
          !
@@ -192 +197 @@
          END DO
          !                                         ! fill sf with slf_i and control print
+
          CALL fld_fill( sf, slf_i, cn_dir, 'sbc_blk_core', 'flux formulation for ocean surface boundary condition', 'namsbc_core' )
          !
@@ -258 +264 @@
       REAL(wp), DIMENSION(:,:), POINTER ::   zq_zu             ! air spec. hum.  at wind speed height
       !!---------------------------------------------------------------------
+      REAL(wp) ::   ztheta               ! local variable, wind direction,V.Alari 27.06.2016
+
       !
       IF( nn_timing == 1 )  CALL timing_start('blk_oce_core')
@@ -285 +293 @@
       END DO
 #endif
+        ! Start 27.06.2016 (V.Alari)
+           ! Wind vector at T points
+      IF ( ln_cdgw ) THEN
+         ! Use neutral 10-m wind speed from wave model if available
+         DO jj = 2, jpjm1
+            DO ji = fs_2, fs_jpim1   ! vect. opt.
+               ! Local wind direction from non-neutral 10-m wind vector
+               ! (relative to grid and no current)
+               ! since we do not have the directional information from the wave
+               ! model
+               ztheta = ATAN2(sf(jp_wndi)%fnow(ji,jj,1),sf(jp_wndj)%fnow(ji,jj,1))
+               ! Wind vector magnitude from 10-m neutral wind speed from wave
+               ! model
+               zwnd_i(ji,jj) = wspd_wavepar(ji,jj) * SIN(ztheta)
+               zwnd_j(ji,jj) = wspd_wavepar(ji,jj) * COS(ztheta)
+               ! Correct for surface current, 0.0 <= rn_rrelwind <= 1.0
+               zwnd_i(ji,jj) = (zwnd_i(ji,jj) - 0.5*rn_vfac*(pu(ji-1,jj  )+pu(ji,jj)))
+               zwnd_j(ji,jj) = (zwnd_j(ji,jj) - 0.5*rn_vfac*(pv(ji  ,jj-1)+pv(ji,jj)))
+            END DO
+         END DO
+      ELSE ! If no neutral wind from wave model
       DO jj = 2, jpjm1
          DO ji = fs_2, fs_jpim1   ! vect. opt.
@@ -291 +320 @@
          END DO
       END DO
+    ENDIF ! End (V.Alari)
+
       CALL lbc_lnk( zwnd_i(:,:) , 'T', -1. )
       CALL lbc_lnk( zwnd_j(:,:) , 'T', -1. )
@@ -320 +351 @@
     
       ! ... tau module, i and j component
+! Start 27.06.2016 (V.Alari)
+     IF ( ln_cdgw ) THEN
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               !zztmp = rhoa * wspd_wavepar(ji,jj) * Cd_n10(ji,jj)
+               zztmp = rhoa * wndm(ji,jj) * cdn_wave(ji,jj)
+               !taum(ji,jj) = tauoc_wavepar(ji,jj)
+               !taum(ji,jj) = zztmp * wspd_wavepar(ji,jj)
+               taum  (ji,jj) = zztmp * wndm  (ji,jj)
+               zwnd_i(ji,jj) = zztmp * zwnd_i(ji,jj)
+               zwnd_j(ji,jj) = zztmp * zwnd_j(ji,jj)
+            ENDDO
+         ENDDO
+      ELSE
       DO jj = 1, jpj
          DO ji = 1, jpi
@@ -328 +373 @@
          END DO
       END DO
-
+     ENDIF ! End, (V.Alari)
       ! ... add the HF tau contribution to the wind stress module?
       IF( lhftau ) THEN
@@ -349 +394 @@
       CALL lbc_lnk( vtau(:,:), 'V', -1. )
 
-    
+        ! START: 27.06.2015 (V. Alari). Added tauoc
+
+         IF (ln_tauoc) THEN
+         utau(:,:) = utau(:,:)*tauoc_wavepar(:,:)
+         vtau(:,:) = vtau(:,:)*tauoc_wavepar(:,:)
+         taum(:,:) = taum(:,:)*tauoc_wavepar(:,:)
+         ENDIF
+
+        ! END  (V.Alari)
       !  Turbulent fluxes over ocean
       ! -----------------------------