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 for NEMO/trunk/src/OCE/ISF – NEMO

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/ISF
Files:
5 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 
Note: See TracChangeset for help on using the changeset viewer.