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

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/ISF/isfload.F90

@@ -44 +44 @@
       !!--------------------------------------------------------------------
       INTEGER,                      INTENT(in   ) ::   Kmm        ! ocean time level index      
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(  out) ::   pisfload   ! ice shelf load
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(  out) ::   pisfload   ! ice shelf load
       !!----------------------------------------------------------------------
       !
@@ -73 +73 @@
       !!--------------------------------------------------------------------
       INTEGER,                      INTENT(in   ) ::   Kmm     ! ocean time level index      
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(  out) ::   pload   ! ice shelf load
+      REAL(dp), DIMENSION(jpi,jpj), INTENT(  out) ::   pload   ! ice shelf load
       !
       INTEGER  :: ji, jj, jk
       INTEGER  :: ikt
-      REAL(wp), DIMENSION(jpi,jpj)      :: zrhdtop_isf ! water density    displaced by the ice shelf (at the interface)
-      REAL(wp), DIMENSION(jpi,jpj,jpts) :: zts_top     ! water properties displaced by the ice shelf   
-      REAL(wp), DIMENSION(jpi,jpj,jpk)  :: zrhd        ! water density    displaced by the ice shelf
+      REAL(dp), DIMENSION(jpi,jpj)      :: zrhdtop_isf ! water density    displaced by the ice shelf (at the interface)
+      REAL(dp), DIMENSION(jpi,jpj,jpts) :: zts_top     ! water properties displaced by the ice shelf   
+      REAL(dp), DIMENSION(jpi,jpj,jpk)  :: zrhd        ! water density    displaced by the ice shelf
       !!----------------------------------------------------------------------
       !
@@ -86 +86 @@
       !
       DO jk = 1, jpk                   !- compute density of the water displaced by the ice shelf 
-         CALL eos( zts_top(:,:,:), CASTWP(gdept(:,:,jk,Kmm)), zrhd(:,:,jk) )
+         CALL eos( zts_top(:,:,:), gdept(:,:,jk,Kmm), zrhd(:,:,jk) )
 !!st ==>> CALL eos( zts_top(:,:,:), gdept_0(:,:,jk), zrhd(:,:,jk) )
       END DO
       !
       !                                !- compute rhd at the ice/oce interface (ice shelf side)
-      CALL eos( zts_top , risfdep, zrhdtop_isf )
+      CALL eos( zts_top , CASTDP(risfdep), zrhdtop_isf )
       !
       !                                !- Surface value + ice shelf gradient