Index: /branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/iodef.xml
===================================================================
--- /branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/iodef.xml	(revision 3462)
+++ /branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/iodef.xml	(revision 3463)
@@ -221,5 +221,5 @@
      <group id="ptrc_T" axis_ref="deptht" grid_ref="grid_T">
        <field id="DIC"      description="Dissolved inorganic Concentration"        unit="mmol/m3" />
-       <field id="Alkalini" description="Total Alkalinity Concentration"           unit="mmol/m3"    />
+       <field id="Alkalini" description="Total Alkalinity Concentration"           unit="mmol/m3" />
        <field id="O2"       description="Oxygen Concentration"                     unit="mmol/m3" />
        <field id="CaCO3"    description="Calcite Concentration"                    unit="mmol/m3" />
@@ -240,6 +240,6 @@
        <field id="GSi"      description="Sinking biogenic Silicate Concentration"  unit="mmol/m3" />
        <field id="NFe"      description="Nano iron Concentration"                  unit="mmol/m3" />
-       <field id="NCHL"     description="Nano chlorophyl Concentration"            unit="mg/m3" />
-       <field id="DCHL"     description="Diatoms chlorophyl Concentration"         unit="mg/m3" />
+       <field id="NCHL"     description="Nano chlorophyl Concentration"            unit="mg/m3"   />
+       <field id="DCHL"     description="Diatoms chlorophyl Concentration"         unit="mg/m3"   />
        <field id="NO3"      description="Nitrates Concentration"                   unit="mmol/m3" />
        <field id="NH4"      description="Ammonium Concentration"                   unit="mmol/m3" />
@@ -249,41 +249,57 @@
 
      <group id="diad_T" axis_ref="none" grid_ref="grid_T">
