""", "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=''):