def process_albedos(data, root):
i0met = geti0met()
excentricity = nc.getvar(root,'excentricity')[:]
solarangle = nc.getvar(root,'solarangle')[:]
atmosphericalbedo = nc.getvar(root,'atmosphericalbedo')[:]
t_earth = nc.getvar(root,'t_earth')[:]
t_sat = nc.getvar(root,'t_sat')[:]
say("Calculating observed albedo, apparent albedo, effective albedo and cloud albedo... ")
observedalbedo = geo.getalbedo(data, i0met , excentricity, solarangle)
v_albedo = nc.clonevar(root, 'data', 'observedalbedo')
v_albedo[:] = observedalbedo
nc.sync(root)
apparentalbedo = geo.getapparentalbedo(observedalbedo,atmosphericalbedo, t_earth, t_sat)
v_albedo = nc.clonevar(root, 'data', 'apparentalbedo')
v_albedo[:] = apparentalbedo
nc.sync(root)
effectivealbedo = geo.geteffectivealbedo(solarangle)
v_albedo = nc.clonevar(root, 'data', 'effectivealbedo')
v_albedo[:] = effectivealbedo
nc.sync(root)
cloudalbedo = geo.getcloudalbedo(effectivealbedo,atmosphericalbedo,t_earth,t_sat)
v_albedo = nc.clonevar(root, 'data', 'cloudalbedo')
v_albedo[:] = cloudalbedo
nc.sync(root)
v_albedo = None
def process_atmospheric_irradiance(lat, lon, data, root):
i0met = geti0met()
dc = nc.getvar(root,'dc')[:]
satellitalzenithangle = nc.getvar(root,'satellitalzenithangle')[:]
excentricity = nc.getvar(root,'excentricity')[:]
say("Calculating atmospheric irradiance... ")
atmosphericradiance = geo.getatmosphericradiance(1367.0, i0met,dc, satellitalzenithangle)
atmosphericalbedo = geo.getalbedo(atmosphericradiance, i0met, excentricity, satellitalzenithangle)
satellitalelevation = geo.getelevation(satellitalzenithangle)
v_atmosphericalbedo = nc.getvar(root, 'atmosphericalbedo', 'f4', ('timing','northing','easting',),4)
v_atmosphericalbedo[:] = atmosphericalbedo
v_satellitalelevation = nc.getvar(root, 'satellitalelevation', 'f4', ('northing','easting',),4)
v_satellitalelevation[:] = satellitalelevation
nc.sync(root)
def process_atmospheric_irradiance(root):
i0met = geti0met()
dc = nc.getvar(root,'dc')[:]
satellitalzenithangle = nc.getvar(root,'satellitalzenithangle')[:]
excentricity = nc.getvar(root,'excentricity')[:]
say("Calculating atmospheric irradiance... ")
atmosphericradiance = geo.getatmosphericradiance(1367.0, i0met,dc, satellitalzenithangle)
atmosphericalbedo = geo.getalbedo(atmosphericradiance, i0met, excentricity, satellitalzenithangle)
satellitalelevation = geo.getelevation(satellitalzenithangle)
v_atmosphericalbedo = nc.clonevar(root, 'data', 'atmosphericalbedo')
v_atmosphericalbedo[:] = atmosphericalbedo
v_satellitalelevation = nc.clonevar(root, 'lon', 'satellitalelevation')
v_satellitalelevation[:] = satellitalelevation
nc.sync(root)
v_atmosphericalbedo, v_satellitalelevation = None, None
def process_albedos(lat, lon, data, root):
i0met = geti0met()
excentricity = nc.getvar(root,'excentricity')[:]
solarangle = nc.getvar(root,'solarangle')[:]
atmosphericalbedo = nc.getvar(root,'atmosphericalbedo')[:]
t_earth = nc.getvar(root,'t_earth')[:]
t_sat = nc.getvar(root,'t_sat')[:]
say("Calculating observed albedo, apparent albedo, effective albedo and cloud albedo... ")
observedalbedo = geo.getalbedo(data, i0met , excentricity, solarangle)
apparentalbedo = geo.getapparentalbedo(observedalbedo,atmosphericalbedo, t_earth, t_sat)
effectivealbedo = geo.geteffectivealbedo(solarangle)
cloudalbedo = geo.getcloudalbedo(effectivealbedo,atmosphericalbedo,t_earth,t_sat)
v_observedalbedo = nc.getvar(root, 'observedalbedo', 'f4', ('timing','northing','easting',),4)
v_observedalbedo[:] = observedalbedo
v_apparentalbedo = nc.getvar(root, 'apparentalbedo', 'f4', ('timing','northing','easting',),4)
v_apparentalbedo[:] = apparentalbedo
v_effectivealbedo = nc.getvar(root, 'effectivealbedo', 'f4', ('timing','northing','easting',),4)
v_effectivealbedo[:] = effectivealbedo
v_cloudalbedo = nc.getvar(root, 'cloudalbedo', 'f4', ('timing','northing','easting',),4)
v_cloudalbedo[:] = cloudalbedo
nc.sync(root)