Ignore:
Timestamp:
11/21/18 16:16:06 (6 years ago)
Author:
jisesh
Message:

devel/unstructured : laplacian operator

File:
1 edited

Legend:

Unmodified
Added
Removed

  • codes/icosagcm/devel/Python/src/kernels_transport.jin

    r687 r784  
    310310    END_BLOCK 
    311311END_BLOCK 
     312 
     313KERNEL(scalar_laplacian) 
     314 
     315  FORALL_CELLS_EXT() 
     316    ON_EDGES 
     317       grad(EDGE) = SIGN*(b(CELL2)-b(CELL1)) 
     318    END_BLOCK 
     319  END_BLOCK 
     320 
     321  FORALL_CELLS_EXT() 
     322    ON_PRIMAL 
     323      div_ij=0. 
     324      FORALL_EDGES 
     325        div_ij = div_ij + SIGN*LE_DE*grad(EDGE) 
     326      END_BLOCK 
     327      divu(CELL) = div_ij / AI 
     328    END_BLOCK 
     329  END_BLOCK 
     330 
     331END_BLOCK 
Note: See TracChangeset for help on using the changeset viewer.