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 13222 – NEMO

Changeset 13222


Ignore:
Timestamp:
2020-07-02T15:04:19+02:00 (4 years ago)
Author:
orioltp
Message:

dev_r12512_HPC-04_mcastril_Mixed_Precision_implementation: bugfix in mpp_lnk_icb_generic.h90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12512_HPC-04_mcastril_Mixed_Precision_implementation/src/OCE/LBC/mpp_lnk_icb_generic.h90

    r12632 r13222  
    6565         ! 
    6666      ELSE                                        !* closed 
    67          IF( .NOT. cd_type == 'F' )   pt2d(  1-kexti   :nn_hls   ,:) = 0._PRECISION    ! east except at F-point 
    68                                       pt2d(jpi-nn_hls+1:jpi+kexti,:) = 0._PRECISION    ! west 
     67# if defined SINGLE_PRECISION 
     68         IF( .NOT. cd_type == 'F' )   pt2d(  1-kexti   :nn_hls   ,:) = 0._sp    ! east except at F-point 
     69                                      pt2d(jpi-nn_hls+1:jpi+kexti,:) = 0._sp    ! west 
     70# else 
     71         IF( .NOT. cd_type == 'F' )   pt2d(  1-kexti   :nn_hls   ,:) = 0._dp    ! east except at F-point 
     72                                      pt2d(jpi-nn_hls+1:jpi+kexti,:) = 0._dp    ! west 
     73# endif 
    6974      ENDIF 
    7075      !                                      ! North-South boundaries 
     
    7378         pt2d(:,  jpj  :jpj+kextj) = pt2d(:,     2     :2+kextj)       ! south 
    7479      ELSE                                        !* closed 
    75          IF( .NOT. cd_type == 'F' )   pt2d(:,  1-kextj   :nn_hls   ) = 0._PRECISION    ! north except at F-point 
    76                                       pt2d(:,jpj-nn_hls+1:jpj+kextj) = 0._PRECISION    ! south 
     80# if defined SINGLE_PRECISION 
     81         IF( .NOT. cd_type == 'F' )   pt2d(:,  1-kextj   :nn_hls   ) = 0._sp    ! north except at F-point 
     82                                      pt2d(:,jpj-nn_hls+1:jpj+kextj) = 0._sp    ! south 
     83# else 
     84         IF( .NOT. cd_type == 'F' )   pt2d(:,  1-kextj   :nn_hls   ) = 0._dp    ! north except at F-point 
     85                                      pt2d(:,jpj-nn_hls+1:jpj+kextj) = 0._dp    ! south 
     86# endif 
    7787      ENDIF 
    7888      ! 
Note: See TracChangeset for help on using the changeset viewer.