Changeset 10673


Ignore:
Timestamp:
2019-02-13T14:27:25+01:00 (16 months ago)
Author:
mathiot
Message:

creation of branch for ticket #2238 suggestion 1

Location:
NEMO/branches/2019/fix_ticket2238_solution1
Files:
1 edited
1 copied

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/fix_ticket2238_solution1/src/OCE/ICB/icbutl.F90

    r10570 r10673  
    223223      ! 
    224224      ! find position in this processor. Prevent near edge problems (see #1389) 
    225       ! 
    226       IF    ( ii < mig( 1 ) ) THEN   ;   ii = 1 
    227       ELSEIF( ii > mig(jpi) ) THEN   ;   ii = jpi 
    228       ELSE                           ;   ii = mi1(ii) 
    229       ENDIF 
    230       IF    ( ij < mjg( 1 ) ) THEN   ;   ij = 1 
    231       ELSEIF( ij > mjg(jpj) ) THEN   ;   ij = jpj 
    232       ELSE                           ;   ij  = mj1(ij) 
    233       ENDIF 
    234       ! 
    235       IF( ii == jpi )   ii = ii-1       
    236       IF( ij == jpj )   ij = ij-1 
     225      ! (PM) will be useless if extra halo is used in NEMO 
     226      ! 
     227      IF    ( ii <= mig(1)-1 ) THEN   ;   ii = 0 
     228      ELSEIF( ii  > mig(jpi) ) THEN   ;   ii = jpi 
     229      ELSE                            ;   ii = mi1(ii) 
     230      ENDIF 
     231      IF    ( ij <= mjg(1)-1 ) THEN   ;   ij = 0 
     232      ELSEIF( ij  > mjg(jpj) ) THEN   ;   ij = jpj 
     233      ELSE                            ;   ij = mj1(ij) 
     234      ENDIF 
     235      ! 
    237236      ! 
    238237      icb_utl_bilin_h = ( pfld(ii,ij  ) * (1.-zi) + pfld(ii+1,ij  ) * zi ) * (1.-zj)   & 
Note: See TracChangeset for help on using the changeset viewer.