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 10321 for NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z – NEMO

Ignore:
Timestamp:
2018-11-16T15:59:30+01:00 (5 years ago)
Author:
davestorkey
Message:

UKMO/NEMO4_beta_mirror: Update to version 10279 of the trunk.

Location:
NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z
Files:
2 deleted
23 edited
2 copied

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zagg.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    2525 
    2626   !!---------------------------------------------------------------------- 
    27    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     27   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    2828   !! $Id$ 
    29    !! Software governed by the CeCILL licence     (./LICENSE) 
     29   !! Software governed by the CeCILL license (see ./LICENSE) 
    3030   !!---------------------------------------------------------------------- 
    3131CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zbio.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    1515   USE p4zsink         !  vertical flux of particulate matter due to sinking 
    1616   USE p4zopt          !  optical model 
    17    USE p4zice          !  Co-limitations of differents nutrients 
     17   USE p4zlim          !  Co-limitations of differents nutrients 
    1818   USE p4zprod         !  Growth rate of the 2 phyto groups 
    1919   USE p4zmort         !  Mortality terms for phytoplankton 
    2020   USE p4zmicro        !  Sources and sinks of microzooplankton 
    2121   USE p4zmeso         !  Sources and sinks of mesozooplankton 
    22    USE p5zice          !  Co-limitations of differents nutrients 
     22   USE p5zlim          !  Co-limitations of differents nutrients 
    2323   USE p5zprod         !  Growth rate of the 2 phyto groups 
    2424   USE p5zmort         !  Mortality terms for phytoplankton 
     
    3939 
    4040   !!---------------------------------------------------------------------- 
    41    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     41   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    4242   !! $Id$ 
    43    !! Software governed by the CeCILL licence (./LICENSE) 
     43   !! Software governed by the CeCILL license (see ./LICENSE) 
    4444   !!---------------------------------------------------------------------- 
    4545CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zche.F90

    r9950 r10321  
    3737   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   tempis   ! In situ temperature 
    3838 
    39    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   akb3       !: ??? 
    40    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   akw3       !: ??? 
    41    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   akf3       !: ??? 
    42    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   aks3       !: ??? 
    43    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   ak1p3      !: ??? 
    44    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   ak2p3      !: ??? 
    45    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   ak3p3      !: ??? 
    46    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   aksi3      !: ??? 
    47    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   borat      !: ??? 
    48    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   fluorid    !: ??? 
    49    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sulfat     !: ??? 
     39   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   akb3       !: ??? 
     40   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   akw3       !: ??? 
     41   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   akf3       !: ??? 
     42   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   aks3       !: ??? 
     43   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   ak1p3      !: ??? 
     44   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   ak2p3      !: ??? 
     45   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   ak3p3      !: ??? 
     46   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   aksi3      !: ??? 
     47   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   borat      !: ??? 
     48   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   fluorid    !: ??? 
     49   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sulfat     !: ??? 
    5050 
    5151   !!* Variable for chemistry of the CO2 cycle 
     
    131131 
    132132   !!---------------------------------------------------------------------- 
    133    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     133   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    134134   !! $Id$ 
    135    !! Software governed by the CeCILL licence     (./LICENSE) 
     135   !! Software governed by the CeCILL license (see ./LICENSE) 
    136136   !!---------------------------------------------------------------------- 
    137137CONTAINS 
     
    233233        END DO 
    234234 
    235  
    236  
    237235      ! CHEMICAL CONSTANTS - DEEP OCEAN 
    238236      ! ------------------------------- 
     
    449447      ! 
    450448   END SUBROUTINE p4z_che 
    451  
    452449 
    453450   SUBROUTINE ahini_for_at(p_hini) 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zfechem.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    3636 
    3737   !!---------------------------------------------------------------------- 
    38    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     38   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    3939   !! $Id$ 
    40    !! Software governed by the CeCILL licence     (./LICENSE) 
     40   !! Software governed by the CeCILL license (see ./LICENSE) 
    4141   !!---------------------------------------------------------------------- 
    4242CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zflx.F90

    r9950 r10321  
    5353 
    5454   !!---------------------------------------------------------------------- 
    55    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     55   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    5656   !! $Id$ 
    57    !! Software governed by the CeCILL licence     (./LICENSE) 
     57   !! Software governed by the CeCILL license (see ./LICENSE) 
    5858   !!---------------------------------------------------------------------- 
    5959CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zint.F90

    r9950 r10321  
    2020 
    2121   !!---------------------------------------------------------------------- 
    22    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     22   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    2323   !! $Id$ 
    24    !! Software governed by the CeCILL licence     (./LICENSE) 
     24   !! Software governed by the CeCILL license (see ./LICENSE) 
    2525   !!---------------------------------------------------------------------- 
    2626CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zligand.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    2727 
    2828   !!---------------------------------------------------------------------- 
    29    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     29   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    3030   !! $Id$ 
    31    !! Software governed by the CeCILL licence     (./LICENSE) 
     31   !! Software governed by the CeCILL license (see ./LICENSE) 
    3232   !!---------------------------------------------------------------------- 
    3333CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zlim.F90

    r10320 r10321  
    6969   !!---------------------------------------------------------------------- 
    7070   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    71    !! $Id: p4zlim.F90 10069 2018-08-28 14:12:24Z nicolasmartin $  
     71   !! $Id$ 
    7272   !! Software governed by the CeCILL license (see ./LICENSE) 
    7373   !!---------------------------------------------------------------------- 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zlys.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    3636  
    3737   !!---------------------------------------------------------------------- 
    38    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     38   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    3939   !! $Id$ 
    40    !! Software governed by the CeCILL licence (./LICENSE) 
     40   !! Software governed by the CeCILL license (see ./LICENSE) 
    4141   !!---------------------------------------------------------------------- 
    4242 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zmeso.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    4444 
    4545   !!---------------------------------------------------------------------- 
    46    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     46   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    4747   !! $Id$ 
    48    !! Software governed by the CeCILL licence (./LICENSE) 
     48   !! Software governed by the CeCILL license (see ./LICENSE) 
    4949   !!---------------------------------------------------------------------- 
    5050CONTAINS 
     
    108108               zdenom    = zfoodlim / ( xkgraz2 + zfoodlim ) 
    109109               zdenom2   = zdenom / ( zfood + rtrn ) 
    110                zgraze2   = grazrat2 * xstep * tgfunc2(ji,jj,jk) * trb(ji,jj,jk,jpmes)  
     110               zgraze2   = grazrat2 * xstep * tgfunc2(ji,jj,jk) * trb(ji,jj,jk,jpmes) * (1. - nitrfac(ji,jj,jk))  
    111111 
    112112               zgrazd    = zgraze2  * xprefc   * zcompadi  * zdenom2  
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zmicro.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    1414   USE trc             ! passive tracers common variables  
    1515   USE sms_pisces      ! PISCES Source Minus Sink variables 
    16    USE p4zice          ! Co-limitations 
     16   USE p4zlim          ! Co-limitations 
    1717   USE p4zprod         ! production 
    1818   USE iom             ! I/O manager 
     
    4242 
    4343   !!---------------------------------------------------------------------- 
    44    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     44   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    4545   !! $Id$ 
    46    !! Software governed by the CeCILL licence (./LICENSE) 
     46   !! Software governed by the CeCILL license (see ./LICENSE) 
    4747   !!---------------------------------------------------------------------- 
    4848CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zmort.F90

    r9950 r10321  
    1414   USE sms_pisces      ! PISCES Source Minus Sink variables 
    1515   USE p4zprod         ! Primary productivity  
    16    USE p4zice          ! Phytoplankton limitation terms 
     16   USE p4zlim          ! Phytoplankton limitation terms 
    1717   USE prtctl_trc      ! print control for debugging 
    1818 
     
    3030 
    3131   !!---------------------------------------------------------------------- 
    32    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     32   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    3333   !! $Id$ 
    34    !! Software governed by the CeCILL licence (./LICENSE) 
     34   !! Software governed by the CeCILL license (see ./LICENSE) 
    3535   !!---------------------------------------------------------------------- 
    3636CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zopt.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    4343    
    4444   !!---------------------------------------------------------------------- 
    45    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     45   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    4646   !! $Id$ 
    47    !! Software governed by the CeCILL licence     (./LICENSE) 
     47   !! Software governed by the CeCILL license (see ./LICENSE) 
    4848   !!---------------------------------------------------------------------- 
    4949CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zpoc.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    3838 
    3939   !!---------------------------------------------------------------------- 
    40    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     40   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    4141   !! $Id$ 
    42    !! Software governed by the CeCILL licence     (./LICENSE) 
     42   !! Software governed by the CeCILL license (see ./LICENSE) 
    4343   !!---------------------------------------------------------------------- 
    4444CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zprod.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    1515   USE trc             ! passive tracers common variables  
    1616   USE sms_pisces      ! PISCES Source Minus Sink variables 
    17    USE p4zice          ! Co-limitations of differents nutrients 
     17   USE p4zlim          ! Co-limitations of differents nutrients 
    1818   USE prtctl_trc      ! print control for debugging 
    1919   USE iom             ! I/O manager 
     
    4949 
    5050   !!---------------------------------------------------------------------- 
    51    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     51   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    5252   !! $Id$ 
    53    !! Software governed by the CeCILL licence     (./LICENSE) 
     53   !! Software governed by the CeCILL license (see ./LICENSE) 
    5454   !!---------------------------------------------------------------------- 
    5555CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zrem.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    1717   USE p4zche          !  chemical model 
    1818   USE p4zprod         !  Growth rate of the 2 phyto groups 
    19    USE p4zice 
     19   USE p4zlim 
    2020   USE prtctl_trc      !  print control for debugging 
    2121   USE iom             !  I/O manager 
     
    4343 
    4444   !!---------------------------------------------------------------------- 
    45    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     45   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    4646   !! $Id$ 
    47    !! Software governed by the CeCILL licence     (./LICENSE) 
     47   !! Software governed by the CeCILL license (see ./LICENSE) 
    4848   !!---------------------------------------------------------------------- 
    4949CONTAINS 
     
    116116                  zammonic = zremik * nitrfac(ji,jj,jk) * trb(ji,jj,jk,jpdoc) 
    117117                  denitr(ji,jj,jk)  = zammonic * ( 1. - nitrfac2(ji,jj,jk) ) 
    118                   zoxyrem           = zammonic *        nitrfac2(ji,jj,jk) 
     118                  denitr(ji,jj,jk)  = MIN( ( trb(ji,jj,jk,jpno3) - rtrn ) / rdenit, denitr(ji,jj,jk) ) 
     119                  zoxyrem           = zammonic - denitr(ji,jj,jk) 
    119120                  ! 
    120121                  zolimi (ji,jj,jk) = MAX( 0.e0, zolimi (ji,jj,jk) ) 
     
    189190               &         / ( 1.+ emoy(ji,jj,jk) ) * ( 1. + fr_i(ji,jj) * emoy(ji,jj,jk) )  
    190191               zdenitnh4 = nitrif * xstep * trb(ji,jj,jk,jpnh4) * nitrfac(ji,jj,jk) 
     192               zdenitnh4 = MIN(  ( trb(ji,jj,jk,jpno3) - rtrn ) / rdenita, zdenitnh4 )  
    191193               ! Update of the tracers trends 
    192194               ! ---------------------------- 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zsbc.F90

    r9950 r10321  
    3737   REAL(wp), PUBLIC ::   concfediaz   !: Fe half-saturation Cste for diazotrophs  
    3838   REAL(wp)         ::   hratio       !: Fe:3He ratio assumed for vent iron supply 
     39   REAL(wp)         ::   distcoast    !: Distance off the coast for Iron from sediments 
    3940   REAL(wp), PUBLIC ::   fep_rats     !: Fep/Fer ratio from sed  sources 
    4041   REAL(wp), PUBLIC ::   fep_rath     !: Fep/Fer ratio from hydro sources 
     
    7475   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   hydrofe          !: Hydrothermal vent supply of iron 
    7576 
    76    REAL(wp), PUBLIC ::   rivalkinput, rivdicinput, nitdepinput, sumdepsi 
    77    REAL(wp), PUBLIC ::   rivdininput, rivdipinput, rivdsiinput 
     77   REAL(wp), PUBLIC :: sedsilfrac, sedcalfrac 
     78   REAL(wp), PUBLIC :: rivalkinput, rivdicinput 
     79   REAL(wp), PUBLIC :: rivdininput, rivdipinput, rivdsiinput 
    7880 
    7981   !! * Substitutions 
    8082#  include "vectopt_loop_substitute.h90" 
    8183   !!---------------------------------------------------------------------- 
    82    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     84   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    8385   !! $Id$ 
    84    !! Software governed by the CeCILL licence (./LICENSE) 
     86   !! Software governed by the CeCILL license (see ./LICENSE) 
    8587   !!---------------------------------------------------------------------- 
    8688CONTAINS 
     
    109111         IF( kt == nit000 .OR. ( kt /= nit000 .AND. ntimes_dust > 1 ) ) THEN 
    110112            CALL fld_read( kt, 1, sf_dust ) 
    111             IF( nn_ice_tr == -1 .AND. .NOT.ln_ironice ) THEN   ;   dust(:,:) = sf_dust(1)%fnow(:,:,1) 
    112             ELSE                                               ;   dust(:,:) = sf_dust(1)%fnow(:,:,1) * ( 1.-fr_i(:,:) ) 
     113            IF( nn_ice_tr == -1 .AND. .NOT.ln_ironice ) THEN   ;   dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) ) 
     114            ELSE                                               ;   dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) * ( 1.-fr_i(:,:) ) ) 
    113115            ENDIF 
    114116         ENDIF 
     
    173175             zcoef = rno3 * 14E6 * ryyss 
    174176             CALL fld_read( kt, 1, sf_ndepo ) 
    175              nitdep(:,:) = sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1)  
     177             nitdep(:,:) = MAX( rtrn, sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1) ) 
    176178         ENDIF 
    177179         IF( .NOT.ln_linssh ) THEN 
    178180           zcoef = rno3 * 14E6 * ryyss 
    179            nitdep(:,:) = sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1)  
     181           nitdep(:,:) = MAX( rtrn, sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1) ) 
    180182         ENDIF 
    181183      ENDIF 
     
    205207      INTEGER  :: ik50                !  last level where depth less than 50 m 
    206208      INTEGER  :: isrow             ! index for ORCA1 starting row 
    207       REAL(wp) :: zexpide, zdenitide, zmaskt 
     209      REAL(wp) :: zexpide, zdenitide, zmaskt, zsurfc, zsurfp,ze3t, ze3t2, zcslp 
    208210      REAL(wp) :: ztimes_dust, ztimes_riv, ztimes_ndep  
    209211      REAL(wp), DIMENSION(nbtimes) :: zsteps                 ! times records 
    210212      REAL(wp), DIMENSION(:), ALLOCATABLE :: rivinput 
    211       REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: zdust, zndepo, zriver, zcmask 
     213      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: zriver, zcmask 
    212214      ! 
    213215      CHARACTER(len=100) ::  cn_dir          ! Root directory for location of ssr files 
     
    220222        &                sn_riverdip, sn_riverdop, sn_riverdsi, sn_ndepo, sn_ironsed, sn_hydrofe, & 
    221223        &                ln_dust, ln_solub, ln_river, ln_ndepo, ln_ironsed, ln_ironice, ln_hydrofe,    & 
    222         &                sedfeinput, dustsolub, icefeinput, wdust, mfrac, nitrfix, diazolight, concfediaz, & 
     224        &                sedfeinput, distcoast, dustsolub, icefeinput, wdust, mfrac, nitrfix, diazolight, concfediaz, & 
    223225        &                hratio, fep_rats, fep_rath, lgw_rath 
    224226      !!---------------------------------------------------------------------- 
     
    248250         WRITE(numout,*) '      fe input from hydrothermal vents         ln_hydrofe  = ', ln_hydrofe 
    249251         WRITE(numout,*) '      coastal release of iron                  sedfeinput  = ', sedfeinput 
     252         WRITE(numout,*) '      distance off the coast                   distcoast   = ', distcoast 
    250253         WRITE(numout,*) '      solubility of the dust                   dustsolub   = ', dustsolub 
    251254         WRITE(numout,*) '      Mineral Fe content of the dust           mfrac       = ', mfrac 
     
    307310            CALL iom_open (  TRIM( sn_dust%clname ) , numdust ) 
    308311            CALL iom_gettime( numdust, zsteps, kntime=ntimes_dust)  ! get number of record in file 
    309             ALLOCATE( zdust(jpi,jpj,ntimes_dust) ) 
    310             DO jm = 1, ntimes_dust 
    311                CALL iom_get( numdust, jpdom_data, TRIM( sn_dust%clvar ), zdust(:,:,jm), jm ) 
    312             END DO 
    313             CALL iom_close( numdust ) 
    314             ztimes_dust = 1._wp / REAL(ntimes_dust, wp)  
    315             sumdepsi = 0.e0 
    316             DO jm = 1, ntimes_dust 
    317                sumdepsi = sumdepsi + glob_sum( zdust(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_dust ) 
    318             END DO 
    319             sumdepsi = sumdepsi / ( nyear_len(1) * rday ) * 12. * 8.8 * 0.075 * mfrac / 28.1  
    320             DEALLOCATE( zdust) 
    321          ENDIF 
    322       ELSE 
    323          sumdepsi  = 0._wp 
     312         END IF 
    324313      END IF 
    325314 
     
    419408            CALL iom_open ( TRIM( sn_ndepo%clname ), numdepo ) 
    420409            CALL iom_gettime( numdepo, zsteps, kntime=ntimes_ndep) 
    421             ALLOCATE( zndepo(jpi,jpj,ntimes_ndep) ) 
    422             DO jm = 1, ntimes_ndep 
    423                CALL iom_get( numdepo, jpdom_data, TRIM( sn_ndepo%clvar ), zndepo(:,:,jm), jm ) 
    424             END DO 
    425             CALL iom_close( numdepo ) 
    426             ztimes_ndep = 1._wp / REAL(ntimes_ndep, wp)  
    427             nitdepinput = 0._wp 
    428             DO jm = 1, ntimes_ndep 
    429               nitdepinput = nitdepinput + glob_sum( zndepo(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_ndep ) 
    430             ENDDO 
    431             nitdepinput = nitdepinput / rno3 / 14E6  
    432             DEALLOCATE( zndepo) 
    433          ENDIF 
    434       ELSE 
    435          nitdepinput = 0._wp 
     410         ENDIF 
    436411      ENDIF 
    437412 
     
    459434            DO jj = 2, jpjm1 
    460435               DO ji = fs_2, fs_jpim1 
    461                   IF( tmask(ji,jj,jk) /= 0. ) THEN 
    462                      zmaskt = tmask(ji+1,jj,jk) * tmask(ji-1,jj,jk) * tmask(ji,jj+1,jk)    & 
    463                         &                       * tmask(ji,jj-1,jk) * tmask(ji,jj,jk+1) 
    464                      IF( zmaskt == 0. )   zcmask(ji,jj,jk ) = MAX( 0.1, zcmask(ji,jj,jk) )  
    465                   END IF 
     436                  ze3t   = e3t_0(ji,jj,jk) 
     437                  zsurfc =  e1u(ji,jj) * ( 1. - umask(ji  ,jj  ,jk) )   & 
     438                          + e1u(ji,jj) * ( 1. - umask(ji-1,jj  ,jk) )   & 
     439                          + e2v(ji,jj) * ( 1. - vmask(ji  ,jj  ,jk) )   & 
     440                          + e2v(ji,jj) * ( 1. - vmask(ji  ,jj-1,jk) ) 
     441                  zsurfp = zsurfc * ze3t / e1e2t(ji,jj) 
     442                  ! estimation of the coastal slope : 5 km off the coast 
     443                  ze3t2 = ze3t * ze3t 
     444                  zcslp = SQRT( ( distcoast*distcoast + ze3t2 ) / ze3t2 ) 
     445                  ! 
     446                  zcmask(ji,jj,jk) = zcmask(ji,jj,jk) + zcslp * zsurfp 
    466447               END DO 
    467448            END DO 
     
    519500         WRITE(numout,*) '    DIC Supply : ', rivdicinput*1E3*12./1E12     ,' TgC/yr' 
    520501         WRITE(numout,*)  
    521          WRITE(numout,*) '    Total input of elements from atmospheric supply' 
    522          WRITE(numout,*) '    ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    523          WRITE(numout,*) '    N Supply   : ', nitdepinput*rno3*1E3/1E12*14.,' TgN/yr' 
    524          WRITE(numout,*)  
    525       ENDIF 
     502      ENDIF 
     503      ! 
     504      sedsilfrac = 0.03     ! percentage of silica loss in the sediments 
     505      sedcalfrac = 0.6      ! percentage of calcite loss in the sediments 
    526506      ! 
    527507   END SUBROUTINE p4z_sbc_init 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zsed.F90

    r9950 r10321  
    1414   USE trc             !  passive tracers common variables  
    1515   USE sms_pisces      !  PISCES Source Minus Sink variables 
    16    USE p4zice          !  Co-limitations of differents nutrients 
     16   USE p4zlim          !  Co-limitations of differents nutrients 
    1717   USE p4zsbc          !  External source of nutrients  
    1818   USE p4zint          !  interpolation and computation of various fields 
     19   USE sed             !  Sediment module 
    1920   USE iom             !  I/O manager 
    2021   USE prtctl_trc      !  print control for debugging 
     
    2930   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:  ) :: sdenit     !: Nitrate reduction in the sediments 
    3031   REAL(wp) :: r1_rday                  !: inverse of rday 
     32   LOGICAL, SAVE :: lk_sed 
    3133 
    3234   !!---------------------------------------------------------------------- 
    33    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     35   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    3436   !! $Id$ 
    35    !! Software governed by the CeCILL licence (./LICENSE) 
     37   !! Software governed by the CeCILL license (see ./LICENSE) 
    3638   !!---------------------------------------------------------------------- 
    3739CONTAINS 
     
    4951      ! 
    5052      INTEGER, INTENT(in) ::   kt, knt ! ocean time step 
    51       INTEGER  ::   ji, jj, jk, ikt 
    52       REAL(wp) ::   zsumsedsi, zsumsedpo4, zsumsedcal 
    53       REAL(wp) ::   zrivalk, zrivsil, zrivno3 
     53      INTEGER  ::  ji, jj, jk, ikt 
     54      REAL(wp) ::  zrivalk, zrivsil, zrivno3 
    5455      REAL(wp) ::  zwflux, zfminus, zfplus 
    5556      REAL(wp) ::  zlim, zfact, zfactcal 
     
    6263      ! 
    6364      CHARACTER (len=25) :: charout 
    64       REAL(wp), DIMENSION(jpi,jpj    ) :: zwork1, zwork2, zwork3 
    65       REAL(wp), DIMENSION(jpi,jpj    ) :: zdenit2d, zbureff 
     65      REAL(wp), DIMENSION(jpi,jpj    ) :: zdenit2d, zbureff, zwork 
    6666      REAL(wp), DIMENSION(jpi,jpj    ) :: zwsbio3, zwsbio4, zwscal 
    6767      REAL(wp), DIMENSION(jpi,jpj    ) :: zsedcal, zsedsi, zsedc 
     
    7373      IF( ln_timing )  CALL timing_start('p4z_sed') 
    7474      ! 
     75      IF( kt == nittrc000 .AND. knt == 1 )   THEN 
     76          r1_rday  = 1. / rday 
     77          IF (ln_sediment .AND. ln_sed_2way) THEN 
     78             lk_sed = .TRUE. 
     79          ELSE 
     80             lk_sed = .FALSE. 
     81          ENDIF 
     82      ENDIF 
     83      ! 
    7584      IF( kt == nittrc000 .AND. knt == 1 )   r1_rday  = 1. / rday 
    7685      ! 
     
    8291      zdenit2d(:,:) = 0.e0 
    8392      zbureff (:,:) = 0.e0 
    84       zwork1  (:,:) = 0.e0 
    85       zwork2  (:,:) = 0.e0 
    86       zwork3  (:,:) = 0.e0 
     93      zwork   (:,:) = 0.e0 
    8794      zsedsi  (:,:) = 0.e0 
    8895      zsedcal (:,:) = 0.e0 
     
    189196      ENDIF 
    190197 
    191       ! Add the external input of iron from sediment mobilization 
    192       ! ------------------------------------------------------ 
    193       IF( ln_ironsed ) THEN 
    194                          tra(:,:,:,jpfer) = tra(:,:,:,jpfer) + ironsed(:,:,:) * rfact2 
    195          IF( ln_ligand ) tra(:,:,:,jpfep) = tra(:,:,:,jpfep) + ( ironsed(:,:,:) * fep_rats ) * rfact2 
    196          ! 
    197          IF( lk_iomput .AND. knt == nrdttrc .AND. iom_use( "Ironsed" ) )   & 
    198             &   CALL iom_put( "Ironsed", ironsed(:,:,:) * 1.e+3 * tmask(:,:,:) ) ! iron inputs from sediments 
    199       ENDIF 
    200  
    201198      ! Add the external input of iron from hydrothermal vents 
    202199      ! ------------------------------------------------------ 
     
    235232 
    236233      IF( .NOT.lk_sed ) THEN 
     234! 
     235         ! Add the external input of iron from sediment mobilization 
     236         ! ------------------------------------------------------ 
     237         IF( ln_ironsed ) THEN 
     238                            tra(:,:,:,jpfer) = tra(:,:,:,jpfer) + ironsed(:,:,:) * rfact2 
     239            IF( ln_ligand ) tra(:,:,:,jpfep) = tra(:,:,:,jpfep) + ( ironsed(:,:,:) * fep_rats ) * rfact2 
     240            ! 
     241            IF( lk_iomput .AND. knt == nrdttrc .AND. iom_use( "Ironsed" ) )   & 
     242               &   CALL iom_put( "Ironsed", ironsed(:,:,:) * 1.e+3 * tmask(:,:,:) ) ! iron inputs from sediments 
     243         ENDIF 
     244 
    237245         ! Computation of the sediment denitrification proportion: The metamodel from midlleburg (2006) is being used 
    238246         ! Computation of the fraction of organic matter that is permanently buried from Dunne's model 
     
    255263                   &     + trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) ) * 1E6 
    256264                 zbureff(ji,jj) = 0.013 + 0.53 * zflx**2 / ( 7.0 + zflx )**2 
    257                 ENDIF 
    258               END DO 
    259            END DO  
    260  
    261            ! Loss of biogenic silicon, Caco3 organic carbon in the sediments.  
    262            ! First, the total loss is computed. 
    263            ! The factor for calcite comes from the alkalinity effect 
    264            ! ------------------------------------------------------------- 
    265            DO jj = 1, jpj 
    266               DO ji = 1, jpi 
    267                  IF( tmask(ji,jj,1) == 1 ) THEN 
    268                     ikt = mbkt(ji,jj)  
    269                     zwork1(ji,jj) = trb(ji,jj,ikt,jpgsi) * zwsbio4(ji,jj) 
    270                     zwork2(ji,jj) = trb(ji,jj,ikt,jpgoc) * zwsbio4(ji,jj) + trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj)  
    271                     ! For calcite, burial efficiency is made a function of saturation 
    272                     zfactcal      = MIN( excess(ji,jj,ikt), 0.2 ) 
    273                     zfactcal      = MIN( 1., 1.3 * ( 0.2 - zfactcal ) / ( 0.4 - zfactcal ) ) 
    274                     zwork3(ji,jj) = trb(ji,jj,ikt,jpcal) * zwscal(ji,jj) * 2.e0 * zfactcal 
    275                 ENDIF 
    276             END DO 
    277          END DO 
    278          zsumsedsi  = glob_sum( zwork1(:,:) * e1e2t(:,:) ) * r1_rday 
    279          zsumsedpo4 = glob_sum( zwork2(:,:) * e1e2t(:,:) ) * r1_rday 
    280          zsumsedcal = glob_sum( zwork3(:,:) * e1e2t(:,:) ) * r1_rday 
     265              ENDIF 
     266            END DO 
     267         END DO  
    281268         ! 
    282269      ENDIF 
     
    285272      ! Thus, the amount of silica lost in the sediments equal the supply at the surface (dust+rivers) 
    286273      ! ------------------------------------------------------ 
    287       IF( .NOT.lk_sed )  zrivsil =  1._wp - ( sumdepsi + rivdsiinput * r1_ryyss ) / ( zsumsedsi + rtrn ) 
     274      IF( .NOT.lk_sed )  zrivsil = 1._wp - sedsilfrac 
    288275 
    289276      DO jj = 1, jpj 
     
    312299               zfactcal = MIN( excess(ji,jj,ikt), 0.2 ) 
    313300               zfactcal = MIN( 1., 1.3 * ( 0.2 - zfactcal ) / ( 0.4 - zfactcal ) ) 
    314                zrivalk  =  1._wp - ( rivalkinput * r1_ryyss ) * zfactcal / ( zsumsedcal + rtrn ) 
     301               zrivalk  = sedcalfrac * zfactcal 
    315302               tra(ji,jj,ikt,jptal) =  tra(ji,jj,ikt,jptal) + zcaloss * zrivalk * 2.0 
    316303               tra(ji,jj,ikt,jpdic) =  tra(ji,jj,ikt,jpdic) + zcaloss * zrivalk 
     
    492479            IF( iom_use("Nfix"   ) ) CALL iom_put( "Nfix", nitrpot(:,:,:) * nitrfix * rno3 * zfact * tmask(:,:,:) )  ! nitrogen fixation  
    493480            IF( iom_use("INTNFIX") ) THEN   ! nitrogen fixation rate in ocean ( vertically integrated ) 
    494                zwork1(:,:) = 0. 
     481               zwork(:,:) = 0. 
    495482               DO jk = 1, jpkm1 
    496                  zwork1(:,:) = zwork1(:,:) + nitrpot(:,:,jk) * nitrfix * rno3 * zfact * e3t_n(:,:,jk) * tmask(:,:,jk) 
     483                 zwork(:,:) = zwork(:,:) + nitrpot(:,:,jk) * nitrfix * rno3 * zfact * e3t_n(:,:,jk) * tmask(:,:,jk) 
    497484               ENDDO 
    498                CALL iom_put( "INTNFIX" , zwork1 )  
     485               CALL iom_put( "INTNFIX" , zwork )  
    499486            ENDIF 
    500487            IF( iom_use("SedCal" ) ) CALL iom_put( "SedCal", zsedcal(:,:) * zfact ) 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zsink.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    3939 
    4040   !!---------------------------------------------------------------------- 
    41    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     41   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    4242   !! $Id$ 
    43    !! Software governed by the CeCILL licence     (./LICENSE) 
     43   !! Software governed by the CeCILL license (see ./LICENSE) 
    4444   !!---------------------------------------------------------------------- 
    4545CONTAINS 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p4zsms.F90

    • Property svn:keywords set to Id
    r9950 r10321  
    4040 
    4141   !!---------------------------------------------------------------------- 
    42    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     42   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    4343   !! $Id$ 
    44    !! Software governed by the CeCILL licence (./LICENSE) 
     44   !! Software governed by the CeCILL license (see ./LICENSE) 
    4545   !!---------------------------------------------------------------------- 
    4646CONTAINS 
     
    101101      ENDIF 
    102102      ! 
     103      IF( ll_sbc ) CALL p4z_sbc( kt )   ! external sources of nutrients  
     104      ! 
     105#if ! defined key_sed_off 
    103106      CALL p4z_che              ! computation of chemical constants 
    104107      CALL p4z_int( kt )        ! computation of various rates for biogeochemistry 
    105108      ! 
    106       IF( ll_sbc ) CALL p4z_sbc( kt )   ! external sources of nutrients  
    107  
    108109      DO jnt = 1, nrdttrc          ! Potential time splitting if requested 
    109110         ! 
     
    149150         END DO 
    150151      END IF 
    151       ! 
    152       IF( lk_sed ) THEN  
     152#endif 
     153      ! 
     154      IF( ln_sediment ) THEN  
    153155         ! 
    154156         CALL sed_model( kt )     !  Main program of Sediment model 
     157         ! 
     158         IF( ln_top_euler ) THEN 
     159            DO jn = jp_pcs0, jp_pcs1 
     160               trn(:,:,:,jn) = trb(:,:,:,jn) 
     161            END DO 
     162         ENDIF 
    155163         ! 
    156164      ENDIF 
     
    352360      IF(lwp)  WRITE(numout,*) 
    353361 
    354       IF( cn_cfg == "orca" .AND. .NOT. lk_c1d ) THEN      ! ORCA configuration (not 1D) ! 
    355          !                                                ! --------------------------- ! 
    356          ! set total alkalinity, phosphate, nitrate & silicate 
    357          zarea          = 1._wp / glob_sum( cvol(:,:,:) ) * 1e6               
    358  
    359          zalksumn = glob_sum( trn(:,:,:,jptal) * cvol(:,:,:)  ) * zarea 
    360          zpo4sumn = glob_sum( trn(:,:,:,jppo4) * cvol(:,:,:)  ) * zarea * po4r 
    361          zno3sumn = glob_sum( trn(:,:,:,jpno3) * cvol(:,:,:)  ) * zarea * rno3 
    362          zsilsumn = glob_sum( trn(:,:,:,jpsil) * cvol(:,:,:)  ) * zarea 
     362      IF( cn_cfg == "ORCA" .OR. cn_cfg == "orca") THEN 
     363         IF( .NOT. lk_c1d ) THEN      ! ORCA configuration (not 1D) ! 
     364            !                                                ! --------------------------- ! 
     365            ! set total alkalinity, phosphate, nitrate & silicate 
     366            zarea          = 1._wp / glob_sum( cvol(:,:,:) ) * 1e6               
     367 
     368            zalksumn = glob_sum( trn(:,:,:,jptal) * cvol(:,:,:)  ) * zarea 
     369            zpo4sumn = glob_sum( trn(:,:,:,jppo4) * cvol(:,:,:)  ) * zarea * po4r 
     370            zno3sumn = glob_sum( trn(:,:,:,jpno3) * cvol(:,:,:)  ) * zarea * rno3 
     371            zsilsumn = glob_sum( trn(:,:,:,jpsil) * cvol(:,:,:)  ) * zarea 
    363372  
    364          IF(lwp) WRITE(numout,*) '       TALKN mean : ', zalksumn 
    365          trn(:,:,:,jptal) = trn(:,:,:,jptal) * alkmean / zalksumn 
    366  
    367          IF(lwp) WRITE(numout,*) '       PO4N  mean : ', zpo4sumn 
    368          trn(:,:,:,jppo4) = trn(:,:,:,jppo4) * po4mean / zpo4sumn 
    369  
    370          IF(lwp) WRITE(numout,*) '       NO3N  mean : ', zno3sumn 
    371          trn(:,:,:,jpno3) = trn(:,:,:,jpno3) * no3mean / zno3sumn 
    372  
    373          IF(lwp) WRITE(numout,*) '       SiO3N mean : ', zsilsumn 
    374          trn(:,:,:,jpsil) = MIN( 400.e-6,trn(:,:,:,jpsil) * silmean / zsilsumn ) 
    375          ! 
    376          ! 
    377          IF( .NOT. ln_top_euler ) THEN 
    378             zalksumb = glob_sum( trb(:,:,:,jptal) * cvol(:,:,:)  ) * zarea 
    379             zpo4sumb = glob_sum( trb(:,:,:,jppo4) * cvol(:,:,:)  ) * zarea * po4r 
    380             zno3sumb = glob_sum( trb(:,:,:,jpno3) * cvol(:,:,:)  ) * zarea * rno3 
    381             zsilsumb = glob_sum( trb(:,:,:,jpsil) * cvol(:,:,:)  ) * zarea 
     373            IF(lwp) WRITE(numout,*) '       TALKN mean : ', zalksumn 
     374            trn(:,:,:,jptal) = trn(:,:,:,jptal) * alkmean / zalksumn 
     375 
     376            IF(lwp) WRITE(numout,*) '       PO4N  mean : ', zpo4sumn 
     377            trn(:,:,:,jppo4) = trn(:,:,:,jppo4) * po4mean / zpo4sumn 
     378 
     379            IF(lwp) WRITE(numout,*) '       NO3N  mean : ', zno3sumn 
     380            trn(:,:,:,jpno3) = trn(:,:,:,jpno3) * no3mean / zno3sumn 
     381 
     382            IF(lwp) WRITE(numout,*) '       SiO3N mean : ', zsilsumn 
     383            trn(:,:,:,jpsil) = MIN( 400.e-6,trn(:,:,:,jpsil) * silmean / zsilsumn ) 
     384            ! 
     385            ! 
     386            IF( .NOT. ln_top_euler ) THEN 
     387               zalksumb = glob_sum( trb(:,:,:,jptal) * cvol(:,:,:)  ) * zarea 
     388               zpo4sumb = glob_sum( trb(:,:,:,jppo4) * cvol(:,:,:)  ) * zarea * po4r 
     389               zno3sumb = glob_sum( trb(:,:,:,jpno3) * cvol(:,:,:)  ) * zarea * rno3 
     390               zsilsumb = glob_sum( trb(:,:,:,jpsil) * cvol(:,:,:)  ) * zarea 
    382391  
    383             IF(lwp) WRITE(numout,*) ' ' 
    384             IF(lwp) WRITE(numout,*) '       TALKB mean : ', zalksumb 
    385             trb(:,:,:,jptal) = trb(:,:,:,jptal) * alkmean / zalksumb 
    386  
    387             IF(lwp) WRITE(numout,*) '       PO4B  mean : ', zpo4sumb 
    388             trb(:,:,:,jppo4) = trb(:,:,:,jppo4) * po4mean / zpo4sumb 
    389  
    390             IF(lwp) WRITE(numout,*) '       NO3B  mean : ', zno3sumb 
    391             trb(:,:,:,jpno3) = trb(:,:,:,jpno3) * no3mean / zno3sumb 
    392  
    393             IF(lwp) WRITE(numout,*) '       SiO3B mean : ', zsilsumb 
    394             trb(:,:,:,jpsil) = MIN( 400.e-6,trb(:,:,:,jpsil) * silmean / zsilsumb ) 
     392               IF(lwp) WRITE(numout,*) ' ' 
     393               IF(lwp) WRITE(numout,*) '       TALKB mean : ', zalksumb 
     394               trb(:,:,:,jptal) = trb(:,:,:,jptal) * alkmean / zalksumb 
     395 
     396               IF(lwp) WRITE(numout,*) '       PO4B  mean : ', zpo4sumb 
     397               trb(:,:,:,jppo4) = trb(:,:,:,jppo4) * po4mean / zpo4sumb 
     398 
     399               IF(lwp) WRITE(numout,*) '       NO3B  mean : ', zno3sumb 
     400               trb(:,:,:,jpno3) = trb(:,:,:,jpno3) * no3mean / zno3sumb 
     401 
     402               IF(lwp) WRITE(numout,*) '       SiO3B mean : ', zsilsumb 
     403               trb(:,:,:,jpsil) = MIN( 400.e-6,trb(:,:,:,jpsil) * silmean / zsilsumb ) 
     404           ENDIF 
    395405        ENDIF 
    396406        ! 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p5zlim.F90

    r10320 r10321  
    123123   !!---------------------------------------------------------------------- 
    124124   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    125    !! $Id: p5zlim.F90 10070 2018-08-28 14:30:54Z nicolasmartin $  
     125   !! $Id$ 
    126126   !! Software governed by the CeCILL license (see ./LICENSE) 
    127127   !!---------------------------------------------------------------------- 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p5zmeso.F90

    • Property svn:executable deleted
    • Property svn:keywords set to Id
    r9950 r10321  
    5151 
    5252   !!---------------------------------------------------------------------- 
    53    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     53   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    5454   !! $Id$ 
    55    !! Software governed by the CeCILL licence (./LICENSE) 
     55   !! Software governed by the CeCILL license (see ./LICENSE) 
    5656   !!---------------------------------------------------------------------- 
    5757 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p5zmicro.F90

    • Property svn:executable deleted
    • Property svn:keywords set to Id
    r9950 r10321  
    1515   USE trc             !  passive tracers common variables  
    1616   USE sms_pisces      !  PISCES Source Minus Sink variables 
    17    USE p5zice          !  Phytoplankton limitation terms 
     17   USE p5zlim          !  Phytoplankton limitation terms 
    1818   USE iom             !  I/O manager 
    1919   USE prtctl_trc      !  print control for debugging 
     
    5151 
    5252   !!---------------------------------------------------------------------- 
    53    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     53   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    5454   !! $Id$ 
    55    !! Software governed by the CeCILL licence (./LICENSE) 
     55   !! Software governed by the CeCILL license (see ./LICENSE) 
    5656   !!---------------------------------------------------------------------- 
    5757 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p5zmort.F90

    • Property svn:executable deleted
    • Property svn:keywords set to Id
    r9950 r10321  
    1414   USE trc             !  passive tracers common variables  
    1515   USE sms_pisces      !  PISCES Source Minus Sink variables 
    16    USE p5zice          !  Phytoplankton limitation terms 
     16   USE p5zlim          !  Phytoplankton limitation terms 
    1717   USE prtctl_trc      !  print control for debugging 
    1818 
     
    3333 
    3434   !!---------------------------------------------------------------------- 
    35    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     35   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    3636   !! $Id$ 
    37    !! Software governed by the CeCILL licence (./LICENSE) 
     37   !! Software governed by the CeCILL license (see ./LICENSE) 
    3838   !!---------------------------------------------------------------------- 
    3939 
  • NEMO/branches/UKMO/NEMO4_beta_mirror/src/TOP/PISCES/P4Z/p5zprod.F90

    • Property svn:executable deleted
    • Property svn:keywords set to Id
    r9950 r10321  
    1616   USE trc             !  passive tracers common variables  
    1717   USE sms_pisces      !  PISCES Source Minus Sink variables 
    18    USE p5zice          !  Co-limitations of differents nutrients 
     18   USE p5zlim          !  Co-limitations of differents nutrients 
    1919   USE prtctl_trc      !  print control for debugging 
    2020   USE iom             !  I/O manager 
     
    5353 
    5454   !!---------------------------------------------------------------------- 
    55    !! NEMO/TOP 3.3 , NEMO Consortium (2018) 
     55   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    5656   !! $Id$ 
    57    !! Software governed by the CeCILL licence     (./LICENSE) 
     57   !! Software governed by the CeCILL license (see ./LICENSE) 
    5858   !!---------------------------------------------------------------------- 
    5959CONTAINS 
Note: See TracChangeset for help on using the changeset viewer.