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