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/TRA/eosbn2.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/TRA/eosbn2.F90

    r12511 r13540  
    180180   !! * Substitutions 
    181181#  include "do_loop_substitute.h90" 
     182#  include "domzgr_substitute.h90" 
    182183   !!---------------------------------------------------------------------- 
    183184   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    … …  
    237238      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    238239         ! 
    239          DO_3D_11_11( 1, jpkm1 ) 
     240         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    240241            ! 
    241242            zh  = pdep(ji,jj,jk) * r1_Z0                                  ! depth 
    … …  
    273274      CASE( np_seos )                !==  simplified EOS  ==! 
    274275         ! 
    275          DO_3D_11_11( 1, jpkm1 ) 
     276         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    276277            zt  = pts  (ji,jj,jk,jp_tem) - 10._wp 
    277278            zs  = pts  (ji,jj,jk,jp_sal) - 35._wp 
    … …  
    337338            END DO 
    338339            ! 
    339             DO_3D_11_11( 1, jpkm1 ) 
     340            DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    340341               ! 
    341342               ! compute density (2*nn_sto_eos) times: 
    … …  
    387388         ! Non-stochastic equation of state 
    388389         ELSE 
    389             DO_3D_11_11( 1, jpkm1 ) 
     390            DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    390391               ! 
    391392               zh  = pdep(ji,jj,jk) * r1_Z0                                  ! depth 
    … …  
    425426      CASE( np_seos )                !==  simplified EOS  ==! 
    426427         ! 
    427          DO_3D_11_11( 1, jpkm1 ) 
     428         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    428429            zt  = pts  (ji,jj,jk,jp_tem) - 10._wp 
    429430            zs  = pts  (ji,jj,jk,jp_sal) - 35._wp 
    … …  
    479480      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    480481         ! 
    481          DO_2D_11_11 
     482         DO_2D( 1, 1, 1, 1 ) 
    482483            ! 
    483484            zh  = pdep(ji,jj) * r1_Z0                                  ! depth 
    … …  
    514515      CASE( np_seos )                !==  simplified EOS  ==! 
    515516         ! 
    516          DO_2D_11_11 
     517         DO_2D( 1, 1, 1, 1 ) 
    517518            ! 
    518519            zt    = pts  (ji,jj,jp_tem)  - 10._wp 
    … …  
    562563      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    563564         ! 
    564          DO_3D_11_11( 1, jpkm1 ) 
     565         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    565566            ! 
    566567            zh  = gdept(ji,jj,jk,Kmm) * r1_Z0                                ! depth 
    … …  
    615616      CASE( np_seos )                  !==  simplified EOS  ==! 
    616617         ! 
    617          DO_3D_11_11( 1, jpkm1 ) 
     618         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    618619            zt  = pts (ji,jj,jk,jp_tem) - 10._wp   ! pot. temperature anomaly (t-T0) 
    619620            zs  = pts (ji,jj,jk,jp_sal) - 35._wp   ! abs. salinity anomaly (s-S0) 
    … …  
    669670      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    670671         ! 
    671          DO_2D_11_11 
     672         DO_2D( 1, 1, 1, 1 ) 
    672673            ! 
    673674            zh  = pdep(ji,jj) * r1_Z0                                  ! depth 
    … …  
    722723      CASE( np_seos )                  !==  simplified EOS  ==! 
    723724         ! 
    724          DO_2D_11_11 
     725         DO_2D( 1, 1, 1, 1 ) 
    725726            ! 
    726727            zt    = pts  (ji,jj,jp_tem) - 10._wp   ! pot. temperature anomaly (t-T0) 
    … …  
    872873      IF( ln_timing )   CALL timing_start('bn2') 
    873874      ! 
    874       DO_3D_11_11( 2, jpkm1 ) 
     875      DO_3D( 1, 1, 1, 1, 2, jpkm1 )      ! interior points only (2=< jk =< jpkm1 ); surface and bottom value set to zero one for all in istate.F90 
    875876         zrw =   ( gdepw(ji,jj,jk  ,Kmm) - gdept(ji,jj,jk,Kmm) )   & 
    876877            &  / ( gdept(ji,jj,jk-1,Kmm) - gdept(ji,jj,jk,Kmm) )  
    … …  
    920921      z1_T0   = 1._wp/40._wp 
    921922      ! 
    922       DO_2D_11_11 
     923      DO_2D( 1, 1, 1, 1 ) 
    923924         ! 
    924925         zt  = ctmp   (ji,jj) * z1_T0 
    … …  
    973974         ! 
    974975         z1_S0 = 1._wp / 35.16504_wp 
    975          DO_2D_11_11 
     976         DO_2D( 1, 1, 1, 1 ) 
    976977            zs= SQRT( ABS( psal(ji,jj) ) * z1_S0 )           ! square root salinity 
    977978            ptf(ji,jj) = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & 
    … …  
    10801081      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    10811082         ! 
    1082          DO_3D_11_11( 1, jpkm1 ) 
     1083         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    10831084            ! 
    10841085            zh  = gdept(ji,jj,jk,Kmm) * r1_Z0                                ! depth 
    … …  
    11391140      CASE( np_seos )                !==  Vallis (2006) simplified EOS  ==! 
    11401141         ! 
    1141          DO_3D_11_11( 1, jpkm1 ) 
     1142         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    11421143            zt  = pts(ji,jj,jk,jp_tem) - 10._wp  ! temperature anomaly (t-T0) 
    11431144            zs = pts (ji,jj,jk,jp_sal) - 35._wp  ! abs. salinity anomaly (s-S0) 
Note: See TracChangeset for help on using the changeset viewer.