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

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/SBC/sbcssm.F90

@@ -70 +70 @@
          ssu_m(:,:) = uu(:,:,1,Kbb)
          ssv_m(:,:) = vv(:,:,1,Kbb)
-         IF( l_useCT )  THEN    ;   sst_m(:,:) = eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) )
+         IF( l_useCT )  THEN    ;   sst_m(:,:) =eos_pt_from_ct( CASTDP(zts(:,:,jp_tem)), CASTDP(zts(:,:,jp_sal)) )
          ELSE                   ;   sst_m(:,:) = zts(:,:,jp_tem)
          ENDIF
@@ -93 +93 @@
             ssu_m(:,:) = zcoef * uu(:,:,1,Kbb)
             ssv_m(:,:) = zcoef * vv(:,:,1,Kbb)
-            IF( l_useCT )  THEN    ;   sst_m(:,:) = zcoef * eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) )
+            IF( l_useCT )  THEN    ;   sst_m(:,:) =zcoef * eos_pt_from_ct( CASTDP(zts(:,:,jp_tem)), CASTDP(zts(:,:,jp_sal)) )
             ELSE                    ;   sst_m(:,:) = zcoef * zts(:,:,jp_tem)
             ENDIF
@@ -121 +121 @@
          ssu_m(:,:) = ssu_m(:,:) + uu(:,:,1,Kbb)
          ssv_m(:,:) = ssv_m(:,:) + vv(:,:,1,Kbb)
-         IF( l_useCT )  THEN     ;   sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) )
+         IF( l_useCT )  THEN     ;   sst_m(:,:) =sst_m(:,:) + eos_pt_from_ct( CASTDP(zts(:,:,jp_tem)), CASTDP(zts(:,:,jp_sal)) )
          ELSE                    ;   sst_m(:,:) = sst_m(:,:) + zts(:,:,jp_tem)
          ENDIF
@@ -190 +190 @@
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in) ::   Kbb, Kmm   ! ocean time level indices
-      REAL(wp) ::   zcoef, zf_sbc   ! local scalar
+      REAL(wp)  :: zcoef! local scalar
+      REAL(dp)  :: zf_sbc! local scalar
       !!----------------------------------------------------------------------
       !
@@ -242 +243 @@
          ssu_m(:,:) = uu(:,:,1,Kbb)
          ssv_m(:,:) = vv(:,:,1,Kbb)
-         IF( l_useCT )  THEN    ;   sst_m(:,:) = eos_pt_from_ct( CASTWP(ts(:,:,1,jp_tem,Kmm)), CASTWP(ts(:,:,1,jp_sal,Kmm)) )
+         IF( l_useCT )  THEN    ;   sst_m(:,:) = eos_pt_from_ct( ts(:,:,1,jp_tem,Kmm), ts(:,:,1,jp_sal,Kmm) )
          ELSE                   ;   sst_m(:,:) = ts(:,:,1,jp_tem,Kmm)
          ENDIF