Last change
on this file since 614 was
614,
checked in by dubos, 7 years ago
|
devel : added unstructured-mesh kernels
|
File size:
1.0 KB
|
Rev | Line | |
---|
[614] | 1 | !-------------------------------------------------------------------------- |
---|
| 2 | !---------------------------- pvort_only ---------------------------------- |
---|
| 3 | !$OMP DO SCHEDULE(STATIC) |
---|
| 4 | DO ij = 1, dual_num |
---|
| 5 | DO l = 1, llm |
---|
| 6 | etav = 0.d0 |
---|
| 7 | DO iedge = 1, dual_deg(ij) |
---|
| 8 | edge = dual_edge(iedge,ij) |
---|
| 9 | etav = etav + dual_ne(iedge,ij)*u(l,edge) |
---|
| 10 | END DO |
---|
| 11 | hv=0. |
---|
| 12 | DO ivertex = 1, dual_deg(ij) |
---|
| 13 | vertex = dual_vertex(ivertex,ij) |
---|
| 14 | hv = hv + Riv2(ivertex,ij)*rhodz(l,vertex) |
---|
| 15 | END DO |
---|
| 16 | qv(l,ij) = (etav + fv(ij)*Av(ij) )/(hv*Av(ij)) |
---|
| 17 | END DO |
---|
| 18 | END DO |
---|
| 19 | !$OMP END DO |
---|
| 20 | !$OMP DO SCHEDULE(STATIC) |
---|
| 21 | DO edge = 1, edge_num |
---|
| 22 | ij_up = up(edge) |
---|
| 23 | ij_down = down(edge) |
---|
| 24 | DO l = 1, llm |
---|
| 25 | qu(l,edge)=0.5d0*(qv(l,ij_down)+qv(l,ij_up)) |
---|
| 26 | END DO |
---|
| 27 | END DO |
---|
| 28 | !$OMP END DO |
---|
| 29 | !---------------------------- pvort_only ---------------------------------- |
---|
| 30 | !-------------------------------------------------------------------------- |
---|
Note: See
TracBrowser
for help on using the repository browser.