Changes between Version 11 and Version 12 of DevelopmentActivities/AccelerationSpinup

Timestamp:

2011-07-22T09:30:20+02:00 (13 years ago)

Author:

dsolyga

Comment:

--

DevelopmentActivities/AccelerationSpinup

@@ -151 +151 @@
 }}}
 
-This means we smooth our curves. Finally, we got :
+
+This means we smooth our curves except the last one with the global stock. Finally, we got :
 
 [[Image(error_LITTER_STR_AB_inferior_to_0.2%.gif)]] [[Image(error_LITTER_STR_BE_inferior_to_0.2%.gif)]]
-[[Image(error_LITTER_MET_AB_inferior_to_0.2%.gif)]]  [[Image(error_LITTER_MET_BE_inferior_to_0.5%.gif)]]
-[[Image(error_CARBON_ACTIVE_inferior_to_0.5%.gif)]]   [[Image(error_CARBON_SLOW_inferior_to_0.5%.gif)]]
-[[Image(error_CARBON_PASSIVE_inferior_to_2%.gif)]]  
-
+[[Image(error_LITTER_MET_AB_inferior_to_0.2%.gif)]] [[Image(error_LITTER_MET_BE_inferior_to_0.5%.gif)]]
+[[Image(error_CARBON_ACTIVE_inferior_to_0.5%.gif)]] [[[Image(error_CARBON_SLOW_inferior_to_0.5%.gif)]]
+[[Image(error_CARBON_PASSIVE_inferior_to_2%.gif)]]  [[Image(error_stock.gif)]]
+
+The vertical axis correspond to the error in percentage. Notice that the carbon stock could be considered
+at equilibrium before the passive pool is at equilibrium (we got more than 1000 years difference!). [[BR]]
+
+For this reason, we have chosen in a first time step to implement a stopping condition based on the maximum
+of the relative errors for each pool. It is still not satisfactory. The next step will be to implement a stopping
+condition based on the carbon flux.   
 
 === Comparison with Forcesoil ===
@@ -207 +214 @@
         
 The output files can be found at : /dmnfs12/cont003/p529sol/IGCM_OUT/OL2/SpinUp_classic_setup_nofire
+
+Forcesoil is very performing. Let's take the example of the passive pool.[[BR]]
+
+ORCHIDEE+forcesoil runs only 500Y to find a error of 0.05%.[[BR]]
+
+The algebraic method have to run 1500Y to have a similar error (and 500Y more to find the exact result) 
+and ORCHIDEE alone got the same value as forcesoil after 6200Y. [[BR]]
+
+Our analytic method is 3 to 4 more faster than ORCHIDEE alone but less performing than forcesoil. [[BR]]
+
+The analytic method found the same results after 500Y than ORCHIDEE+forcesoil in 51Y.[[BR]]
+(see /dmnfs12/cont003/p529sol/IGCM_OUT/OL2/50Y_ORC_FORCESOIL_2/)
+
+
+=== Another test ===
+
+A another test was made using the same point and the pft 6. The results obtained are globally the same : they
+are not detailed but could be found here : /dmnfs12/cont003/p529sol/IGCM_OUT/OL2/spinup_7000Y_pft6.tar
+
+
+
+== Summary ==
+
+=== Work done ===
+
+- Implementation of the method of PASIM in ORCHIDEE : adaptation to 7 pools + flux to the fire[[BR]]
+
+- Implementation of a stopping condition based on the maximum of all the relative errors calculated for each pool[[BR]]
+
+- Modification of the script and job of OOL_SEC_STO : if the stopping condition chosen by the user is reached before the end of the run, the job will stop[[BR]]
+
+- Tests made on the temperate regions for pfts 10 and 6
+
+=== To be done ===
+
+- New tests on a boreal region, tropical region  and global[[BR]]
+
+- Implementation of a new stopping condition based on the evaluation of the flux and not on the stock