Changeset 587 for codes/icosagcm/devel/src/kernels
- Timestamp:
- 10/16/17 15:22:59 (7 years ago)
- Location:
- codes/icosagcm/devel/src/kernels
- Files:
-
- 2 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 -
codes/icosagcm/devel/src/kernels/wind_centered.k90
r585 r587 34 34 uy = uy + ue_le*(.5*(xyz_v(ij+z_rup,2)+xyz_v(ij+z_rdown,2))-cy) 35 35 uz = uz + ue_le*(.5*(xyz_v(ij+z_rup,3)+xyz_v(ij+z_rdown,3))-cz) 36 ucenter(ij,l,1)=ux 37 ucenter(ij,l,2)=uy 38 ucenter(ij,l,3)=uz 36 fac = scale*(1./Ai(ij)) 37 ucenter(ij,l,1)=ux*fac 38 ucenter(ij,l,2)=uy*fac 39 ucenter(ij,l,3)=uz*fac 39 40 END DO 40 41 END DO
Note: See TracChangeset
for help on using the changeset viewer.