| 245 | |
| 246 | === Suggestions, remarques === |
| 247 | |
| 248 | * Simplification du code : on peut éliminer quelques paramètres redondants comme maint_resp_slope, pheno_gdd_cri et senescence_temp. Par exemple, maint_resp_slope est défini comme ceci : |
| 249 | |
| 250 | {{{ |
| 251 | REAL(r_std), ALLOCATABLE, SAVE, DIMENSION(:,:) :: maint_resp_slope !! slope of maintenance respiration coefficient |
| 252 | !! (1/K, 1/K^2, 1/K^3), used in the code |
| 253 | |
| 254 | maint_resp_slope(:,1) = maint_resp_slope_c(:) |
| 255 | maint_resp_slope(:,2) = maint_resp_slope_b(:) |
| 256 | maint_resp_slope(:,3) = maint_resp_slope_a(:) |
| 257 | }}} |
| 258 | On l'utilise dans l'équation suivante dan stomate_npp et stomate_resp : |
| 259 | |
| 260 | {{{ |
| 261 | slope(:) = maint_resp_slope(j,1) + tl(:) * maint_resp_slope(j,2) + & |
| 262 | tl(:)*tl(:) * maint_resp_slope(j,3) |
| 263 | }}} |
| 264 | alors que l'on pourrait utiliser directement maint_resp_slope_a, maint_resp_slope_b et maint_resp_slope_c : |
| 265 | {{{ |
| 266 | slope(:) = maint_resp_slope_c(j) + tl(:) * maint_resp_slope_b(j) + & |
| 267 | tl(:)*tl(:) * maint_resp_slope_a(j) |
| 268 | }}} |
| 269 | * Documentation : Martial a converti en latex les documents techniques que Fabienne et lui avaient écrits. Ils sont disponibles dans sa version personnelle sur svn. Il faudrait qu'ils soient mis à jour. |