     float corrected = getcorrectedelevation(solarelevation);
     float opticalpath = getopticalpath(corrected, dem,
                                        HEIGHT);
-    float opticaldepth = getopticaldepth(opticalpath);  
+    float opticaldepth = getopticaldepth(opticalpath);
     return gethorizontalirradiance(EXT_RAD,
                                    excentricity, zenithangle) *
            getbeamtransmission(linke, opticalpath,

{
    float anglerelation = pow(0.5f / cos(satellitalzenithangle * DEG2RAD),
                              0.8f);
    return ((i0met[0] * diffuseclearsky * anglerelation) / 
            (PI * EXT_RAD[0]));
}


__device__ void getalbedo(float *albedo, float radiance,
float *i0met, float *excentricity, float zenithangle)
{
    albedo[i_dxyt] = ((PI * radiance) /
            (i0met[0] * excentricity[i_dt] *
            cos(zenithangle * DEG2RAD)));
}
