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

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

@@ -154 +154 @@
          &                 surfdata%mi,   surfdata%mj,   & 
          &                 surfdata%rlam, surfdata%rphi, &
-         &                 CASTWP(glamt),        CASTWP(gphit),        &
+         &                 glamt,        gphit,        &
          &                 tmask(:,:,1), surfdata%nqc,  &
          &                 iosdsobs,     ilansobs,     &
@@ -279 +279 @@
       INTEGER, DIMENSION(imaxavtypes), OPTIONAL :: &
          & kdailyavtypes                          ! Types for daily averages
-      REAL(wp), INTENT(IN), DIMENSION(kpi,kpj,kpk,profdata%nvar) :: &
+      REAL(dp), INTENT(IN), DIMENSION(kpi,kpj,kpk,profdata%nvar) :: &
          & zmask
-      REAL(wp), INTENT(IN), DIMENSION(kpi,kpj,profdata%nvar) :: &
+      REAL(dp), INTENT(IN), DIMENSION(kpi,kpj,profdata%nvar) :: &
          & pglam, &
          & pgphi
@@ -609 +609 @@
       INTEGER :: iskip
       INTEGER :: idaystp
-      REAL(KIND=wp) :: zminstp
-      REAL(KIND=wp) :: zhoustp
-      REAL(KIND=wp) :: zobsstp 
+      REAL(KIND=dp) :: zminstp
+      REAL(KIND=dp) :: zhoustp
+      REAL(KIND=dp) :: zobsstp 
       INTEGER, DIMENSION(12) ::   imonth_len    !: length in days of the months of the current year
  
@@ -887 +887 @@
       INTEGER , INTENT(in   )                     ::   kpi    , kpj      ! Number of grid points in (i,j)
       INTEGER , INTENT(in   ), DIMENSION(kobsno)  ::   kobsi  , kobsj    ! Observation (i,j) coordinates
-      REAL(wp), INTENT(in   ), DIMENSION(kobsno)  ::   pobslam, pobsphi  ! Observation (lon,lat) coordinates
-      REAL(wp), INTENT(in   ), DIMENSION(kpi,kpj) ::   plam   , pphi     ! Model (lon,lat) coordinates
-      REAL(wp), INTENT(in   ), DIMENSION(kpi,kpj) ::   pmask             ! Land mask array
+      REAL(dp), INTENT(in   ), DIMENSION(kobsno)  ::   pobslam, pobsphi  ! Observation (lon,lat) coordinates
+      REAL(dp), INTENT(in   ), DIMENSION(kpi,kpj) ::   plam   , pphi     ! Model (lon,lat) coordinates
+      REAL(dp), INTENT(in   ), DIMENSION(kpi,kpj) ::   pmask             ! Land mask array
       INTEGER , INTENT(inout), DIMENSION(kobsno)  ::   kobsqc            ! Observation quality control
       INTEGER , INTENT(inout)                     ::   kosdobs           ! Observations outside space domain
@@ -899 +899 @@
       INTEGER , INTENT(in   )                     ::   kqc_cutoff        ! Cutoff QC value
       !
-      REAL(KIND=wp), DIMENSION(2,2,kobsno) ::   zgmsk          ! Grid mask
-      REAL(KIND=wp), DIMENSION(2,2,kobsno) ::   zbmsk          ! Boundary mask
-      REAL(KIND=wp), DIMENSION(jpi,jpj)    ::   zbdymask
-      REAL(KIND=wp), DIMENSION(2,2,kobsno) ::   zglam, zgphi   ! Model Lon/lat at grid points
+      REAL(KIND=dp), DIMENSION(2,2,kobsno) ::   zgmsk          ! Grid mask
+      REAL(KIND=dp), DIMENSION(2,2,kobsno) ::   zbmsk          ! Boundary mask
+      REAL(KIND=dp), DIMENSION(jpi,jpj)    ::   zbdymask
+      REAL(KIND=dp), DIMENSION(2,2,kobsno) ::   zglam, zgphi   ! Model Lon/lat at grid points
       INTEGER      , DIMENSION(2,2,kobsno) ::   igrdi, igrdj   ! Grid i,j
       LOGICAL ::   lgridobs           ! Is observation on a model grid point.
@@ -1096 +1096 @@
       INTEGER, DIMENSION(kobsno), INTENT(IN) :: &
          & kobsk              ! Observation k coordinate
-      REAL(KIND=wp), DIMENSION(kprofno), INTENT(IN) :: &
+      REAL(KIND=dp), DIMENSION(kprofno), INTENT(IN) :: &
          & pobslam, &         ! Observation (lon,lat) coordinates
          & pobsphi
-      REAL(KIND=wp), DIMENSION(kobsno), INTENT(INOUT) :: &
+      REAL(KIND=dp), DIMENSION(kobsno), INTENT(INOUT) :: &
          & pobsdep            ! Observation depths  
-      REAL(KIND=wp), DIMENSION(kpi,kpj), INTENT(IN) :: &
+      REAL(KIND=dp), DIMENSION(kpi,kpj), INTENT(IN) :: &
          & plam, pphi         ! Model (lon,lat) coordinates
-      REAL(KIND=wp), DIMENSION(kpk), INTENT(IN) :: &
+      REAL(KIND=dp), DIMENSION(kpk), INTENT(IN) :: &
          & pdep               ! Model depth coordinates
-      REAL(KIND=wp), DIMENSION(kpi,kpj,kpk), INTENT(IN) :: &
+      REAL(KIND=dp), DIMENSION(kpi,kpj,kpk), INTENT(IN) :: &
          & pmask              ! Land mask array
       INTEGER, DIMENSION(kprofno), INTENT(INOUT) :: &
@@ -1121 +1121 @@
 
       !! * Local declarations
-      REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: &
+      REAL(KIND=dp), DIMENSION(2,2,kpk,kprofno) :: &
          & zgmsk              ! Grid mask
-      REAL(KIND=wp), DIMENSION(2,2,kprofno) :: &
+      REAL(KIND=dp), DIMENSION(2,2,kprofno) :: &
          & zbmsk              ! Boundary mask
-      REAL(KIND=wp), DIMENSION(jpi,jpj) :: zbdymask
-      REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: &
+      REAL(KIND=dp), DIMENSION(jpi,jpj) :: zbdymask
+      REAL(KIND=dp), DIMENSION(2,2,kpk,kprofno) :: &
          & zgdepw
-      REAL(KIND=wp), DIMENSION(2,2,kprofno) :: &
+      REAL(KIND=dp), DIMENSION(2,2,kprofno) :: &
          & zglam, &           ! Model longitude at grid points
          & zgphi              ! Model latitude at grid points
-      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zdepw
+      REAL(dp), DIMENSION(jpi,jpj,jpk) :: zdepw
       INTEGER, DIMENSION(2,2,kprofno) :: &
          & igrdi, &           ! Grid i,j