Ignore:
Timestamp:
2015-06-29T17:19:38+02:00 (5 years ago)
Author:
clevy
Message:

eORCA1 hard coded indexes, see ticket #1554

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

    r5385 r5506  
    402402      IF( cp_cfg == "orca" .AND. jp_cfg == 1 ) THEN   ! ORCA R1 configuration 
    403403         !                                                 ! Increased lateral friction near of some straits 
    404          ! This dirty section will be suppressed by simplification process: all this will come back in input files 
    405          ! Currently these hard-wired indices relate to the original (pre-v3.6) configuration  
    406          ! which had a grid-size of 362x292. 
    407          ! This grid has been extended southwards for use with the under ice-shelf options (isf) introduced in v3.6.  
    408          ! The original domain can still be used optionally if the isf code is not activated.  
    409          ! An adjustment (isrow) is made to the hard-wired indices if the extended domain (362x332) is being used. 
    410          !  
    411          IF    ( jpjglo == 292 ) THEN  ;  isrow = 0  ! Using pre-v3.6 files or adjusted start row from isf-extended grid 
    412          ELSEIF( jpjglo == 332 ) THEN  ;  isrow = 40 ! Using full isf­extended domain.  
    413          ENDIF      
    414  
     404         ! This dirty section will be suppressed by simplification process: 
     405         ! all this will come back in input files 
     406         ! Currently these hard-wired indices relate to configuration with 
     407         ! extend grid (jpjglo=332) 
     408         ! 
     409         isrow = 332 - jpjglo 
     410         ! 
    415411         IF(lwp) WRITE(numout,*) 
    416412         IF(lwp) WRITE(numout,*) '   orca_r1: increase friction near the following straits : ' 
    417413         IF(lwp) WRITE(numout,*) '      Gibraltar ' 
    418414         ii0 = 282           ;   ii1 = 283        ! Gibraltar Strait  
    419          ij0 = 201 + isrow   ;   ij1 = 201 + isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     415         ij0 = 201 + isrow   ;   ij1 = 241 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    420416 
    421417         IF(lwp) WRITE(numout,*) '      Bhosporus ' 
    422418         ii0 = 314           ;   ii1 = 315        ! Bhosporus Strait  
    423          ij0 = 208 + isrow   ;   ij1 = 208 + isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     419         ij0 = 208 + isrow   ;   ij1 = 248 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    424420 
    425421         IF(lwp) WRITE(numout,*) '      Makassar (Top) ' 
    426422         ii0 =  48           ;   ii1 =  48        ! Makassar Strait (Top)  
    427          ij0 = 149 + isrow   ;   ij1 = 150 + isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
     423         ij0 = 149 + isrow   ;   ij1 = 190 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
    428424 
    429425         IF(lwp) WRITE(numout,*) '      Lombok ' 
    430426         ii0 =  44           ;   ii1 =  44        ! Lombok Strait  
    431          ij0 = 124 + isrow   ;   ij1 = 125 + isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     427         ij0 = 124 + isrow   ;   ij1 = 165 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    432428 
    433429         IF(lwp) WRITE(numout,*) '      Ombai ' 
    434430         ii0 =  53           ;   ii1 =  53        ! Ombai Strait  
    435          ij0 = 124 + isrow   ;   ij1 = 125 + isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     431         ij0 = 124 + isrow   ;   ij1 = 165 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    436432 
    437433         IF(lwp) WRITE(numout,*) '      Timor Passage ' 
    438434         ii0 =  56           ;   ii1 =  56        ! Timor Passage  
    439          ij0 = 124 + isrow   ;   ij1 = 125 + isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     435         ij0 = 124 + isrow   ;   ij1 = 165 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    440436 
    441437         IF(lwp) WRITE(numout,*) '      West Halmahera ' 
    442438         ii0 =  58           ;   ii1 =  58        ! West Halmahera Strait  
    443          ij0 = 141 + isrow   ;   ij1 = 142 + isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
     439         ij0 = 141 + isrow   ;   ij1 = 182 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
    444440 
    445441         IF(lwp) WRITE(numout,*) '      East Halmahera ' 
    446442         ii0 =  55           ;   ii1 =  55        ! East Halmahera Strait  
    447          ij0 = 141 + isrow   ;   ij1 = 142 + isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
     443         ij0 = 141 + isrow   ;   ij1 = 182 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
    448444         ! 
    449445      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.