Changeset 426


Ignore:
Timestamp:
2006-03-21T16:54:15+01:00 (15 years ago)
Author:
opalod
Message:

nemo_v1_bugfix_035 : CT : take into account the case with no split along the i/j-direction (nbondi/nbondj = 2)

Location:
trunk/NEMO
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/prtctl.F90

    r387 r426  
    145145               njctle = nlejtl(jn) + overlap * MIN( 1, nlcjtl(jn) - nlejtl(jn)) 
    146146               ! Do not take into account the bound of the domain 
    147                IF( ibonitl(jn) == -1 ) nictls = MAX(2, nictls) 
    148                IF( ibonitl(jn) ==  1 ) nictle = MIN(nictle, nleitl(jn) - 1) 
    149                IF( ibonjtl(jn) == -1 ) njctls = MAX(2, njctls) 
    150                IF( ibonjtl(jn) ==  1 ) njctle = MIN(njctle, nlejtl(jn) - 1) 
     147               IF( ibonitl(jn) == -1 .OR. ibonitl(jn) == 2 ) nictls = MAX(2, nictls) 
     148               IF( ibonjtl(jn) == -1 .OR. ibonjtl(jn) == 2 ) njctls = MAX(2, njctls) 
     149               IF( ibonitl(jn) ==  1 .OR. ibonitl(jn) == 2 ) nictle = MIN(nictle, nleitl(jn) - 1) 
     150               IF( ibonjtl(jn) ==  1 .OR. ibonjtl(jn) == 2 ) njctle = MIN(njctle, nlejtl(jn) - 1) 
    151151            ELSE 
    152152               nictls = MAX( 1, nimpptl(jn) + nlditl(jn) - 1 - overlap ) 
     
    155155               njctle = njmpptl(jn) + nlejtl(jn) - 1 + overlap * MIN( 1, nlcjtl(jn) - nlejtl(jn) )  
    156156               ! Do not take into account the bound of the domain 
    157                IF( ibonitl(jn) == -1 ) nictls = MAX(2, nictls) 
    158                IF( ibonjtl(jn) == -1 ) njctls = MAX(2, njctls) 
    159                IF( ibonitl(jn) ==  1 ) nictle = MIN(nictle, nimpptl(jn) + nleitl(jn) - 2) 
    160                IF( ibonjtl(jn) ==  1 ) njctle = MIN(njctle, njmpptl(jn) + nlejtl(jn) - 2) 
     157               IF( ibonitl(jn) == -1 .OR. ibonitl(jn) == 2 ) nictls = MAX(2, nictls) 
     158               IF( ibonjtl(jn) == -1 .OR. ibonjtl(jn) == 2 ) njctls = MAX(2, njctls) 
     159               IF( ibonitl(jn) ==  1 .OR. ibonitl(jn) == 2 ) nictle = MIN(nictle, nimpptl(jn) + nleitl(jn) - 2) 
     160               IF( ibonjtl(jn) ==  1 .OR. ibonjtl(jn) == 2 ) njctle = MIN(njctle, njmpptl(jn) + nlejtl(jn) - 2) 
    161161            ENDIF 
    162162         ENDIF 
  • trunk/NEMO/TOP_SRC/prtctl_trc.F90

    r331 r426  
    129129               njctle = nlejtl(js) + overlap * MIN( 1, nlcjtl(js) - nlejtl(js)) 
    130130               ! Do not take into account the bound of the domain 
    131                IF( ibonitl(js) == -1 ) nictls = MAX(2, nictls) 
    132                IF( ibonitl(js) ==  1 ) nictle = MIN(nictle, nleitl(js) - 1) 
    133                IF( ibonjtl(js) == -1 ) njctls = MAX(2, njctls) 
    134                IF( ibonjtl(js) ==  1 ) njctle = MIN(njctle, nlejtl(js) - 1) 
     131               IF( ibonitl(js) == -1 .OR. ibonitl(js) == 2 ) nictls = MAX(2, nictls) 
     132               IF( ibonitl(js) ==  1 .OR. ibonitl(js) == 2 ) nictle = MIN(nictle, nleitl(js) - 1) 
     133               IF( ibonjtl(js) == -1 .OR. ibonjtl(js) == 2 ) njctls = MAX(2, njctls) 
     134               IF( ibonjtl(js) ==  1 .OR. ibonjtl(js) == 2 ) njctle = MIN(njctle, nlejtl(js) - 1) 
    135135            ELSE 
    136136               nictls = MAX( 1, nimpptl(js) + nlditl(js) - 1 - overlap ) 
     
    139139               njctle = njmpptl(js) + nlejtl(js) - 1 + overlap * MIN( 1, nlcjtl(js) - nlejtl(js) )  
    140140               ! Do not take into account the bound of the domain 
    141                IF( ibonitl(js) == -1 ) nictls = MAX(2, nictls) 
    142                IF( ibonjtl(js) == -1 ) njctls = MAX(2, njctls) 
    143                IF( ibonitl(js) ==  1 ) nictle = MIN(nictle, nimpptl(js) + nleitl(js) - 2) 
    144                IF( ibonjtl(js) ==  1 ) njctle = MIN(njctle, njmpptl(js) + nlejtl(js) - 2) 
     141               IF( ibonitl(js) == -1 .OR. ibonitl(js) == 2 ) nictls = MAX(2, nictls) 
     142               IF( ibonjtl(js) == -1 .OR. ibonjtl(js) == 2 ) njctls = MAX(2, njctls) 
     143               IF( ibonitl(js) ==  1 .OR. ibonitl(js) == 2 ) nictle = MIN(nictle, nimpptl(js) + nleitl(js) - 2) 
     144               IF( ibonjtl(js) ==  1 .OR. ibonjtl(js) == 2 ) njctle = MIN(njctle, njmpptl(js) + nlejtl(js) - 2) 
    145145            ENDIF 
    146146         ENDIF 
Note: See TracChangeset for help on using the changeset viewer.