Changeset 3085 for branches/2011/dev_MERCATOR_INGV_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

Timestamp:

2011-11-14T14:13:32+01:00 (12 years ago)

Author:

cbricaud

Message:

commit changes from dev_INGV_2011

File:

• branches/2011/dev_MERCATOR_INGV_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (6 diffs)

branches/2011/dev_MERCATOR_INGV_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -34 +34 @@
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
    USE prtctl          ! Print control
+   USE sbcwave,ONLY :  cdn_wave !wave module 
 #if defined key_lim3
    USE sbc_ice         ! Surface boundary condition: ice fields
@@ -43 +44 @@
    PUBLIC   sbc_blk_core         ! routine called in sbcmod module
    PUBLIC   blk_ice_core         ! routine called in sbc_ice_lim module
+   PUBLIC   turb_core_2z         ! routine calles in sbcblk_ecmwf module
 
    INTEGER , PARAMETER ::   jpfld   = 9           ! maximum number of files to read 
@@ -682 +684 @@
       !! Neutral Drag Coefficient
       stab    = 0.5 + sign(0.5,dT)    ! stable : stab = 1 ; unstable : stab = 0 
-      Cd_n10  = 1E-3 * ( 2.7/dU10 + 0.142 + dU10/13.09 )    !   L & Y eq. (6a)
+      IF  ( ln_cdgw ) THEN
+        cdn_wave = cdn_wave - rsmall*(tmask(:,:,1)-1)
+        Cd_n10(:,:) =   cdn_wave
+      ELSE
+        Cd_n10  = 1E-3 * ( 2.7/dU10 + 0.142 + dU10/13.09 )    !   L & Y eq. (6a)
+      ENDIF
       sqrt_Cd_n10 = sqrt(Cd_n10)
       Ce_n10  = 1E-3 * ( 34.6 * sqrt_Cd_n10 )               !   L & Y eq. (6b)
@@ -705 +712 @@
          zpsi_m  = psi_m(zeta)
 
-         !! Shifting the wind speed to 10m and neutral stability :
-         U_n10 = dU10*1./(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) !  L & Y eq. (9a)
-
-         !! Updating the neutral 10m transfer coefficients :
-         Cd_n10  = 1E-3 * (2.7/U_n10 + 0.142 + U_n10/13.09)              !  L & Y eq. (6a)
-         sqrt_Cd_n10 = sqrt(Cd_n10)
-         Ce_n10  = 1E-3 * (34.6 * sqrt_Cd_n10)                           !  L & Y eq. (6b)
-         stab    = 0.5 + sign(0.5,zeta)
-         Ch_n10  = 1E-3*sqrt_Cd_n10*(18.*stab + 32.7*(1-stab))           !  L & Y eq. (6c), (6d)
-
-         !! Shifting the neutral  10m transfer coefficients to ( zu , zeta ) :
-         !!
-         xct = 1. + sqrt_Cd_n10/kappa*(log(zu/10) - zpsi_m)
-         Cd  = Cd_n10/(xct*xct) ;  sqrt_Cd = sqrt(Cd)
+         IF  ( ln_cdgw ) THEN
+           sqrt_Cd=kappa/((kappa/sqrt_Cd_n10) - zpsi_m) ; Cd=sqrt_Cd*sqrt_Cd;
+         ELSE
+           !! Shifting the wind speed to 10m and neutral stability :
+           U_n10 = dU10*1./(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) !  L & Y eq. (9a)
+
+           !! Updating the neutral 10m transfer coefficients :
+           Cd_n10  = 1E-3 * (2.7/U_n10 + 0.142 + U_n10/13.09)              !  L & Y eq. (6a)
+           sqrt_Cd_n10 = sqrt(Cd_n10)
+           Ce_n10  = 1E-3 * (34.6 * sqrt_Cd_n10)                           !  L & Y eq. (6b)
+           stab    = 0.5 + sign(0.5,zeta)
+           Ch_n10  = 1E-3*sqrt_Cd_n10*(18.*stab + 32.7*(1-stab))           !  L & Y eq. (6c), (6d)
+
+           !! Shifting the neutral  10m transfer coefficients to ( zu , zeta ) :
+           !!
+           xct = 1. + sqrt_Cd_n10/kappa*(log(zu/10) - zpsi_m)
+           Cd  = Cd_n10/(xct*xct) ;  sqrt_Cd = sqrt(Cd)
+         ENDIF
          !!
          xlogt = log(zu/10.) - zpsi_h
@@ -812 +823 @@
       !! Neutral Drag Coefficient :
       stab = 0.5 + sign(0.5,dT)                 ! stab = 1  if dT > 0  -> STABLE
-      Cd_n10  = 1E-3*( 2.7/dU10 + 0.142 + dU10/13.09 ) 
+      IF  ( ln_cdgw ) THEN
+        cdn_wave = cdn_wave - rsmall*(tmask(:,:,1)-1)
+        Cd_n10(:,:) =   cdn_wave
+      ELSE
+        Cd_n10  = 1E-3*( 2.7/dU10 + 0.142 + dU10/13.09 ) 
+      ENDIF
       sqrt_Cd_n10 = sqrt(Cd_n10)
       Ce_n10  = 1E-3*( 34.6 * sqrt_Cd_n10 )
@@ -853 +869 @@
          stab = 0.5 + sign(0.5,q_zu) ;  q_zu = stab*q_zu
          !!
-         !! Updating the neutral 10m transfer coefficients :
-         Cd_n10  = 1E-3 * (2.7/U_n10 + 0.142 + U_n10/13.09)    ! L & Y eq. (6a)
-         sqrt_Cd_n10 = sqrt(Cd_n10)
-         Ce_n10  = 1E-3 * (34.6 * sqrt_Cd_n10)                 ! L & Y eq. (6b)
-         stab    = 0.5 + sign(0.5,zeta_u)
-         Ch_n10  = 1E-3*sqrt_Cd_n10*(18.*stab + 32.7*(1-stab)) ! L & Y eq. (6c-6d)
-         !!
-         !!
-         !! Shifting the neutral 10m transfer coefficients to (zu,zeta_u) :
-!        xct = 1. + sqrt_Cd_n10/kappa*(log(zu/10.) - psi_m(zeta_u))
-         xct = 1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)
-         Cd = Cd_n10/(xct*xct) ; sqrt_Cd = sqrt(Cd)
-         !!
-!        xlogt = log(zu/10.) - psi_h(zeta_u)
+         IF  ( ln_cdgw ) THEN
+            sqrt_Cd=kappa/((kappa/sqrt_Cd_n10) - zpsi_m) ; Cd=sqrt_Cd*sqrt_Cd;
+         ELSE
+           !! Updating the neutral 10m transfer coefficients :
+           Cd_n10  = 1E-3 * (2.7/U_n10 + 0.142 + U_n10/13.09)    ! L & Y eq. (6a)
+           sqrt_Cd_n10 = sqrt(Cd_n10)
+           Ce_n10  = 1E-3 * (34.6 * sqrt_Cd_n10)                 ! L & Y eq. (6b)
+           stab    = 0.5 + sign(0.5,zeta_u)
+           Ch_n10  = 1E-3*sqrt_Cd_n10*(18.*stab + 32.7*(1-stab)) ! L & Y eq. (6c-6d)
+           !!
+           !!
+           !! Shifting the neutral 10m transfer coefficients to (zu,zeta_u) :
+           xct = 1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)
+           Cd = Cd_n10/(xct*xct) ; sqrt_Cd = sqrt(Cd)
+         ENDIF
+         !!
          xlogt = log(zu/10.) - zpsi_hu
          !!