Changeset 11209 for NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare3p5.F90

Timestamp:

2019-07-03T09:09:17+02:00 (5 years ago)

Author:

laurent

Message:

LB: making more efficient and more centralized use of physical/thermodynamics functions defined into "OCE/SBC/sbcblk_phy.F90".

File:

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

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

@@ -35 +35 @@
    USE sbc_ice         ! Surface boundary condition: ice fields
 #endif
-   !
+   USE sbcblk_phy     !LB: all thermodynamics functions, rho_air, q_sat, etc... #LB
    USE iom             ! I/O manager library
    USE lib_mpp         ! distribued memory computing library
@@ -255 +255 @@
    END SUBROUTINE turb_coare3p5
 
-
-
-   FUNCTION One_on_L( ptha, pqa, pus, pts, pqs )
-      !!------------------------------------------------------------------------
-      !!
-      !! Evaluates the 1./(Monin Obukhov length) from air temperature and
-      !!  specific humidity, and frictional scales u*, t* and q*
-      !!
-      !! Author: L. Brodeau, june 2016 / AeroBulk
-      !!         (https://github.com/brodeau/aerobulk/)
-      !!------------------------------------------------------------------------
-      REAL(wp), DIMENSION(jpi,jpj)             :: One_on_L         !: 1./(Monin Obukhov length) [m^-1]
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: ptha,  &  !: average potetntial air temperature [K]
-         &                                        pqa,   &  !: average specific humidity of air   [kg/kg]
-         &                                      pus, pts, pqs   !: frictional velocity, temperature and humidity
-      !
-      INTEGER  ::   ji, jj         ! dummy loop indices
-      REAL(wp) ::     zqa          ! local scalar
-      !
-      DO jj = 1, jpj
-         DO ji = 1, jpi
-            !
-            zqa = (1. + rctv0*pqa(ji,jj))
-            !
-            One_on_L(ji,jj) =  grav*vkarmn*(pts(ji,jj)*zqa + rctv0*ptha(ji,jj)*pqs(ji,jj)) &
-               &                      / ( pus(ji,jj)*pus(ji,jj) * ptha(ji,jj)*zqa )
-            !
-         END DO
-      END DO
-      !
-   END FUNCTION One_on_L
-
-
    FUNCTION psi_m_coare( pzeta )
       !!----------------------------------------------------------------------------------
@@ -389 +356 @@
    END FUNCTION psi_h_coare
 
-
-   FUNCTION visc_air( ptak )
-      !!---------------------------------------------------------------------
-      !! Air kinetic viscosity (m^2/s) given from temperature in degrees...
-      !!
-      !! Author: L. Brodeau, june 2016 / AeroBulk
-      !!         (https://github.com/brodeau/aerobulk/)
-      !!---------------------------------------------------------------------
-      REAL(wp), DIMENSION(jpi,jpj)             ::   visc_air
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) ::   ptak       ! air temperature   [K]
-      !
-      INTEGER  ::   ji, jj         ! dummy loop indices
-      REAL(wp) ::   ztc, ztc2      ! local scalar
-      !
-      DO jj = 1, jpj
-         DO ji = 1, jpi
-            ztc  = ptak(ji,jj) - rt0   ! air temp, in deg. C
-            ztc2 = ztc*ztc
-            visc_air(ji,jj) = 1.326E-5*(1. + 6.542E-3*ztc + 8.301E-6*ztc2 - 4.84E-9*ztc2*ztc)
-         END DO
-      END DO
-      !
-   END FUNCTION visc_air
-
    !!======================================================================
 END MODULE sbcblk_algo_coare3p5