Changeset 587 for codes/icosagcm/devel/src/kernels/flux_centered.k90
- Timestamp:
- 10/16/17 15:22:59 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/kernels/flux_centered.k90
r585 r587 2 2 !---------------------------- flux_centered ---------------------------------- 3 3 ! NB : here the input data is a flux and has already the factor l_e in it 4 ! Input data is rescaled by factor scale 4 5 ! Perot reconstruction based on Gauss theorem 5 6 ! u = sum( u.edge_normal * edge_length * (edge_midpoint-cell_centroid) ) /cell_area … … 35 36 uy = uy + ue_le*(.5*(xyz_v(ij+z_rup,2)+xyz_v(ij+z_rdown,2))-cy) 36 37 uz = uz + ue_le*(.5*(xyz_v(ij+z_rup,3)+xyz_v(ij+z_rdown,3))-cz) 37 ucenter(ij,l,1)=ux 38 ucenter(ij,l,2)=uy 39 ucenter(ij,l,3)=uz 38 fac = scale*(1./Ai(ij)) 39 ucenter(ij,l,1)=ux*fac 40 ucenter(ij,l,2)=uy*fac 41 ucenter(ij,l,3)=uz*fac 40 42 END DO 41 43 END DO
Note: See TracChangeset
for help on using the changeset viewer.