Changes between Version 17 and Version 18 of DevelopmentActivities/inputRoutine

Timestamp:

2017-07-30T16:39:32+02:00 (7 years ago)

Author:

dgoll

Comment:

--

DevelopmentActivities/inputRoutine

@@ -79 +79 @@
 == Deposition data ==
 We have time series of N and P deposition from Wang et al. (2017) 'Global forest carbon uptake due to nitrogen and phosphorus deposition from 1850 to 2100' ready to use for ORCHIDEE.
-
 === since r4535 ===
 {{{
@@ -89 +88 @@
 curie:/ccc/work/cont003/dsm/p529goll/NP_deposition/final_NEW
 }}}
+
+== use libIGCM to update yearly the fluxes in run.def in case of IMPOSE_NMAP=n == 
+Add this to orchidee_ol.driver
+COMP/orchidee_ol.driver:
+{{{
+##########################################################################
+# NP DEPOSITION FROM FILE - SITE SIMULATION
+##############################################################################################
+    if [ -f NPdeposition.txt ] ; then
+        # If this file exist, then read N and P deposition 
+        NOY=`grep Annee_${year} NPdeposition.txt | awk -F= '{print $2}'`
+        NHX=`grep Annee_${year} NPdeposition.txt | awk -F= '{print $3}'`
+        Pdep=`grep Annee_${year} NPdeposition.txt | awk -F= '{print $4}'`
+
+        if [ X"${NOY}" = X ] ; then
+            # The grep returned empty variable: stop execution
+            IGCM_debug_Exit "The file NPdeposition.txt.txt do not contain the current year."
+            IGCM_debug_Verif_Exit
+        fi
+        IGCM_comp_modifyDefFile nonblocker run.def Ndep_NHX ${NHX}
+        IGCM_comp_modifyDefFile nonblocker run.def Ndep_NOY ${NOY}
+        IGCM_comp_modifyDefFile nonblocker run.def Pdep_p   ${Pdep}
+    fi
+
+##############################################################################################
+}}}
+add to COMP/stomate.card the  link to your deposition file
+{{{
+[ParametersFiles]
+List=   (${SUBMIT_DIR}/PARAM/run.def, .),\
+        (/ccc/work/cont003/dsm/p529goll/AmaFACE_forcing/v01/CO2_1850_2100_AmaFACE.txt, CO2.txt),\
+        (/ccc/work/cont003/dsm/p529goll/AmaFACE_forcing/AmaFACE18502100_NPdeposition.txt, NPdeposition.txt)
+}}}
+the deposition file should store yearly fluxes in g/m2/yr; sequence is Annee_YEAR,NOY,NHX,Pdep:
+{{{
+Annee_1850,0.0715,0.0715,0.0144
+Annee_1851,0.0715,0.0715,0.0144
+Annee_1852,0.0715,0.0715,0.0144
+Annee_1853,0.0715,0.0715,0.0144
+Annee_1854,0.0715,0.0715,0.0144
+Annee_1855,0.0715,0.0715,0.0144
+Annee_1856,0.0715,0.0715,0.0144
+Annee_1857,0.072,0.072,0.0144
+Annee_1858,0.072,0.072,0.0144
+}}}
+
+
 
 = For ENSEMBLE setup =