-       <field id="PH"          description="PH"                                       unit="-"          axis_ref="deptht" />
-       <field id="CO3"         description="Bicarbonates"                             unit="mol/m3"     axis_ref="deptht" />
-       <field id="CO3sat"      description="CO3 saturation"                           unit="mol/m3"     axis_ref="deptht" />
-       <field id="PAR"         description="Photosynthetically Available Radiation"   unit="W/m2"       axis_ref="deptht" />
-       <field id="PPPHY"       description="Primary production of nanophyto"          unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="PPPHY2"      description="Primary production of diatoms"            unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="PPNEWN"      description="New Primary production of nanophyto"      unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="PPNEWD"      description="New Primary production of diatoms"        unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="PBSi"        description="Primary production of Si diatoms"         unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="PFeN"        description="Primary production of nano iron"          unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="PFeD"        description="Primary production of diatoms iron"       unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="PCAL"        description="Calcite production"                       unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="DCAL"        description="Calcite dissolution"                      unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="GRAZ"        description="Grazing by zooplankton"                   unit="mol/m3/s"   axis_ref="deptht" />
-       <field id="Mumax"       description="Maximum growth rate"                      unit="s-1"        axis_ref="deptht" />
-       <field id="MuN"         description="Realized growth rate for nanophyto"       unit="s-1"        axis_ref="deptht" />
-       <field id="MuD"         description="Realized growth rate for diatomes"        unit="s-1"        axis_ref="deptht" />
-       <field id="MuNlight"    description="Light limited growth rate for nanophyto"  unit="s-1"        axis_ref="deptht" />
-       <field id="MuDlight"    description="Light limited growth rate for diatomes"   unit="s-1"        axis_ref="deptht" />
-       <field id="LNnut"       description="Nutrient limitation term in Nanophyto"    unit="-"          axis_ref="deptht" />
-       <field id="LDnut"       description="Nutrient limitation term in Diatoms"      unit="-"          axis_ref="deptht" />
-       <field id="LNFe"        description="Iron limitation term in Nanophyto"        unit="-"          axis_ref="deptht" />
-       <field id="LDFe"        description="Iron limitation term in Diatoms"          unit="-"          axis_ref="deptht" />
-       <field id="LNlight"     description="Light limitation term in Nanophyto"       unit="-"          axis_ref="deptht" />
-       <field id="LDlight"     description="Light limitation term in Diatoms"         unit="-"          axis_ref="deptht" />
-       <field id="Nfix"        description="Nitrogen fixation at surface"             unit="mol/m2/s"                     />
-       <field id="EPC100"      description="Export of carbon particles at 100 m"      unit="mol/m2/s"                     />
-       <field id="EPFE100"     description="Export of biogenic iron at 100 m"         unit="mol/m2/s"                     />
-       <field id="EPSI100"     description="Export of Silicate at 100 m"              unit="mol/m2/s"                     />
-       <field id="EPCAL100"    description="Export of Calcite at 100 m"               unit="mol/m2/s"                     />
-       <field id="Cflx"        description="DIC flux"                                 unit="mol/m2/s"                     />
-       <field id="Oflx"        description="Oxygen flux"                              unit="mol/m2/s"                     />
-       <field id="Kg"          description="Gas transfer"                             unit="mol/m2/s/uatm"                />
-       <field id="Dpco2"       description="Delta CO2"                                unit="uatm"                         />
-       <field id="Dpo2"        description="Delta O2"                                 unit="uatm"                         />
-       <field id="Heup"        description="Euphotic layer depth"                     unit="m"                            />
-       <field id="Irondep"     description="Iron deposition"                          unit="mol/m2/s"                     />
+       <field id="PH"          description="PH"                                      unit="-"          axis_ref="deptht" />
+       <field id="CO3"         description="Bicarbonates"                            unit="mol/m3"     axis_ref="deptht" />
+       <field id="CO3sat"      description="CO3 saturation"                          unit="mol/m3"     axis_ref="deptht" />
+       <field id="PAR"         description="Photosynthetically Available Radiation"  unit="W/m2"       axis_ref="deptht" />
+       <field id="PPPHY"       description="Primary production of nanophyto"         unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="PPPHY2"      description="Primary production of diatoms"           unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="PPNEWN"      description="New Primary production of nanophyto"     unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="PPNEWD"      description="New Primary production of diatoms"       unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="PBSi"        description="Primary production of Si diatoms"        unit="molSi/m3/s" axis_ref="deptht" />
+       <field id="PFeN"        description="Primary production of nano iron"         unit="molFe/m3/s" axis_ref="deptht" />
+       <field id="PFeD"        description="Primary production of diatoms iron"      unit="molFe/m3/s" axis_ref="deptht" />
+       <field id="xfracal"     description="Calcifying fraction"                     unit="-"          axis_ref="deptht" />
+       <field id="PCAL"        description="Calcite production"                      unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="DCAL"        description="Calcite dissolution"                     unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="GRAZ1"       description="Grazing by microzooplankton"             unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="GRAZ2"       description="Grazing by mesozooplankton"              unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="REMIN"       description="Oxic remineralization of OM"             unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="DENIT"       description="Anoxic remineralization of OM"           unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="Nfix"        description="Nitrogen fixation"                       unit="molN/m3/s"  axis_ref="deptht" />
+       <field id="Mumax"       description="Maximum growth rate"                     unit="s-1"        axis_ref="deptht" />
+       <field id="MuN"         description="Realized growth rate for nanophyto"      unit="s-1"        axis_ref="deptht" />
+       <field id="MuD"         description="Realized growth rate for diatomes"       unit="s-1"        axis_ref="deptht" />
+       <field id="LNnut"       description="Nutrient limitation term in Nanophyto"   unit="-"          axis_ref="deptht" />
+       <field id="LDnut"       description="Nutrient limitation term in Diatoms"     unit="-"          axis_ref="deptht" />
+       <field id="LNFe"        description="Iron limitation term in Nanophyto"       unit="-"          axis_ref="deptht" />
+       <field id="LDFe"        description="Iron limitation term in Diatoms"         unit="-"          axis_ref="deptht" />
+       <field id="LNlight"     description="Light limitation term in Nanophyto"      unit="-"          axis_ref="deptht" />
+       <field id="LDlight"     description="Light limitation term in Diatoms"        unit="-"          axis_ref="deptht" />
+       <field id="Fe2"         description="Iron II concentration"                   unit="nmol/L"     axis_ref="deptht" />
+       <field id="Fe3"         description="Iron III concentration"                  unit="nmol/L"     axis_ref="deptht" />
+       <field id="FeL1"        description="Complexed Iron concentration with L1"    unit="nmol/L"     axis_ref="deptht" />
+       <field id="FeL2"        description="Complexed Iron concentration with L2"    unit="nmol/L"     axis_ref="deptht" />
+       <field id="FeP"         description="Precipitated Iron III"                   unit="nmol/L"     axis_ref="deptht" />
+       <field id="TL1"         description="Total L1 concentration"                  unit="nmol/L"     axis_ref="deptht" />
+       <field id="TL2"         description="Total L2 concentration"                  unit="nmol/L"     axis_ref="deptht" />
+       <field id="pdust"       description="dust concentration"                      unit="g/L"                          />
+       <field id="Totlig"      description="Total ligand concentation"               unit="nmol/L"     axis_ref="deptht" />
+       <field id="Biron"       description="Bioavailable iron"                       unit="nmol/L"     axis_ref="deptht" />
+       <field id="Sdenit"      description="Nitrate reduction in the sediments"      unit="molN/m2/s"                    />
+       <field id="Ironice"     description="Iron input/uptake due to sea ice"        unit="molFe/m2/s"                   />
+       <field id="HYDR"        description="Iron input from hydrothemal vents"       unit="molFe/m2/s"                   />
+       <field id="EPC100"      description="Export of carbon particles at 100 m"     unit="mol/m2/s"                     />
+       <field id="EPFE100"     description="Export of biogenic iron at 100 m"        unit="mol/m2/s"                     />
+       <field id="EPSI100"     description="Export of Silicate at 100 m"             unit="mol/m2/s"                     />
+       <field id="EPCAL100"    description="Export of Calcite at 100 m"              unit="mol/m2/s"                     />
+       <field id="Cflx"        description="DIC flux"                                unit="mol/m2/s"                     />
+       <field id="Oflx"        description="Oxygen flux"                             unit="mol/m2/s"                     />
+       <field id="Kg"          description="Gas transfer"                            unit="mol/m2/s/uatm"                />
+       <field id="Dpco2"       description="Delta CO2"                               unit="uatm"                         />
+       <field id="Dpo2"        description="Delta O2"                                unit="uatm"                         />
+       <field id="Heup"        description="Euphotic layer depth"                    unit="m"                            />
+       <field id="Irondep"     description="Iron deposition from dust"               unit="mol/m2/s"                     />
+       <field id="Ironsed"     description="Iron deposition from sediment"           unit="mol/m2/s"  axis_ref="deptht"  />
      </group>
        
