Changes between Version 2 and Version 3 of UgCellArrays


Ignore:
Timestamp:
11/09/06 14:33:55 (13 years ago)
Author:
tomgross
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • UgCellArrays

    v2 v3  
    1111A simple cell is defined by the indices which call out which nodes, edges, faces comprise the cell.  In this strawman CDL case the simple cell for the grid will be a triangle defined by three nodes and having some variables specified at edges and the central face  of the triangle.  Below in the grid meta data I call this "tri_face_edged" .
    1212
    13 The data set will consist of a grid of cells.  The index array for the grid will dimensioned to the number of cells.  The data set will also consist of a collection of nodes, defined by their spatial location.  The data set will also contain a collection of edges.  Edges are defined by two nodes.  Edges must have a direction associated with them, so it is necessary to define an edge independent of the two different cells which it may be part of.     
     13The data set will consist of a grid of cells.  The index array for the grid will dimensioned to the number of cells.  The data set will also consist of a collection of nodes, defined by their spatial location.  The data set will also contain a collection of edges.  Edges are defined by two nodes.  Edges must have a direction associated with them, so it is necessary to define an edge independent of the two different cells which it may be part of.   So, grid(icell,:) provides the indices to all the other data:
     14 * grid(icell, :) = [ 4 5 8 2 3 4 10 ]
     15 * The triangle is defined by lon([4 5 8]), lat([4 5 8])
     16 * The face or area average temperature is at:  temp(2)
     17 * The Heat flux to the inside of the triangle is given by the sum( Ftemp([3 4 10]).*length_edge([3 4 10]) ) with something clever done to get the sign of the inside to outside of the triangle correct.  That info is in there, but more later.....(sorry)
    1418
    1519I find it interesting that there is enough information specified in the grid and edge arrays to define the connectivity of the irregular grid.  But it all makes no sense without data which needs to be pointed to by indices of the grid.  So there must be locations for the nodes, longitude and latitude.  A variable which is specified at the nodes will be dimensioned by the number of nodes, like the zeta variable below.  Edge variables are not defined at the center of an edge.  They are defined for the whole line segment and you would expect them to be objects like normal fluxes, Uf and Ftemp below.  If you want to locate a variable like flux at a point, I think you are not thinking about flux anymore.  To have any transport a flux must be applied to a length or area, not a single point.