Changeset 658 for codes/icosagcm/devel/src/kernels_unst/caldyn_vert.k90
- Timestamp:
- 12/30/17 02:00:38 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/kernels_unst/caldyn_vert.k90
r624 r658 4 4 !$OMP DO SCHEDULE(STATIC) 5 5 DO ij = 1, primal_num 6 !DIR$ SIMD 6 7 DO l = 2, llm 7 8 dtheta_rhodz(l,ij,iq) = dtheta_rhodz(l,ij,iq) + 0.5*(theta(l,ij,iq)+theta(l-1,ij,iq))*wflux(l,ij) … … 11 12 !$OMP DO SCHEDULE(STATIC) 12 13 DO ij = 1, primal_num 14 !DIR$ SIMD 13 15 DO l = 1, llm-1 14 16 dtheta_rhodz(l,ij,iq) = dtheta_rhodz(l,ij,iq) - 0.5*(theta(l,ij,iq)+theta(l+1,ij,iq))*wflux(l+1,ij) … … 23 25 ij_left = left(edge) 24 26 ij_right = right(edge) 27 !DIR$ SIMD 25 28 DO l = 2, llm 26 29 wwuu(l,edge) = .25*(wflux(l,ij_left)+wflux(l,ij_right))*(u(l,edge)+u(l-1,edge)) ! Fu … … 34 37 ij_left = left(edge) 35 38 ij_right = right(edge) 39 !DIR$ SIMD 36 40 DO l = 1, llm 37 41 dFu_deta = wwuu(l+1,edge)-wwuu(l,edge) ! d/deta (F*u) … … 46 50 ij_left = left(edge) 47 51 ij_right = right(edge) 52 !DIR$ SIMD 48 53 DO l = 2, llm 49 54 wwuu(l,edge) = .5*(wflux(l,ij_left)+wflux(l,ij_right))*(u(l,edge)-u(l-1,edge)) … … 57 62 ij_left = left(edge) 58 63 ij_right = right(edge) 64 !DIR$ SIMD 59 65 DO l = 1, llm 60 66 du(l,edge) = du(l,edge) - (wwuu(l,edge)+wwuu(l+1,edge)) / (rhodz(l,ij_left)+rhodz(l,ij_right))
Note: See TracChangeset
for help on using the changeset viewer.