surface_conductance = ice_column.get_conductance_at_z() # should return U_surface
SC = -temp_surface*surface_conductance # freezing temp is 0C
SC *= time_span_in_sec # energy over 24hrs
SC /= 1000 # J to kJ
return SC, surface_conductance
if __name__ == "__main__":
import getFiledata as gfd
from setEnvironment import data_path
icecols = gfd.importColumns("{0}Semsvann observasjoner 2012-2013.csv".format(data_path))
icecol = icecols[6]
eb = temp_surface_from_eb(utm33_x=130513, utm33_y=6802070, ice_column=icecol, temp_atm=5,
prec=0.01,prec_snow=0.1,time_span_in_sec=24*60*60,albedo_prim=0.37)
S, s_inn, albedo, albedo_prim, age_factor_tau\
= get_short_wave(utm33_x=130513, utm33_y=6802070, day_no=180, temp_atm= -1.5, cloud_cover=0.5,
snow_depth=1., snow_density=0.3, prec_snow=0.01,
time_hour=0, time_span_in_sec=24*60*60, temp_surface=-2.,
age_factor_tau=None, albedo_prim=0.30, albedo_method="ueb")
a = 1
if __name__ == "__main__":
import getFiledata as gfd
from setEnvironment import data_path
# some tests of the functions in the IceColumn file
#date = datetime.datetime.strptime("2011-10-05", "%Y-%m-%d")
#icecol = IceColumn(date, 0)
# icecol.add_layer_at_index(0, 0.1, 'new_snow')
#icecol.update_slush_level()
#icecol.merge_and_remove_excess_layers()
#icecol.merge_snow_layers_and_compress(-5)
icecols = gfd.importColumns("{0}test ice columns.csv".format(data_path))
icecol = icecols[4]
icecol.update_column_temperatures(-15)
icecol.update_total_column_height()
layer = icecol.get_layer_at_z(0.9)
conductance = icecol.get_conductance_at_z()
height = icecol.get_depth_at_conductance(const.U_surface)
R13 = 0.3/0.11+0.18/1.12+0.28/2.24
U11 = 0.11/0.30
U12 = 1 / (1/0.3667 + 0.18/1.12)
U13 = 1 / (1/0.346 + 0.28/2.24)
R13_2 = 1/U13
a = 1
surface_conductance = ice_column.get_conductance_at_z(
) # should return U_surface
SC = -temp_surface * surface_conductance # freezing temp is 0C
SC *= time_span_in_sec # energy over 24hrs
SC /= 1000 # J to kJ
return SC, surface_conductance
if __name__ == "__main__":
import getFiledata as gfd
from setenvironment import data_path
icecols = gfd.importColumns(
"{0}Semsvann observasjoner 2012-2013.csv".format(data_path))
icecol = icecols[6]
eb = temp_surface_from_eb(utm33_x=130513,
utm33_y=6802070,
ice_column=icecol,
temp_atm=5,
prec=0.01,
prec_snow=0.1,
time_span_in_sec=24 * 60 * 60,
albedo_prim=0.37)
S, s_inn, albedo, albedo_prim, age_factor_tau\
= get_short_wave(utm33_x=130513, utm33_y=6802070, day_no=180, temp_atm= -1.5, cloud_cover=0.5,
snow_depth=1., snow_density=0.3, prec_snow=0.01,
time_hour=0, time_span_in_sec=24*60*60, temp_surface=-2.,
if __name__ == "__main__":
import getFiledata as gfd
from setEnvironment import data_path
# some tests of the functions in the IceColumn file
#date = datetime.datetime.strptime("2011-10-05", "%Y-%m-%d")
#icecol = IceColumn(date, 0)
# icecol.add_layer_at_index(0, 0.1, 'new_snow')
#icecol.update_slush_level()
#icecol.merge_and_remove_excess_layers()
#icecol.merge_snow_layers_and_compress(-5)
icecols = gfd.importColumns("{0}test ice columns.csv".format(data_path))
icecol = icecols[4]
icecol.update_column_temperatures(-15)
icecol.update_total_column_height()
layer = icecol.get_layer_at_z(0.9)
conductance = icecol.get_conductance_at_z()
height = icecol.get_depth_at_conductance(const.U_surface)
R13 = 0.3 / 0.11 + 0.18 / 1.12 + 0.28 / 2.24
U11 = 0.11 / 0.30
U12 = 1 / (1 / 0.3667 + 0.18 / 1.12)
U13 = 1 / (1 / 0.346 + 0.28 / 2.24)
R13_2 = 1 / U13
a = 1
