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

    float corrected = getcorrectedelevation(solarelevation);
    float opticalpath = getopticalpath(corrected, dem,
                                       HEIGHT);
    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]));
}

