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 15574 for NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/lib_fortran_generic.h90 – NEMO

Ignore:
Timestamp:
2021-12-03T20:32:50+01:00 (3 years ago)
Author:
techene
Message:

#2605 #2715 trunk merged into dev_r14318_RK3_stage1

Location:
NEMO/branches/2021/dev_r14318_RK3_stage1
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2021/dev_r14318_RK3_stage1

    • Property svn:externals

      •  

        old new  
        99 
        1010# SETTE 
        11 ^/utils/CI/sette@14244        sette 
         11^/utils/CI/sette@HEAD        sette 
         12 

  • NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/lib_fortran_generic.h90

    r13226 r15574  
    3838      !!----------------------------------------------------------------------- 
    3939      ! 
    40       REAL(wp)                              ::   FUNCTION_GLOB_OP   ! global sum 
     40      REAL(wp) ::   FUNCTION_GLOB_OP   ! global sum 
    4141      !! 
    4242      COMPLEX(dp)::   ctmp 
    4343      REAL(wp)   ::   ztmp 
    44       INTEGER    ::   ji, jj, jk   ! dummy loop indices 
    45       INTEGER    ::   ipi, ipj, ipk    ! dimensions 
     44      INTEGER    ::   ji, jj, jk           ! dummy loop indices 
     45      INTEGER    ::   ipi,ipj, ipk         ! dimensions 
     46      INTEGER    ::   iis, iie, ijs, ije   ! loop start and end 
    4647      !!----------------------------------------------------------------------- 
    4748      ! 
    … …  
    5051      ipk = K_SIZE(ptab)   ! 3rd dimension 
    5152      ! 
     53      IF( ipi == jpi .AND. ipj == jpj ) THEN   ! do 2D loop only over the inner domain (-> avoid to use undefined values) 
     54         iis = Nis0   ;   iie = Nie0 
     55         ijs = Njs0   ;   ije = Nje0 
     56      ELSE 
     57         iis = 1   ;   iie = jpi 
     58         ijs = 1   ;   ije = jpj 
     59      ENDIF 
     60      ! 
    5261      ctmp = CMPLX( 0.e0, 0.e0, dp )   ! warning ctmp is cumulated 
    53     
    5462      DO jk = 1, ipk 
    55         DO jj = 1, ipj 
    56           DO ji = 1, ipi 
     63        DO jj = ijs, ije 
     64          DO ji = iis, iie 
    5765             ztmp =  ARRAY_IN(ji,jj,jk) * MASK_ARRAY(ji,jj) 
    5866             CALL DDPDD( CMPLX( ztmp, 0.e0, dp ), ctmp ) 
Note: See TracChangeset for help on using the changeset viewer.