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 13295 for NEMO/trunk/src/TOP/PISCES/SED – NEMO

Ignore:
Timestamp:
2020-07-10T20:24:21+02:00 (4 years ago)
Author:
acc
Message:

Replace do-loop macros in the trunk with alternative forms with greater flexibility for extra halo applications. This alters a lot of routines but does not change any behaviour or results. do_loop_substitute.h90 is greatly simplified by this change. SETTE results are identical to those with the previous revision

Location:
NEMO/trunk/src/TOP/PISCES/SED
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/src/TOP/PISCES/SED/sedchem.F90

    r12839 r13295  
    138138         CALL sed_chem_cst 
    139139      ELSE 
    140          DO_2D_11_11 
     140         DO_2D( 1, 1, 1, 1 ) 
    141141            ikt = mbkt(ji,jj)  
    142142            IF ( tmask(ji,jj,ikt) == 1 ) THEN 

  • NEMO/trunk/src/TOP/PISCES/SED/seddta.F90

    r13237 r13295  
    9696      !    ----------------------------------------------------------- 
    9797      IF (ln_sediment_offline) THEN 
    98          DO_2D_11_11 
     98         DO_2D( 1, 1, 1, 1 ) 
    9999            ikt = mbkt(ji,jj) 
    100100            zwsbio4(ji,jj) = wsbio2 / rday 
    … …  
    102102         END_2D 
    103103      ELSE 
    104          DO_2D_11_11 
     104         DO_2D( 1, 1, 1, 1 ) 
    105105            ikt = mbkt(ji,jj) 
    106106            zdep = e3t(ji,jj,ikt,Kmm) / rDt_trc 
    … …  
    111111 
    112112      trc_data(:,:,:) = 0. 
    113       DO_2D_11_11 
     113      DO_2D( 1, 1, 1, 1 ) 
    114114         ikt = mbkt(ji,jj) 
    115115         IF ( tmask(ji,jj,ikt) == 1 ) THEN 

  • NEMO/trunk/src/TOP/PISCES/SED/sedini.F90

    r12489 r13295  
    135135      ! Determination of sediments number of points and allocate global variables 
    136136      epkbot(:,:) = 0. 
    137       DO_2D_11_11 
     137      DO_2D( 1, 1, 1, 1 ) 
    138138         ikt = mbkt(ji,jj)  
    139139         IF( tmask(ji,jj,ikt) == 1 ) epkbot(ji,jj) = e3t_1d(ikt) 
    … …  
    247247      ! Computation of 1D array of sediments points 
    248248      indoce = 0 
    249       DO_2D_11_11 
     249      DO_2D( 1, 1, 1, 1 ) 
    250250         IF (  epkbot(ji,jj) > 0. ) THEN 
    251251            indoce          = indoce + 1 

  • NEMO/trunk/src/TOP/PISCES/SED/sedsfc.F90

    r12377 r13295  
    4848 
    4949 
    50       DO_2D_11_11 
     50      DO_2D( 1, 1, 1, 1 ) 
    5151         ikt = mbkt(ji,jj) 
    5252         IF ( tmask(ji,jj,ikt) == 1 ) THEN 

  • NEMO/trunk/src/TOP/PISCES/SED/trcdmp_sed.F90

    r13286 r13295  
    9393               CALL trc_dta( kt, Kmm, sf_trcdta(jl), rf_trfac(jl), ztrcdta )   ! read tracer data at nit000 
    9494               ! 
    95                DO_2D_11_11 
     95               DO_2D( 1, 1, 1, 1 ) 
    9696                  ikt = mbkt(ji,jj) 
    9797                  tr(ji,jj,ikt,jn,Kbb) = ztrcdta(ji,jj,ikt) + ( tr(ji,jj,ikt,jn,Kbb) -  ztrcdta(ji,jj,ikt) )     & 
Note: See TracChangeset for help on using the changeset viewer.