Changeset 15540 for NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/USR/usrdef_hgr.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/USR/usrdef_hgr.F90 (modified) (1 diff)

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/USR/usrdef_hgr.F90

@@ -61 +61 @@
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: plamt
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: plamf
-      REAL(wp), DIMENSION(:,:), INTENT(out) :: plamu, plamv   ! longitude outputs                     [degrees]
+      REAL(dp), DIMENSION(:,:), INTENT(out) :: plamu, plamv   ! longitude outputs                     [degrees]
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: pphit
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: pphif
-      REAL(wp), DIMENSION(:,:), INTENT(out) :: pphiu, pphiv   ! latitude outputs                      [degrees]
+      REAL(dp), DIMENSION(:,:), INTENT(out) :: pphiu, pphiv   ! latitude outputs                      [degrees]
       INTEGER                 , INTENT(out) ::   kff                          ! =1 Coriolis parameter computed here, =0 otherwise
-      REAL(wp), DIMENSION(:,:), INTENT(out) ::   pff_f, pff_t                 ! Coriolis factor at f-point                [1/s]
+      REAL(dp), DIMENSION(:,:), INTENT(out) ::   pff_f, pff_t                 ! Coriolis factor at f-point                [1/s]
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: pe1t
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: pe1u
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: pe1f
-      REAL(wp), DIMENSION(:,:), INTENT(out) :: pe1v       ! i-scale factors                             [m]
+      REAL(dp), DIMENSION(:,:), INTENT(out) :: pe1v       ! i-scale factors                             [m]
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: pe2t
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: pe2v
       REAL(dp), DIMENSION(:,:), INTENT(out)  :: pe2f
-      REAL(wp), DIMENSION(:,:), INTENT(out) :: pe2u       ! j-scale factors                             [m]
+      REAL(dp), DIMENSION(:,:), INTENT(out) :: pe2u       ! j-scale factors                             [m]
       INTEGER                 , INTENT(out) ::   ke1e2u_v                     ! =1 u- & v-surfaces computed here, =0 otherwise 
-      REAL(wp), DIMENSION(:,:), INTENT(out) ::   pe1e2u, pe1e2v               ! u- & v-surfaces (if reduction in strait)   [m2]
+      REAL(dp), DIMENSION(:,:), INTENT(out) ::   pe1e2u, pe1e2v               ! u- & v-surfaces (if reduction in strait)   [m2]
       !
       INTEGER  ::   ji, jj               ! dummy loop indices
-      REAL(wp) ::   zlam1, zlam0, zcos_alpha, zim1 , zjm1 , ze1  , ze1deg, zf0 ! local scalars
-      REAL(wp) ::   zphi1, zphi0, zsin_alpha, zim05, zjm05, zbeta, znorme      !   -      -
+      REAL(dp) ::   zlam1, zlam0, zcos_alpha, zim1 , zjm1 , ze1  , ze1deg, zf0 ! local scalars
+      REAL(dp) ::   zphi1, zphi0, zsin_alpha, zim05, zjm05, zbeta, znorme      !   -      -
       !!-------------------------------------------------------------------------------
       !