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 13540 for NEMO/branches/2020/r12377_ticket2386/src/OCE/DIA/diaptr.F90 – NEMO

Ignore:
Timestamp:
2020-09-29T12:41:06+02:00 (4 years ago)
Author:
andmirek
Message:

Ticket #2386: update to latest trunk

Location:
NEMO/branches/2020/r12377_ticket2386
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r12377_ticket2386

    • Property svn:externals
      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
        88 
        99# SETTE 
        10 ^/utils/CI/sette@HEAD         sette 
         10^/utils/CI/sette@13507        sette 
  • NEMO/branches/2020/r12377_ticket2386/src/OCE/DIA/diaptr.F90

    r12511 r13540  
    6060 
    6161   LOGICAL ::   ll_init = .TRUE.        !: tracers  trend flag (set from namelist in trdini) 
     62    
    6263   !! * Substitutions 
    6364#  include "do_loop_substitute.h90" 
     65#  include "domzgr_substitute.h90" 
    6466   !!---------------------------------------------------------------------- 
    6567   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
     
    117119            zmask(:,:,:) = 0._wp 
    118120            zts(:,:,:,:) = 0._wp 
    119             DO_3D_10_11( 1, jpkm1 ) 
     121            DO_3D( 1, 0, 1, 1, 1, jpkm1 ) 
    120122               zvfc = e1v(ji,jj) * e3v(ji,jj,jk,Kmm) 
    121123               zmask(ji,jj,jk)      = vmask(ji,jj,jk)      * zvfc 
     
    188190         zts(:,:,:,:) = 0._wp 
    189191         IF( iom_use( 'zotem' ) .OR. iom_use( 'zosal' ) .OR. iom_use( 'zosrf' )  ) THEN    ! i-mean i-k-surface  
    190             DO_3D_11_11( 1, jpkm1 ) 
     192            DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    191193               zsfc = e1t(ji,jj) * e3t(ji,jj,jk,Kmm) 
    192194               zmask(ji,jj,jk)      = tmask(ji,jj,jk)      * zsfc 
     
    278280         IF( iom_use( 'sopstvtr' ) .OR. iom_use( 'sophtvtr' ) ) THEN 
    279281            zts(:,:,:,:) = 0._wp 
    280             DO_3D_10_11( 1, jpkm1 ) 
     282            DO_3D( 1, 0, 1, 1, 1, jpkm1 ) 
    281283               zvfc = e1v(ji,jj) * e3v(ji,jj,jk,Kmm) 
    282284               zts(ji,jj,jk,jp_tem) = (ts(ji,jj,jk,jp_tem,Kmm)+ts(ji,jj+1,jk,jp_tem,Kmm)) * 0.5 * zvfc  !Tracers averaged onto V grid 
     
    353355         btmsk(:,:,1) = tmask_i(:,:)                  
    354356         CALL iom_open( 'subbasins', inum,  ldstop = .FALSE.  ) 
    355          CALL iom_get( inum, jpdom_data, 'atlmsk', btmsk(:,:,2) )   ! Atlantic basin 
    356          CALL iom_get( inum, jpdom_data, 'pacmsk', btmsk(:,:,3) )   ! Pacific  basin 
    357          CALL iom_get( inum, jpdom_data, 'indmsk', btmsk(:,:,4) )   ! Indian   basin 
     357         CALL iom_get( inum, jpdom_global, 'atlmsk', btmsk(:,:,2) )   ! Atlantic basin 
     358         CALL iom_get( inum, jpdom_global, 'pacmsk', btmsk(:,:,3) )   ! Pacific  basin 
     359         CALL iom_get( inum, jpdom_global, 'indmsk', btmsk(:,:,4) )   ! Indian   basin 
    358360         CALL iom_close( inum ) 
    359361         btmsk(:,:,5) = MAX ( btmsk(:,:,3), btmsk(:,:,4) )          ! Indo-Pacific basin 
     
    503505      ijpj = jpj 
    504506      p_fval(:) = 0._wp 
    505       DO_3D_00_00( 1, jpkm1 ) 
     507      DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    506508         p_fval(jj) = p_fval(jj) + pvflx(ji,jj,jk) * pmsk(ji,jj) * tmask_i(ji,jj) 
    507509      END_3D 
     
    536538      ijpj = jpj 
    537539      p_fval(:) = 0._wp 
    538       DO_2D_00_00 
     540      DO_2D( 0, 0, 0, 0 ) 
    539541         p_fval(jj) = p_fval(jj) + pvflx(ji,jj) * pmsk(ji,jj) * tmask_i(ji,jj) 
    540542      END_2D 
     
    565567      p_fval(:,:) = 0._wp 
    566568      DO jc = 1, jpnj ! looping over all processors in j axis 
    567          DO_2D_00_00 
     569         DO_2D( 0, 0, 0, 0 ) 
    568570            p_fval(ji,jj) = p_fval(ji,jj-1) + pva(ji,jj) * tmask_i(ji,jj) 
    569571         END_2D 
    570          CALL lbc_lnk( 'diaptr', p_fval, 'U', -1. ) 
     572         CALL lbc_lnk( 'diaptr', p_fval, 'U', -1.0_wp ) 
    571573      END DO 
    572574      !  
     
    604606      p_fval(:,:) = 0._wp 
    605607      ! 
    606       DO_3D_00_00( 1, jpkm1 ) 
     608      DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 
    607609         p_fval(jj,jk) = p_fval(jj,jk) + pta(ji,jj,jk) * pmsk(ji,jj) * tmask_i(ji,jj) 
    608610      END_3D 
Note: See TracChangeset for help on using the changeset viewer.