New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 13138 for NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/USR – NEMO

Ignore:
Timestamp:
2020-06-22T11:13:03+02:00 (4 years ago)
Author:
smasson
Message:

Extra_Halo: minor bugfixes and cleaning, see #2366

Location:
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/USR
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/USR/usrdef_fmask.F90

    r12960 r13138  
    7676            ! 
    7777            IF(lwp) WRITE(numout,*) '      Bab el Mandeb ' 
    78             ij0 =  87 + nn_hls       ;   ij1 =  88  + nn_hls          ! Bab el Mandeb : partial slip (pfmsk=1) 
    79             ii0 = 160  + nn_hls - 1  ;   ii1 = 160 + nn_hls - 1 
     78            ij0 =  87 + nn_hls       ;   ij1 = 88  + nn_hls          ! Bab el Mandeb : partial slip (pfmsk=1) 
     79            ii0 = 160 + nn_hls - 1   ;   ii1 = 160 + nn_hls - 1 
    8080            pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) =  1._wp 
    8181            ij0 =  88 + nn_hls       ;   ij1 =  88 + nn_hls 
    … …  
    9898!!gm    ! Currently these hard-wired indices relate to configuration with extend grid (jpjglo=332) 
    9999            ! 
    100             isrow = 332 - Nj0glo - 2 
     100            isrow = 332 - (Nj0glo + 1)   ! was 332 - jpjglo -> jpjglo_old_version = Nj0glo + 1 
    101101            ! 
    102102            IF(lwp) WRITE(numout,*) 
    103103            IF(lwp) WRITE(numout,*) '   orca_r1: increase friction near the following straits : ' 
    104104            IF(lwp) WRITE(numout,*) '      Gibraltar ' 
    105             ii0 = 282           ;   ii1 = 283        ! Gibraltar Strait  
    106             ij0 = 241 - isrow   ;   ij1 = 241 - isrow   ;   pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     105            ii0 = 282 + nn_hls - 1       ;   ii1 = 283 + nn_hls - 1        ! Gibraltar Strait  
     106            ij0 = 241 + nn_hls - isrow   ;   ij1 = 241 + nn_hls - isrow 
     107            pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    107108            ! 
    108109            IF(lwp) WRITE(numout,*) '      Bhosporus ' 
    109             ii0 = 314           ;   ii1 = 315        ! Bhosporus Strait  
    110             ij0 = 248 - isrow   ;   ij1 = 248 - isrow   ;   pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     110            ii0 = 314 + nn_hls - 1       ;   ii1 = 315 + nn_hls - 1        ! Bhosporus Strait  
     111            ij0 = 248 + nn_hls - isrow   ;   ij1 = 248 + nn_hls - isrow 
     112            pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    111113            ! 
    112114            IF(lwp) WRITE(numout,*) '      Makassar (Top) ' 
    113             ii0 =  48           ;   ii1 =  48        ! Makassar Strait (Top)  
    114             ij0 = 189 - isrow   ;   ij1 = 190 - isrow   ;   pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
     115            ii0 =  48 + nn_hls - 1       ;   ii1 =  48 + nn_hls - 1        ! Makassar Strait (Top)  
     116            ij0 = 189 + nn_hls - isrow   ;   ij1 = 190 + nn_hls - isrow 
     117            pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
    115118            ! 
    116119            IF(lwp) WRITE(numout,*) '      Lombok ' 
    117             ii0 =  44           ;   ii1 =  44        ! Lombok Strait  
    118             ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     120            ii0 =  44 + nn_hls - 1       ;   ii1 =  44 + nn_hls - 1        ! Lombok Strait  
     121            ij0 = 164 + nn_hls - isrow   ;   ij1 = 165 + nn_hls - isrow 
     122            pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    119123            ! 
    120124            IF(lwp) WRITE(numout,*) '      Ombai ' 
    121             ii0 =  53           ;   ii1 =  53        ! Ombai Strait  
    122             ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     125            ii0 =  53 + nn_hls - 1       ;   ii1 =  53 + nn_hls - 1        ! Ombai Strait  
     126            ij0 = 164 + nn_hls - isrow   ;   ij1 = 165 + nn_hls - isrow 
     127            pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    123128            ! 
    124129            IF(lwp) WRITE(numout,*) '      Timor Passage ' 
    125             ii0 =  56           ;   ii1 =  56        ! Timor Passage  
    126             ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     130            ii0 =  56 + nn_hls - 1       ;   ii1 =  56 + nn_hls - 1        ! Timor Passage  
     131            ij0 = 164 + nn_hls - isrow   ;   ij1 = 165 + nn_hls - isrow 
     132            pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    127133            ! 
    128134            IF(lwp) WRITE(numout,*) '      West Halmahera ' 
    129             ii0 =  58           ;   ii1 =  58        ! West Halmahera Strait  
    130             ij0 = 181 - isrow   ;   ij1 = 182 - isrow   ;   pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
     135            ii0 =  58 + nn_hls - 1       ;   ii1 =  58 + nn_hls - 1        ! West Halmahera Strait  
     136            ij0 = 181 + nn_hls - isrow   ;   ij1 = 182 + nn_hls - isrow 
     137            pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
    131138            ! 
    132139            IF(lwp) WRITE(numout,*) '      East Halmahera ' 
    133             ii0 =  55           ;   ii1 =  55        ! East Halmahera Strait  
    134             ij0 = 181 - isrow   ;   ij1 = 182 - isrow   ;   pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
     140            ii0 =  55 + nn_hls - 1       ;   ii1 =  55 + nn_hls - 1        ! East Halmahera Strait  
     141            ij0 = 181 + nn_hls - isrow   ;   ij1 = 182 + nn_hls - isrow 
     142            pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
    135143            ! 
    136144         CASE DEFAULT 

  • NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/USR/usrdef_zgr.F90

    r12939 r13138  
    199199      ! 
    200200      z2d(:,:) = REAL( jpkm1 , wp )                              ! flat bottom 
    201       z2d(mi0(  1 ):mi1(jpiglo),mj0(Njs0):mj1( Njs0 )) = 0._wp   ! line   number Njs0 at 0 
    202       z2d(mi0(Nis0):mi1( Nis0 ),mj0(  1 ):mj1(jpjglo)) = 0._wp   ! column number Nis0 at 0 
    203       ! 
    204       CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) 
    205201      ! 
    206202      k_bot(:,:) = NINT( z2d(:,:) )          ! =jpkm1 over the ocean point, =0 elsewhere 
Note: See TracChangeset for help on using the changeset viewer.