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 13710 for NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/TOP/PISCES/P4Z/p4zsed.F90 – NEMO

Ignore:
Timestamp:
2020-11-02T10:56:42+01:00 (4 years ago)
Author:
emanuelaclementi
Message:

branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves: merge with trunk@13708, see #2155 and #2339

Location:
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves

    • 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@13559        sette 

  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/TOP/PISCES/P4Z/p4zsed.F90

    r12377 r13710  
    1818   USE sed             !  Sediment module 
    1919   USE iom             !  I/O manager 
    20    USE prtctl_trc      !  print control for debugging 
     20   USE prtctl          !  print control for debugging 
    2121 
    2222   IMPLICIT NONE 
    … …  
    3939   !! * Substitutions 
    4040#  include "do_loop_substitute.h90" 
     41#  include "domzgr_substitute.h90" 
    4142   !!---------------------------------------------------------------------- 
    4243   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    … …  
    9394         ! OA: Warning, the following part is necessary to avoid CFL problems above the sediments 
    9495         ! -------------------------------------------------------------------- 
    95          DO_2D_11_11 
     96         DO_2D( 1, 1, 1, 1 ) 
    9697            ikt  = mbkt(ji,jj) 
    9798            zdep = e3t(ji,jj,ikt,Kmm) / xstep 
    … …  
    103104         ! Computation of the fraction of organic matter that is permanently buried from Dunne's model 
    104105         ! ------------------------------------------------------- 
    105          DO_2D_11_11 
     106         DO_2D( 1, 1, 1, 1 ) 
    106107           IF( tmask(ji,jj,1) == 1 ) THEN 
    107108              ikt = mbkt(ji,jj) 
    … …  
    129130      IF( .NOT.lk_sed )  zrivsil = 1._wp - sedsilfrac 
    130131 
    131       DO_2D_11_11 
     132      DO_2D( 1, 1, 1, 1 ) 
    132133         ikt  = mbkt(ji,jj) 
    133134         zdep = xstep / e3t(ji,jj,ikt,Kmm)  
    … …  
    141142      ! 
    142143      IF( .NOT.lk_sed ) THEN 
    143          DO_2D_11_11 
     144         DO_2D( 1, 1, 1, 1 ) 
    144145            ikt  = mbkt(ji,jj) 
    145146            zdep = xstep / e3t(ji,jj,ikt,Kmm)  
    … …  
    159160      ENDIF 
    160161      ! 
    161       DO_2D_11_11 
     162      DO_2D( 1, 1, 1, 1 ) 
    162163         ikt  = mbkt(ji,jj) 
    163164         zdep = xstep / e3t(ji,jj,ikt,Kmm)  
    … …  
    171172      ! 
    172173      IF( ln_p5z ) THEN 
    173          DO_2D_11_11 
     174         DO_2D( 1, 1, 1, 1 ) 
    174175            ikt  = mbkt(ji,jj) 
    175176            zdep = xstep / e3t(ji,jj,ikt,Kmm)  
    … …  
    186187         ! The 0.5 factor in zpdenit is to avoid negative NO3 concentration after 
    187188         ! denitrification in the sediments. Not very clever, but simpliest option. 
    188          DO_2D_11_11 
     189         DO_2D( 1, 1, 1, 1 ) 
    189190            ikt  = mbkt(ji,jj) 
    190191            zdep = xstep / e3t(ji,jj,ikt,Kmm)  
    … …  
    223224      ENDDO 
    224225      IF( ln_p4z ) THEN 
    225          DO_3D_11_11( 1, jpkm1 ) 
     226         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    226227            !                      ! Potential nitrogen fixation dependant on temperature and iron 
    227228            ztemp = ts(ji,jj,jk,jp_tem,Kmm) 
    … …  
    239240         END_3D 
    240241      ELSE       ! p5z 
    241          DO_3D_11_11( 1, jpkm1 ) 
     242         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    242243            !                      ! Potential nitrogen fixation dependant on temperature and iron 
    243244            ztemp = ts(ji,jj,jk,jp_tem,Kmm) 
    … …  
    260261      ! ---------------------------------------- 
    261262      IF( ln_p4z ) THEN 
    262          DO_3D_11_11( 1, jpkm1 ) 
     263         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    263264            zfact = nitrpot(ji,jj,jk) * nitrfix 
    264265            tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) + zfact / 3.0 
    … …  
    277278         END_3D 
    278279      ELSE    ! p5z 
    279          DO_3D_11_11( 1, jpkm1 ) 
     280         DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
    280281            zfact = nitrpot(ji,jj,jk) * nitrfix 
    281282            tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) + zfact / 3.0 
    … …  
    312313      ENDIF 
    313314      ! 
    314       IF(sn_cfctl%l_prttrc) THEN  ! print mean trends (USEd for debugging) 
     315      IF(sn_cfctl%l_prttrc) THEN  ! print mean trneds (USEd for debugging) 
    315316         WRITE(charout, fmt="('sed ')") 
    316          CALL prt_ctl_trc_info(charout) 
    317          CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 
     317         CALL prt_ctl_info( charout, cdcomp = 'top' ) 
     318         CALL prt_ctl(tab4d_1=tr(:,:,:,:,Krhs), mask1=tmask, clinfo=ctrcnm) 
    318319      ENDIF 
    319320      ! 
    … …  
    365366      lk_sed = ln_sediment .AND. ln_sed_2way  
    366367      ! 
     368      nitrpot(:,:,jpk) = 0._wp   ! define last level for iom_put 
     369      ! 
    367370   END SUBROUTINE p4z_sed_init 
    368371 
Note: See TracChangeset for help on using the changeset viewer.