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 10416 for NEMO/trunk/src/TOP/PISCES/P4Z/p4zligand.F90 – NEMO

Ignore:
Timestamp:
2018-12-19T12:45:43+01:00 (5 years ago)
Author:
aumont
Message:

Remove LFe from standard PISCES code

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/TOP/PISCES/P4Z/p4zligand.F90

    r10362 r10416  
    2525   REAL(wp), PUBLIC ::  rlig     !: Remin ligand production 
    2626   REAL(wp), PUBLIC ::  prlgw    !: Photochemical of weak ligand 
    27    REAL(wp), PUBLIC ::  rfep     !: Dissolution rate of FeP 
    2827 
    2928   !!---------------------------------------------------------------------- 
     
    4342      ! 
    4443      INTEGER  ::   ji, jj, jk 
    45       REAL(wp) ::   zlgwp, zlgwpr, zlgwr, zlablgw, zrfepa, zfepr 
     44      REAL(wp) ::   zlgwp, zlgwpr, zlgwr, zlablgw 
    4645      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zligrem, zligpr, zrligprod 
    4746      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   zw3d 
     
    7170               zligpr(ji,jj,jk)    = zlgwpr 
    7271               zrligprod(ji,jj,jk) = zlgwp 
    73                ! 
    74                ! ---------------------------------------------------------- 
    75                ! Dissolution of nanoparticle Fe 
    76                ! ---------------------------------------------------------- 
    77                ! dissolution rate is maximal in the presence of light and  
    78                ! lower in the aphotici zone 
    79                ! ! 25 Wm-2 constant 
    80                zrfepa = rfep * ( 1. - EXP( -1. * etot(ji,jj,jk) / 25. ) ) * (1.- fr_i(ji,jj)) 
    81                zrfepa = MAX( (zrfepa / 10.0), zrfepa ) ! min of 10 days lifetime 
    82                zfepr  = rfep * xstep * trb(ji,jj,jk,jpfep) 
    83                tra(ji,jj,jk,jpfep) = tra(ji,jj,jk,jpfep) - zfepr 
    84                tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zfepr 
    8572               ! 
    8673            END DO 
     
    130117      INTEGER ::   ios   ! Local integer  
    131118      ! 
    132       NAMELIST/nampislig/ rlgw, prlgw, rlgs, rfep, rlig 
     119      NAMELIST/nampislig/ rlgw, prlgw, rlgs, rlig 
    133120      !!---------------------------------------------------------------------- 
    134121      ! 
     
    148135      IF(lwp) THEN                         ! control print 
    149136         WRITE(numout,*) '   Namelist : nampislig' 
    150          WRITE(numout,*) '      Dissolution rate of FeP                      rfep  =', rfep 
    151137         WRITE(numout,*) '      Lifetime (years) of weak ligands             rlgw  =', rlgw 
    152138         WRITE(numout,*) '      Remin ligand production per unit C           rlig  =', rlig 
Note: See TracChangeset for help on using the changeset viewer.