示例#1
0
""", "splat_coords")

    op2.par_loop(kernel, crds_latlon.cell_set,
                 crds_latlon.dat(op2.RW, crds_latlon.cell_node_map()))
    return Mesh(crds_latlon)


if write_latlon:
    outfile_ll = File('{0}_ll.pvd'.format(fname))
    field_output_ll = []
    mesh_ll = get_latlon_mesh(mesh)
    for f in field_output:
        FSpwG = functionspaceimpl.WithGeometry
        field_ll = Function(FSpwG(f.function_space(), mesh_ll),
                            val=f.topological,
                            name=f.name() + '_ll')
        field_output_ll.append(field_ll)
else:
    outfile_ll = None
    field_output_ll = None

nc_temp = zeros((nc_h5_dumpfreq, len(nc_diag)))

# Create checkpoint
chkpt = DumbCheckpoint("chkpt_{0}".format(fname), mode=FILE_CREATE)
if h5_safe:
    chkpt_bckp = DumbCheckpoint("chkpt_bckp_{0}".format(fname),
                                mode=FILE_CREATE)


# Nc file functions
}
""", "splat_coords")

    op2.par_loop(kernel, crds_latlon.cell_set,
                 crds_latlon.dat(op2.RW, crds_latlon.cell_node_map()))
    return Mesh(crds_latlon)


if write_latlon:
    outfile_ll = File('{0}_ll.pvd'.format(fname))
    field_output_ll = []
    mesh_ll = get_latlon_mesh(mesh)
    for f in field_output:
        FSpwG = functionspaceimpl.WithGeometry
        field_ll = Function(FSpwG(f.function_space(), mesh_ll),
                            val=f.topological, name=f.name() + '_ll')
        field_output_ll.append(field_ll)
else:
    outfile_ll = None
    field_output_ll = None

nc_temp = zeros((nc_h5_dumpfreq, len(nc_diag)))

# Create checkpoint
chkpt = DumbCheckpoint("chkpt_{0}".format(fname), mode=FILE_CREATE)
if h5_safe:
    chkpt_bckp = DumbCheckpoint("chkpt_bckp_{0}".format(fname),
                                mode=FILE_CREATE)

# Nc file functions
def create_nc_file(pfx=''):