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

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/OBS/obs_oper.F90

@@ -108 +108 @@
       INTEGER       , INTENT(in   ) ::   kdaystp         ! Number of time steps per day
       INTEGER       , INTENT(in   ) ::   kvar            ! Number of variables in prodatqc
-      REAL(KIND=wp) , INTENT(in   ), DIMENSION(kpi,kpj,kpk) ::   pvar             ! Model field
-      REAL(KIND=wp) , INTENT(in   ), DIMENSION(kpi,kpj,kpk) ::   pmask            ! Land-sea mask
-      REAL(KIND=wp) , INTENT(in   ), DIMENSION(kpi,kpj)     ::   plam             ! Model longitude
-      REAL(KIND=wp) , INTENT(in   ), DIMENSION(kpi,kpj)     ::   pphi             ! Model latitudes
-      REAL(KIND=wp) , INTENT(in   ), DIMENSION(kpi,kpj,kpk) ::   pgdept, pgdepw   ! depth of T and W levels 
+      REAL(KIND=dp) , INTENT(in   ), DIMENSION(kpi,kpj,kpk) ::   pvar             ! Model field
+      REAL(KIND=dp) , INTENT(in   ), DIMENSION(kpi,kpj,kpk) ::   pmask            ! Land-sea mask
+      REAL(KIND=dp) , INTENT(in   ), DIMENSION(kpi,kpj)     ::   plam             ! Model longitude
+      REAL(KIND=dp) , INTENT(in   ), DIMENSION(kpi,kpj)     ::   pphi             ! Model latitudes
+      REAL(KIND=dp) , INTENT(in   ), DIMENSION(kpi,kpj,kpk) ::   pgdept, pgdepw   ! depth of T and W levels 
       INTEGER, DIMENSION(imaxavtypes), OPTIONAL ::   kdailyavtypes             ! Types for daily averages
 
@@ -135 +135 @@
       INTEGER, ALLOCATABLE, DIMENSION(:) :: iv_indic
 
-      REAL(KIND=wp) :: zlam
-      REAL(KIND=wp) :: zphi
-      REAL(KIND=wp) :: zdaystp
-      REAL(KIND=wp), DIMENSION(kpk) :: &
+      REAL(KIND=dp) :: zlam
+      REAL(KIND=dp) :: zphi
+      REAL(KIND=dp) :: zdaystp
+      REAL(KIND=dp), DIMENSION(kpk) :: &
          & zobsk,  &
          & zobs2k
-      REAL(KIND=wp), DIMENSION(2,2,1) :: &
+      REAL(KIND=dp), DIMENSION(2,2,1) :: &
          & zweig1, &
          & zweig
-      REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: &
+      REAL(dp), DIMENSION(:,:,:,:), ALLOCATABLE :: &
          & zmask,  &
          & zint,   &
@@ -150 +150 @@
          & zgdept, & 
          & zgdepw
-      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: &
+      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE :: &
          & zglam,  &
          & zgphi
-      REAL(KIND=wp), DIMENSION(1) :: zmsk
-      REAL(KIND=wp), DIMENSION(:,:,:), ALLOCATABLE :: interp_corner
+      REAL(KIND=dp), DIMENSION(1) :: zmsk
+      REAL(KIND=dp), DIMENSION(:,:,:), ALLOCATABLE :: interp_corner
 
       LOGICAL :: ld_dailyav
@@ -492 +492 @@
       INTEGER, INTENT(IN) :: kdaystp   ! Number of time steps per day
       INTEGER, INTENT(IN) :: k2dint    ! Horizontal interpolation type (see header)
-      REAL(wp), INTENT(IN), DIMENSION(kpi,kpj) :: &
+      REAL(dp), INTENT(IN), DIMENSION(kpi,kpj) :: &
          & psurf,  &                   ! Model surface field
          & psurfmask                   ! Land-sea mask
       LOGICAL, INTENT(IN) :: ldnightav ! Logical for averaging night-time data
-      REAL(KIND=wp), INTENT(IN) :: &
+      REAL(KIND=dp), INTENT(IN) :: &
          & plamscl, &                  ! Diameter in metres of obs footprint in E/W, N/S directions
          & pphiscl                     ! This is the full width (rather than half-width)
@@ -520 +520 @@
          & icount_night,      &
          & imask_night
-      REAL(wp) :: zlam
-      REAL(wp) :: zphi
-      REAL(wp), DIMENSION(1) :: zext, zobsmask
-      REAL(wp) :: zdaystp
-      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: &
+      REAL(dp) :: zlam
+      REAL(dp) :: zphi
+      REAL(dp), DIMENSION(1) :: zext, zobsmask
+      REAL(dp) :: zdaystp
+      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE :: &
          & zweig,  &
          & zmask,  &
@@ -535 +535 @@
          & zgphif
 
-      REAL(wp), DIMENSION(:,:), SAVE, ALLOCATABLE :: &
+      REAL(dp), DIMENSION(:,:), SAVE, ALLOCATABLE :: &
          & zintmp,  &
          & zouttmp, &
@@ -655 +655 @@
 
       CALL obs_int_comm_2d( imaxifp, imaxjfp, isurf, kpi, kpj, &
-         &                  igrdi, igrdj, CASTWP(glamt), zglam )
+         &                  igrdi, igrdj, glamt, zglam )
       CALL obs_int_comm_2d( imaxifp, imaxjfp, isurf, kpi, kpj, &
-         &                  igrdi, igrdj, CASTWP(gphit), zgphi )
+         &                  igrdi, igrdj, gphit, zgphi )
       CALL obs_int_comm_2d( imaxifp, imaxjfp, isurf, kpi, kpj, &
          &                  igrdi, igrdj, psurfmask, zmask )
@@ -663 +663 @@
          &                  igrdi, igrdj, psurf, zsurf )
       CALL obs_int_comm_2d( imaxifp+1, imaxjfp+1, isurf, kpi, kpj, &
-         &                  igrdip1, igrdjp1, CASTWP(glamf), zglamf )
+         &                  igrdip1, igrdjp1, glamf, zglamf )
       CALL obs_int_comm_2d( imaxifp+1, imaxjfp+1, isurf, kpi, kpj, &
-         &                  igrdip1, igrdjp1, CASTWP(gphif), zgphif )
+         &                  igrdip1, igrdjp1, gphif, zgphif )
 
       ! At the end of the day get interpolated means