Subversion Repositories public iLand

Rev

Rev 485 | Rev 536 | Go to most recent revision | Show entire file | Regard 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 ...
74 
        utilizable_rad = calculateUtilizablePAR(i); // utilizable radiation of the month times ... (MJ/m2)
75 
        epsilon = calculateEpsilon(i); // ... photosynthetic efficiency ...
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