Changeset 14219 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/ISF/isfcpl.F90

Timestamp:

2020-12-18T18:52:57+01:00 (4 years ago)

Author:

mcastril

Message:

Add Mixed Precision support by Oriol Tintó

File:

• NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/ISF/isfcpl.F90 (modified) (6 diffs)

NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/ISF/isfcpl.F90

@@ -50 +50 @@
 #  include "do_loop_substitute.h90"
 #  include "domzgr_substitute.h90"
+#  include "single_precision_substitute.h90"
+
    !!----------------------------------------------------------------------
    !! NEMO/OCE 4.0 , NEMO Consortium (2018)
@@ -90 +92 @@
       IF(lwp) WRITE(numout,*) ' isfcpl_init:', id
       IF (id == 0) THEN
-         IF(lwp) WRITE(numout,*) ' isfcpl_init: restart variables for ice sheet coupling are missing, skip coupling for this leg '
+         IF(lwp) WRITE(numout,*) ' isfcpl_init: restart variables for ice sheet coupling are missing, skip coupling for this leg ' 
          IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~'
          IF(lwp) WRITE(numout,*) ''
@@ -214 +216 @@
       ssh(:,:,Kbb) = ssh(:,:,Kmm)
       !
-      IF ( ln_isfdebug ) CALL debug('isfcpl_ssh: sshn',ssh(:,:,Kmm))
+      IF ( ln_isfdebug ) CALL debug('isfcpl_ssh: sshn',CASTWP(ssh(:,:,Kmm)))
       !
       ! recompute the vertical scale factor, depth and water thickness
@@ -586 +588 @@
       ENDDO
       !
-      ! global
+      ! global 
       CALL mpp_sum('isfcpl',nisfl  )
       !
@@ -592 +594 @@
       ALLOCATE(zisfpts(nisfl(narea)))
       !
-      zisfpts(:) = isfcons(0,0,0,-HUGE(1.0), -HUGE(1.0), -HUGE(1.0), -HUGE(1.0), -HUGE(1.0), 0)
+      zisfpts(:) = isfcons(0,0,0,-HUGE(1.0_wp), -HUGE(1.0_wp), -HUGE(1.0_wp), -HUGE(1.0_wp), -HUGE(1.0_wp), 0)
       !
       ! start computing the correction and fill zisfpts
@@ -665 +667 @@
             ELSE
                iig  =0   ; ijg  =0   ; ik   =0
-               zdvol=-HUGE(1.0) ; zdsal=-HUGE(1.0) ; zdtem=-HUGE(1.0)
-               zlat =-HUGE(1.0) ; zlon =-HUGE(1.0)
+               zdvol=-HUGE(1.0_wp) ; zdsal=-HUGE(1.0_wp) ; zdtem=-HUGE(1.0)
+               zlat =-HUGE(1.0_wp) ; zlon =-HUGE(1.0_wp)
                ingb = 0
             END IF