Changeset 15540 for NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/SBC/sbcblk_skin_coare.F90

Timestamp:

2021-11-26T12:27:56+01:00 (3 years ago)

Author:

sparonuz

Message:

Mixed precision version, tested up to 30 years on ORCA2.

File:

• NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/SBC/sbcblk_skin_coare.F90 (modified) (5 diffs)

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/SBC/sbcblk_skin_coare.F90

@@ -36 +36 @@
 
    !! Cool-skin related parameters:
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: dT_cs !: dT due to cool-skin effect
+   REAL(dp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: dT_cs !: dT due to cool-skin effect
    !                                                      ! => temperature difference between air-sea interface (z=0)
    !                                                      !    and right below viscous layer (z=delta)
 
    !! Warm-layer related parameters:
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: dT_wl !: dT due to warm-layer effect
+   REAL(dp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: dT_wl !: dT due to warm-layer effect
    !                                                      ! => difference between "almost surface (right below
    !                                                      !    viscous layer, z=delta) and depth of bulk SST (z=gdept_1d(1))
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: Hz_wl !: depth (aka thickness) of warm-layer [m]
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: Qnt_ac !: time integral / accumulated heat stored by the warm layer
+   REAL(dp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: Hz_wl !: depth (aka thickness) of warm-layer [m]
+   REAL(dp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: Qnt_ac !: time integral / accumulated heat stored by the warm layer
    !                                                      !         Qxdt => [J/m^2] (reset to zero every midnight)
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: Tau_ac !: time integral / accumulated momentum
+   REAL(dp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC :: Tau_ac !: time integral / accumulated momentum
    !                                                      !         Tauxdt => [N.s/m^2] (reset to zero every midnight)
 
-   REAL(wp), PARAMETER, PUBLIC :: Hwl_max = 20._wp    !: maximum depth of warm layer (adjustable)
+   REAL(dp), PARAMETER, PUBLIC :: Hwl_max = 20._wp    !: maximum depth of warm layer (adjustable)
    !
-   REAL(wp), PARAMETER :: rich   = 0.65_wp   !: critical Richardson number
+   REAL(dp), PARAMETER :: rich   = 0.65_wp   !: critical Richardson number
    !
-   REAL(wp), PARAMETER :: zfr0   = 0.5_wp    !: initial value of solar flux absorption
+   REAL(dp), PARAMETER :: zfr0   = 0.5_wp    !: initial value of solar flux absorption
    !
    !!----------------------------------------------------------------------
@@ -81 +81 @@
       !!------------------------------------------------------------------
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: pQsw   ! net solar a.k.a shortwave radiation into the ocean (after albedo) [W/m^2]
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: pQnsol ! non-solar heat flux to the ocean [W/m^2]
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: pustar  ! friction velocity, temperature and humidity (u*,t*,q*)
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: pSST ! bulk SST [K]
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: pQlat  ! latent heat flux [W/m^2]
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(in) :: pQnsol ! non-solar heat flux to the ocean [W/m^2]
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(in) :: pustar  ! friction velocity, temperature and humidity (u*,t*,q*)
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(in) :: pSST ! bulk SST [K]
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(in) :: pQlat  ! latent heat flux [W/m^2]
       !!---------------------------------------------------------------------
       INTEGER  :: ji, jj, jc
-      REAL(wp) :: zQabs, zdlt, zfr, zalfa, zqlat, zus
+      REAL(dp) :: zQabs, zdlt, zfr, zalfa, zqlat, zus
       !!---------------------------------------------------------------------
       DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
@@ -131 +131 @@
       !!---------------------------------------------------------------------
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in)  :: pQsw     ! surface net solar radiation into the ocean [W/m^2]     => >= 0 !
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in)  :: pQnsol   ! surface net non-solar heat flux into the ocean [W/m^2] => normally < 0 !
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in)  :: pTau     ! wind stress [N/m^2]
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in)  :: pSST     ! bulk SST at depth gdept_1d(1) [K]
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(in)  :: pQnsol   ! surface net non-solar heat flux into the ocean [W/m^2] => normally < 0 !
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(in)  :: pTau     ! wind stress [N/m^2]
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(in)  :: pSST     ! bulk SST at depth gdept_1d(1) [K]
       INTEGER ,                     INTENT(in)  :: iwait    ! if /= 0 then wait before updating accumulated fluxes
       !!
       INTEGER :: ji,jj
       !
-      REAL(wp) :: zdTwl, zHwl, zQabs, zfr
-      REAL(wp) :: zqac, ztac
-      REAL(wp) :: zalfa, zcd1, zcd2, flg
-      !!---------------------------------------------------------------------
-
-      REAL(wp) :: ztime, znoon, zmidn
+      REAL(dp) :: zdTwl, zHwl, zQabs, zfr
+      REAL(dp) :: zqac, ztac
+      REAL(dp) :: zalfa, zcd1, zcd2, flg
+      !!---------------------------------------------------------------------
+
+      REAL(dp) :: ztime, znoon, zmidn
       INTEGER  :: jl
 
@@ -271 +271 @@
       !! L. Brodeau, october 2019
       !!---------------------------------------------------------------------
-      REAL(wp)                :: delta_skin_layer
-      REAL(wp), INTENT(in)    :: palpha   ! thermal expansion coefficient of sea-water (SST accurate enough!)
-      REAL(wp), INTENT(in)    :: pQd ! < 0 !!! part of the net heat flux actually absorbed in the WL [W/m^2]
+      REAL(dp)                :: delta_skin_layer
+      REAL(dp), INTENT(in)    :: palpha   ! thermal expansion coefficient of sea-water (SST accurate enough!)
+      REAL(dp), INTENT(in)    :: pQd ! < 0 !!! part of the net heat flux actually absorbed in the WL [W/m^2]
       !                              !  => term "Q + Rs*fs" in eq.6 of Fairall et al. 1996
-      REAL(wp), INTENT(in)    :: pQlat    ! latent heat flux [W/m^2]
-      REAL(wp), INTENT(in)    :: pustar_a ! friction velocity in the air (u*) [m/s]
-      !!---------------------------------------------------------------------
-      REAL(wp) :: zusw, zusw2, zlamb, zQd, ztf, ztmp
+      REAL(dp), INTENT(in)    :: pQlat    ! latent heat flux [W/m^2]
+      REAL(dp), INTENT(in)    :: pustar_a ! friction velocity in the air (u*) [m/s]
+      !!---------------------------------------------------------------------
+      REAL(dp) :: zusw, zusw2, zlamb, zQd, ztf, ztmp
       !!---------------------------------------------------------------------
 
@@ -302 +302 @@
       !! Fraction of solar heat flux absorbed inside warm layer
       !!---------------------------------------------------------------------
-      REAL(wp)             :: frac_solar_abs
-      REAL(wp), INTENT(in) :: pHwl   ! thickness of warm-layer [m]
+      REAL(dp)             :: frac_solar_abs
+      REAL(dp), INTENT(in) :: pHwl   ! thickness of warm-layer [m]
       !!---------------------------------------------------------------------
       frac_solar_abs = 1._wp - ( 0.28*0.014  *(1._wp - EXP(-pHwl/0.014)) &