@@ -507,5 +523,5 @@
           <field ref="CO3"      />
           <field ref="CO3sat"   />
-          <field ref="PAR"      />
+          <field ref="PAR"   />
           <field ref="PPPHY"    />
           <field ref="PPPHY2"   />
@@ -515,18 +531,9 @@
           <field ref="PFeN"     />
           <field ref="PFeD"     />
+          <field ref="xfracal"  />
           <field ref="PCAL"     />
           <field ref="DCAL"     />
-          <field ref="GRAZ"     />
-          <field ref="Mumax"    />
-          <field ref="MuN"      />
-          <field ref="MuD"      />
-          <field ref="MuNlight" />
-          <field ref="MuDlight" />
-          <field ref="LNnut"    />
-          <field ref="LDnut"    />
-          <field ref="LNFe"     />
-          <field ref="LDFe"     />
-          <field ref="LNlight"  />
-          <field ref="LDlight"  />
+          <field ref="GRAZ1"     />
+          <field ref="GRAZ2"     />
           <field ref="EPC100"   />
           <field ref="EPFE100"  />
@@ -540,5 +547,25 @@
           <field ref="Heup"     />
           <field ref="Irondep"  />
+          <field ref="Ironsed"  />
+          <field ref="Ironice"  />
           <field ref="Nfix"     />
+          <field ref="MuN"      />
+          <field ref="MuD"      />
+          <field ref="LNnut"      />
+          <field ref="LDnut"      />
+          <field ref="LNFe"      />
+          <field ref="LDFe"      />
+          <field ref="LNlight"      />
+          <field ref="LDlight"      />
+          <field ref="pdust"      />
+          <field ref="Fe2"      />
+          <field ref="Fe3"      />
+          <field ref="FeL1"      />
+          <field ref="FeL2"      />
+          <field ref="FeP"      />
+          <field ref="TL1"      />
+          <field ref="TL2"      />
+          <field ref="Sdenit"      />
+          <field ref="Totlig" />
        </file>
 
Index: /branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_pisces
===================================================================
--- /branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_pisces	(revision 3462)
+++ /branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_pisces	(revision 3463)
@@ -40,26 +40,40 @@
    ferat3     =  10.E-6   ! Fe/C in zooplankton 
    wsbio2     =  30.      ! Big particles sinking speed
+   niter1max  =  1        ! Maximum number of iterations for POC
+   niter2max  =  1        ! Maximum number of iterations for GOC
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
 &nampislim     !   parameters for nutrient limitations
 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
-   conc0      =  1.e-6    ! Phosphate half saturation
-   conc1      =  8E-6     ! Phosphate half saturation for diatoms
-   conc2      =  1E-9     ! Iron half saturation for phyto
-   conc2m     =  3E-9     ! Max iron half saturation for phyto
-   conc3      =  3E-9     ! Iron half saturation for diatoms
-   conc3m     =  8E-9     ! Maxi iron half saturation for diatoms
+   concnno3   =  1.e-6    ! Nitrate half saturation of nanophytoplankton
+   concdno3   =  3.E-6   ! Phosphate half saturation for diatoms
+   concnnh4   =  1.E-7    ! NH4 half saturation for phyto
+   concdnh4   =  3.E-7   ! NH4 half saturation for diatoms
+   concnfer   =  1.E-9     ! Iron half saturation for phyto
+   concdfer   =  3.E-9   ! Iron half saturation for diatoms
+   concbfe    =  1.E-11   ! Half-saturation for Fe limitation of Bacteria
+   concbnh4   =  2.5E-8   ! NH4 half saturation for phyto
+   concbno3   =  2.5E-7   ! Phosphate half saturation for diatoms
    xsizedia   =  1.E-6    ! Minimum size criteria for diatoms
    xsizephy   =  1.E-6    ! Minimum size criteria for phyto
