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

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/TRD/trdglo.F90

@@ -38 +38 @@
 
    !                     !!! Variables used for diagnostics
-   REAL(wp) ::   tvolt    ! volume of the whole ocean computed at t-points
-   REAL(wp) ::   tvolu    ! volume of the whole ocean computed at u-points
-   REAL(wp) ::   tvolv    ! volume of the whole ocean computed at v-points
-   REAL(wp) ::   rpktrd   ! potential to kinetic energy conversion
-   REAL(wp) ::   peke     ! conversion potential energy - kinetic energy trend
+   REAL(dp) ::   tvolt    ! volume of the whole ocean computed at t-points
+   REAL(dp) ::   tvolu    ! volume of the whole ocean computed at u-points
+   REAL(dp) ::   tvolv    ! volume of the whole ocean computed at v-points
+   REAL(dp) ::   rpktrd   ! potential to kinetic energy conversion
+   REAL(dp) ::   peke     ! conversion potential energy - kinetic energy trend
 
    !                     !!! domain averaged trends
-   REAL(wp), DIMENSION(jptot_tra) ::   tmo, smo   ! temperature and salinity trends 
-   REAL(wp), DIMENSION(jptot_tra) ::   t2 , s2    ! T^2 and S^2 trends 
-   REAL(wp), DIMENSION(jptot_dyn) ::   umo, vmo   ! momentum trends 
-   REAL(wp), DIMENSION(jptot_dyn) ::   hke        ! kinetic energy trends (u^2+v^2) 
+   REAL(dp), DIMENSION(jptot_tra) ::   tmo, smo   ! temperature and salinity trends 
+   REAL(dp), DIMENSION(jptot_tra) ::   t2 , s2    ! T^2 and S^2 trends 
+   REAL(dp), DIMENSION(jptot_dyn) ::   umo, vmo   ! momentum trends 
+   REAL(dp), DIMENSION(jptot_dyn) ::   hke        ! kinetic energy trends (u^2+v^2) 
 
    !! * Substitutions
@@ -78 +78 @@
       INTEGER ::   ji, jj, jk      ! dummy loop indices
       INTEGER ::   ikbu, ikbv      ! local integers
-      REAL(wp)::   zvm, zvt, zvs, z1_2rho0   ! local scalars
-      REAL(wp), DIMENSION(jpi,jpj)  :: ztswu, ztswv, z2dx, z2dy   ! 2D workspace 
+      REAL(dp)::   zvm, zvt, zvs, z1_2rho0   ! local scalars
+      REAL(dp), DIMENSION(jpi,jpj)  :: ztswu, ztswv, z2dx, z2dy   ! 2D workspace 
       !!----------------------------------------------------------------------
       !
@@ -186 +186 @@
       !
       INTEGER  ::   ji, jj, jk   ! dummy loop indices
-      REAL(wp) ::   zcof         ! local scalar
-      REAL(wp), DIMENSION(jpi,jpj,jpk)  ::  zkx, zky, zkz, zkepe  
+      REAL(dp) ::   zcof         ! local scalar
+      REAL(dp), DIMENSION(jpi,jpj,jpk)  ::  zkx, zky, zkz, zkepe  
       !!----------------------------------------------------------------------
 
@@ -203 +203 @@
          zkepe(:,:,:) = 0._wp
    
-         CALL eos( CASTWP(ts(:,:,:,:,Kmm)), rhd, CASTWP(rhop) )       ! now potential density
+         CALL eos( ts(:,:,:,:,Kmm), rhd, rhop )       ! now potential density
 
          zcof = 0.5_wp / rho0             ! Density flux at w-point