Changes between Version 12 and Version 13 of DevelopmentActivities/AccelerationSpinup
- Timestamp:
- 2011-07-22T10:33:59+02:00 (13 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
DevelopmentActivities/AccelerationSpinup
v12 v13 2 2 3 3 This page presents a work which consists to implement a numerical method that 4 can accelerate the convergence of the SPINUP. ORCHIDEE usually use a iterativ 5 and specific method based on forcesoil, consisting by running the module 6 stomate_soilcarbon.f90 stand alone between 5000 and 10000 times.Before 7 forcesoil, you have to run ORCHIDEE a given length just to have the litter pools 4 can accelerate the convergence of the SPINUP. [[BR]] 5 6 7 ORCHIDEE usually use a iterative and specific method based on forcesoil, consisting by running the module 8 stomate_soilcarbon.f90 stand alone between 5000 and 10000 times.[[BR]] 9 10 Before forcesoil, you have to run ORCHIDEE a given length just to have the litter pools 8 11 at equilibrium. The problem with the FORCESOIL job is that you got some biaises 9 and takes longer at the global scale. 12 and takes longer at the global scale.[[BR]] 13 10 14 That's why we implement a new method for accelerating the spinup. 11 This method is based used by the model PASIM(see paper and description below). 15 This method is based used by the model PASIM(send a email to get the original papers or look 16 at the presentation available at the end of the page). 12 17 13 18 == Description of the method used by PASIM == 14 19 15 PASIM is a model of prairial management developped by the INRA institute. 20 PASIM is a model of prairial management developped by the INRA institute.[[BR]] 21 16 22 The soil module of PASIM is based on the CENTURY model like ORCHIDEE. It 17 23 consists on 5 carbon which are linked by fluxes : litter structural, 18 litter metabolic, carbon active, carbon slow and carbon passive. (See flow chart below).24 litter metabolic, carbon active, carbon slow and carbon passive. 19 25 we can represent symbolically the fluxes between the pools by a matrix of size (5,5). 20 26 21 27 Each line of the matrix represents a pool and each row the fluxes received (or leaving) by 22 the pool. We can consider the following equation (for more details see papers below).28 the pool. All these pools are set by the same differential equation.(see presentation at the end of the page). 23 29 24 30 Using the approximation of the derivative in time and two sequences, we obtain a equation … … 92 98 93 99 To be precise, these values are taken from the last restart file of the run. 94 The writing of the restart file happens at the same instant than the computation of the solution by the matrix method. 100 The writing of the restart file happens at the same instant than the computation of the solution by the matrix method.[[BR]] 101 95 102 In order to compare the results, I have to take the values for the carbon pools found in the restart files (variables 96 103 LITTER_xx_AB, LITTER_xx_BE ans CARBON_xx). … … 117 124 118 125 119 We obtained the following time series :126 We obtained the following time series (in red the results obtained by the matrix method) : 120 127 121 128 [[Image(comp_LITTER_STR_AB_LITTER_STR_AB_star.gif)]] [[Image(comp_LITTER_STR_BE_LITTER_STR_BE_star.gif)]] … … 239 246 === Work done === 240 247 241 - Implementation of the method of PASIM in ORCHIDEE : adaptation to 7 pools + f lux to the fire[[BR]]248 - Implementation of the method of PASIM in ORCHIDEE : adaptation to 7 pools + fire flux [[BR]] 242 249 243 250 - Implementation of a stopping condition based on the maximum of all the relative errors calculated for each pool[[BR]] … … 245 252 - 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]] 246 253 247 - Tests made on the temperate regions for pfts 10 and 6 254 - Tests made on the temperate regions for pfts 10 and 6 [[BR]] 255 256 - Creation of a branche on ths svn server called Spinup_analytic based on the externalized version 257 258 248 259 249 260 === To be done === … … 252 263 253 264 - Implementation of a new stopping condition based on the evaluation of the flux and not on the stock 265 266 267 === PS === 268 269 - If you want to know more about this work, look at the presentation below and/or send a email to Didier.Solyga.at.lsce.ipsl.fr