-   concnnh4   =  1.E-7    ! NH4 half saturation for phyto
-   concdnh4   =  8.E-7    ! NH4 half saturation for diatoms
+   xsizern    =  3.0      ! Size ratio for nanophytoplankton
+   xsizerd    =  3.0      ! Size ratio for diatoms
    xksi1      =  2.E-6    ! half saturation constant for Si uptake
-   xksi2      =  3.33E-6  ! half saturation constant for Si/C
+   xksi2      =  20E-6  ! half saturation constant for Si/C
    xkdoc      =  417.E-6  ! half-saturation constant of DOC remineralization
-   concfebac  =  1.E-11   ! Half-saturation for Fe limitation of Bacteria
    qnfelim    =  7.E-6    ! Optimal quota of phyto
    qdfelim    =  7.E-6    ! Optimal quota of diatoms
-   caco3r     =  0.16     ! mean rain ratio
-/
+   caco3r     =  0.3      ! mean rain ratio
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisopt     !   parameters for optics
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!              !  file name        ! frequency (hours) ! variable   ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation !
+!              !                   !  (if <0  months)  !   name     !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !
+   sn_par     = 'par.orca'       ,     24            , 'fr_par'     ,  .true.      , .true. ,   'yearly'  , ''       , ''
+   cn_dir      = './'      !  root directory for the location of the dynamical files
+   ln_varpar   =  .true.   ! boolean for PAR variable
+   parlux      =  0.43      ! Fraction of shortwave as PAR
+/ 
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
 &nampisprod     !   parameters for phytoplankton growth
@@ -69,18 +83,19 @@
    excret     =  0.05     ! excretion ratio of phytoplankton
    excret2    =  0.05     ! excretion ratio of diatoms
-   ln_newprod =  .true.  ! Enable new parame. of production (T/F) 
+   ln_newprod =  .true.   ! Enable new parame. of production (T/F) 
    bresp      =  0.00333  ! Basal respiration rate
    chlcnm     =  0.033    ! Minimum Chl/C in nanophytoplankton
    chlcdm     =  0.05     ! Minimum Chl/C in diatoms
-   chlcmin    =  0.0033   ! Maximum Chl/c in phytoplankton
+   chlcmin    =  0.004    ! Maximum Chl/c in phytoplankton
    fecnm      =  40E-6    ! Maximum Fe/C in nanophytoplankton
    fecdm      =  40E-6    ! Minimum Fe/C in diatoms
-   grosip     =  0.151    ! mean Si/C ratio
+   grosip     =  0.159    ! mean Si/C ratio
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
 &nampismort     !   parameters for phytoplankton sinks
 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
-   wchl       =  0.001    ! quadratic mortality of phytoplankton
-   wchld      =  0.02     ! maximum quadratic mortality of diatoms
+   wchl       =  0.01    ! quadratic mortality of phytoplankton
+   wchld      =  0.01     ! maximum quadratic mortality of diatoms
+   wchldm     =  0.03     ! maximum quadratic mortality of diatoms
    mprat      =  0.01     ! phytoplankton mortality rate
    mprat2     =  0.01     ! Diatoms mortality rate
@@ -90,6 +105,6 @@
 &nampismes     !   parameters for mesozooplankton
 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
-   part2      =  0.75    ! part of calcite not dissolved in mesozoo guts
-   grazrat2   =  0.7      ! maximal mesozoo grazing rate
+   part2      =  0.75     ! part of calcite not dissolved in mesozoo guts
+   grazrat2   =  0.75     ! maximal mesozoo grazing rate
    resrat2    =  0.005    ! exsudation rate of mesozooplankton
    mzrat2     =  0.03     ! mesozooplankton mortality rate
@@ -102,5 +117,5 @@
    xthresh2phy = 1E-8     ! nanophyto feeding threshold for mesozooplankton 
    xthresh2poc = 1E-8     ! poc feeding threshold for mesozooplankton 
-   xthresh2   =  2E-7    ! Food threshold for grazing
+   xthresh2   =  3E-7     ! Food threshold for grazing
    xkgraz2    =  20.E-6   ! half sturation constant for meso grazing
    epsher2    =  0.3      ! Efficicency of Mesozoo growth
@@ -115,5 +130,5 @@
    grazrat    =  3.0      ! maximal zoo grazing rate
    resrat     =  0.03     ! exsudation rate of zooplankton
