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

NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/ZDF/zdfmfc.F90

@@ -43 +43 @@
    PUBLIC   zdf_mfc_init    ! routine called in nemo module
    !
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  edmfa, edmfb, edmfc   !: diagonal term of the matrix.
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::  edmftra               !: y term for matrix inversion
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  edmfm               !: y term for matrix inversion
+   REAL(dp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  edmfa, edmfb, edmfc   !: diagonal term of the matrix.
+   REAL(dp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::  edmftra               !: y term for matrix inversion
+   REAL(dp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  edmfm               !: y term for matrix inversion
    !
    !! ** Namelist  namzdf_edmf  **
-   REAL(wp) ::   rn_cemf           ! entrain of T/S
-   REAL(wp) ::   rn_cwmf           ! detrain of T/S
-   REAL(wp) ::   rn_cent           ! entrain of the convective mass flux
-   REAL(wp) ::   rn_cdet           ! detrain of the convective mass flux
-   REAL(wp) ::   rn_cap            ! Factor of computation for convective area (negative => area constant)
-   REAL(wp) ::   App_max           ! Maximum of the convective area
+   REAL(dp) ::   rn_cemf           ! entrain of T/S
+   REAL(dp) ::   rn_cwmf           ! detrain of T/S
+   REAL(dp) ::   rn_cent           ! entrain of the convective mass flux
+   REAL(dp) ::   rn_cdet           ! detrain of the convective mass flux
+   REAL(dp) ::   rn_cap            ! Factor of computation for convective area (negative => area constant)
+   REAL(dp) ::   App_max           ! Maximum of the convective area
    LOGICAL, PUBLIC, SAVE  ::   ln_edmfuv         !: EDMF flag for velocity  !
    !
@@ -97 +97 @@
       INTEGER                                  , INTENT(in)    :: Kmm, Krhs ! time level indices
       REAL(dp), DIMENSION(jpi,jpj,jpk,jpts,jpt), INTENT(inout) :: pts       ! active tracers and RHS of tracer equation
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk,2) ::   ztsp         ! T/S of the plume
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk,2) ::   ztse         ! T/S at W point
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk) :: zrwp          !
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk) :: zrwp2         !
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk) :: zapp          !
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk) :: zedmf         !
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk) :: zepsT, zepsW  !
-      !
-      REAL(wp), DIMENSION(A2D(nn_hls)) :: zustar, zustar2   !
-      REAL(wp), DIMENSION(A2D(nn_hls)) :: zuws, zvws, zsws, zfnet          !
-      REAL(wp), DIMENSION(A2D(nn_hls)) :: zfbuo, zrautbm1
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk,2) ::   ztsp         ! T/S of the plume
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk,2) ::   ztse         ! T/S at W point
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk) :: zrwp          !
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk) :: zrwp2         !
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk) :: zapp          !
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk) :: zedmf         !
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk) :: zepsT, zepsW  !
+      !
+      REAL(dp), DIMENSION(A2D(nn_hls)) :: zustar, zustar2   !
+      REAL(dp), DIMENSION(A2D(nn_hls)) :: zuws, zvws, zsws, zfnet          !
+      REAL(dp), DIMENSION(A2D(nn_hls)) :: zfbuo, zrautbm1
       REAL(dp), DIMENSION(A2D(nn_hls)) :: zrautb, zraupl
-      REAL(wp), DIMENSION(A2D(nn_hls)) :: zwpsurf            !
-      REAL(wp), DIMENSION(A2D(nn_hls)) :: zop0 , zsp0 !
-      REAL(wp), DIMENSION(A2D(nn_hls)) :: zrwp_0, zrwp2_0  !
-      REAL(wp), DIMENSION(A2D(nn_hls)) :: zapp0           !
-      REAL(wp), DIMENSION(A2D(nn_hls)) :: zphp, zph, zphpm1, zphm1, zNHydro
-      REAL(wp), DIMENSION(A2D(nn_hls)) :: zhcmo          !
-      !
-      REAL(wp), DIMENSION(A2D(nn_hls),jpk)   ::   zn2    ! N^2
-      REAL(wp), DIMENSION(A2D(nn_hls),2  ) ::   zab, zabm1, zabp ! alpha and beta
+      REAL(dp), DIMENSION(A2D(nn_hls)) :: zwpsurf            !
+      REAL(dp), DIMENSION(A2D(nn_hls)) :: zop0 , zsp0 !
+      REAL(dp), DIMENSION(A2D(nn_hls)) :: zrwp_0, zrwp2_0  !
+      REAL(dp), DIMENSION(A2D(nn_hls)) :: zapp0           !
+      REAL(dp), DIMENSION(A2D(nn_hls)) :: zphp, zph, zphpm1, zphm1, zNHydro
+      REAL(dp), DIMENSION(A2D(nn_hls)) :: zhcmo          !
+      !
+      REAL(dp), DIMENSION(A2D(nn_hls),jpk)   ::   zn2    ! N^2
+      REAL(dp), DIMENSION(A2D(nn_hls),2  ) ::   zab, zabm1, zabp ! alpha and beta
      
-      REAL(wp), PARAMETER :: zepsilon = 1.e-30                 ! local small value
-
-      REAL(wp) :: zrho, zrhop
-      REAL(wp) :: zcnh, znum, zden, zcoef1, zcoef2
-      REAL(wp) :: zca, zcb, zcd, zrw, zxl, zcdet, zctre
-      REAL(wp) :: zaw, zbw, zxw
-      REAL(wp) :: alpha
+      REAL(dp), PARAMETER :: zepsilon = 1.e-30                 ! local small value
+
+      REAL(dp) :: zrho, zrhop
+      REAL(dp) :: zcnh, znum, zden, zcoef1, zcoef2
+      REAL(dp) :: zca, zcb, zcd, zrw, zxl, zcdet, zctre
+      REAL(dp) :: zaw, zbw, zxw
+      REAL(dp) :: alpha
      !
       INTEGER, INTENT(in   )    ::   kt   ! ocean time-step index      !
@@ -214 +214 @@
          ! Compute the buoyancy acceleration on T-points at jk-1
          zrautbm1(:,:) = zrautb(:,:)
-         CALL eos( CASTWP(pts (:,:,jk  ,:,Kmm)) ,  zrautb(:,:)   )
-         CALL eos( CASTWP(ztsp(:,:,jk-1,:    )) ,  zraupl(:,:)   )
+         CALL eos( pts (:,:,jk  ,:,Kmm) ,  zrautb(:,:)   )
+         CALL eos( ztsp(:,:,jk-1,:    ) ,  zraupl(:,:)   )
 
          DO_2D( 0, 0, 0, 0 )
@@ -399 +399 @@
 
       REAL(dp), DIMENSION(A2D(nn_hls),jpk), INTENT(inout) ::  zdiagi, zdiagd, zdiags  ! inout: tridaig. terms
-      REAL(wp)                            , INTENT(in   ) ::   p2dt                   ! tracer time-step
+      REAL(dp)                            , INTENT(in   ) ::   p2dt                   ! tracer time-step
       INTEGER                             , INTENT(in   ) ::   Kaa                    ! ocean time level indices
 
@@ -445 +445 @@
       INTEGER ::   jk    ! dummy loop indices
       INTEGER ::   ios   ! Local integer output status for namelist read
-      REAL(wp)::   zcr   ! local scalar
+      REAL(dp)::   zcr   ! local scalar
       !!
       NAMELIST/namzdf_mfc/ ln_edmfuv, rn_cemf, rn_cwmf, rn_cent, rn_cdet, rn_cap, App_max