Subversion Repositories public iLand

Rev

Rev 557 | Rev 564 | Go to most recent revision | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 557 Rev 562
Line 174... Line 174...
174 
    int iday=0;
174 
    int iday=0;
175 
    // calculate the water-factor for each month (see Adair et al 2008)
175 
    // calculate the water-factor for each month (see Adair et al 2008)
176 
    double fw_month[12];
176 
    double fw_month[12];
177 
    double ratio;
177 
    double ratio;
178 
    for (int m=0;m<12;m++) {
178 
    for (int m=0;m<12;m++) {
179 
        if (mRU->waterCycle()->potentialEvapotranspiration()[m]>0.)
179 
        if (mRU->waterCycle()->referenceEvapotranspiration()[m]>0.)
180 
            ratio = mRU->climate()->precipitationMonth()[m] /  mRU->waterCycle()->potentialEvapotranspiration()[m];
180 
            ratio = mRU->climate()->precipitationMonth()[m] /  mRU->waterCycle()->referenceEvapotranspiration()[m];
181 
        else
181 
        else
182 
            ratio = 0;
182 
            ratio = 0;
183 
        fw_month[m] = 1. / (1. + 30.*exp(-8.5*ratio));
183 
        fw_month[m] = 1. / (1. + 30.*exp(-8.5*ratio));
184 
        qDebug() <<"month"<< m << "PET" << mRU->waterCycle()->potentialEvapotranspiration()[m] << "prec" <<mRU->climate()->precipitationMonth()[m];
184 
        qDebug() <<"month"<< m << "PET" << mRU->waterCycle()->referenceEvapotranspiration()[m] << "prec" <<mRU->climate()->precipitationMonth()[m];
185 
    }
185 
    }
186 

            186
            

          

          

            187
            
    for (const ClimateDay *day=mRU->climate()->begin(); day!=mRU->climate()->end(); ++day, ++iday)


            187
            
    for (const ClimateDay *day=mRU->climate()->begin(); day!=mRU->climate()->end(); ++day, ++iday)


          

          

            188
            
    {


            188
            
    {


          

          

            189
            
        ft = exp(308.56*(1./56.02-1./((273.+day->temperature)-227.13)));  // empirical variable Q10 model of Lloyd and Taylor (1994), see also Adair et al. (2008)


            189
            
        ft = exp(308.56*(1./56.02-1./((273.+day->temperature)-227.13)));  // empirical variable Q10 model of Lloyd and Taylor (1994), see also Adair et al. (2008)