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 12962 for NEMO/branches/UKMO/dev_r12866_HPC-02_Daley_Tiling_trial_extra_halo/src/OCE/DOM/domain.F90 – NEMO

Ignore:
Timestamp:
2020-05-22T16:18:18+02:00 (4 years ago)
Author:
hadcv
Message:

Update with [12960] fixes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/dev_r12866_HPC-02_Daley_Tiling_trial_extra_halo/src/OCE/DOM/domain.F90

    r12958 r12962  
    217217      !!              - mig0, mjg0: local  domain indices ==> global domain, excluding halos, indices 
    218218      !!              - mi0 , mi1 : global domain indices ==> local  domain indices 
    219       !!              - mj0,, mj1   (global point not in the local domain ==> mi0>mi1 and/or mj0>mj1) 
     219      !!              - mj0 , mj1   (if global point not in the local domain ==> mi0>mi1 and/or mj0>mj1) 
    220220      !!---------------------------------------------------------------------- 
    221221      INTEGER ::   ji, jj   ! dummy loop argument 
     
    230230      !                              ! local domain indices ==> global domain, excluding halos, indices 
    231231      ! 
     232      mig0(:) = mig(:) - nn_hls 
     233      mjg0(:) = mjg(:) - nn_hls   
    232234      ! WARNING: to keep compatibility with the trunk that was including periodocity into the input data,  
    233235      ! we must define mig0 and mjg0 as bellow. 
    234236      ! Once we decide to forget trunk compatibility, we must simply define mig0 and mjg0 as: 
    235 !!$      mig0(:) = mig(:) - nn_hls 
    236 !!$      mjg0(:) = mjg(:) - nn_hls   
    237       mig0(:) = mig(:) - nn_hls + COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 
    238       mjg0(:) = mjg(:) - nn_hls + COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) 
     237      mig0_oldcmp(:) = mig0(:) + COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 
     238      mjg0_oldcmp(:) = mjg0(:) + COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) 
    239239      ! 
    240240      !                              ! global domain, including halos, indices ==> local domain indices 
Note: See TracChangeset for help on using the changeset viewer.