Ignore:
Timestamp:
2015-06-05T13:34:40+02:00 (5 years ago)
Author:
cetlod
Message:

NEMOGCM_dev_r5204_CNRS_PISCES_dcy : merge in dev_r5171_CNRS_LIM3_seaicebgc

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

    r5348 r5367  
    436436            END DO 
    437437         END DO 
    438          IF( cp_cfg == 'orca' ) THEN 
    439             IF( jp_cfg == 2 ) THEN 
    440                ii0 = 176   ;   ii1 =  176        ! Southern Island : Kerguelen 
    441                ij0 =  37   ;   ij1 =   37  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    442                ! 
    443                ii0 = 119   ;   ii1 =  119        ! South Georgia 
    444                ij0 =  29   ;   ij1 =   29  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    445                ! 
    446                ii0 = 111   ;   ii1 =  111        ! Falklands 
    447                ij0 =  35   ;   ij1 =   35  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    448                ! 
    449                ii0 = 168   ;   ii1 =  168        ! Crozet 
    450                ij0 =  40   ;   ij1 =   40  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    451                ! 
    452                ii0 = 119   ;   ii1 =  119        ! South Orkney 
    453                ij0 =  28   ;   ij1 =   28  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    454                ! 
    455                ii0 = 140   ;   ii1 =  140        ! Bouvet Island 
    456                ij0 =  33   ;   ij1 =   33  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    457                ! 
    458                ii0 = 178   ;   ii1 =  178        ! Prince edwards 
    459                ij0 =  34   ;   ij1 =   34  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    460                ! 
    461                ii0 =  43   ;   ii1 =   43        ! Balleny islands 
    462                ij0 =  21   ;   ij1 =   21  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    463                ! 
    464             ELSE IF( jp_cfg == 1 ) THEN 
    465                ! 
    466                ! This dirty section will be suppressed by simplification process: all this will come back in input files 
    467                ! Currently these hard-wired indices relate to the original (pre-v3.6) configuration  
    468                ! which had a grid-size of 362x292. 
    469                ! This grid has been extended southwards for use with the under ice-shelf options (isf) introduced in v3.6.  
    470                ! The original domain can still be used optionally if the isf code is not activated.  
    471                ! An adjustment (isrow) is made to the hard-wired indices if the extended domain (362x332) is being used. 
    472                !  
    473                IF    ( jpjglo == 292 ) THEN  ;  isrow = 0  ! Using pre-v3.6 files or adjusted start row from isf-extended grid 
    474                ELSEIF( jpjglo == 332 ) THEN  ;  isrow = 40 ! Using full isf­extended domain.  
    475                ENDIF      
    476  
    477                ii0 = 357             ;   ii1 =  357        ! Southern Island : Kerguelen 
    478                ij0 =  75  + isrow   ;   ij1 =  76   + isrow   
    479                zcmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1), 1:jpk ) = 0.3_wp 
    480                ! 
    481                ii0 = 243             ;   ii1 =  243        ! South Georgia 
    482                ij0 =  57  + isrow   ;   ij1 =   59  + isrow  
    483                zcmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1), 1:jpk ) = 0.3_wp 
    484                ! 
    485                ii0 = 227             ;   ii1 =  227        ! Falklands 
    486                ij0 =  71  + isrow   ;   ij1 =   73  + isrow   
    487                zcmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1), 1:jpk ) = 0.3_wp 
    488                ! 
    489                ii0 = 326             ;   ii1 =  327        ! Crozet 
    490                ij0 =  79  + isrow   ;   ij1 =   79  + isrow  
    491                zcmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1), 1:jpk ) = 0.3_wp 
    492                !  
    493                ii0 = 243             ;   ii1 =  243        ! South Orkney 
    494                ij0 =  56  + isrow   ;   ij1 =   56  + isrow   
    495                zcmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1), 1:jpk ) = 0.3_wp 
    496                ! 
    497                ii0 = 283             ;   ii1 =  286        ! Bouvet Island 
    498                ij0 =  66  + isrow   ;   ij1 =   66  + isrow   
    499                zcmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1), 1:jpk ) = 0.3_wp 
    500                ! 
    501                ii0 = 361             ;   ii1 =  361        ! Prince edwards 
    502                ij0 =  67  + isrow   ;   ij1 =   68  + isrow  
    503                zcmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1), 1:jpk ) = 0.3_wp 
    504                ! 
    505                ii0 =  91             ;   ii1 =   92        ! Balleny islands 
    506                ij0 =  42  + isrow   ;   ij1 =   43  + isrow  
    507                zcmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1), 1:jpk ) = 0.3_wp 
    508                ! 
    509             ELSE IF( jp_cfg == 05 ) THEN 
    510                ! 
    511                CALL ctl_warn( ' Coastal supply of iron modifications for southern Islands in ORCA_R05' ) 
    512                ! 
    513             ELSE IF( jp_cfg == 025 ) THEN 
    514                ! 
    515                CALL ctl_warn( ' Coastal supply of iron modifications for southern Islands in ORCA_R025' ) 
    516                ! 
    517             ENDIF 
    518             ! 
    519          ENDIF 
     438         ! 
    520439         CALL lbc_lnk( zcmask , 'T', 1. )      ! lateral boundary conditions on cmask   (sign unchanged) 
     440         ! 
    521441         DO jk = 1, jpk 
    522442            DO jj = 1, jpj 
Note: See TracChangeset for help on using the changeset viewer.