-   mzrat      =  0.001    ! zooplankton mortality rate
+   mzrat      =  0.004    ! zooplankton mortality rate
    xpref2c    =  0.1      ! Microzoo preference for POM
    xpref2p    =  1.       ! Microzoo preference for Nanophyto
@@ -122,5 +137,5 @@
    xthreshphy =  1.E-8    ! Nanophyto feeding threshold for microzooplankton 
    xthreshpoc =  1.E-8    ! POC feeding threshold for microzooplankton 
-   xthresh    =  2.E-7    ! Food threshold for feeding
+   xthresh    =  3.E-7    ! Food threshold for feeding
    xkgraz     =  20.E-6   ! half sturation constant for grazing
    epsher     =  0.3      ! Efficiency of microzoo growth
@@ -129,15 +144,22 @@
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisfer     !   parameters for iron chemistry
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+   ln_fechem =  .false.   ! complex iron chemistry ( T/F )
+   ln_ligvar =  .true.   ! variable ligand concentration
+   xlam1     =  0.005     ! scavenging rate of Iron
+   xlamdust  =  100.0     ! Scavenging rate of dust
+   ligand    =  0.6E-9    ! Ligands concentration 
+/  
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
 &nampisrem     !   parameters for remineralization
 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
-   xremik    =  0.25      ! remineralization rate of DOC
+   xremik    =  0.35      ! remineralization rate of DOC
    xremip    =  0.025     ! remineralisation rate of POC
    nitrif    =  0.05      ! NH4 nitrification rate
    xsirem    =  0.003     ! remineralization rate of Si
-   xsiremlab =  0.025     ! fast remineralization rate of Si
-   xsilab    =  0.31      ! Fraction of labile biogenic silica
-   xlam1     =  0.005     ! scavenging rate of Iron
+   xsiremlab =  0.03      ! fast remineralization rate of Si
+   xsilab    =  0.5       ! Fraction of labile biogenic silica
    oxymin    =  1.E-6     ! Half-saturation constant for anoxia
-   ligand    =  0.6E-9    ! Ligands concentration
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
@@ -148,25 +170,37 @@
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
-&nampissed     !   parameters for inputs deposition
+&nampissbc     !   parameters for inputs deposition
 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
 !              !  file name        ! frequency (hours) ! variable   ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation !
 !              !                   !  (if <0  months)  !   name     !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !
    sn_dust     = 'dust.orca'       ,     -1            , 'dust'     ,  .true.      , .true. ,   'yearly'  , ''       , ''
-   sn_riverdic = 'river.orca'      ,    -12            , 'riverdic' ,  .false.     , .true. ,   'yearly'  , ''       , ''
-   sn_riverdoc = 'river.orca'      ,    -12            , 'riverdoc' ,  .false.     , .true. ,   'yearly'  , ''       , ''
+   sn_solub    = 'solubility.orca' ,    -12            , 'solubility1' ,  .false.      , .true. ,   'yearly'  , ''       , ''
+   sn_riverdic = 'river.orca'      ,    120            , 'riverdic' ,  .true.     , .true. ,   'yearly'  , ''       , ''
+   sn_riverdoc = 'river.orca'      ,    120            , 'riverdoc' ,  .true.     , .true. ,   'yearly'  , ''       , ''
+   sn_riverdin = 'river.orca'      ,    120            , 'riverdin' ,  .true.     , .true. ,   'yearly'  , ''       , ''
+   sn_riverdon = 'river.orca'      ,    120            , 'riverdon' ,  .true.     , .true. ,   'yearly'  , ''       , ''
+   sn_riverdip = 'river.orca'      ,    120            , 'riverdip' ,  .true.     , .true. ,   'yearly'  , ''       , ''
+   sn_riverdop = 'river.orca'      ,    120            , 'riverdop' ,  .true.     , .true. ,   'yearly'  , ''       , ''
+   sn_riverdsi = 'river.orca'      ,    120            , 'riverdsi' ,  .true.     , .true. ,   'yearly'  , ''       , ''
    sn_ndepo    = 'ndeposition.orca',    -12            , 'ndep'     ,  .false.     , .true. ,   'yearly'  , ''       , ''
    sn_ironsed  = 'bathy.orca'      ,    -12            , 'bathy'    ,  .false.     , .true. ,   'yearly'  , ''       , ''
+   sn_hydrofe  = 'hydrofe'         ,    -12            , 'epsdb'    ,  .false.     , .true. ,   'yearly'  , ''       , ''
 !
    cn_dir      = './'      !  root directory for the location of the dynamical files
    ln_dust     =  .true.   ! boolean for dust input from the atmosphere
