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

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/DYN/dynadv_ubs.F90

@@ -27 +27 @@
    PRIVATE
 
-   REAL(wp), PARAMETER :: gamma1 = 1._wp/3._wp  ! =1/4 quick      ; =1/3  3rd order UBS
-   REAL(wp), PARAMETER :: gamma2 = 1._wp/32._wp ! =0   2nd order  ; =1/32 4th order centred
+   REAL(dp), PARAMETER :: gamma1 = 1._wp/3._wp  ! =1/4 quick      ; =1/3  3rd order UBS
+   REAL(dp), PARAMETER :: gamma2 = 1._wp/32._wp ! =0   2nd order  ; =1/32 4th order centred
 
    PUBLIC   dyn_adv_ubs   ! routine called by step.F90
@@ -75 +75 @@
       !
       INTEGER  ::   ji, jj, jk   ! dummy loop indices
-      REAL(wp) ::   zui, zvj, zfuj, zfvi, zl_u, zl_v   ! local scalars
+      REAL(dp) ::   zui, zvj, zfuj, zfvi, zl_u, zl_v   ! local scalars
       REAL(dp), DIMENSION(A2D(nn_hls),jpk)   :: zfu_t, zfu_uw
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk)   :: zfu_f, zfu
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk)   :: zfu_f, zfu
       REAL(dp), DIMENSION(A2D(nn_hls),jpk)    :: zfv_t, zfv_vw
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk)   :: zfv_f, zfv, zfw
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk,2) ::   zlu_uu, zlu_uv
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk,2) ::   zlv_vv, zlv_vu
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk)   :: zfv_f, zfv, zfw
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk,2) ::   zlu_uu, zlu_uv
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk,2) ::   zlv_vv, zlv_vu
       !!----------------------------------------------------------------------
       !
@@ -250 +250 @@
       ENDIF
       !                                         ! Control print
-      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=CASTWP(puu(:,:,:,Krhs)), clinfo1=' ubs2 adv - Ua: ', mask1=umask,   &
-         &                                  tab3d_2=CASTWP(pvv(:,:,:,Krhs)), clinfo2=           ' Va: ', mask2=vmask, clinfo3='dyn' )
+      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=puu(:,:,:,Krhs), clinfo1=' ubs2 adv - Ua: ', mask1=umask,   &
+         &                                  tab3d_2=pvv(:,:,:,Krhs), clinfo2=           ' Va: ', mask2=vmask, clinfo3='dyn' )
       !
    END SUBROUTINE dyn_adv_ubs