33 | | During our validations on fluxnet sites, we modify the herbivory model and the formulation for the variable |
34 | | "black_carbon". |
| 33 | During our validations on fluxnet sites, we modified the herbivory model and the formulation for the variable |
| 34 | "black_carbon". In '''stomate_season.f90''', we replace : |
| 35 | {{{ |
| 36 | DO j = 2,nvm |
| 37 | ! |
| 38 | IF ( natural(j) ) THEN |
| 39 | ! |
| 40 | WHERE ( nlflong_nat(:) .GT. zero ) |
| 41 | consumption(:) = hvc1 * nlflong_nat(:) ** hvc2 |
| 42 | herbivores(:,j) = one_year * green_age(:) * nlflong_nat(:) / consumption(:) |
| 43 | ELSEWHERE |
| 44 | herbivores(:,j) = 100000. |
| 45 | ENDWHERE |
| 46 | ! |
| 47 | ELSE |
| 48 | ! |
| 49 | herbivores(:,j) = 100000. |
| 50 | ! |
| 51 | ENDIF |
| 52 | ! |
| 53 | ENDDO |
| 54 | }}} |
| 55 | by : |
| 56 | {{{ |
| 57 | DO j = 2,nvm |
| 58 | ! |
| 59 | IF ( natural(j) ) THEN |
| 60 | ! |
| 61 | WHERE ( nlflong_nat(:,j) .GT. zero ) |
| 62 | consumption(:) = hvc1 * nlflong_nat(:,j) ** hvc2 |
| 63 | herbivores(:,j) = one_year * green_age(:,j) * nlflong_nat(:,j) / consumption(:) |
| 64 | ELSEWHERE |
| 65 | herbivores(:,j) = 100000. |
| 66 | ENDWHERE |
| 67 | ! |
| 68 | ELSE |
| 69 | ! |
| 70 | herbivores(:,j) = 100000. |
| 71 | ! |
| 72 | ENDIF |
| 73 | ! |
| 74 | ENDDO |
| 75 | }}} |
| 76 | in order to make herbivores more pft-specific. |