42 | | 3. Le calcul de la concentration mose fait dans |
| 42 | 3. Le calcul de la concentration mensuelle se fait dans ATM_Update. |
| 43 | {{{ |
| 44 | co2_ppm = co2_ppm_prec + ( fco2_ff + t_oce_co2_flx + net_co2_flux_monthly + fco2_luc ) / 2.12 |
| 45 | }}} |
| 46 | 2.12 est le nombre de PetaGC / ppm. |
| 47 | |
| 48 | Attention aux signes des flux : convention atmosphérique. Ils doivent être comptés positivement si ils forment une source de carbone pour l'atmosphère. |
| 49 | Pour l'océan, on doit le vérifier : |
| 50 | {{{ |
| 51 | ! Conversion in GtC/yr ; negative for outgoing from ocean |
| 52 | t_oce_co2_flx = (-1.) * t_oce_co2_flx * 12. / 1.e15 |
| 53 | }}} |
| 54 | Pour la surface, c'est OK : |
| 55 | {{{ |
| 56 | ! CO2 flux in gC/m**2/sec (positive towards the atmosphere) is sum of: |
| 57 | co2_flux(:,:) = resp_hetero(:,:) + resp_maint(:,:) + resp_growth(:,:) & |
| 58 | & + (co2_fire(:,:)-co2_to_bm_dgvm(:,:))*veget_cov_max(:,:)/one_day & |
| 59 | & - gpp(:,:) |
| 60 | }}} |
| 61 | A voir ensuite: les signes des flux totaux de déforestation. La NBP est IPCC est positive pour la surface. |
| 62 | {{{ |
| 63 | nbp(:)=(SUM((gpp_daily-(resp_maint+resp_growth+resp_hetero)-co2_fire) & |
| 64 | & *veget_max,dim=2)-cflux_prod_total-harvest_above)/1e3/one_day*contfrac |
| 65 | }}} |
| 66 | |