-   ln_river    =  .false.   ! boolean for river input of nutrients
+   ln_solub    =  .true.   ! boolean for variable solubility of atm. Iron
+   ln_river    =  .true.   ! boolean for river input of nutrients
    ln_ndepo    =  .true.   ! boolean for atmospheric deposition of N
    ln_ironsed  =  .true.   ! boolean for Fe input from sediments
+   ln_ironice  =  .true.   ! boolean for Fe input from sea ice
+   ln_hydrofe  =  .false.   ! boolean for from hydrothermal vents
    sedfeinput  =  1E-9     ! Coastal release of Iron
    dustsolub   =  0.02     ! Solubility of the dust
-   wdust       =  2.0      ! Dust sinking speed
+   wdust       =  2.0      ! Dust sinking speed 
+   icefeinput  =  10E-9    ! Iron concentration in sea ice
    nitrfix     =  1E-7     ! Nitrogen fixation rate
    diazolight  =  50.      ! Diazotrophs sensitivity to light (W/m2)
    concfediaz  =  1.E-10   ! Diazotrophs half-saturation Cste for Iron
+   hratio      =  9.E-5    ! Fe to 3He ratio assumed for vent iron supply 
 /
 !'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
@@ -175,4 +209,5 @@
    xkr_eta      = 1.17    ! Sinking  exponent
    xkr_zeta     = 2.28    ! N content exponent
+   xkr_ncontent = 5.7E-6  ! N content factor
    xkr_mass_min = 0.0002  ! Minimum mass for Aggregates
    xkr_mass_max = 1.      ! Maximum mass for Aggregates
Index: /branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/iodef.xml
===================================================================
--- /branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/iodef.xml	(revision 3462)
+++ /branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/iodef.xml	(revision 3463)
@@ -19,5 +19,5 @@
      <group id="ptrc_T" axis_ref="deptht" grid_ref="grid_T">
        <field id="DIC"      description="Dissolved inorganic Concentration"        unit="mmol/m3" />
-       <field id="Alkalini" description="Total Alkalinity Concentration"           unit="mmol/m3"    />
+       <field id="Alkalini" description="Total Alkalinity Concentration"           unit="mmol/m3" />
        <field id="O2"       description="Oxygen Concentration"                     unit="mmol/m3" />
        <field id="CaCO3"    description="Calcite Concentration"                    unit="mmol/m3" />
@@ -38,6 +38,6 @@
        <field id="GSi"      description="Sinking biogenic Silicate Concentration"  unit="mmol/m3" />
        <field id="NFe"      description="Nano iron Concentration"                  unit="mmol/m3" />
-       <field id="NCHL"     description="Nano chlorophyl Concentration"            unit="mg/m3" />
-       <field id="DCHL"     description="Diatoms chlorophyl Concentration"         unit="mg/m3" />
+       <field id="NCHL"     description="Nano chlorophyl Concentration"            unit="mg/m3"   />
+       <field id="DCHL"     description="Diatoms chlorophyl Concentration"         unit="mg/m3"   />
        <field id="NO3"      description="Nitrates Concentration"                   unit="mmol/m3" />
        <field id="NH4"      description="Ammonium Concentration"                   unit="mmol/m3" />
@@ -47,57 +47,57 @@
 
      <group id="diad_T" axis_ref="none" grid_ref="grid_T">
