Changeset 5315


Ignore:
Timestamp:
2015-05-29T12:28:06+02:00 (5 years ago)
Author:
cetlod
Message:

dev_r5204_CNRS_PISCES_dcy : An adjustment (isrow) is made to the hard-wired

Location:
branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/OPA_SRC/DIA/diafwb.F90

    r4990 r5315  
    5151      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
    5252      !! 
    53       INTEGER :: inum             ! temporary logical unit 
    54       INTEGER :: ji, jj, jk, jt   ! dummy loop indices 
    55       INTEGER :: ii0, ii1, ij0, ij1 
    56       REAL(wp) ::   zarea, zvol, zwei 
    57       REAL(wp) ::  ztemi(4), ztemo(4), zsali(4), zsalo(4), zflxi(4), zflxo(4) 
    58       REAL(wp) ::  zt, zs, zu   
    59       REAL(wp) ::  zsm0, zfwfnew 
     53      INTEGER  :: inum             ! temporary logical unit 
     54      INTEGER  :: ji, jj, jk, jt   ! dummy loop indices 
     55      INTEGER  :: ii0, ii1, ij0, ij1 
     56      INTEGER  :: isrow         ! index for ORCA1 starting row 
     57      REAL(wp) :: zarea, zvol, zwei 
     58      REAL(wp) :: ztemi(4), ztemo(4), zsali(4), zsalo(4), zflxi(4), zflxo(4) 
     59      REAL(wp) :: zt, zs, zu   
     60      REAL(wp) :: zsm0, zfwfnew 
    6061      IF( cp_cfg == "orca" .AND. jp_cfg == 1 .OR. jp_cfg == 2 .OR. jp_cfg == 4 ) THEN 
    6162      !!---------------------------------------------------------------------- 
     
    165166         CASE ( 1 )                                  !  ORCA_R1 configurations 
    166167            !                                        ! ======================= 
    167             ii0 = 283   ;   ii1 = 283 
    168             ij0 = 200   ;   ij1 = 200 
     168            ! This dirty section will be suppressed by simplification process: all this will come back in input files 
     169            ! Currently these hard-wired indices relate to the original (pre-v3.6) configuration  
     170            ! which had a grid-size of 362x292. 
     171            ! This grid has been extended southwards for use with the under ice-shelf options (isf) introduced in v3.6.  
     172            ! The original domain can still be used optionally if the isf code is not activated.  
     173            ! An adjustment (isrow) is made to the hard-wired indices if the extended domain (362x332) is being used. 
     174            !  
     175            IF    ( jpjglo == 292 ) THEN  ;  isrow = 0  ! Using pre-v3.6 files or adjusted start row from isf-extended grid 
     176            ELSEIF( jpjglo == 332 ) THEN  ;  isrow = 40 ! Using full isf­extended domain.  
     177            ENDIF                                       ! Adjust j­indices to account for more southerly starting latitude 
     178            ii0 = 283           ;   ii1 = 283 
     179            ij0 = 201 + isrow   ;   ij1 = 201 + isrow 
    169180            !                                        ! ======================= 
    170181         CASE DEFAULT                                !    ORCA R05 or R025 
     
    212223         CASE ( 1 )                                  !  ORCA_R1 configurations 
    213224            !                                        ! ======================= 
    214             ii0 = 282   ;   ii1 = 282 
    215             ij0 = 200   ;   ij1 = 200 
     225            !  This dirty section will be suppressed by simplification process: 
     226            !  all this will come back in input files 
     227            IF    ( jpjglo == 292 ) THEN  ;  isrow = 0   
     228            ELSEIF( jpjglo == 332 ) THEN  ;  isrow = 39 
     229            ENDIF                                       
     230            ii0 = 282           ;   ii1 = 282 
     231            ij0 = 201 + isrow   ;   ij1 = 201 + isrow 
    216232            !                                        ! ======================= 
    217233         CASE DEFAULT                                !    ORCA R05 or R025 
     
    259275         CASE ( 1 )                                  !  ORCA_R1 configurations 
    260276            !                                        ! ======================= 
    261             ii0 = 331   ;   ii1 = 331 
    262             ij0 = 176   ;   ij1 = 176 
     277            !  This dirty section will be suppressed by simplification process: 
     278            !  all this will come back in input files 
     279            IF    ( jpjglo == 292 ) THEN  ;  isrow = 0   
     280            ELSEIF( jpjglo == 332 ) THEN  ;  isrow = 39 
     281            ENDIF                                       
     282            ii0 = 331           ;   ii1 = 331 
     283            ij0 = 176 + isrow   ;   ij1 = 176 + isrow 
    263284            !                                        ! ======================= 
    264285         CASE DEFAULT                                !    ORCA R05 or R025 
     
    306327         CASE ( 1 )                                  !  ORCA_R1 configurations 
    307328            !                                        ! ======================= 
    308             ii0 = 297   ;   ii1 = 297  
    309             ij0 = 230   ;   ij1 = 230 
     329            !  This dirty section will be suppressed by simplification process: 
     330            !  all this will come back in input files 
     331            IF    ( jpjglo == 292 ) THEN  ;  isrow = 0   
     332            ELSEIF( jpjglo == 332 ) THEN  ;  isrow = 39 
     333            ENDIF                                       
     334            ii0 = 297           ;   ii1 = 297 
     335            ij0 = 230 + isrow   ;   ij1 = 230 + isrow 
    310336            !                                        ! ======================= 
    311337         CASE DEFAULT                                !    ORCA R05 or R025 
  • branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

    r5123 r5315  
    105105      REAL(wp) ::   zlam1, zcos_alpha, zim1 , zjm1 , ze1, ze1deg 
    106106      REAL(wp) ::   zphi1, zsin_alpha, zim05, zjm05 
     107      INTEGER  ::   isrow                ! index for ORCA1 starting row 
     108 
    107109      !!---------------------------------------------------------------------- 
    108110      ! 
     
    159161         IF( cp_cfg == "orca" .AND. jp_cfg == 1 ) THEN    ! ORCA R1 configuration 
    160162            !                                             ! ===================== 
    161  
    162             ii0 = 281   ;   ii1 = 282        ! Gibraltar Strait (e2u = 20 km) 
    163             ij0 = 200   ;   ij1 = 200   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  20.e3 
     163            ! This dirty section will be suppressed by simplification process: all this will come back in input files 
     164            ! Currently these hard-wired indices relate to the original (pre-v3.6) configuration  
     165            ! which had a grid-size of 362x292. 
     166            ! This grid has been extended southwards for use with the under ice-shelf options (isf) introduced in v3.6.  
     167            ! The original domain can still be used optionally if the isf code is not activated.  
     168            ! An adjustment (isrow) is made to the hard-wired indices if the extended domain (362x332) is being used. 
     169            !  
     170            IF    ( jpjglo == 292 ) THEN  ;  isrow = 0  ! Using pre-v3.6 files or adjusted start row from isf-extended grid 
     171            ELSEIF( jpjglo == 332 ) THEN  ;  isrow = 40 ! Using full isf­extended domain.  
     172            ENDIF      
     173            ! 
     174            ii0 = 282           ;   ii1 = 283        ! Gibraltar Strait (e2u = 20 km) 
     175            ij0 = 201 + isrow   ;   ij1 = 201 + isrow   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  20.e3 
    164176            IF(lwp) WRITE(numout,*) 
    165177            IF(lwp) WRITE(numout,*) '             orca_r1: Gibraltar : e2u reduced to 20 km' 
    166178 
    167             ii0 = 314   ;   ii1 = 315        ! Bhosporus Strait (e2u = 10 km) 
    168             ij0 = 208   ;   ij1 = 208   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  10.e3 
     179            ii0 = 314           ;   ii1 = 315        ! Bhosporus Strait (e2u = 10 km) 
     180            ij0 = 208 + isrow   ;   ij1 = 208 + isrow   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  10.e3 
    169181            IF(lwp) WRITE(numout,*) 
    170182            IF(lwp) WRITE(numout,*) '             orca_r1: Bhosporus : e2u reduced to 10 km' 
    171183 
    172             ii0 =  44   ;   ii1 =  44        ! Lombok Strait (e1v = 13 km) 
    173             ij0 = 124   ;   ij1 = 125   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  13.e3 
     184            ii0 =  44           ;   ii1 =  44        ! Lombok Strait (e1v = 13 km) 
     185            ij0 = 124 + isrow   ;   ij1 = 125 + isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  13.e3 
    174186            IF(lwp) WRITE(numout,*) 
    175187            IF(lwp) WRITE(numout,*) '             orca_r1: Lombok : e1v reduced to 10 km' 
    176188 
    177             ii0 =  48   ;   ii1 =  48        ! Sumba Strait (e1v = 8 km) [closed from bathy_11 on] 
    178             ij0 = 124   ;   ij1 = 125   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  8.e3 
     189            ii0 =  48           ;   ii1 =  48        ! Sumba Strait (e1v = 8 km) [closed from bathy_11 on] 
     190            ij0 = 124 + isrow   ;   ij1 = 125 + isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  8.e3 
    179191            IF(lwp) WRITE(numout,*) 
    180192            IF(lwp) WRITE(numout,*) '             orca_r1: Sumba : e1v reduced to 8 km' 
    181193 
    182             ii0 =  53   ;   ii1 =  53        ! Ombai Strait (e1v = 13 km) 
    183             ij0 = 124   ;   ij1 = 125   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 13.e3 
     194            ii0 =  53           ;   ii1 =  53        ! Ombai Strait (e1v = 13 km) 
     195            ij0 = 124 + isrow   ;   ij1 = 125 + isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 13.e3 
    184196            IF(lwp) WRITE(numout,*) 
    185197            IF(lwp) WRITE(numout,*) '             orca_r1: Ombai : e1v reduced to 13 km' 
    186198 
    187             ii0 =  56   ;   ii1 =  56        ! Timor Passage (e1v = 20 km) 
    188             ij0 = 124   ;   ij1 = 125   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 20.e3 
     199            ii0 =  56           ;   ii1 =  56        ! Timor Passage (e1v = 20 km) 
     200            ij0 = 124 + isrow   ;   ij1 = 125 + isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 20.e3 
    189201            IF(lwp) WRITE(numout,*) 
    190202            IF(lwp) WRITE(numout,*) '             orca_r1: Timor Passage : e1v reduced to 20 km' 
    191203 
    192             ii0 =  55   ;   ii1 =  55        ! West Halmahera Strait (e1v = 30 km) 
    193             ij0 = 141   ;   ij1 = 142   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 30.e3 
     204            ii0 =  55           ;   ii1 =  55        ! West Halmahera Strait (e1v = 30 km) 
     205            ij0 = 141 + isrow   ;   ij1 = 142 + isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 30.e3 
    194206            IF(lwp) WRITE(numout,*) 
    195207            IF(lwp) WRITE(numout,*) '             orca_r1: W Halmahera : e1v reduced to 30 km' 
    196208 
    197             ii0 =  58   ;   ii1 =  58        ! East Halmahera Strait (e1v = 50 km) 
    198             ij0 = 141   ;   ij1 = 142   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 50.e3 
     209            ii0 =  58           ;   ii1 =  58        ! East Halmahera Strait (e1v = 50 km) 
     210            ij0 = 141 + isrow   ;   ij1 = 142 + isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 50.e3 
    199211            IF(lwp) WRITE(numout,*) 
    200212            IF(lwp) WRITE(numout,*) '             orca_r1: E Halmahera : e1v reduced to 50 km' 
    201  
    202             ! 
    203  
    204             ! 
    205             ! 
    206213            ! 
    207214            ! 
  • branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

    r5120 r5315  
    134134      INTEGER  ::   ijf, ijl, ij0, ij1       !   -       - 
    135135      INTEGER  ::   ios 
     136      INTEGER  ::   isrow                    ! index for ORCA1 starting row 
    136137      INTEGER , POINTER, DIMENSION(:,:) ::  imsk 
    137138      REAL(wp), POINTER, DIMENSION(:,:) ::  zwf 
     
    401402      IF( cp_cfg == "orca" .AND. jp_cfg == 1 ) THEN   ! ORCA R1 configuration 
    402403         !                                                 ! 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 
    403415         IF(lwp) WRITE(numout,*) 
    404416         IF(lwp) WRITE(numout,*) '   orca_r1: increase friction near the following straits : ' 
    405417         IF(lwp) WRITE(numout,*) '      Gibraltar ' 
    406          ii0 = 283   ;   ii1 = 284        ! Gibraltar Strait  
    407          ij0 = 200   ;   ij1 = 200   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp   
     418         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   
    408420 
    409421         IF(lwp) WRITE(numout,*) '      Bhosporus ' 
    410          ii0 = 314   ;   ii1 = 315        ! Bhosporus Strait  
    411          ij0 = 208   ;   ij1 = 208   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp   
     422         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   
    412424 
    413425         IF(lwp) WRITE(numout,*) '      Makassar (Top) ' 
    414          ii0 =  48   ;   ii1 =  48        ! Makassar Strait (Top)  
    415          ij0 = 149   ;   ij1 = 150   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3._wp   
     426         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   
    416428 
    417429         IF(lwp) WRITE(numout,*) '      Lombok ' 
    418          ii0 =  44   ;   ii1 =  44        ! Lombok Strait  
    419          ij0 = 124   ;   ij1 = 125   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp   
     430         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   
    420432 
    421433         IF(lwp) WRITE(numout,*) '      Ombai ' 
    422          ii0 =  53   ;   ii1 =  53        ! Ombai Strait  
    423          ij0 = 124   ;   ij1 = 125   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp   
     434         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   
    424436 
    425437         IF(lwp) WRITE(numout,*) '      Timor Passage ' 
    426          ii0 =  56   ;   ii1 =  56        ! Timor Passage  
    427          ij0 = 124   ;   ij1 = 125   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp   
     438         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   
    428440 
    429441         IF(lwp) WRITE(numout,*) '      West Halmahera ' 
    430          ii0 =  58   ;   ii1 =  58        ! West Halmahera Strait  
    431          ij0 = 141   ;   ij1 = 142   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3._wp   
     442         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   
    432444 
    433445         IF(lwp) WRITE(numout,*) '      East Halmahera ' 
    434          ii0 =  55   ;   ii1 =  55        ! East Halmahera Strait  
    435          ij0 = 141   ;   ij1 = 142   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3._wp   
     446         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   
    436448         ! 
    437449      ENDIF 
  • branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r5120 r5315  
    10391039      INTEGER ::   ji, jj, jk                                          ! dummy loop indices 
    10401040      INTEGER ::   ij0, ij1, ii0, ii1                                  ! dummy loop indices 
     1041      INTEGER ::   isrow                                               ! index for ORCA1 starting row 
    10411042      !! acc 
    10421043      !! Hmm with the time splitting these "fixes" seem to do more harm than good. Temporarily disabled for 
     
    11221123      IF( cp_cfg == "orca" .AND. jp_cfg == 1 ) THEN    ! ORCA R1 configuration 
    11231124         !                                             ! ===================== 
    1124          ! 
    1125          ii0 = 281   ;   ii1 = 282        ! Gibraltar Strait (e2u was modified) 
    1126          ij0 = 200   ;   ij1 = 200 
     1125         ! This dirty section will be suppressed by simplification process: all this will come back in input files 
     1126         ! Currently these hard-wired indices relate to the original (pre-v3.6) configuration  
     1127         ! which had a grid-size of 362x292. 
     1128         ! This grid has been extended southwards for use with the under ice-shelf options (isf) introduced in v3.6.  
     1129         ! The original domain can still be used optionally if the isf code is not activated.  
     1130         ! An adjustment (isrow) is made to the hard-wired indices if the extended domain (362x332) is being used. 
     1131         !  
     1132         IF    ( jpjglo == 292 ) THEN  ;  isrow = 0  ! Using pre-v3.6 files or adjusted start row from isf-extended grid 
     1133         ELSEIF( jpjglo == 332 ) THEN  ;  isrow = 40 ! Using full isf­extended domain.  
     1134         ENDIF      
     1135         ! 
     1136         ii0 = 282           ;   ii1 = 283        ! Gibraltar Strait (e2u was modified) 
     1137         ij0 = 201 + isrow   ;   ij1 = 201 + isrow 
    11271138         DO jk = 1, jpkm1 
    11281139            DO jj = mj0(ij0), mj1(ij1) 
     
    11441155         END DO 
    11451156         ! 
    1146          ii0 = 314   ;   ii1 = 315        ! Bhosporus Strait (e2u was modified) 
    1147          ij0 = 208   ;   ij1 = 208 
     1157         ii0 = 314           ;   ii1 = 315        ! Bhosporus Strait (e2u was modified) 
     1158         ij0 = 208 + isrow   ;   ij1 = 208 + isrow 
    11481159         DO jk = 1, jpkm1 
    11491160            DO jj = mj0(ij0), mj1(ij1) 
     
    11651176         END DO 
    11661177         ! 
    1167          ii0 =  44   ;   ii1 =  44        ! Lombok Strait (e1v was modified) 
    1168          ij0 = 124   ;   ij1 = 125 
     1178         ii0 =  44           ;   ii1 =  44        ! Lombok Strait (e1v was modified) 
     1179         ij0 = 124 + isrow   ;   ij1 = 125 + isrow 
    11691180         DO jk = 1, jpkm1 
    11701181            DO jj = mj0(ij0), mj1(ij1) 
     
    11811192         END DO 
    11821193         ! 
    1183          ii0 =  48   ;   ii1 =  48        ! Sumba Strait (e1v was modified) [closed from bathy_11 on] 
    1184          ij0 = 124   ;   ij1 = 125 
     1194         ii0 =  48           ;   ii1 =  48        ! Sumba Strait (e1v was modified) [closed from bathy_11 on] 
     1195         ij0 = 124 + isrow   ;   ij1 = 125 + isrow 
    11851196         DO jk = 1, jpkm1 
    11861197            DO jj = mj0(ij0), mj1(ij1) 
     
    11971208         END DO 
    11981209         ! 
    1199          ii0 =  53   ;   ii1 =  53        ! Ombai Strait (e1v was modified) 
    1200          ij0 = 124   ;   ij1 = 125 
     1210         ii0 =  53          ;   ii1 =  53        ! Ombai Strait (e1v was modified) 
     1211         ij0 = 124 + isrow  ;   ij1 = 125  + isrow   
    12011212         DO jk = 1, jpkm1 
    12021213            DO jj = mj0(ij0), mj1(ij1) 
     
    12131224         END DO 
    12141225         ! 
    1215          ii0 =  56   ;   ii1 =  56        ! Timor Passage (e1v was modified) 
    1216          ij0 = 124   ;   ij1 = 125 
     1226         ii0 =  56             ;   ii1 =  56        ! Timor Passage (e1v was modified) 
     1227         ij0 = 124  + isrow    ;   ij1 = 125  + isrow   
    12171228         DO jk = 1, jpkm1 
    12181229            DO jj = mj0(ij0), mj1(ij1) 
     
    12291240         END DO 
    12301241         ! 
    1231          ii0 =  55   ;   ii1 =  55        ! West Halmahera Strait (e1v was modified) 
    1232          ij0 = 141   ;   ij1 = 142 
     1242         ii0 =  55             ;   ii1 =  55        ! West Halmahera Strait (e1v was modified) 
     1243         ij0 = 141  + isrow    ;   ij1 = 142  + isrow   
    12331244         DO jk = 1, jpkm1 
    12341245            DO jj = mj0(ij0), mj1(ij1) 
     
    12451256         END DO 
    12461257         ! 
    1247          ii0 =  58   ;   ii1 =  58        ! East Halmahera Strait (e1v was modified) 
    1248          ij0 = 141   ;   ij1 = 142 
     1258         ii0 =  58             ;   ii1 =  58        ! East Halmahera Strait (e1v was modified) 
     1259         ij0 = 141  + isrow    ;   ij1 = 142  + isrow   
    12491260         DO jk = 1, jpkm1 
    12501261            DO jj = mj0(ij0), mj1(ij1) 
  • branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r5120 r5315  
    505505            CALL iom_close( inum ) 
    506506            mbathy(:,:) = INT( bathy(:,:) ) 
     507            ! 
     508            ! CL :  add Amazon deeper 
     509            IF( cp_cfg == "orca" .AND. jp_cfg == 1 ) THEN    ! ORCA R1 configuration  
     510                 ii0 = 230   ;   ii1 = 245                   ! Amazon area  
     511                 ij0 = 140   ;   ij1 = 155                   ! no ocean shallower than 30 meters 
     512                 DO ji = mi0(ii0), mi1(ii1) 
     513                    DO jj = mj0(ij0), mj1(ij1) 
     514                       IF( bathy(ji,jj) .LE. 30. .AND. bathy(ji,jj) .GT. 0.0 ) bathy(ji,jj) = 30._wp 
     515                    END DO 
     516                 END DO 
     517                 IF(lwp) WRITE(numout,*) 
     518                 IF(lwp) WRITE(numout,*) '             orca_r1: Amazon area not shallower than 30 meters for: ' 
     519                 IF(lwp) WRITE(numout,*) '                    Longitude index ',ii0, ii0 
     520                 IF(lwp) WRITE(numout,*) '                    Latitude index  ',ij0, ij0 
     521            ENDIF 
    507522            !                                                ! ===================== 
    508523            IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN    ! ORCA R2 configuration 
  • branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

    r5226 r5315  
    187187      INTEGER  :: ierr, ierr1, ierr2, ierr3 
    188188      INTEGER  :: ios                 ! Local integer output status for namelist read 
    189       INTEGER  :: ik50               !  last level where depth less than 50 m 
     189      INTEGER  :: ik50                !  last level where depth less than 50 m 
     190      INTEGER  ::   isrow             ! index for ORCA1 starting row 
    190191      REAL(wp) :: zexpide, zdenitide, zmaskt 
    191192      REAL(wp) :: ztimes_dust, ztimes_riv, ztimes_ndep  
     
    463464            ELSE IF( jp_cfg == 1 ) THEN 
    464465               ! 
    465                ii0 = 357   ;   ii1 =  357        ! Southern Island : Kerguelen 
    466                ij0 =  75   ;   ij1 =   76  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp 
    467                ! 
    468                ii0 = 243   ;   ii1 =  243        ! South Georgia 
    469                ij0 =  57   ;   ij1 =   59  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp 
    470                ! 
    471                ii0 = 227   ;   ii1 =  227        ! Falklands 
    472                ij0 =  71   ;   ij1 =   73  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp 
    473                ! 
    474                ii0 = 326   ;   ii1 =  327        ! Crozet 
    475                ij0 =  79   ;   ij1 =   79  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp 
    476                ! 
    477                ii0 = 243   ;   ii1 =  243        ! South Orkney 
    478                ij0 =  56   ;   ij1 =   56  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp 
    479                ! 
    480                ii0 = 283   ;   ii1 =  286        ! Bouvet Island 
    481                ij0 =  66   ;   ij1 =   66  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp 
    482                ! 
    483                ii0 = 361   ;   ii1 =  361        ! Prince edwards 
    484                ij0 =  67   ;   ij1 =   68  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp 
    485                ! 
    486                ii0 =  91   ;   ii1 =   92        ! Balleny islands 
    487                ij0 =  42   ;   ij1 =   43  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp 
     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  + iswrow   ;   ij1 =  76   + iswrow   
     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  + iswrow   ;   ij1 =   59  + iswrow  
     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  + iswrow   ;   ij1 =   73  + iswrow   
     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  + iswrow   ;   ij1 =   79  + iswrow  
     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  + iswrow   ;   ij1 =   56  + iswrow   
     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  + iswrow   ;   ij1 =   66  + iswrow   
     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  + iswrow   ;   ij1 =   68  + iswrow  
     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  + iswrow   ;   ij1 =   43  + iswrow  
     507               zcmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1), 1:jpk ) = 0.3_wp 
    488508               ! 
    489509            ELSE IF( jp_cfg == 05 ) THEN 
Note: See TracChangeset for help on using the changeset viewer.