Changeset 852 for codes/icosagcm/devel/src
- Timestamp:
- 05/05/19 14:58:31 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dynamics/compute_pvort_only.F90
r840 r852 17 17 END SUBROUTINE check_interface 18 18 19 #if BEGIN_DYSL 20 21 KERNEL(pvort_only) 22 FORALL_CELLS_EXT() 23 ON_DUAL 24 etav = 0.d0 25 FORALL_EDGES 26 etav = etav + SIGN*u(EDGE) 27 END_BLOCK 28 hv=0. 29 FORALL_VERTICES 30 hv = hv + RIV2*rhodz(VERTEX) 31 END_BLOCK 32 qv(DUAL_CELL) = (etav + FV*AV )/(hv*AV) 33 END_BLOCK 34 END_BLOCK 35 36 FORALL_CELLS() 37 ON_EDGES 38 qu(EDGE)=0.5d0*(qv(VERTEX1)+qv(VERTEX2)) 39 END_BLOCK 40 END_BLOCK 41 END_BLOCK 42 43 #endif END_DYSL 44 19 45 SUBROUTINE compute_pvort_only_unst(u,rhodz,qu,qv, hv_) 20 46 USE ISO_C_BINDING, only : C_DOUBLE, C_FLOAT … … 92 118 93 119 END MODULE compute_pvort_only_mod 120
Note: See TracChangeset
for help on using the changeset viewer.