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 15540 for NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/TRA/tradmp.F90 – NEMO

Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/TRA/tradmp.F90

    r14986 r15540  
    4949   CHARACTER(LEN=200) , PUBLIC ::   cn_resto    !: name of netcdf file containing restoration coefficient field 
    5050   ! 
    51    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   resto    !: restoring coeff. on T and S (s-1) 
     51   REAL(dp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   resto    !: restoring coeff. on T and S (s-1) 
    5252 
    5353   !! * Substitutions 
     
    9898      INTEGER ::   ji, jj, jk, jn   ! dummy loop indices 
    9999      REAL(dp), DIMENSION(A2D(nn_hls),jpk,jpts)     ::  zts_dta 
    100       REAL(wp), DIMENSION(:,:,:)  , ALLOCATABLE ::  zwrk 
    101       REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE ::  ztrdts 
     100      REAL(dp), DIMENSION(:,:,:)  , ALLOCATABLE ::  zwrk 
     101      REAL(dp), DIMENSION(:,:,:,:), ALLOCATABLE ::  ztrdts 
    102102      !!---------------------------------------------------------------------- 
    103103      ! 
     
    175175      ENDIF 
    176176      !                           ! Control print 
    177 IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=CASTWP(pts(:,:,:,jp_tem,Krhs)), clinfo1=' dmp  - Ta: ', mask1=tmask, tab3d_2=CASTWP(pts(:,:,:,jp_sal,Krhs)), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
     177IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' dmp  - Ta: ', mask1=tmask, tab3d_2=pts(:,:,:,jp_sal,Krhs), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    178178 
    179179      ! 
Note: See TracChangeset for help on using the changeset viewer.