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

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/ASM/asminc.F90

@@ -69 +69 @@
    INTEGER, PUBLIC :: nn_divdmp     !: Apply divergence damping filter nn_divdmp times
 
-   REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE ::   t_bkg   , s_bkg      !: Background temperature and salinity
-   REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE ::   u_bkg   , v_bkg      !: Background u- & v- velocity components
-   REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE ::   t_bkginc, s_bkginc   !: Increment to the background T & S
-   REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE ::   u_bkginc, v_bkginc   !: Increment to the u- & v-components
-   REAL(wp), PUBLIC, DIMENSION(:)    , ALLOCATABLE ::   wgtiau               !: IAU weights for each time step
+   REAL(dp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE ::   t_bkg   , s_bkg      !: Background temperature and salinity
+   REAL(dp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE ::   u_bkg   , v_bkg      !: Background u- & v- velocity components
+   REAL(dp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE ::   t_bkginc, s_bkginc   !: Increment to the background T & S
+   REAL(dp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE ::   u_bkginc, v_bkginc   !: Increment to the u- & v-components
+   REAL(dp), PUBLIC, DIMENSION(:)    , ALLOCATABLE ::   wgtiau               !: IAU weights for each time step
 #if defined key_asminc
    REAL(wp), PUBLIC, DIMENSION(:,:)  , ALLOCATABLE ::   ssh_iau              !: IAU-weighted sea surface height increment
@@ -85 +85 @@
    INTEGER , PUBLIC ::   niaufn      !: Type of IAU weighing function: = 0   Constant weighting
    !                                 !: = 1   Linear hat-like, centred in middle of IAU interval
-   REAL(wp), PUBLIC ::   salfixmin   !: Ensure that the salinity is larger than this  value if (ln_salfix)
-
-   REAL(wp), DIMENSION(:,:), ALLOCATABLE ::   ssh_bkg, ssh_bkginc   ! Background sea surface height and its increment
-   REAL(wp), DIMENSION(:,:), ALLOCATABLE ::   seaice_bkginc         ! Increment to the background sea ice conc
+   REAL(dp), PUBLIC ::   salfixmin   !: Ensure that the salinity is larger than this  value if (ln_salfix)
+
+   REAL(dp), DIMENSION(:,:), ALLOCATABLE ::   ssh_bkg, ssh_bkginc   ! Background sea surface height and its increment
+   REAL(dp), DIMENSION(:,:), ALLOCATABLE ::   seaice_bkginc         ! Increment to the background sea ice conc
 #if defined key_cice && defined key_asminc
    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::   ndaice_da             ! ice increment tendency into CICE
@@ -127 +127 @@
       REAL(KIND=dp) :: ditiaufin_date  ! Date YYYYMMDD.HHMMSS of IAU interval final time step
 
-      REAL(wp) :: znorm        ! Normalization factor for IAU weights
-      REAL(wp) :: ztotwgt      ! Value of time-integrated IAU weights (should be equal to one)
-      REAL(wp) :: z_inc_dateb  ! Start date of interval on which increment is valid
-      REAL(wp) :: z_inc_datef  ! End date of interval on which increment is valid
-      REAL(wp) :: zdate_bkg    ! Date in background state file for DI
-      REAL(wp) :: zdate_inc    ! Time axis in increments file
-      !
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   zhdiv   ! 2D workspace
+      REAL(dp) :: znorm        ! Normalization factor for IAU weights
+      REAL(dp) :: ztotwgt      ! Value of time-integrated IAU weights (should be equal to one)
+      REAL(dp) :: z_inc_dateb  ! Start date of interval on which increment is valid
+      REAL(dp) :: z_inc_datef  ! End date of interval on which increment is valid
+      REAL(dp) :: zdate_bkg    ! Date in background state file for DI
+      REAL(dp) :: zdate_inc    ! Time axis in increments file
+      !
+      REAL(dp), ALLOCATABLE, DIMENSION(:,:) ::   zhdiv   ! 2D workspace
       !!
       NAMELIST/nam_asminc/ ln_bkgwri,                                      &
@@ -520 +520 @@
       INTEGER  :: ji, jj, jk
       INTEGER  :: it
-      REAL(wp) :: zincwgt  ! IAU weight for current time step
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk) :: fzptnz ! 3d freezing point values
+      REAL(dp) :: zincwgt  ! IAU weight for current time step
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk) :: fzptnz ! 3d freezing point values
       !!----------------------------------------------------------------------
       !
@@ -528 +528 @@
       IF( ln_temnofreeze ) THEN
          DO jk = 1, jpkm1
-           CALL eos_fzp( CASTWP(pts(:,:,jk,jp_sal,Kmm)), fzptnz(:,:,jk), CASTWP(gdept(:,:,jk,Kmm)) )
+           CALL eos_fzp( pts(:,:,jk,jp_sal,Kmm), fzptnz(:,:,jk), gdept(:,:,jk,Kmm) )
          END DO
       ENDIF
@@ -613 +613 @@
             pts(:,:,:,:,Kbb) = pts(:,:,:,:,Kmm)                 ! Update before fields
 
-            CALL eos( CASTWP(pts(:,:,:,:,Kbb)), rhd, rhop, gdept_0(:,:,:) )  ! Before potential and in situ densities
+            CALL eos( pts(:,:,:,:,Kbb), rhd, rhop, gdept_0(:,:,:) )  ! Before potential and in situ densities
 !!gm  fabien
 !            CALL eos( pts(:,:,:,:,Kbb), rhd, rhop )                ! Before potential and in situ densities
@@ -654 +654 @@
       INTEGER :: ji, jj, jk
       INTEGER :: it
-      REAL(wp) :: zincwgt  ! IAU weight for current time step
+      REAL(dp) :: zincwgt  ! IAU weight for current time step
       !!----------------------------------------------------------------------
       !
@@ -729 +729 @@
       INTEGER :: it
       INTEGER :: ji, jj, jk
-      REAL(wp) :: zincwgt  ! IAU weight for current time step
+      REAL(dp) :: zincwgt  ! IAU weight for current time step
       !!----------------------------------------------------------------------
       !
@@ -813 +813 @@
       INTEGER, INTENT(IN) :: kt                               ! ocean time-step index
       INTEGER, INTENT(IN) :: Kbb, Kmm                         ! time level indices
-      REAL(wp), DIMENSION(:,:,:), INTENT(inout) ::   phdivn   ! horizontal divergence
+      REAL(dp), DIMENSION(:,:,:), INTENT(inout) ::   phdivn   ! horizontal divergence
       !!
       INTEGER  ::   ji, jj, jk                                ! dummy loop index
-      REAL(wp), DIMENSION(:,:)  , POINTER       ::   ztim     ! local array
+      REAL(dp), DIMENSION(:,:)  , POINTER       ::   ztim     ! local array
       !!----------------------------------------------------------------------
       !
@@ -858 +858 @@
       INTEGER  ::   ji, jj
       INTEGER  ::   it
-      REAL(wp) ::   zincwgt   ! IAU weight for current time step
+      REAL(dp) ::   zincwgt   ! IAU weight for current time step
 #if defined key_si3
       REAL(wp), DIMENSION(A2D(nn_hls)) ::   zofrld, zohicif, zseaicendg, zhicifinc