Changeset 13552


Ignore:
Timestamp:
2020-10-01T12:09:56+02:00 (4 months ago)
Author:
hadcv
Message:

#2365: Remove diagnostic output lbc_lnk calls

Location:
NEMO/branches/2020/dev_r13383_HPC-02_Daley_Tiling/src/OCE
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13383_HPC-02_Daley_Tiling/src/OCE/DIA/diaar5.F90

    r13519 r13552  
    308308   END SUBROUTINE dia_ar5 
    309309 
    310    ! TEMP: These changes and lbc_lnk not necessary if using XIOS (subdomain support, will not output haloes) 
     310   ! TEMP: These changes not necessary if using XIOS (subdomain support, will not output haloes) 
    311311   SUBROUTINE dia_ar5_hst( ktra, cptr, puflx, pvflx ) 
    312312      !!---------------------------------------------------------------------- 
     
    348348      IF( ntile == 0 .OR. ntile == nijtile ) THEN 
    349349         IF( cptr == 'adv' ) THEN 
    350             CALL lbc_lnk( 'diaar5', hstr_adv(:,:,ktra,1), 'U', -1.0_wp ) 
    351350            IF( ktra == jp_tem ) CALL iom_put( 'uadv_heattr' , rho0_rcp * hstr_adv(:,:,ktra,1) )  ! advective heat transport in i-direction 
    352351            IF( ktra == jp_sal ) CALL iom_put( 'uadv_salttr' , rho0     * hstr_adv(:,:,ktra,1) )  ! advective salt transport in i-direction 
    353             CALL lbc_lnk( 'diaar5', hstr_adv(:,:,ktra,2), 'V', -1.0_wp ) 
    354352            IF( ktra == jp_tem ) CALL iom_put( 'vadv_heattr' , rho0_rcp * hstr_adv(:,:,ktra,2) )  ! advective heat transport in j-direction 
    355353            IF( ktra == jp_sal ) CALL iom_put( 'vadv_salttr' , rho0     * hstr_adv(:,:,ktra,2) )  ! advective salt transport in j-direction 
    356354         ENDIF 
    357355         IF( cptr == 'ldf' ) THEN 
    358             CALL lbc_lnk( 'diaar5', hstr_ldf(:,:,ktra,1), 'U', -1.0_wp ) 
    359356            IF( ktra == jp_tem ) CALL iom_put( 'udiff_heattr' , rho0_rcp * hstr_ldf(:,:,ktra,1) ) ! diffusive heat transport in i-direction 
    360357            IF( ktra == jp_sal ) CALL iom_put( 'udiff_salttr' , rho0     * hstr_ldf(:,:,ktra,1) ) ! diffusive salt transport in i-direction 
    361             CALL lbc_lnk( 'diaar5', hstr_ldf(:,:,ktra,2), 'V', -1.0_wp ) 
    362358            IF( ktra == jp_tem ) CALL iom_put( 'vdiff_heattr' , rho0_rcp * hstr_ldf(:,:,ktra,2) ) ! diffusive heat transport in j-direction 
    363359            IF( ktra == jp_sal ) CALL iom_put( 'vdiff_salttr' , rho0     * hstr_ldf(:,:,ktra,2) ) ! diffusive salt transport in j-direction 
  • NEMO/branches/2020/dev_r13383_HPC-02_Daley_Tiling/src/OCE/DIA/diaptr.F90

    r13519 r13552  
    703703            p_fval(ji,jj) = p_fval(ji,jj-1) + pva(ji,jj) * tmask_i(ji,jj) 
    704704         END_2D 
    705          CALL lbc_lnk( 'diaptr', p_fval, 'U', -1.0_wp ) 
    706705      END DO 
    707706      !  
  • NEMO/branches/2020/dev_r13383_HPC-02_Daley_Tiling/src/OCE/LDF/ldftra.F90

    r13515 r13552  
    796796!!gm     to be redesigned....    
    797797      !                                                  !==  eiv stream function: output  ==! 
    798       CALL lbc_lnk_multi( 'ldftra', psi_uw, 'U', -1.0_wp , psi_vw, 'V', -1.0_wp ) 
    799       ! 
    800798!!gm      CALL iom_put( "psi_eiv_uw", psi_uw )                 ! output 
    801799!!gm      CALL iom_put( "psi_eiv_vw", psi_vw ) 
     
    805803      zw3d(:,:,jpk) = 0._wp                                    ! bottom value always 0 
    806804      ! 
    807       DO_3D( 1, 1, 1, 1, 1, jpkm1 )                                  ! e2u e3u u_eiv = -dk[psi_uw] 
     805      DO_3D( 0, 0, 0, 0, 1, jpkm1 )                                  ! e2u e3u u_eiv = -dk[psi_uw] 
    808806         zw3d(ji,jj,jk) = ( psi_uw(ji,jj,jk+1) - psi_uw(ji,jj,jk) ) / ( e2u(ji,jj) * e3u(ji,jj,jk,Kmm) ) 
    809807      END_3D 
    810808      CALL iom_put( "uoce_eiv", zw3d ) 
    811809      ! 
    812       DO_3D( 1, 1, 1, 1, 1, jpkm1 )                                  ! e1v e3v v_eiv = -dk[psi_vw] 
     810      DO_3D( 0, 0, 0, 0, 1, jpkm1 )                                  ! e1v e3v v_eiv = -dk[psi_vw] 
    813811         zw3d(ji,jj,jk) = ( psi_vw(ji,jj,jk+1) - psi_vw(ji,jj,jk) ) / ( e1v(ji,jj) * e3v(ji,jj,jk,Kmm) ) 
    814812      END_3D 
     
    819817            &              + psi_uw(ji,jj,jk) - psi_uw(ji-1,jj  ,jk)  ) / e1e2t(ji,jj) 
    820818      END_3D 
    821       CALL lbc_lnk( 'ldftra', zw3d, 'T', 1.0_wp )      ! lateral boundary condition 
    822819      CALL iom_put( "woce_eiv", zw3d ) 
    823820      ! 
    824821      IF( iom_use('weiv_masstr') ) THEN   ! vertical mass transport & its square value 
    825          DO_2D( 1, 1, 1, 1 ) 
     822         DO_2D( 0, 0, 0, 0 ) 
    826823            zw2d(ji,jj) = rho0 * e1e2t(ji,jj) 
    827824         END_2D 
     
    849846           zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 
    850847        END_3D 
    851         CALL lbc_lnk( 'ldftra', zw2d, 'U', -1.0_wp ) 
    852         CALL lbc_lnk( 'ldftra', zw3d, 'U', -1.0_wp ) 
    853848        CALL iom_put( "ueiv_heattr"  , zztmp * zw2d )                  ! heat transport in i-direction 
    854849        CALL iom_put( "ueiv_heattr3d", zztmp * zw3d )                  ! heat transport in i-direction 
     
    870865         zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 
    871866      END_3D 
    872       CALL lbc_lnk( 'ldftra', zw2d, 'V', -1.0_wp ) 
    873       CALL lbc_lnk( 'ldftra', zw3d, 'V', -1. ) 
    874867      CALL iom_put( "veiv_heattr"  , zztmp * zw2d )                  !  heat transport in j-direction 
    875868      CALL iom_put( "veiv_heattr3d", zztmp * zw3d )                  !  heat transport in j-direction 
     
    886879           zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 
    887880        END_3D 
    888         CALL lbc_lnk( 'ldftra', zw2d, 'U', -1.0_wp ) 
    889         CALL lbc_lnk( 'ldftra', zw3d, 'U', -1.0_wp ) 
    890881        CALL iom_put( "ueiv_salttr", zztmp * zw2d )                  ! salt transport in i-direction 
    891882        CALL iom_put( "ueiv_salttr3d", zztmp * zw3d )                ! salt transport in i-direction 
     
    898889         zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 
    899890      END_3D 
    900       CALL lbc_lnk( 'ldftra', zw2d, 'V', -1.0_wp ) 
    901       CALL lbc_lnk( 'ldftra', zw3d, 'V', -1. ) 
    902891      CALL iom_put( "veiv_salttr"  , zztmp * zw2d )                  !  salt transport in j-direction 
    903892      CALL iom_put( "veiv_salttr3d", zztmp * zw3d )                  !  salt transport in j-direction 
  • NEMO/branches/2020/dev_r13383_HPC-02_Daley_Tiling/src/OCE/TRA/trabbl.F90

    r13551 r13552  
    128128         IF( ntile == 0 .OR. ntile == nijtile ) THEN                       ! Do only on the last tile 
    129129            ! lateral boundary conditions ; just need for outputs 
    130             CALL lbc_lnk_multi( 'trabbl', ahu_bbl, 'U', 1.0_wp , ahv_bbl, 'V', 1.0_wp ) 
    131130            CALL iom_put( "ahu_bbl", ahu_bbl )   ! bbl diffusive flux i-coef 
    132131            CALL iom_put( "ahv_bbl", ahv_bbl )   ! bbl diffusive flux j-coef 
     
    143142         IF( ntile == 0 .OR. ntile == nijtile ) THEN                       ! Do only on the last tile 
    144143            ! lateral boundary conditions ; just need for outputs 
    145             CALL lbc_lnk_multi( 'trabbl', utr_bbl, 'U', 1.0_wp , vtr_bbl, 'V', 1.0_wp ) 
    146144            CALL iom_put( "uoce_bbl", utr_bbl )  ! bbl i-transport 
    147145            CALL iom_put( "voce_bbl", vtr_bbl )  ! bbl j-transport 
  • NEMO/branches/2020/dev_r13383_HPC-02_Daley_Tiling/src/OCE/TRA/traldf_triad.F90

    r13539 r13552  
    1313   USE oce            ! ocean dynamics and active tracers 
    1414   USE dom_oce        ! ocean space and time domain 
    15    ! TEMP: This change not necessary if lbc_lnk is removed from ldf_eiv_dia and XIOS has subdomain support 
     15   ! TEMP: This change not necessary if XIOS has subdomain support 
    1616   USE domain, ONLY : dom_tile 
    1717   USE domutl, ONLY : is_tile 
     
    122122      REAL(wp), DIMENSION(ST_2D(nn_hls)        ) ::   z2d                            ! 2D workspace 
    123123      REAL(wp), DIMENSION(ST_2D(nn_hls)    ,jpk) ::   zdit, zdjt, zftu, zftv, ztfw   ! 3D     - 
    124       ! TEMP: This can be ST_2D(nn_hls) if lbc_lnk is removed from ldf_eiv_dia and XIOS has subdomain support 
     124      ! TEMP: This can be ST_2D(nn_hls) if XIOS has subdomain support 
    125125      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zpsi_uw, zpsi_vw 
    126126      !!---------------------------------------------------------------------- 
     
    215215         ENDIF 
    216216         ! 
    217          ! TEMP: These changes not necessary if lbc_lnk is removed from ldf_eiv_dia and XIOS has subdomain support 
     217         ! TEMP: These changes not necessary if XIOS has subdomain support 
    218218         IF( ntile == 0 .OR. ntile == nijtile )  THEN                ! Do only for the full domain 
    219219            IF( ln_ldfeiv_dia .AND. cdtype == 'TRA' ) THEN 
  • NEMO/branches/2020/dev_r13383_HPC-02_Daley_Tiling/src/OCE/TRA/tramle.F90

    r13516 r13552  
    217217         IF( kt == nit000 .AND. (ntile == 0 .OR. ntile == 1) )  THEN             ! Do only on the first tile and timestep 
    218218            ALLOCATE( zLf_NH(jpi,jpj), zpsiu_mle(jpi,jpj,jpk), zpsiv_mle(jpi,jpj,jpk) ) 
     219            zpsiu_mle(:,:,:) = 0._wp ; zpsiv_mle(:,:,:) = 0._wp 
    219220         ENDIF 
    220221         ! 
    221          DO_2D( 1, 1, 1, 1 ) 
     222         DO_2D( 0, 0, 0, 0 ) 
    222223            zLf_NH(ji,jj) = SQRT( rb_c * zmld(ji,jj) ) * r1_ft(ji,jj)      ! Lf = N H / f 
    223224         END_2D 
    224225         ! 
    225226         ! divide by cross distance to give streamfunction with dimensions m^2/s 
    226          DO_3D( 1, 1, 1, 1, 1, ikmax+1 ) 
     227         DO_3D( 0, 0, 0, 0, 1, ikmax+1 ) 
    227228            zpsiu_mle(ji,jj,jk) = zpsi_uw(ji,jj,jk) * r1_e2u(ji,jj) 
    228229            zpsiv_mle(ji,jj,jk) = zpsi_vw(ji,jj,jk) * r1_e1v(ji,jj) 
Note: See TracChangeset for help on using the changeset viewer.