Rev 485 | Rev 536 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 485 | Rev 513 | ||
---|---|---|---|
Line 69... | Line 69... | ||
69 | clear(); |
69 | clear(); |
70 | double utilizable_rad, epsilon; |
70 | double utilizable_rad, epsilon; |
71 | // conversion from gC to kg Biomass: C/Biomass=0.5
|
71 | // conversion from gC to kg Biomass: C/Biomass=0.5
|
72 | const double gC_to_kg_biomass = 1. / (biomassCFraction * 1000.); |
72 | const double gC_to_kg_biomass = 1. / (biomassCFraction * 1000.); |
73 | for (int i=0;i<12;i++) { |
73 | for (int i=0;i<12;i++) { |
74 | utilizable_rad = calculateUtilizablePAR(i); // utilizable radiation of the month times ... |
- | |
75 | epsilon = calculateEpsilon(i); // ... photosynthetic efficiency ... |
- | |
- | 74 | utilizable_rad = calculateUtilizablePAR(i); // utilizable radiation of the month times ... (MJ/m2) |
|
- | 75 | epsilon = calculateEpsilon(i); // ... photosynthetic efficiency ... (gC/MJ) |
|
76 | mUPAR[i] = utilizable_rad ; |
76 | mUPAR[i] = utilizable_rad ; |
77 | mGPP[i] =utilizable_rad * epsilon * gC_to_kg_biomass; // ... results in GPP of the month kg Biomass/m2 (converted from gC/m2) |
77 | mGPP[i] =utilizable_rad * epsilon * gC_to_kg_biomass; // ... results in GPP of the month kg Biomass/m2 (converted from gC/m2) |
78 | year_raw_gpp += mGPP[i]; // kg Biomass/m2 |
78 | year_raw_gpp += mGPP[i]; // kg Biomass/m2 |
79 | }
|
79 | }
|
80 | 80 |