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 15145 for NEMO/trunk/tests/ISOMIP+/MY_SRC/eosbn2.F90 – NEMO

Ignore:
Timestamp:
2021-07-26T18:16:45+02:00 (3 years ago)
Author:
smasson
Message:

trunk: pass all sette tests in debug with nn_hls = 2

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/ISOMIP+/MY_SRC/eosbn2.F90

    r14995 r15145  
    11851185      z1_T0   = 1._wp/40._wp 
    11861186      ! 
    1187       DO_2D( 1, 1, 1, 1 ) 
     1187      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    11881188         ! 
    11891189         zt  = ctmp   (ji,jj) * z1_T0 
    1190          zs  = SQRT( ABS( psal(ji,jj) + zdeltaS ) * r1_S0 ) 
     1190         zs  = SQRT( ABS( psal(ji,jj) + zdeltaS ) * z1_S0 ) 
    11911191         ztm = tmask(ji,jj,1) 
    11921192         ! 
     
    12491249         ! 
    12501250         z1_S0 = 1._wp / 35.16504_wp 
    1251          DO_2D( 1, 1, 1, 1 ) 
     1251         DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    12521252            zs= SQRT( ABS( psal(ji,jj) ) * z1_S0 )           ! square root salinity 
    12531253            ptf(ji,jj) = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & 
     
    13561356      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    13571357         ! 
    1358          DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
     1358         DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 
    13591359            ! 
    13601360            zh  = gdept(ji,jj,jk,Kmm) * r1_Z0                                ! depth 
     
    14151415      CASE( np_seos )                !==  Vallis (2006) simplified EOS  ==! 
    14161416         ! 
    1417          DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
     1417         DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 
    14181418            zt  = pts(ji,jj,jk,jp_tem) - 10._wp  ! temperature anomaly (t-T0) 
    14191419            zs = pts (ji,jj,jk,jp_sal) - 35._wp  ! abs. salinity anomaly (s-S0) 
Note: See TracChangeset for help on using the changeset viewer.