Changeset 11266 for NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare.F90

Timestamp:

2019-07-15T12:09:30+02:00 (5 years ago)

Author:

laurent

Message:

LB: CSWL param now uses NEMO time step (rdt) and previous-t value of t_skin as first guess.

File:

• NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare.F90 (modified) (7 diffs)

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare.F90

@@ -57 +57 @@
       &                   Cd, Ch, Ce, t_zu, q_zu, U_blk,      &
       &                   Cdn, Chn, Cen,                      &
-      &                   Qsw, rad_lw, slp                   )
+      &                   Qsw, rad_lw, slp, Tsk_b             )
       !!----------------------------------------------------------------------
       !!                      ***  ROUTINE  turb_coare  ***
@@ -81 +81 @@
       !! INPUT/OUTPUT:
       !! -------------
-      !!    *  T_s  : SST or skin temperature                                 [K]
+      !!    *  T_s  : always "bulk SST" as input                              [K]
+      !!              -> unchanged "bulk SST" as output if CSWL not used      [K]
+      !!              -> skin temperature as output if CSWL used              [K]
+      !!
       !!    *  q_s  : SSQ aka saturation specific humidity at temp. T_s       [kg/kg]
       !!              -> doesn't need to be given a value if skin temp computed (in case l_use_skin=True)
@@ -91 +94 @@
       !!    *  rad_lw : downwelling longwave radiation at the surface  (>0)   [W/m^2]
       !!    *  slp    : sea-level pressure                                    [Pa]
+      !!    *  Tsk_b  : estimate of skin temperature at previous time-step    [K]
       !!
       !! OUTPUT :
@@ -122 +126 @@
       REAL(wp), INTENT(in   ), OPTIONAL, DIMENSION(jpi,jpj) ::   rad_lw   !             [W/m^2]
       REAL(wp), INTENT(in   ), OPTIONAL, DIMENSION(jpi,jpj) ::   slp      !             [Pa]
+      REAL(wp), INTENT(in   ), OPTIONAL, DIMENSION(jpi,jpj) ::   Tsk_b    !             [Pa]
       !
       INTEGER :: j_itt
@@ -152 +157 @@
 
       !! Initialization for cool skin:
+      zsst   = T_s    ! save the bulk SST
       IF( l_use_skin ) THEN
-         zsst   = T_s    ! save the bulk SST
-         T_s    = T_s - 0.25                      ! First guess of correction
+         ! First guess for skin temperature:
+         IF( PRESENT(Tsk_b) ) THEN
+            T_s = Tsk_b
+         ELSE
+            T_s = T_s - 0.25     ! sst - 0.25
+         END IF
          q_s    = rdct_qsat_salt*q_sat(MAX(T_s, 200._wp), slp) ! First guess of q_s
       END IF
@@ -202 +212 @@
       q_star = dq_zu*ztmp0
 
-      ! What's need to be done if zt /= zu:
+      ! What needs to be done if zt /= zu:
       IF( .NOT. l_zt_equal_zu ) THEN
          !! First update of values at zu (or zt for wind)
@@ -275 +285 @@
             ztmp1 = U_blk*MAX(rho_air(t_zu, q_zu, slp), 1._wp)     ! rho*U10
             ztmp2 = T_s*T_s
-            ztmp1 = ztmp1 * ( Ce*rLevap*(q_zu - q_s) + Ch*rCp_dry*(t_zu - T_s) ) & ! Total turb. heat flux
-               &       +    (rad_lw - emiss_w*stefan*ztmp2*ztmp2)                  ! Net longwave flux
+            ztmp1 = ztmp1 * ( Ce*L_vap(T_s)*(q_zu - q_s) + Ch*cp_air(q_zu)*(t_zu - T_s) ) & ! Total turb. heat flux
+               &       +      rad_lw - emiss_w*stefan*ztmp2*ztmp2                           ! Net longwave flux
+            !!         => "ztmp1" is the net non-solar surface heat flux !
             !! Updating the values of the skin temperature T_s and q_s :
             CALL CSWL_ECMWF( Qsw, ztmp1, u_star, zsst, T_s ) ! yes ECMWF, because more advanced than COARE (warm-layer added!)