Changeset 15605 for NEMO/branches/CNRS/dev_r14723_tides_under_isf/src/OCE/ISF/isfcavgam.F90

Timestamp:

2021-12-16T14:58:14+01:00 (2 years ago)

Author:

khutchinson

Message:

changes required to add background tidal velocities in isf cavities (not sette tested yet)

File:

• NEMO/branches/CNRS/dev_r14723_tides_under_isf/src/OCE/ISF/isfcavgam.F90 (modified) (6 diffs)

NEMO/branches/CNRS/dev_r14723_tides_under_isf/src/OCE/ISF/isfcavgam.F90

@@ -17 +17 @@
    USE phycst  , ONLY: grav, vkarmn        ! physical constant
    USE eosbn2  , ONLY: eos_rab             ! equation of state
-   USE zdfdrg  , ONLY: rCd0_top, r_ke0_top ! vertical physics: top/bottom drag coef.
+! tipaccs 2d top tidal velocity
+   USE zdfdrg  , ONLY: rCd0_top, rke0_top  ! vertical physics: top/bottom drag coef.
+! end tipaccs 2d top tidal velocity
    USE iom     , ONLY: iom_put             !
    USE lib_mpp , ONLY: ctl_stop
@@ -92 +94 @@
          pgs(:,:) = rn_gammas0
       CASE ( 'vel' ) ! gamma is proportional to u*
-         CALL gammats_vel      (                   zutbl, zvtbl, rCd0_top, r_ke0_top,               pgt, pgs )
+! tipaccs 2d top tidal velocity
+         CALL gammats_vel      (                   zutbl, zvtbl, rCd0_top, rke0_top,               pgt, pgs )
       CASE ( 'vel_stab' ) ! gamma depends of stability of boundary layer and u*
-         CALL gammats_vel_stab (Kmm, pttbl, pstbl, zutbl, zvtbl, rCd0_top, r_ke0_top, pqoce, pqfwf, pgt, pgs )
+         CALL gammats_vel_stab (Kmm, pttbl, pstbl, zutbl, zvtbl, rCd0_top, rke0_top, pqoce, pqfwf, pgt, pgs )
+! end tipaccs 2d top tidal velocity
       CASE DEFAULT
          CALL ctl_stop('STOP','method to compute gamma (cn_gammablk) is unknown (should not see this)')
@@ -131 +135 @@
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: putbl, pvtbl ! velocity in the losch top boundary layer
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pCd          ! drag coefficient
-      REAL(wp),                     INTENT(in   ) :: pke2         ! background velocity
+! tipaccs (2d tidal velocity)
+      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pke2         ! background velocity
+! end tipaccs (2d tidal velocity)
       !!---------------------------------------------------------------------
       REAL(wp), DIMENSION(jpi,jpj) :: zustar
@@ -137 +143 @@
       !
       ! compute ustar (AD15 eq. 27)
-      zustar(:,:) = SQRT( pCd(:,:) * ( putbl(:,:) * putbl(:,:) + pvtbl(:,:) * pvtbl(:,:) + pke2 ) ) * mskisf_cav(:,:)
+! tipaccs (2d tidal velocity)
+      zustar(:,:) = SQRT( pCd(:,:) * ( putbl(:,:) * putbl(:,:) + pvtbl(:,:) * pvtbl(:,:) + pke2(:,:) ) ) * mskisf_cav(:,:)
+! end tipaccs (2d tidal velocity)
       !
       ! Compute gammats
@@ -161 +169 @@
       !!-------------------------- IN  -------------------------------------
       INTEGER                                     :: Kmm            ! ocean time level index
-      REAL(wp),                     INTENT(in   ) :: pke2           ! background velocity squared
+! tipaccs (2d tidal velocity)
+      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pke2           ! background velocity squared
+! end tipaccs (2d tidal velocity)
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pqoce, pqfwf   ! surface heat flux and fwf flux
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pCd            ! drag coeficient
@@ -186 +196 @@
       !
       ! compute ustar
-      zustar(:,:) = SQRT( pCd * ( putbl(:,:) * putbl(:,:) + pvtbl(:,:) * pvtbl(:,:) + pke2 ) )
+! tipaccs (2d tidal velocity)
+      zustar(:,:) = SQRT( pCd * ( putbl(:,:) * putbl(:,:) + pvtbl(:,:) * pvtbl(:,:) + pke2(:,:) ) )
+! end tipaccs (2d tidal velocity)
       !
       ! output ustar