surface = [] f = open(file, 'r') for line in f: line = line.strip() columns = line.split() # If data point is from open ocean (1) or from a lead (2) if columns[0] == '1' or columns[0] == '2': # If data point is listed as 'valid' if columns[1] == '1': lat.append(float(columns[5])) lon.append(float(columns[6])) ssha.append((float(columns[7]) - float(columns[8]))) surface.append(int(columns[0])) tracker_type = funct.mode_points(lat, lon, month) for point in range(len(tracker_type)): if tracker_type[point] == 1 and surface[point] == '1': ssha[point] += funct.apply_offset(month, 'SAR') # Find the boundaries iedge_ocean_ice = [] for it in range(len(surface)): # Find ice boundaries if surface[it:it + len_bound] == ocean_to_ice: iedge_ocean_ice.append(it + len_bound//2) elif surface[it:it + len_bound] == ice_to_ocean: iedge_ocean_ice.append(it + len_bound//2) for step in iedge_ocean_ice:
surface_type = [] f = open(file, 'r') for line in f: line = line.strip() columns = line.split() # If data point is from open ocean (1) or from a lead (2) if columns[0] == '1' or columns[0] == '2': # If data point is listed as 'valid' if columns[1] == '1': lat.append(float(columns[5])) lon.append(float(columns[6])) ssha.append(float(columns[7]) - float(columns[8])) # Generate a list of retracker modes for this track surface_type = funct.mode_points(lat, lon, month) iedge_SAR = [] iedge_SARIn = [] for it in range(len(surface_type)): if surface_type[it:it + Len_SAR] == SAR_edge: iedge_SAR.append(it + Len_SAR//2) elif surface_type[it:it + Len_SAR] == SAR_edge2: iedge_SAR.append(it + Len_SAR//2) elif surface_type[it:it + Len_SARIn] == SARIn_edge: iedge_SARIn.append(it + Len_SARIn//2) elif surface_type[it:it + Len_SARIn] == SARIn_edge2: iedge_SARIn.append(it + Len_SARIn//2) for step in iedge_SAR:
line.strip() columns = line.split() lat_asc_filt.append(float(columns[0])) ssha_asc_filt.append(float(columns[1])) output_ssh.close() os.system('rm ../OUTPUT_ssh.dat') lat = list(lat_desc_filt) + list(lat_asc_filt) ssha = list(ssha_desc_filt) + list(ssha_asc_filt) lon = list(lon_desc) + list(lon_asc) surface = list(surface_desc) + list(surface_asc) if len(lat) == len(lon): # Generate a list of retracker modes for this track tracker_type = funct.mode_points(lat, lon, month) # Find the boundaries iedge_LRM_SAR = [] for it in range(len(tracker_type)): # Find retracker boundaries if tracker_type[it:it + len_bound] == LRM_to_SAR: iedge_LRM_SAR.append(it + len_bound // 2) elif tracker_type[it:it + len_bound] == SAR_to_LRM: iedge_LRM_SAR.append(it + len_bound // 2) # For the LRM-SAR boundary for step in iedge_LRM_SAR: # If the boundary is ALL OCEAN if surface[step - len_bound // 2:step + len_bound // 2] == [1] * len_bound: