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 15053 – NEMO

Changeset 15053


Ignore:
Timestamp:
2021-06-24T17:39:38+02:00 (3 years ago)
Author:
clem
Message:

nn_hls=2: solve a number of floating points because some arrays are not defined on the first halo. Looks like this is not the end.

Location:
NEMO/trunk/src/OCE
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/src/OCE/ISF/isfcavmlt.F90

    r15004 r15053  
    213213      ! compute upward heat flux zhtflx and upward water flux zwflx 
    214214      ! Resolution of a 3d equation from equation 24, 25 and 26 (note conduction through the ice has been added to Eq 24) 
    215       DO_2D( 1, 1, 1, 1 ) 
     215      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    216216         ! 
    217217         ! compute coeficient to solve the 2nd order equation 

  • NEMO/trunk/src/OCE/ISF/isfcpl.F90

    r15004 r15053  
    374374      ! ----------------------------------------------------------------------------------------- 
    375375      ! case we open a cell but no neigbour cells available to get an estimate of T and S 
    376       DO_3D( 1, 1, 1, 1, 1,jpk-1 ) 
     376      DO_3D( 0, 0, 0, 0, 1,jpk-1 ) 
    377377         IF (tmask(ji,jj,jk) == 1._wp .AND. ts(ji,jj,jk,2,Kmm) == 0._wp)              & 
    378378            &   CALL ctl_stop('STOP', 'failing to fill all new weet cell,     & 

  • NEMO/trunk/src/OCE/ISF/isfdiags.F90

    r15004 r15053  
    101101      zvar3d(:,:,:) = 0._wp 
    102102      ! 
    103       DO_2D( 1, 1, 1, 1 ) 
     103      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    104104         ikt = ktop(ji,jj) 
    105105         ikb = kbot(ji,jj) 

  • NEMO/trunk/src/OCE/ISF/isfload.F90

    r14064 r15053  
    9494      !                                !- Surface value + ice shelf gradient 
    9595      pload(:,:) = 0._wp                      ! compute pressure due to ice shelf load  
    96       DO_2D( 1, 1, 1, 1 ) 
     96      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    9797         ikt = mikt(ji,jj) 
    9898         ! 

  • NEMO/trunk/src/OCE/ISF/isftbl.F90

    r14834 r15053  
    7878         ! compute tbl property at T point 
    7979         pvarout(1,:) = 0._wp 
    80          DO_2D( 0, 1, 1, 1 ) 
     80         DO_2D( nn_hls-1, nn_hls, nn_hls, nn_hls ) 
    8181            pvarout(ji,jj) = 0.5_wp * (zvarout(ji,jj) + zvarout(ji-1,jj)) 
    8282         END_2D 
    … …  
    9999         ! pvarout is an averaging of wet point 
    100100         pvarout(:,1) = 0._wp 
    101          DO_2D( 1, 1, 0, 1 ) 
     101         DO_2D( nn_hls, nn_hls, nn_hls-1, nn_hls ) 
    102102            pvarout(ji,jj) = 0.5_wp * (zvarout(ji,jj) + zvarout(ji,jj-1)) 
    103103         END_2D 
    … …  
    138138      ! 
    139139      ! compute tbl top.bottom level and thickness 
    140       DO_2D( 1, 1, 1, 1 ) 
     140      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    141141         ! 
    142142         ! tbl top/bottom indices initialisation 

  • NEMO/trunk/src/OCE/TRA/eosbn2.F90

    r14834 r15053  
    10781078      z1_T0   = 1._wp/40._wp 
    10791079      ! 
    1080       DO_2D( 1, 1, 1, 1 ) 
     1080      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    10811081         ! 
    10821082         zt  = ctmp   (ji,jj) * z1_T0 
    … …  
    11421142         ! 
    11431143         z1_S0 = 1._wp / 35.16504_wp 
    1144          DO_2D( 1, 1, 1, 1 ) 
     1144         DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    11451145            zs= SQRT( ABS( psal(ji,jj) ) * z1_S0 )           ! square root salinity 
    11461146            ptf(ji,jj) = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & 
    … …  
    12491249      CASE( np_teos10, np_eos80 )                !==  polynomial TEOS-10 / EOS-80 ==! 
    12501250         ! 
    1251          DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
     1251         DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 
    12521252            ! 
    12531253            zh  = gdept(ji,jj,jk,Kmm) * r1_Z0                                ! depth 
    … …  
    13081308      CASE( np_seos )                !==  Vallis (2006) simplified EOS  ==! 
    13091309         ! 
    1310          DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
     1310         DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 
    13111311            zt  = pts(ji,jj,jk,jp_tem) - 10._wp  ! temperature anomaly (t-T0) 
    13121312            zs = pts (ji,jj,jk,jp_sal) - 35._wp  ! abs. salinity anomaly (s-S0) 

  • NEMO/trunk/src/OCE/TRA/trabbl.F90

    r14834 r15053  
    339339      ENDIF 
    340340      !                                        !* bottom variables (T, S, alpha, beta, depth, velocity) 
    341       DO_2D( 1, 1, 1, 1 ) 
     341      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 
    342342         ik = mbkt(ji,jj)                             ! bottom T-level index 
    343343         zts (ji,jj,jp_tem) = ts(ji,jj,ik,jp_tem,Kbb) ! bottom before T and S 
Note: See TracChangeset for help on using the changeset viewer.