~florian-rathgeber/fluidity/femtools-python-ctypes

« back to all changes in this revision

Viewing changes to tools/petsc_readnsolve.F90

  • Committer: Colin Cotter
  • Date: 2011-06-06 20:05:17 UTC
  • mfrom: (3470.2.26 element-numbering)
  • Revision ID: colin.cotter@imperial.ac.uk-20110606200517-088m0syfzs0diwqs
This merge from element-numbering branch implements trace spaces as a 
mesh_type.
This required to change the interface of allocate_element to
allocate_element(element, ele_num, ngi, type, stat)
where ele_num is the element numbering type. Therefore I had to change the 
interface in various places.

Added trace spaces as a mesh type to mesh_options.rn*, and included it
as an option in shallow water options. This hasn't been included in
Fluidity yet (could do if there was a demand for it) so no changes to
fluidity_options.rn*.

Added Test_Trace_Space.F90 to tools, which tests allocating the mesh
in populate state from options, and then does some numbering checks on
the resulting field.

Bon appetit!

Show diffs side-by-side

added added

removed removed

Lines of Context:
342
342
    integer i, n, istate, stat
343
343
    type(halo_type), pointer :: my_halo
344
344
    integer nstates, universal_nodes, components
 
345
    type(ele_numbering_type), pointer :: ele_num
345
346
 
346
347
    ewrite(1,*) "Opening flml file ", trim(flml)
347
348
    call load_options(flml)
443
444
    else
444
445
    
445
446
      ! allocate a dummy shape and mesh:
446
 
      call allocate(shape, 1, 1, 1, 1)
 
447
       ele_num => find_element_numbering(&
 
448
         &vertices = 1, dimension = 1, degree = 1)
 
449
      call allocate(shape,ele_num,1)
447
450
      call allocate(mesh, n, 1, shape, "Mesh")
448
451
      
449
452
      ! setup trivial petsc numbering