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/diahth.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/diahth.F90

    r12511 r13540  
    4242   !! * Substitutions 
    4343#  include "do_loop_substitute.h90" 
     44#  include "domzgr_substitute.h90" 
    4445   !!---------------------------------------------------------------------- 
    4546   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    … …  
    129130            zdepinv(:,:) = 0._wp   
    130131            zmaxdzT(:,:) = 0._wp   
    131             DO_2D_11_11 
     132            DO_2D( 1, 1, 1, 1 ) 
    132133               zztmp = gdepw(ji,jj,mbkt(ji,jj)+1,Kmm)  
    133134               hth     (ji,jj) = zztmp 
    … …  
    138139            END_2D 
    139140            IF( nla10 > 1 ) THEN  
    140                DO_2D_11_11 
     141               DO_2D( 1, 1, 1, 1 ) 
    141142                  zztmp = gdepw(ji,jj,mbkt(ji,jj)+1,Kmm)  
    142143                  zrho0_3(ji,jj) = zztmp 
    … …  
    147148            ! Preliminary computation 
    148149            ! computation of zdelr = (dr/dT)(T,S,10m)*(-0.2 degC) 
    149             DO_2D_11_11 
     150            DO_2D( 1, 1, 1, 1 ) 
    150151               IF( tmask(ji,jj,nla10) == 1. ) THEN 
    151152                  zu  =  1779.50 + 11.250 * ts(ji,jj,nla10,jp_tem,Kmm) - 3.80   * ts(ji,jj,nla10,jp_sal,Kmm)  & 
    … …  
    169170            ! MLD: rho = rho(1) + zrho1                                     ! 
    170171            ! ------------------------------------------------------------- ! 
    171             DO_3DS_11_11( jpkm1, 2, -1 ) 
     172            DO_3DS( 1, 1, 1, 1, jpkm1, 2, -1 )   ! loop from bottom to 2 
    172173               ! 
    173174               zzdep = gdepw(ji,jj,jk,Kmm) 
    … …  
    206207            ! depth of temperature inversion                                ! 
    207208            ! ------------------------------------------------------------- ! 
    208             DO_3DS_11_11( jpkm1, nlb10, -1 ) 
     209            DO_3DS( 1, 1, 1, 1, jpkm1, nlb10, -1 )   ! loop from bottom to nlb10 
    209210               ! 
    210211               zzdep = gdepw(ji,jj,jk,Kmm) * tmask(ji,jj,1) 
    … …  
    304305      ! --------------------------------------- ! 
    305306      iktem(:,:) = 1 
    306       DO_3D_11_11( 1, jpkm1 ) 
     307      DO_3D( 1, 1, 1, 1, 1, jpkm1 )   ! beware temperature is not always decreasing with depth => loop from top to bottom 
    307308         zztmp = ts(ji,jj,jk,jp_tem,Kmm) 
    308309         IF( zztmp >= ptem )   iktem(ji,jj) = jk 
    … …  
    312313      !  Depth of ptem isotherm         ! 
    313314      ! ------------------------------- ! 
    314       DO_2D_11_11 
     315      DO_2D( 1, 1, 1, 1 ) 
    315316         ! 
    316317         zzdep = gdepw(ji,jj,mbkt(ji,jj)+1,Kmm)       ! depth of the ocean bottom 
    … …  
    350351      ! 
    351352      ilevel(:,:) = 1 
    352       DO_3D_11_11( 2, jpkm1 ) 
     353      DO_3D( 1, 1, 1, 1, 2, jpkm1 ) 
    353354         IF( ( gdept(ji,jj,jk,Kmm) < pdep ) .AND. ( tmask(ji,jj,jk) == 1 ) ) THEN 
    354355             ilevel(ji,jj) = jk 
    … …  
    358359      END_3D 
    359360      ! 
    360       DO_2D_11_11 
     361      DO_2D( 1, 1, 1, 1 ) 
    361362         ik = ilevel(ji,jj) 
    362363         zthick(ji,jj) = pdep - zthick(ji,jj)   !   remaining thickness to reach depht pdep 
    363          phtc(ji,jj)   = phtc(ji,jj) + pt(ji,jj,ik+1) * MIN( e3t(ji,jj,ik+1,Kmm), zthick(ji,jj) ) & 
     364         phtc(ji,jj)   = phtc(ji,jj)    & 
     365            &           + pt (ji,jj,ik+1) * MIN( e3t(ji,jj,ik+1,Kmm), zthick(ji,jj) ) & 
    364366                                                       * tmask(ji,jj,ik+1) 
    365367      END_2D 
Note: See TracChangeset for help on using the changeset viewer.