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

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/DIU/diu_bulk.F90

@@ -26 +26 @@
 
    ! Parameters
-   REAL(wp), PRIVATE, PARAMETER :: pp_alpha = 2.0e-4_wp
-   REAL(wp), PRIVATE, PARAMETER :: pp_veltol = 0._wp
-   REAL(wp), PRIVATE, PARAMETER :: pp_min_fvel = 1.e-10_wp 
+   REAL(dp), PRIVATE, PARAMETER :: pp_alpha = 2.0e-4_wp
+   REAL(dp), PRIVATE, PARAMETER :: pp_veltol = 0._wp
+   REAL(dp), PRIVATE, PARAMETER :: pp_min_fvel = 1.e-10_wp 
    
    ! Key variables
-   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: x_dsst     ! Delta SST
-   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: x_solfrac  ! Fraction of 
+   REAL(dp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: x_dsst     ! Delta SST
+   REAL(dp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: x_solfrac  ! Fraction of 
    !                                                           ! absorbed radiation
 
@@ -96 +96 @@
       REAL(wp), DIMENSION(jpi,jpj)          , INTENT(in) ::   psolflux       ! solar flux (Watts)
       REAL(wp), DIMENSION(jpi,jpj)          , INTENT(in) ::   pqflux         ! heat (non-solar) flux (Watts)
-      REAL(wp), DIMENSION(jpi,jpj)          , INTENT(in) ::   ptauflux       ! wind stress  (kg/ m s^2)
+      REAL(dp), DIMENSION(jpi,jpj)          , INTENT(in) ::   ptauflux       ! wind stress  (kg/ m s^2)
       REAL(dp), DIMENSION(jpi,jpj)          , INTENT(in) ::   prho           ! water density  (kg/m^3)
       REAL(dp)                              , INTENT(in) ::   p_rdt          ! time-step
@@ -108 +108 @@
       INTEGER :: ji,jj
       LOGICAL  :: ll_calcfrac
-      REAL(wp), DIMENSION(jpi,jpj) :: z_fvel              ! friction velocity     
-      REAL(wp), DIMENSION(jpi,jpj) :: zthick, zcoolthick, zmu, zla
-      REAL(wp), DIMENSION(jpi,jpj) :: z_abflux            ! absorbed flux           
-      REAL(wp), DIMENSION(jpi,jpj) :: z_fla               ! Langmuir function value 
+      REAL(dp), DIMENSION(jpi,jpj) :: z_fvel              ! friction velocity     
+      REAL(dp), DIMENSION(jpi,jpj) :: zthick, zcoolthick, zmu, zla
+      REAL(dp), DIMENSION(jpi,jpj) :: z_abflux            ! absorbed flux           
+      REAL(dp), DIMENSION(jpi,jpj) :: z_fla               ! Langmuir function value 
       !!----------------------------------------------------------------------
 
@@ -177 +177 @@
       
       ! Function definition
-      REAL(wp), DIMENSION(jpi,jpj) :: t_imp
+      REAL(dp), DIMENSION(jpi,jpj) :: t_imp
       ! Dummy variables
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(IN) :: p_dsst     ! Delta SST
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(IN) :: p_dsst     ! Delta SST
       REAL(dp), INTENT(IN)                     :: p_rdt      ! Time-step
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(IN) :: p_abflux   ! Heat forcing
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(IN) :: p_fvel     ! Friction velocity
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(IN) :: p_fla      ! Langmuir number
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(IN) :: pmu        ! Structure parameter
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(IN) :: pthick     ! Layer thickness
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(IN) :: p_abflux   ! Heat forcing
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(IN) :: p_fvel     ! Friction velocity
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(IN) :: p_fla      ! Langmuir number
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(IN) :: pmu        ! Structure parameter
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(IN) :: pthick     ! Layer thickness
       REAL(dp), DIMENSION(jpi,jpj), INTENT(IN) :: prho       ! Water density
    
       ! Local variables
-      REAL(wp) :: z_olength          ! Obukhov length
-      REAL(wp) :: z_sigma, z_sigma2
-      REAL(wp) :: z_term1, z_term2     
-      REAL(wp) :: z_stabfunc         ! stability function value
-      REAL(wp) :: z_fvel      
+      REAL(dp) :: z_olength          ! Obukhov length
+      REAL(dp) :: z_sigma, z_sigma2
+      REAL(dp) :: z_term1, z_term2     
+      REAL(dp) :: z_stabfunc         ! stability function value
+      REAL(dp) :: z_fvel      
       
       CHARACTER(200) :: warn_string