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

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

@@ -82 +82 @@
       INTEGER, DIMENSION(knumtypes) :: &
          & ibiastypes             ! Array of the bias types in each file
-      REAL(wp), DIMENSION(jpi,jpj,knumtypes) :: & 
+      REAL(dp), DIMENSION(jpi,jpj,knumtypes) :: & 
          & z_sstbias              ! Array to store the SST bias values
-      REAL(wp), DIMENSION(jpi,jpj) :: & 
+      REAL(dp), DIMENSION(jpi,jpj) :: & 
          & z_sstbias_2d           ! Array to store the SST bias values   
-      REAL(wp), DIMENSION(1) :: &
+      REAL(dp), DIMENSION(1) :: &
          & zext, &
          & zobsmask
-      REAL(wp), DIMENSION(2,2,1) :: &
+      REAL(dp), DIMENSION(2,2,1) :: &
          & zweig
-      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: &
+      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE :: &
          & zmask, &
          & zglam, &
          & zgphi
-      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: &
+      REAL(dp), DIMENSION(:,:,:), ALLOCATABLE :: &
          & zmask_tmp, &
          & zglam_tmp, &
          & zgphi_tmp   
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::  zbias   
-      REAL(wp) :: zlam
-      REAL(wp) :: zphi
+      REAL(dp), ALLOCATABLE, DIMENSION(:,:,:) ::  zbias   
+      REAL(dp) :: zlam
+      REAL(dp) :: zphi
       INTEGER, DIMENSION(:,:,:), ALLOCATABLE :: &
          & igrdi, &
@@ -171 +171 @@
       END DO
       CALL obs_int_comm_2d( 2, 2, sstdata%nsurf, jpi, jpj, &
-         &                  igrdi, igrdj, CASTWP(glamt), zglam )
+         &                  igrdi, igrdj, glamt, zglam )
       CALL obs_int_comm_2d( 2, 2, sstdata%nsurf, jpi, jpj, &
-         &                  igrdi, igrdj, CASTWP(gphit), zgphi )
+         &                  igrdi, igrdj, gphit, zgphi )
       CALL obs_int_comm_2d( 2, 2, sstdata%nsurf, jpi, jpj, &
          &                  igrdi, igrdj, tmask(:,:,1), zmask )