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 10170 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests – NEMO

Ignore:
Timestamp:
2018-10-03T16:49:50+02:00 (5 years ago)
Author:
smasson
Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: action 2a: add report calls of lbc_lnk, see #2133

Location:
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests
Files:
16 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/BENCH/MY_SRC/diawri.F90

    r10166 r10170  
    181181            END DO 
    182182         END DO 
    183          CALL lbc_lnk( z2d, 'T', 1. ) 
     183         CALL lbc_lnk( 'diawri', z2d, 'T', 1. ) 
    184184         CALL iom_put( "taubot", z2d )            
    185185      ENDIF 
     
    237237            END DO 
    238238         END DO 
    239          CALL lbc_lnk( z2d, 'T', 1. ) 
     239         CALL lbc_lnk( 'diawri', z2d, 'T', 1. ) 
    240240         CALL iom_put( "sstgrad2",  z2d )          ! square of module of sst gradient 
    241241         z2d(:,:) = SQRT( z2d(:,:) ) 
     
    281281            END DO 
    282282         END DO 
    283          CALL lbc_lnk( z3d, 'T', 1. ) 
     283         CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 
    284284         CALL iom_put( "eken", z3d )                 ! kinetic energy 
    285285      ENDIF 
     
    307307            END DO 
    308308         END DO 
    309          CALL lbc_lnk( z2d, 'U', -1. ) 
     309         CALL lbc_lnk( 'diawri', z2d, 'U', -1. ) 
    310310         CALL iom_put( "u_heattr", 0.5*rcp * z2d )    ! heat transport in i-direction 
    311311      ENDIF 
     
    320320            END DO 
    321321         END DO 
    322          CALL lbc_lnk( z2d, 'U', -1. ) 
     322         CALL lbc_lnk( 'diawri', z2d, 'U', -1. ) 
    323323         CALL iom_put( "u_salttr", 0.5 * z2d )        ! heat transport in i-direction 
    324324      ENDIF 
     
    342342            END DO 
    343343         END DO 
    344          CALL lbc_lnk( z2d, 'V', -1. ) 
     344         CALL lbc_lnk( 'diawri', z2d, 'V', -1. ) 
    345345         CALL iom_put( "v_heattr", 0.5*rcp * z2d )    !  heat transport in j-direction 
    346346      ENDIF 
     
    355355            END DO 
    356356         END DO 
    357          CALL lbc_lnk( z2d, 'V', -1. ) 
     357         CALL lbc_lnk( 'diawri', z2d, 'V', -1. ) 
    358358         CALL iom_put( "v_salttr", 0.5 * z2d )        !  heat transport in j-direction 
    359359      ENDIF 
     
    368368            END DO 
    369369         END DO 
    370          CALL lbc_lnk( z2d, 'T', -1. ) 
     370         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    371371         CALL iom_put( "tosmint", rau0 * z2d )        ! Vertical integral of temperature 
    372372      ENDIF 
     
    380380            END DO 
    381381         END DO 
    382          CALL lbc_lnk( z2d, 'T', -1. ) 
     382         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    383383         CALL iom_put( "somint", rau0 * z2d )         ! Vertical integral of salinity 
    384384      ENDIF 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/BENCH/MY_SRC/usrdef_istate.F90

    r10166 r10170  
    8080      ENDDO 
    8181      ! 
    82       CALL lbc_lnk(pssh, 'T',  1. )            ! apply boundary conditions 
    83       CALL lbc_lnk( pts, 'T',  1. )            ! apply boundary conditions 
    84       CALL lbc_lnk(  pu, 'U', -1. )            ! apply boundary conditions 
    85       CALL lbc_lnk(  pv, 'V', -1. )            ! apply boundary conditions 
     82      CALL lbc_lnk('usrdef_istate', pssh, 'T',  1. )            ! apply boundary conditions 
     83      CALL lbc_lnk( 'usrdef_istate', pts, 'T',  1. )            ! apply boundary conditions 
     84      CALL lbc_lnk(  'usrdef_istate', pu, 'U', -1. )            ! apply boundary conditions 
     85      CALL lbc_lnk(  'usrdef_istate', pv, 'V', -1. )            ! apply boundary conditions 
    8686       
    8787   END SUBROUTINE usr_def_istate 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/BENCH/MY_SRC/usrdef_zgr.F90

    r10166 r10170  
    194194      z2d(:,:) = REAL( jpkm1 , wp )          ! flat bottom 
    195195      ! 
    196       CALL lbc_lnk( z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) 
     196      CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) 
    197197      ! 
    198198      k_bot(:,:) = INT( z2d(:,:) )           ! =jpkm1 over the ocean point, =0 elsewhere 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/diawri.F90

    r10074 r10170  
    185185            END DO 
    186186         END DO 
    187          CALL lbc_lnk( z2d, 'T', 1. ) 
     187         CALL lbc_lnk( 'diawri', z2d, 'T', 1. ) 
    188188         CALL iom_put( "taubot", z2d )            
    189189      ENDIF 
     
    244244            END DO 
    245245         END DO 
    246          CALL lbc_lnk( z3d, 'T', 1. ) 
     246         CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 
    247247         CALL iom_put( "salgrad2",  z3d )          ! square of module of sal gradient 
    248248         z3d(:,:,:) = SQRT( z3d(:,:,:) ) 
     
    260260            END DO 
    261261         END DO 
    262          CALL lbc_lnk( z2d, 'T', 1. ) 
     262         CALL lbc_lnk( 'diawri', z2d, 'T', 1. ) 
    263263         CALL iom_put( "sstgrad2",  z2d )          ! square of module of sst gradient 
    264264         z2d(:,:) = SQRT( z2d(:,:) ) 
     
    314314            END DO 
    315315         END DO 
    316          CALL lbc_lnk( z3d, 'T', 1. ) 
     316         CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 
    317317         CALL iom_put( "eken", z3d )                 ! kinetic energy 
    318318      ENDIF 
     
    335335         END DO 
    336336          
    337          CALL lbc_lnk( z3d, 'T', 1. ) 
     337         CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 
    338338         CALL iom_put( "ke", z3d ) ! kinetic energy 
    339339 
     
    363363            END DO 
    364364         END DO 
    365          CALL lbc_lnk( z3d, 'F', 1. ) 
     365         CALL lbc_lnk( 'diawri', z3d, 'F', 1. ) 
    366366         CALL iom_put( "relvor", z3d )                  ! relative vorticity 
    367367 
     
    387387            END DO 
    388388         END DO 
    389          CALL lbc_lnk( z3d, 'F', 1. ) 
     389         CALL lbc_lnk( 'diawri', z3d, 'F', 1. ) 
    390390         CALL iom_put( "potvor", z3d )                  ! potential vorticity 
    391391 
     
    413413            END DO 
    414414         END DO 
    415          CALL lbc_lnk( z2d, 'U', -1. ) 
     415         CALL lbc_lnk( 'diawri', z2d, 'U', -1. ) 
    416416         CALL iom_put( "u_heattr", 0.5*rcp * z2d )    ! heat transport in i-direction 
    417417      ENDIF 
     
    426426            END DO 
    427427         END DO 
    428          CALL lbc_lnk( z2d, 'U', -1. ) 
     428         CALL lbc_lnk( 'diawri', z2d, 'U', -1. ) 
    429429         CALL iom_put( "u_salttr", 0.5 * z2d )        ! heat transport in i-direction 
    430430      ENDIF 
     
    448448            END DO 
    449449         END DO 
    450          CALL lbc_lnk( z2d, 'V', -1. ) 
     450         CALL lbc_lnk( 'diawri', z2d, 'V', -1. ) 
    451451         CALL iom_put( "v_heattr", 0.5*rcp * z2d )    !  heat transport in j-direction 
    452452      ENDIF 
     
    461461            END DO 
    462462         END DO 
    463          CALL lbc_lnk( z2d, 'V', -1. ) 
     463         CALL lbc_lnk( 'diawri', z2d, 'V', -1. ) 
    464464         CALL iom_put( "v_salttr", 0.5 * z2d )        !  heat transport in j-direction 
    465465      ENDIF 
     
    474474            END DO 
    475475         END DO 
    476          CALL lbc_lnk( z2d, 'T', -1. ) 
     476         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    477477         CALL iom_put( "tosmint", rau0 * z2d )        ! Vertical integral of temperature 
    478478      ENDIF 
     
    486486            END DO 
    487487         END DO 
    488          CALL lbc_lnk( z2d, 'T', -1. ) 
     488         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    489489         CALL iom_put( "somint", rau0 * z2d )         ! Vertical integral of salinity 
    490490      ENDIF 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/domvvl.F90

    r10074 r10170  
    385385         !                       ! d - thickness diffusion transport: boundary conditions 
    386386         !                             (stored for tracer advction and continuity equation) 
    387          CALL lbc_lnk_multi( un_td , 'U' , -1._wp, vn_td , 'V' , -1._wp) 
     387         CALL lbc_lnk_multi( 'domvvl', un_td , 'U' , -1._wp, vn_td , 'V' , -1._wp) 
    388388 
    389389         ! 4 - Time stepping of baroclinic scale factors 
     
    396396            z2dt = 2.0_wp * rdt 
    397397         ENDIF 
    398          CALL lbc_lnk( tilde_e3t_a(:,:,:), 'T', 1._wp ) 
     398         CALL lbc_lnk( 'domvvl', tilde_e3t_a(:,:,:), 'T', 1._wp ) 
    399399         tilde_e3t_a(:,:,:) = tilde_e3t_b(:,:,:) + z2dt * tmask(:,:,:) * tilde_e3t_a(:,:,:) 
    400400 
     
    688688            END DO 
    689689         END DO 
    690          CALL lbc_lnk( pe3_out(:,:,:), 'U', 1._wp ) 
     690         CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'U', 1._wp ) 
    691691         pe3_out(:,:,:) = pe3_out(:,:,:) + e3u_0(:,:,:) 
    692692         ! 
     
    701701            END DO 
    702702         END DO 
    703          CALL lbc_lnk( pe3_out(:,:,:), 'V', 1._wp ) 
     703         CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'V', 1._wp ) 
    704704         pe3_out(:,:,:) = pe3_out(:,:,:) + e3v_0(:,:,:) 
    705705         ! 
     
    715715            END DO 
    716716         END DO 
    717          CALL lbc_lnk( pe3_out(:,:,:), 'F', 1._wp ) 
     717         CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'F', 1._wp ) 
    718718         pe3_out(:,:,:) = pe3_out(:,:,:) + e3f_0(:,:,:) 
    719719         ! 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/trazdf.F90

    r10074 r10170  
    9292         END DO 
    9393!!gm this should be moved in trdtra.F90 and done on all trends 
    94          CALL lbc_lnk_multi( ztrdt, 'T', 1. , ztrds, 'T', 1. ) 
     94         CALL lbc_lnk_multi( 'trazdf', ztrdt, 'T', 1. , ztrds, 'T', 1. ) 
    9595!!gm 
    9696         CALL trd_tra( kt, 'TRA', jp_tem, jptra_zdf, ztrdt ) 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/usrdef_istate.F90

    r10074 r10170  
    169169               END DO 
    170170            END DO 
    171             CALL lbc_lnk( pssh, 'T',  1. ) 
     171            CALL lbc_lnk( 'usrdef_istate', pssh, 'T',  1. ) 
    172172         END DO 
    173173          
     
    293293         pssh(:,:) = pssh(:,:) + ( 0.1  * zrandom(:,:) - 0.05 ) 
    294294      END IF 
    295       CALL lbc_lnk( pssh, 'T',  1. ) 
    296       CALL lbc_lnk(  pts, 'T',  1. ) 
    297       CALL lbc_lnk(   pu, 'U', -1. ) 
    298       CALL lbc_lnk(   pv, 'V', -1. ) 
     295      CALL lbc_lnk( 'usrdef_istate', pssh, 'T',  1. ) 
     296      CALL lbc_lnk(  'usrdef_istate', pts, 'T',  1. ) 
     297      CALL lbc_lnk(   'usrdef_istate', pu, 'U', -1. ) 
     298      CALL lbc_lnk(   'usrdef_istate', pv, 'V', -1. ) 
    299299 
    300300   END SUBROUTINE usr_def_istate 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/usrdef_zgr.F90

    r10074 r10170  
    204204      END SELECT 
    205205      ! 
    206       CALL lbc_lnk( z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) 
     206      CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) 
    207207      ! 
    208208      k_bot(:,:) = INT( z2d(:,:) )           ! =jpkm1 over the ocean point, =0 elsewhere 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/ISOMIP/MY_SRC/usrdef_zgr.F90

    r10074 r10170  
    8989      ! the ocean basin surrounded by land (1 grid-point) set through lbc_lnk call as jperio=0  
    9090      z2d(:,:) = 1._wp                    ! surface ocean is the 1st level 
    91       CALL lbc_lnk( z2d, 'T', 1. )        ! closed basin since jperio = 0 (see userdef_nam.F90) 
     91      CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )        ! closed basin since jperio = 0 (see userdef_nam.F90) 
    9292      zmsk(:,:) = NINT( z2d(:,:) ) 
    9393      ! 
     
    177177            END DO 
    178178         END DO 
    179          CALL lbc_lnk( pe3v , 'V', 1._wp )   ;   CALL lbc_lnk( pe3vw, 'V', 1._wp ) 
    180          CALL lbc_lnk( pe3f , 'F', 1._wp ) 
     179         CALL lbc_lnk( 'usrdef_zgr', pe3v , 'V', 1._wp )   ;   CALL lbc_lnk( 'usrdef_zgr', pe3vw, 'V', 1._wp ) 
     180         CALL lbc_lnk( 'usrdef_zgr', pe3f , 'F', 1._wp ) 
    181181         DO jk = 1, jpk 
    182182            ! set to z-scale factor if zero (i.e. along closed boundaries) because of lbclnk 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/LOCK_EXCHANGE/MY_SRC/usrdef_zgr.F90

    r10074 r10170  
    8888      ! the ocean basin surrounded by land (1 grid-point) set through lbc_lnk call as jperio=0  
    8989      z2d(:,:) = 1._wp                    ! surface ocean is the 1st level 
    90       CALL lbc_lnk( z2d, 'T', 1. )        ! closed basin since jperio = 0 (see userdef_nam.F90) 
     90      CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )        ! closed basin since jperio = 0 (see userdef_nam.F90) 
    9191      k_top(:,:) = NINT( z2d(:,:) ) 
    9292      ! 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/OVERFLOW/MY_SRC/usrdef_zgr.F90

    r10074 r10170  
    9393         zhu(ji,:) = 0.5_wp * ( zht(ji,:) + zht(ji+1,:) ) 
    9494      END DO 
    95       CALL lbc_lnk( zhu, 'U', 1. )     ! boundary condition: this mask the surrouding grid-points 
     95      CALL lbc_lnk( 'usrdef_zgr', zhu, 'U', 1. )     ! boundary condition: this mask the surrouding grid-points 
    9696      !                                ! ==>>>  set by hand non-zero value on first/last columns & rows  
    9797      DO ji = mi0(1), mi1(1)              ! first row of global domain only 
     
    112112      ! the ocean basin surrounded by land (1 grid-point) set through lbc_lnk call as jperio=0  
    113113      z2d(:,:) = 1._wp                    ! surface ocean is the 1st level 
    114       CALL lbc_lnk( z2d, 'T', 1. )        ! closed basin since jperio = 0 (see userdef_nam.F90) 
     114      CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )        ! closed basin since jperio = 0 (see userdef_nam.F90) 
    115115      k_top(:,:) = NINT( z2d(:,:) ) 
    116116      ! 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/VORTEX/MY_SRC/domvvl.F90

    r10074 r10170  
    408408         !                       ! d - thickness diffusion transport: boundary conditions 
    409409         !                             (stored for tracer advction and continuity equation) 
    410          CALL lbc_lnk_multi( un_td , 'U' , -1._wp, vn_td , 'V' , -1._wp) 
     410         CALL lbc_lnk_multi( 'domvvl', un_td , 'U' , -1._wp, vn_td , 'V' , -1._wp) 
    411411 
    412412         ! 4 - Time stepping of baroclinic scale factors 
     
    419419            z2dt = 2.0_wp * rdt 
    420420         ENDIF 
    421          CALL lbc_lnk( tilde_e3t_a(:,:,:), 'T', 1._wp ) 
     421         CALL lbc_lnk( 'domvvl', tilde_e3t_a(:,:,:), 'T', 1._wp ) 
    422422         tilde_e3t_a(:,:,:) = tilde_e3t_b(:,:,:) + z2dt * tmask(:,:,:) * tilde_e3t_a(:,:,:) 
    423423 
     
    711711            END DO 
    712712         END DO 
    713          CALL lbc_lnk( pe3_out(:,:,:), 'U', 1._wp ) 
     713         CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'U', 1._wp ) 
    714714         pe3_out(:,:,:) = pe3_out(:,:,:) + e3u_0(:,:,:) 
    715715         ! 
     
    724724            END DO 
    725725         END DO 
    726          CALL lbc_lnk( pe3_out(:,:,:), 'V', 1._wp ) 
     726         CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'V', 1._wp ) 
    727727         pe3_out(:,:,:) = pe3_out(:,:,:) + e3v_0(:,:,:) 
    728728         ! 
     
    738738            END DO 
    739739         END DO 
    740          CALL lbc_lnk( pe3_out(:,:,:), 'F', 1._wp ) 
     740         CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'F', 1._wp ) 
    741741         pe3_out(:,:,:) = pe3_out(:,:,:) + e3f_0(:,:,:) 
    742742         ! 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/VORTEX/MY_SRC/usrdef_istate.F90

    r10074 r10170  
    136136      END DO 
    137137 
    138       CALL lbc_lnk( pu, 'U', -1. ) 
    139       CALL lbc_lnk( pv, 'V', -1. ) 
     138      CALL lbc_lnk( 'usrdef_istate', pu, 'U', -1. ) 
     139      CALL lbc_lnk( 'usrdef_istate', pv, 'V', -1. ) 
    140140      !    
    141141   END SUBROUTINE usr_def_istate 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/VORTEX/MY_SRC/usrdef_zgr.F90

    r10074 r10170  
    192192      z2d(:,:) = REAL( jpkm1 , wp )          ! flat bottom 
    193193      ! 
    194       CALL lbc_lnk( z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) 
     194      CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) 
    195195      ! 
    196196      k_bot(:,:) = INT( z2d(:,:) )           ! =jpkm1 over the ocean point, =0 elsewhere 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/WAD/MY_SRC/bdyini.F90

    r10074 r10170  
    11331133         END DO 
    11341134      END DO 
    1135       CALL lbc_lnk_multi( bdyumask, 'U', 1. , bdyvmask, 'V', 1. )   ! Lateral boundary cond.  
     1135      CALL lbc_lnk_multi( 'bdyini', bdyumask, 'U', 1. , bdyvmask, 'V', 1. )   ! Lateral boundary cond.  
    11361136 
    11371137      ! bdy masks are now set to zero on boundary points: 
     
    11691169 
    11701170      ! Lateral boundary conditions 
    1171       CALL lbc_lnk( zfmask, 'F', 1. )  
    1172       CALL lbc_lnk_multi( bdyumask, 'U', 1. , bdyvmask, 'V', 1., bdytmask, 'T', 1. ) 
     1171      CALL lbc_lnk( 'bdyini', zfmask, 'F', 1. )  
     1172      CALL lbc_lnk_multi( 'bdyini', bdyumask, 'U', 1. , bdyvmask, 'V', 1., bdytmask, 'T', 1. ) 
    11731173      DO ib_bdy = 1, nb_bdy       ! Indices and directions of rim velocity components 
    11741174 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/WAD/MY_SRC/usrdef_zgr.F90

    r10074 r10170  
    234234         zhu(ji,:) = 0.5_wp * ( zht(ji,:) + zht(ji+1,:) ) 
    235235      END DO 
    236       CALL lbc_lnk( zhu, 'U', 1. )     ! boundary condition: this mask the surrounding grid-points 
     236      CALL lbc_lnk( 'usrdef_zgr', zhu, 'U', 1. )     ! boundary condition: this mask the surrounding grid-points 
    237237      !                                ! ==>>>  set by hand non-zero value on first/last columns & rows  
    238238      DO ji = mi0(1), mi1(1)              ! first row of global domain only 
     
    247247         zhv(:,jj) = 0.5_wp * ( zht(:,jj) + zht(:,jj+1) ) 
    248248      END DO 
    249       CALL lbc_lnk( zhv, 'V', 1. )     ! boundary condition: this mask the surrounding grid-points 
     249      CALL lbc_lnk( 'usrdef_zgr', zhv, 'V', 1. )     ! boundary condition: this mask the surrounding grid-points 
    250250      DO jj = mj0(1), mj1(1)   ! first  row of global domain only 
    251251         zhv(:,jj) = zht(:,jj) 
     
    272272 
    273273 
    274       CALL lbc_lnk( z2d, 'T', 1. )        ! closed basin since jperio = 0 (see userdef_nam.F90) 
     274      CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )        ! closed basin since jperio = 0 (see userdef_nam.F90) 
    275275      k_top(:,:) = NINT( z2d(:,:) ) 
    276276      ! 
     
    311311           END DO       
    312312         END DO       
    313          CALL lbc_lnk( pdept, 'T', 1. ) 
    314          CALL lbc_lnk( pdepw, 'T', 1. ) 
    315          CALL lbc_lnk( pe3t , 'T', 1. ) 
    316          CALL lbc_lnk( pe3w , 'T', 1. ) 
    317          CALL lbc_lnk( pe3u , 'U', 1. ) 
    318          CALL lbc_lnk( pe3uw, 'U', 1. ) 
    319          CALL lbc_lnk( pe3f , 'F', 1. ) 
    320          CALL lbc_lnk( pe3v , 'V', 1. ) 
    321          CALL lbc_lnk( pe3vw, 'V', 1. ) 
     313         CALL lbc_lnk( 'usrdef_zgr', pdept, 'T', 1. ) 
     314         CALL lbc_lnk( 'usrdef_zgr', pdepw, 'T', 1. ) 
     315         CALL lbc_lnk( 'usrdef_zgr', pe3t , 'T', 1. ) 
     316         CALL lbc_lnk( 'usrdef_zgr', pe3w , 'T', 1. ) 
     317         CALL lbc_lnk( 'usrdef_zgr', pe3u , 'U', 1. ) 
     318         CALL lbc_lnk( 'usrdef_zgr', pe3uw, 'U', 1. ) 
     319         CALL lbc_lnk( 'usrdef_zgr', pe3f , 'F', 1. ) 
     320         CALL lbc_lnk( 'usrdef_zgr', pe3v , 'V', 1. ) 
     321         CALL lbc_lnk( 'usrdef_zgr', pe3vw, 'V', 1. ) 
    322322         WHERE( pe3t (:,:,:) == 0._wp )   pe3t (:,:,:) = 1._wp 
    323323         WHERE( pe3u (:,:,:) == 0._wp )   pe3u (:,:,:) = 1._wp 
Note: See TracChangeset for help on using the changeset viewer.