Changeset 587 for codes/icosagcm/devel/src/kernels/flux_centered.k90

Timestamp:

10/16/17 15:22:59 (7 years ago)

Author:

dubos

Message:

devel : multiplicative factor in flux/wind reconstructions

File:

• codes/icosagcm/devel/src/kernels/flux_centered.k90 (modified) (2 diffs)

codes/icosagcm/devel/src/kernels/flux_centered.k90

@@ -2 +2 @@
    !---------------------------- flux_centered ----------------------------------
    ! NB : here the input data is a flux and has already the factor l_e in it
+   ! Input data is rescaled by factor scale
    ! Perot reconstruction based on Gauss theorem
    ! 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)
          uz = uz + ue_le*(.5*(xyz_v(ij+z_rup,3)+xyz_v(ij+z_rdown,3))-cz)
-         ucenter(ij,l,1)=ux
-         ucenter(ij,l,2)=uy
-         ucenter(ij,l,3)=uz
+         fac = scale*(1./Ai(ij))
+         ucenter(ij,l,1)=ux*fac
+         ucenter(ij,l,2)=uy*fac
+         ucenter(ij,l,3)=uz*fac
       END DO
    END DO