New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 14219 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/DYN/sshwzv.F90 – NEMO

Ignore:
Timestamp:
2020-12-18T18:52:57+01:00 (4 years ago)
Author:
mcastril
Message:

Add Mixed Precision support by Oriol Tintó

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/DYN/sshwzv.F90

    r14200 r14219  
    5353#  include "do_loop_substitute.h90" 
    5454#  include "domzgr_substitute.h90" 
     55#  include "single_precision_substitute.h90" 
     56 
    5557   !!---------------------------------------------------------------------- 
    5658   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
     
    7678      INTEGER                         , INTENT(in   ) ::   kt             ! time step 
    7779      INTEGER                         , INTENT(in   ) ::   Kbb, Kmm, Kaa  ! time level index 
    78       REAL(wp), DIMENSION(jpi,jpj,jpt), INTENT(inout) ::   pssh           ! sea-surface height 
     80      REAL(dp), DIMENSION(jpi,jpj,jpt), INTENT(inout) ::   pssh           ! sea-surface height 
    7981      !  
    8082      INTEGER  ::   jk      ! dummy loop index 
     
    127129      !                                           !------------------------------! 
    128130      ! 
    129       IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab2d_1=pssh(:,:,Kaa), clinfo1=' pssh(:,:,Kaa)  - : ', mask1=tmask ) 
     131IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab2d_1=CASTWP(pssh(:,:,Kaa)), clinfo1=' pssh(:,:,Kaa)  - : ', mask1=tmask ) 
    130132      ! 
    131133      IF( ln_timing )   CALL timing_stop('ssh_nxt') 
     
    278280      INTEGER                         , INTENT(in   ) ::   kt             ! ocean time-step index 
    279281      INTEGER                         , INTENT(in   ) ::   Kbb, Kmm, Kaa  ! ocean time level indices 
    280       REAL(wp), DIMENSION(jpi,jpj,jpt)          , TARGET, INTENT(inout) ::   pssh           ! SSH field 
    281       REAL(wp), DIMENSION(jpi,jpj    ), OPTIONAL, TARGET, INTENT(  out) ::   pssh_f         ! filtered SSH field 
    282       ! 
    283       REAL(wp) ::   zcoef   ! local scalar 
    284       REAL(wp), POINTER, DIMENSION(:,:) ::   zssh   ! pointer for filtered SSH  
     282      REAL(dp), DIMENSION(jpi,jpj,jpt)          , TARGET, INTENT(inout) ::   pssh           ! SSH field 
     283      REAL(dp), DIMENSION(jpi,jpj    ), OPTIONAL, TARGET, INTENT(  out) ::   pssh_f         ! filtered SSH field 
     284      ! 
     285      REAL(dp) ::   zcoef   ! local scalar 
     286      REAL(dp), POINTER, DIMENSION(:,:) ::   zssh   ! pointer for filtered SSH  
    285287      !!---------------------------------------------------------------------- 
    286288      ! 
     
    314316      ENDIF 
    315317      ! 
    316       IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab2d_1=pssh(:,:,Kmm), clinfo1=' atf  - pssh(:,:,Kmm): ', mask1=tmask ) 
     318      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab2d_1=CASTWP(pssh(:,:,Kmm)), clinfo1=' pssh(:,:,Kmm)  - : ', mask1=tmask ) 
    317319      ! 
    318320      IF( ln_timing )   CALL timing_stop('ssh_atf') 
Note: See TracChangeset for help on using the changeset viewer.