本文整理汇总了Python中anuga.file.netcdf.NetCDFFile.variables[lat_name][:]方法的典型用法代码示例。如果您正苦于以下问题：Python NetCDFFile.variables[lat_name][:]方法的具体用法？Python NetCDFFile.variables[lat_name][:]怎么用？Python NetCDFFile.variables[lat_name][:]使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类anuga.file.netcdf.NetCDFFile的用法示例。
在下文中一共展示了NetCDFFile.variables[lat_name][:]方法的2个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: most2nc
# 需要导入模块: from anuga.file.netcdf import NetCDFFile [as 别名]
# 或者: from anuga.file.netcdf.NetCDFFile import variables[lat_name][:] [as 别名]
def most2nc(input_file, output_file, inverted_bathymetry=False, verbose=True):
"""Convert a MOST file to NetCDF format.
input_file the input file to convert
output_file the name of the oputput NetCDF file to create
inverted_bathymetry ??
verbose True if the function is to be verbose
"""
# variable names
long_name = 'LON'
lat_name = 'LAT'
elev_name = 'ELEVATION'
# set up bathymetry
if inverted_bathymetry:
up = -1.
else:
up = +1.
# read data from the MOST file
in_file = open(input_file,'r')
if verbose: log.critical('reading header')
nx_ny_str = in_file.readline()
nx_str,ny_str = nx_ny_str.split()
nx = int(nx_str)
ny = int(ny_str)
h1_list=[]
for i in range(nx):
h1_list.append(float(in_file.readline()))
h2_list=[]
for j in range(ny):
h2_list.append(float(in_file.readline()))
h2_list.reverse()
if verbose: log.critical('reading depths')
in_depth_list = in_file.readlines()
in_file.close()
out_depth_list = [[]]
if verbose: log.critical('processing depths')
k=1
for in_line in in_depth_list:
for string in in_line.split():
#j = k/nx
out_depth_list[(k-1)/nx].append(float(string)*up)
if k==nx*ny:
break
if k-(k/nx)*nx ==0:
out_depth_list.append([])
k+=1
in_file.close()
out_depth_list.reverse()
depth_list = out_depth_list
# write the NetCDF file
if verbose: log.critical('writing results')
out_file = NetCDFFile(output_file, netcdf_mode_w)
out_file.createDimension(long_name,nx)
out_file.createVariable(long_name,'d',(long_name,))
out_file.variables[long_name].point_spacing='uneven'
out_file.variables[long_name].units='degrees_east'
out_file.variables[long_name][:] = h1_list
out_file.createDimension(lat_name,ny)
out_file.createVariable(lat_name,'d',(lat_name,))
out_file.variables[lat_name].point_spacing='uneven'
out_file.variables[lat_name].units='degrees_north'
out_file.variables[lat_name][:] = h2_list
out_file.createVariable(elev_name,'d',(lat_name,long_name))
out_file.variables[elev_name].point_spacing='uneven'
out_file.variables[elev_name].units='meters'
out_file.variables[elev_name][:] = depth_list
out_file.close()示例2: setUp
# 需要导入模块: from anuga.file.netcdf import NetCDFFile [as 别名]
# 或者: from anuga.file.netcdf.NetCDFFile import variables[lat_name][:] [as 别名]
def setUp(self):
import time
self.verbose = Test_File_Conversion.verbose
# Create basic mesh
points, vertices, boundary = rectangular(2, 2)
# Create shallow water domain
domain = Domain(points, vertices, boundary)
domain.default_order = 2
# Set some field values
domain.set_quantity('elevation', lambda x,y: -x)
domain.set_quantity('friction', 0.03)
######################
# Boundary conditions
B = Transmissive_boundary(domain)
domain.set_boundary( {'left': B, 'right': B, 'top': B, 'bottom': B})
######################
#Initial condition - with jumps
bed = domain.quantities['elevation'].vertex_values
stage = num.zeros(bed.shape, num.float)
h = 0.3
for i in range(stage.shape[0]):
if i % 2 == 0:
stage[i,:] = bed[i,:] + h
else:
stage[i,:] = bed[i,:]
domain.set_quantity('stage', stage)
domain.distribute_to_vertices_and_edges()
self.initial_stage = copy.copy(domain.quantities['stage'].vertex_values)
self.domain = domain
C = domain.get_vertex_coordinates()
self.X = C[:,0:6:2].copy()
self.Y = C[:,1:6:2].copy()
self.F = bed
#Write A testfile (not realistic. Values aren't realistic)
self.test_MOST_file = 'most_small'
longitudes = [150.66667, 150.83334, 151., 151.16667]
latitudes = [-34.5, -34.33333, -34.16667, -34]
long_name = 'LON'
lat_name = 'LAT'
nx = 4
ny = 4
six = 6
for ext in ['_ha.nc', '_ua.nc', '_va.nc', '_e.nc']:
fid = NetCDFFile(self.test_MOST_file + ext, netcdf_mode_w)
fid.createDimension(long_name,nx)
fid.createVariable(long_name,netcdf_float,(long_name,))
fid.variables[long_name].point_spacing='uneven'
fid.variables[long_name].units='degrees_east'
fid.variables[long_name][:] = longitudes
fid.createDimension(lat_name,ny)
fid.createVariable(lat_name,netcdf_float,(lat_name,))
fid.variables[lat_name].point_spacing='uneven'
fid.variables[lat_name].units='degrees_north'
fid.variables[lat_name][:] = latitudes
fid.createDimension('TIME',six)
fid.createVariable('TIME',netcdf_float,('TIME',))
fid.variables['TIME'].point_spacing='uneven'
fid.variables['TIME'].units='seconds'
fid.variables['TIME'][:] = [0.0, 0.1, 0.6, 1.1, 1.6, 2.1]
name = ext[1:3].upper()
if name == 'E.': name = 'ELEVATION'
fid.createVariable(name,netcdf_float,('TIME', lat_name, long_name))
fid.variables[name].units='CENTIMETERS'
fid.variables[name].missing_value=-1.e+034
fid.variables[name][:] = [[[0.3400644, 0, -46.63519, -6.50198],
[-0.1214216, 0, 0, 0],
[0, 0, 0, 0],
[0, 0, 0, 0]],
[[0.3400644, 2.291054e-005, -23.33335, -6.50198],
[-0.1213987, 4.581959e-005, -1.594838e-007, 1.421085e-012],
[2.291054e-005, 4.582107e-005, 4.581715e-005, 1.854517e-009],
[0, 2.291054e-005, 2.291054e-005, 0]],
[[0.3400644, 0.0001374632, -23.31503, -6.50198],
#.........这里部分代码省略.........