Changeset 658 for codes/icosagcm/devel/src/kernels_unst/caldyn_vert_NH.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_NH.k90
r614 r658 9 9 eta_dot(l,ij) = wflux_ij / mass(l,ij) 10 10 wcov(l,ij) = w_ij*(geopot(l+1,ij)-geopot(l,ij)) 11 !DIR$ SIMD 11 12 DO l = 2, llm-1 12 13 w_ij = .5*( W(l,ij)+W(l+1,ij) )/mass(l,ij) … … 29 30 ij_left = left(edge) 30 31 ij_right = right(edge) 32 !DIR$ SIMD 31 33 DO l = 1, llm 32 34 du(l,edge) = du(l,edge) - .5*(wcov(l,ij_left)+wcov(l,ij_right))*1.*(eta_dot(l,ij_right)-eta_dot(l,ij_left)) … … 38 40 !$OMP DO SCHEDULE(STATIC) 39 41 DO ij = 1, primal_num 42 !DIR$ SIMD 40 43 DO l = 2, llm 41 44 dPhi(l,ij)=dPhi(l,ij)-wflux(l,ij)*(geopot(l+1,ij)-geopot(l-1,ij))/(mass(l-1,ij)+mass(l,ij)) … … 49 52 l=1 50 53 dW(l,ij) = dW(l,ij) - W_etadot(l,ij) 54 !DIR$ SIMD 51 55 DO l = 2, llm 52 56 dW(l,ij) = dW(l,ij) + W_etadot(l-1,ij) - W_etadot(l,ij)
Note: See TracChangeset
for help on using the changeset viewer.