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

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/TRA/traqsr.F90

@@ -61 +61 @@
    !
    INTEGER  ::   nqsr    ! user choice of the type of light penetration
-   REAL(wp) ::   xsi0r   ! inverse of rn_si0
-   REAL(wp) ::   xsi1r   ! inverse of rn_si1
+   REAL(dp) ::   xsi0r   ! inverse of rn_si0
+   REAL(dp) ::   xsi1r   ! inverse of rn_si1
    !
-   REAL(wp) , PUBLIC, DIMENSION(3,61)   ::   rkrgb    ! tabulated attenuation coefficients for RGB absorption
+   REAL(dp) , PUBLIC, DIMENSION(3,61)   ::   rkrgb    ! tabulated attenuation coefficients for RGB absorption
    TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf_chl   ! structure of input Chl (file informations, fields read)
 
@@ -112 +112 @@
       REAL(wp) ::   zchl, zcoef, z1_2        ! local scalars
       REAL(wp) ::   zc0 , zc1 , zc2 , zc3    !    -         -
-      REAL(wp) ::   zzc0, zzc1, zzc2, zzc3   !    -         -
+      REAL(dp) ::   zzc0, zzc1, zzc2, zzc3   !    -         -
       REAL(wp) ::   zz0 , zz1 , ze3t, zlui   !    -         -
       REAL(wp) ::   zCb, zCmax, zpsi, zpsimax, zrdpsi, zCze
-      REAL(wp) ::   zlogc, zlogze, zlogCtot, zlogCze
+      REAL(wp)  :: zlogc
+      REAL(dp)  :: zlogze, zlogCtot, zlogCze
       REAL(wp), ALLOCATABLE, DIMENSION(:,:)   :: ze0, ze1, ze2, ze3
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrdt, zetot, ztmp3d
+      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:)  :: ztmp3d
+      REAL(dp), ALLOCATABLE, DIMENSION(:,:,:)  :: ztrdt, zetot
       !!----------------------------------------------------------------------
       !
@@ -317 +319 @@
       ENDIF
       !                       ! print mean trends (used for debugging)
-      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=CASTWP(pts(:,:,:,jp_tem,Krhs)), clinfo1=' qsr  - Ta: ', mask1=tmask, clinfo3='tra-ta' )
+      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' qsr  - Ta: ', mask1=tmask, clinfo3='tra-ta' )
       !
       IF( ln_timing )   CALL timing_stop('tra_qsr')
@@ -343 +345 @@
       INTEGER  ::   ji, jj, jk                  ! dummy loop indices
       INTEGER  ::   ios, irgb, ierror, ioptio   ! local integer
-      REAL(wp) ::   zz0, zc0 , zc1, zcoef      ! local scalars
-      REAL(wp) ::   zz1, zc2 , zc3, zchl       !   -      -
+      REAL(dp) ::   zz0, zc0 , zc1, zcoef      ! local scalars
+      REAL(dp) ::   zz1, zc2 , zc3, zchl       !   -      -
       !
       CHARACTER(len=100) ::   cn_dir   ! Root directory for location of ssr files
@@ -424 +426 @@
          IF(lwp)  WRITE(numout,*) '   ==>>>   2 bands light penetration'
          !
-         nksr = trc_oce_ext_lev( rn_si1, 100._wp )    ! level of light extinction
+         nksr =trc_oce_ext_lev( CASTDP(rn_si1), 100._wp )    ! level of light extinction
          IF(lwp) WRITE(numout,*) '        level of light extinction = ', nksr, ' ref depth = ', gdepw_1d(nksr+1), ' m'
          !