Changeset 15540 for NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/LDF/ldfslp.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/LDF/ldfslp.F90 (modified) (7 diffs)

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/LDF/ldfslp.F90

@@ -56 +56 @@
    
    !                                                     !! Classic operator (Madec)
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)     ::   uslp, wslpi          !: i_slope at U- and W-points
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)     ::   vslp, wslpj          !: j-slope at V- and W-points
+   REAL(dp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)     ::   uslp, wslpi          !: i_slope at U- and W-points
+   REAL(dp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)     ::   vslp, wslpj          !: j-slope at V- and W-points
    !                                                     !! triad operator (Griffies)
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)     ::   wslp2                !: wslp**2 from Griffies quarter cells
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:,:) ::   triadi_g, triadj_g   !: skew flux  slopes relative to geopotentials
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:,:) ::   triadi  , triadj     !: isoneutral slopes relative to model-coordinate
+   REAL(dp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)     ::   wslp2                !: wslp**2 from Griffies quarter cells
+   REAL(dp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:,:) ::   triadi_g, triadj_g   !: skew flux  slopes relative to geopotentials
+   REAL(dp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:,:) ::   triadi  , triadj     !: isoneutral slopes relative to model-coordinate
    !                                                     !! both operators
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)     ::   ah_wslp2             !: ah * slope^2 at w-point
@@ -68 +68 @@
    !                                                     !! Madec operator
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   omlmask           ! mask of the surface mixed layer at T-pt
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   uslpml, wslpiml   ! i_slope at U- and W-points just below the mixed layer
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   vslpml, wslpjml   ! j_slope at V- and W-points just below the mixed layer
-
-   REAL(wp) ::   repsln = 1.e-25_wp       ! tiny value used as minium of di(rho), dj(rho) and dk(rho)
+   REAL(dp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   uslpml, wslpiml   ! i_slope at U- and W-points just below the mixed layer
+   REAL(dp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   vslpml, wslpjml   ! j_slope at V- and W-points just below the mixed layer
+
+   REAL(dp) ::   repsln = 1.e-25_wp       ! tiny value used as minium of di(rho), dj(rho) and dk(rho)
 
    !! * Substitutions
@@ -110 +110 @@
       INTEGER , INTENT(in)                   ::   kt    ! ocean time-step index
       INTEGER , INTENT(in)                   ::   Kbb, Kmm   ! ocean time level indices
-      REAL(wp), INTENT(in), DIMENSION(:,:,:) ::   prd   ! in situ density
-      REAL(wp), INTENT(in), DIMENSION(:,:,:) ::   pn2   ! Brunt-Vaisala frequency (locally ref.)
+      REAL(dp), INTENT(in), DIMENSION(:,:,:) ::   prd   ! in situ density
+      REAL(dp), INTENT(in), DIMENSION(:,:,:) ::   pn2   ! Brunt-Vaisala frequency (locally ref.)
       !!
       INTEGER  ::   ji , jj , jk    ! dummy loop indices
@@ -122 +122 @@
       REAL(wp) ::   zdepu, zdepv                   !   -      -
       REAL(wp), DIMENSION(jpi,jpj)     ::  zslpml_hmlpu, zslpml_hmlpv
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::  zgru, zwz, zdzr
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::  zgrv, zww
+      REAL(wp), DIMENSION(jpi,jpj,jpk)  :: zgru, zdzr
+      REAL(dp), DIMENSION(jpi,jpj,jpk)  :: zwz
+      REAL(wp), DIMENSION(jpi,jpj,jpk)  :: zgrv
+      REAL(dp), DIMENSION(jpi,jpj,jpk)  :: zww
       !!----------------------------------------------------------------------
       !
@@ -349 +351 @@
       INTEGER  ::   ji, jj, jk, jl, ip, jp, kp  ! dummy loop indices
       INTEGER  ::   iku, ikv                    ! local integer
-      REAL(wp) ::   zfacti, zfactj              ! local scalars
-      REAL(wp) ::   znot_thru_surface           ! local scalars
-      REAL(wp) ::   zdit, zdis, zdkt, zbu, zbti, zisw
-      REAL(wp) ::   zdjt, zdjs, zdks, zbv, zbtj, zjsw
-      REAL(wp) ::   zdxrho_raw, zti_coord, zti_raw, zti_lim, zti_g_raw, zti_g_lim
-      REAL(wp) ::   zdyrho_raw, ztj_coord, ztj_raw, ztj_lim, ztj_g_raw, ztj_g_lim
-      REAL(wp) ::   zdzrho_raw
-      REAL(wp) ::   zbeta0, ze3_e1, ze3_e2
-      REAL(wp), DIMENSION(jpi,jpj)     ::   z1_mlbw
-      REAL(wp), DIMENSION(jpi,jpj,jpk,0:1) ::   zdxrho , zdyrho, zdzrho     ! Horizontal and vertical density gradients
-      REAL(wp), DIMENSION(jpi,jpj,0:1,0:1) ::   zti_mlb, ztj_mlb            ! for Griffies operator only
+      REAL(dp) ::   zfacti, zfactj              ! local scalars
+      REAL(dp) ::   znot_thru_surface           ! local scalars
+      REAL(dp) ::   zdit, zdis, zdkt, zbu, zbti, zisw
+      REAL(dp) ::   zdjt, zdjs, zdks, zbv, zbtj, zjsw
+      REAL(dp) ::   zdxrho_raw, zti_coord, zti_raw, zti_lim, zti_g_raw, zti_g_lim
+      REAL(dp) ::   zdyrho_raw, ztj_coord, ztj_raw, ztj_lim, ztj_g_raw, ztj_g_lim
+      REAL(dp) ::   zdzrho_raw
+      REAL(dp) ::   zbeta0, ze3_e1, ze3_e2
+      REAL(dp), DIMENSION(jpi,jpj)     ::   z1_mlbw
+      REAL(dp), DIMENSION(jpi,jpj,jpk,0:1) ::   zdxrho , zdyrho, zdzrho     ! Horizontal and vertical density gradients
+      REAL(dp), DIMENSION(jpi,jpj,0:1,0:1) ::   zti_mlb, ztj_mlb            ! for Griffies operator only
       !!----------------------------------------------------------------------
       !
@@ -573 +575 @@
       !!                omlmask         :  mixed layer mask
       !!----------------------------------------------------------------------
-      REAL(wp), DIMENSION(:,:,:), INTENT(in) ::   prd            ! in situ density
-      REAL(wp), DIMENSION(:,:,:), INTENT(in) ::   pn2            ! Brunt-Vaisala frequency (locally ref.)
+      REAL(dp), DIMENSION(:,:,:), INTENT(in) ::   prd            ! in situ density
+      REAL(dp), DIMENSION(:,:,:), INTENT(in) ::   pn2            ! Brunt-Vaisala frequency (locally ref.)
       REAL(wp), DIMENSION(:,:,:), INTENT(in) ::   p_gru, p_grv   ! i- & j-gradient of density (u- & v-pts)
       REAL(wp), DIMENSION(:,:,:), INTENT(in) ::   p_dzr          ! z-gradient of density      (T-point)
@@ -582 +584 @@
       INTEGER  ::   iku, ikv, ik, ikm1             ! local integers
       REAL(wp) ::   zeps, zm1_g, zm1_2g, z1_slpmax ! local scalars
-      REAL(wp) ::   zci, zfi, zau, zbu, zai, zbi   !   -      -
-      REAL(wp) ::   zcj, zfj, zav, zbv, zaj, zbj   !   -      -
-      REAL(wp) ::   zck, zfk,      zbw             !   -      -
+      REAL(wp)  :: zci, zau, zbu, zai, zbi!   -      -
+      REAL(dp)  :: zfi                    !   -      -
+      REAL(wp)  :: zcj, zav, zbv, zaj, zbj!   -      -
+      REAL(dp)  :: zfj                    !   -      -
+      REAL(wp)  :: zbw                    !   -      -
+      REAL(dp)  :: zck, zfk               !   -      -
       !!----------------------------------------------------------------------
       !