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

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

@@ -38 +38 @@
    !! ---------------------------------------------------------------------
 
-   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) ::   wdmask   !: u- and v- limiter 
+   REAL(dp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) ::   wdmask   !: u- and v- limiter 
    !                                                           !  (can include negative depths)
-   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) ::   wdramp, wdrampu, wdrampv !: for hpg limiting
+   REAL(dp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) ::   wdramp, wdrampu, wdrampv !: for hpg limiting
 
    LOGICAL,  PUBLIC  ::   ln_wd_il    !: Wetting/drying il activation switch (T:on,F:off)
@@ -46 +46 @@
    REAL(wp), PUBLIC  ::   rn_wdmin0   !: depth at which wetting/drying starts
    REAL(wp), PUBLIC  ::   rn_wdmin1   !: minimum water depth on dried cells
-   REAL(wp), PUBLIC  ::   r_rn_wdmin1 !: 1/minimum water depth on dried cells 
+   REAL(dp), PUBLIC  ::   r_rn_wdmin1 !: 1/minimum water depth on dried cells 
    REAL(wp), PUBLIC  ::   rn_wdmin2   !: tolerance of minimum water depth on dried cells
    REAL(wp), PUBLIC  ::   rn_wd_sbcdep   !: Depth at which to taper sbc fluxes
@@ -55 +55 @@
                                       !: where the flow is from wet points on less than half the barotropic sub-steps  
    LOGICAL,  PUBLIC  ::  ln_wd_dl_rmp !: use a ramp for the dl flux limiter between 2 rn_wdmin1 and rn_wdmin1 (rather than a cut-off at rn_wdmin1)      
-   REAL(wp), PUBLIC  ::   ssh_ref     !: height of z=0 with respect to the geoid; 
+   REAL(dp), PUBLIC  ::   ssh_ref     !: height of z=0 with respect to the geoid; 
 
    LOGICAL,  PUBLIC  ::   ll_wd = .FALSE. !: Wetting/drying activation switch (ln_wd_il or ln_wd_dl) <- default def if wad_init not called
@@ -134 +134 @@
       !!----------------------------------------------------------------------
       REAL(dp), DIMENSION(:,:)            , INTENT(inout) ::   psshb1
-      REAL(wp), DIMENSION(:,:)            , INTENT(in   ) ::   psshemp
-      REAL(wp)                            , INTENT(in   ) ::   z2dt
+      REAL(dp), DIMENSION(:,:)            , INTENT(in   ) ::   psshemp
+      REAL(dp)                            , INTENT(in   ) ::   z2dt
       INTEGER                             , INTENT(in   ) ::   Kmm       ! time level index
       REAL(dp), DIMENSION(jpi,jpj,jpk,jpt), INTENT(inout) ::   puu, pvv  ! velocity arrays
@@ -141 +141 @@
       INTEGER  ::   ji, jj, jk, jk1     ! dummy loop indices
       INTEGER  ::   jflag               ! local scalar
-      REAL(wp) ::   zcoef, zdep1, zdep2 ! local scalars
-      REAL(wp) ::   zzflxp, zzflxn      ! local scalars
-      REAL(wp) ::   zdepwd              ! local scalar, always wet cell depth
-      REAL(wp) ::   ztmp                ! local scalars
-      REAL(wp),  DIMENSION(jpi,jpj) ::   zwdlmtu, zwdlmtv   ! W/D flux limiters
-      REAL(wp),  DIMENSION(jpi,jpj) ::   zflxp  ,  zflxn    ! local 2D workspace
-      REAL(wp),  DIMENSION(jpi,jpj) ::   zflxu  ,  zflxv    ! local 2D workspace
-      REAL(wp),  DIMENSION(jpi,jpj) ::   zflxu1 , zflxv1    ! local 2D workspace
+      REAL(dp) ::   zcoef, zdep1, zdep2 ! local scalars
+      REAL(dp) ::   zzflxp, zzflxn      ! local scalars
+      REAL(dp) ::   zdepwd              ! local scalar, always wet cell depth
+      REAL(dp) ::   ztmp                ! local scalars
+      REAL(dp),  DIMENSION(jpi,jpj) ::   zwdlmtu, zwdlmtv   ! W/D flux limiters
+      REAL(dp),  DIMENSION(jpi,jpj) ::   zflxp  ,  zflxn    ! local 2D workspace
+      REAL(dp),  DIMENSION(jpi,jpj) ::   zflxu  ,  zflxv    ! local 2D workspace
+      REAL(dp),  DIMENSION(jpi,jpj) ::   zflxu1 , zflxv1    ! local 2D workspace
       !!----------------------------------------------------------------------
       IF( ln_timing )   CALL timing_start('wad_lmt')      !
@@ -282 +282 @@
       !! ** Action  : - calculate flux limiter and W/D flag
       !!----------------------------------------------------------------------
-      REAL(wp)                , INTENT(in   ) ::   rDt_e    ! ocean time-step index
+      REAL(dp)                , INTENT(in   ) ::   rDt_e    ! ocean time-step index
       REAL(dp), DIMENSION(:,:), INTENT(inout) ::   sshn_e, zssh_frc
-      REAL(wp), DIMENSION(:,:), INTENT(inout) ::   zflxu,  zflxv  
+      REAL(dp), DIMENSION(:,:), INTENT(inout) ::   zflxu,  zflxv  
       !
       INTEGER  ::   ji, jj, jk, jk1     ! dummy loop indices
       INTEGER  ::   jflag               ! local integer
-      REAL(wp) ::   z2dt
-      REAL(wp) ::   zcoef, zdep1, zdep2 ! local scalars
-      REAL(wp) ::   zzflxp, zzflxn      ! local scalars
-      REAL(wp) ::   zdepwd              ! local scalar, always wet cell depth
-      REAL(wp) ::   ztmp                ! local scalars
-      REAL(wp), DIMENSION(jpi,jpj) ::   zwdlmtu, zwdlmtv         !: W/D flux limiters
-      REAL(wp), DIMENSION(jpi,jpj) ::   zflxp,  zflxn            ! local 2D workspace
-      REAL(wp), DIMENSION(jpi,jpj) ::   zflxu1, zflxv1           ! local 2D workspace
+      REAL(dp) ::   z2dt
+      REAL(dp) ::   zcoef, zdep1, zdep2 ! local scalars
+      REAL(dp) ::   zzflxp, zzflxn      ! local scalars
+      REAL(dp) ::   zdepwd              ! local scalar, always wet cell depth
+      REAL(dp) ::   ztmp                ! local scalars
+      REAL(dp), DIMENSION(jpi,jpj) ::   zwdlmtu, zwdlmtv         !: W/D flux limiters
+      REAL(dp), DIMENSION(jpi,jpj) ::   zflxp,  zflxn            ! local 2D workspace
+      REAL(dp), DIMENSION(jpi,jpj) ::   zflxu1, zflxv1           ! local 2D workspace
       !!----------------------------------------------------------------------
       IF( ln_timing )   CALL timing_start('wad_lmt_bt')      !