-       <field id="PH"          description="PH"                                       unit="-"          axis_ref="deptht" />
-       <field id="CO3"         description="Bicarbonates"                             unit="mol/m3"     axis_ref="deptht" />
-       <field id="CO3sat"      description="CO3 saturation"                           unit="mol/m3"     axis_ref="deptht" />
-       <field id="PAR"         description="Photosynthetically Available Radiation"   unit="W/m2"       axis_ref="deptht" />
-       <field id="PPPHY"       description="Primary production of nanophyto"          unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="PPPHY2"      description="Primary production of diatoms"            unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="PPNEWN"      description="New Primary production of nanophyto"      unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="PPNEWD"      description="New Primary production of diatoms"        unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="PBSi"        description="Primary production of Si diatoms"         unit="molSi/m3/s"   axis_ref="deptht" />
-       <field id="PFeN"        description="Primary production of nano iron"          unit="molFe/m3/s"   axis_ref="deptht" />
-       <field id="PFeD"        description="Primary production of diatoms iron"       unit="molFe/m3/s"   axis_ref="deptht" />
-       <field id="xfracal"     description="Calcifying fraction"                      unit="-"          axis_ref="deptht" />
-       <field id="PCAL"        description="Calcite production"                       unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="DCAL"        description="Calcite dissolution"                      unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="GRAZ1"       description="Grazing by microzooplankton"              unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="GRAZ2"       description="Grazing by mesozooplankton"              unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="REMIN"       description="Oxic remineralization of OM"             unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="DENIT"       description="Anoxic remineralization of OM"           unit="molC/m3/s"   axis_ref="deptht" />
-       <field id="Nfix"        description="Nitrogen fixation"                         unit="molN/m3/s"   axis_ref="deptht" />
-       <field id="Mumax"       description="Maximum growth rate"                      unit="s-1"        axis_ref="deptht" />
-       <field id="MuN"         description="Realized growth rate for nanophyto"       unit="s-1"        axis_ref="deptht" />
-       <field id="MuD"         description="Realized growth rate for diatomes"        unit="s-1"        axis_ref="deptht" />
-       <field id="LNnut"       description="Nutrient limitation term in Nanophyto"    unit="-"          axis_ref="deptht" />
-       <field id="LDnut"       description="Nutrient limitation term in Diatoms"      unit="-"          axis_ref="deptht" />
-       <field id="LNFe"        description="Iron limitation term in Nanophyto"        unit="-"          axis_ref="deptht" />
-       <field id="LDFe"        description="Iron limitation term in Diatoms"          unit="-"          axis_ref="deptht" />
-       <field id="LNlight"     description="Light limitation term in Nanophyto"       unit="-"          axis_ref="deptht" />
-       <field id="LDlight"     description="Light limitation term in Diatoms"         unit="-"          axis_ref="deptht" />
-       <field id="Fe2"        description="Iron II concentration"                      unit="nmol/L"   axis_ref="deptht" />
-       <field id="Fe3"        description="Iron III concentration"                      unit="nmol/L"   axis_ref="deptht" />
-       <field id="FeL1"        description="Complexed Iron concentration with L1"       unit="nmol/L"   axis_ref="deptht" />
-       <field id="FeL2"        description="Complexed Iron concentration with L2"       unit="nmol/L"   axis_ref="deptht" />
-       <field id="FeP"        description="Precipitated Iron III"                      unit="nmol/L"   axis_ref="deptht" />
-       <field id="TL1"        description="Total L1 concentration"                      unit="nmol/L"   axis_ref="deptht" />
-       <field id="TL2"        description="Total L2 concentration"                      unit="nmol/L"   axis_ref="deptht" />
-       <field id="pdust"        description="dust concentration"                      unit="g/L"    />
-       <field id="Totlig"        description="Total ligand concentation"                unit="nmol/L"   axis_ref="deptht" />
-       <field id="Biron"        description="Bioavailable iron"                       unit="nmol/L"   axis_ref="deptht" />
-       <field id="Sdenit"        description="Nitrate reduction in the sediments"     unit="molN/m2/s"      />
-       <field id="Ironice"     description="Iron input/uptake due to sea ice"             unit="molFe/m2/s"      />
-       <field id="HYDR"       description="Iron input from hydrothemal vents"          unit="molFe/m2/s"      />
-       <field id="EPC100"      description="Export of carbon particles at 100 m"      unit="mol/m2/s"                     />
-       <field id="EPFE100"     description="Export of biogenic iron at 100 m"         unit="mol/m2/s"                     />
-       <field id="EPSI100"     description="Export of Silicate at 100 m"              unit="mol/m2/s"                     />
-       <field id="EPCAL100"    description="Export of Calcite at 100 m"               unit="mol/m2/s"                     />
-       <field id="Cflx"        description="DIC flux"                                 unit="mol/m2/s"                     />
-       <field id="Oflx"        description="Oxygen flux"                              unit="mol/m2/s"                     />
-       <field id="Kg"          description="Gas transfer"                             unit="mol/m2/s/uatm"                />
-       <field id="Dpco2"       description="Delta CO2"                                unit="uatm"                         />
-       <field id="Dpo2"        description="Delta O2"                                 unit="uatm"                         />
-       <field id="Heup"        description="Euphotic layer depth"                     unit="m"                            />
-       <field id="Irondep"     description="Iron deposition from dust"                 unit="mol/m2/s"                     />
-       <field id="Ironsed"     description="Iron deposition from sediment"             unit="mol/m2/s"  axis_ref="deptht"  />
+       <field id="PH"          description="PH"                                      unit="-"          axis_ref="deptht" />
+       <field id="CO3"         description="Bicarbonates"                            unit="mol/m3"     axis_ref="deptht" />
+       <field id="CO3sat"      description="CO3 saturation"                          unit="mol/m3"     axis_ref="deptht" />
+       <field id="PAR"         description="Photosynthetically Available Radiation"  unit="W/m2"       axis_ref="deptht" />
+       <field id="PPPHY"       description="Primary production of nanophyto"         unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="PPPHY2"      description="Primary production of diatoms"           unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="PPNEWN"      description="New Primary production of nanophyto"     unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="PPNEWD"      description="New Primary production of diatoms"       unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="PBSi"        description="Primary production of Si diatoms"        unit="molSi/m3/s" axis_ref="deptht" />
+       <field id="PFeN"        description="Primary production of nano iron"         unit="molFe/m3/s" axis_ref="deptht" />
+       <field id="PFeD"        description="Primary production of diatoms iron"      unit="molFe/m3/s" axis_ref="deptht" />
+       <field id="xfracal"     description="Calcifying fraction"                     unit="-"          axis_ref="deptht" />
+       <field id="PCAL"        description="Calcite production"                      unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="DCAL"        description="Calcite dissolution"                     unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="GRAZ1"       description="Grazing by microzooplankton"             unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="GRAZ2"       description="Grazing by mesozooplankton"              unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="REMIN"       description="Oxic remineralization of OM"             unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="DENIT"       description="Anoxic remineralization of OM"           unit="molC/m3/s"  axis_ref="deptht" />
+       <field id="Nfix"        description="Nitrogen fixation"                       unit="molN/m3/s"  axis_ref="deptht" />
+       <field id="Mumax"       description="Maximum growth rate"                     unit="s-1"        axis_ref="deptht" />
+       <field id="MuN"         description="Realized growth rate for nanophyto"      unit="s-1"        axis_ref="deptht" />
+       <field id="MuD"         description="Realized growth rate for diatomes"       unit="s-1"        axis_ref="deptht" />
+       <field id="LNnut"       description="Nutrient limitation term in Nanophyto"   unit="-"          axis_ref="deptht" />
+       <field id="LDnut"       description="Nutrient limitation term in Diatoms"     unit="-"          axis_ref="deptht" />
+       <field id="LNFe"        description="Iron limitation term in Nanophyto"       unit="-"          axis_ref="deptht" />
+       <field id="LDFe"        description="Iron limitation term in Diatoms"         unit="-"          axis_ref="deptht" />
+       <field id="LNlight"     description="Light limitation term in Nanophyto"      unit="-"          axis_ref="deptht" />
+       <field id="LDlight"     description="Light limitation term in Diatoms"        unit="-"          axis_ref="deptht" />
+       <field id="Fe2"         description="Iron II concentration"                   unit="nmol/L"     axis_ref="deptht" />
+       <field id="Fe3"         description="Iron III concentration"                  unit="nmol/L"     axis_ref="deptht" />
+       <field id="FeL1"        description="Complexed Iron concentration with L1"    unit="nmol/L"     axis_ref="deptht" />
+       <field id="FeL2"        description="Complexed Iron concentration with L2"    unit="nmol/L"     axis_ref="deptht" />
+       <field id="FeP"         description="Precipitated Iron III"                   unit="nmol/L"     axis_ref="deptht" />
+       <field id="TL1"         description="Total L1 concentration"                  unit="nmol/L"     axis_ref="deptht" />
+       <field id="TL2"         description="Total L2 concentration"                  unit="nmol/L"     axis_ref="deptht" />
+       <field id="pdust"       description="dust concentration"                      unit="g/L"                          />
+       <field id="Totlig"      description="Total ligand concentation"               unit="nmol/L"     axis_ref="deptht" />
+       <field id="Biron"       description="Bioavailable iron"                       unit="nmol/L"     axis_ref="deptht" />
+       <field id="Sdenit"      description="Nitrate reduction in the sediments"      unit="molN/m2/s"                    />
+       <field id="Ironice"     description="Iron input/uptake due to sea ice"        unit="molFe/m2/s"                   />
+       <field id="HYDR"        description="Iron input from hydrothemal vents"       unit="molFe/m2/s"                   />
+       <field id="EPC100"      description="Export of carbon particles at 100 m"     unit="mol/m2/s"                     />
+       <field id="EPFE100"     description="Export of biogenic iron at 100 m"        unit="mol/m2/s"                     />
+       <field id="EPSI100"     description="Export of Silicate at 100 m"             unit="mol/m2/s"                     />
+       <field id="EPCAL100"    description="Export of Calcite at 100 m"              unit="mol/m2/s"                     />
+       <field id="Cflx"        description="DIC flux"                                unit="mol/m2/s"                     />
+       <field id="Oflx"        description="Oxygen flux"                             unit="mol/m2/s"                     />
+       <field id="Kg"          description="Gas transfer"                            unit="mol/m2/s/uatm"                />
+       <field id="Dpco2"       description="Delta CO2"                               unit="uatm"                         />
+       <field id="Dpo2"        description="Delta O2"                                unit="uatm"                         />
+       <field id="Heup"        description="Euphotic layer depth"                    unit="m"                            />
+       <field id="Irondep"     description="Iron deposition from dust"               unit="mol/m2/s"                     />
+       <field id="Ironsed"     description="Iron deposition from sediment"           unit="mol/m2/s"  axis_ref="deptht"  />
      </group>         
 
