     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,

__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)));
}

    excentricity[i_dt] = getexcentricity(gamma);
    bc = getbeamirradiance(EXT_RAD, excentricity, solarangle,
                           solarelevation[i_dxyt], linke, dem, HEIGHT);

    dc = getdiffuseirradiance(EXT_RAD, excentricity,
                              solarelevation[i_dxyt], linke);
    getglobalirradiance(gc, bc, dc);
