Changeset 799 for codes

Timestamp:

01/14/19 10:54:58 (5 years ago)

Author:

dubos

Message:

devel/Python : shift positions in Cartesian meshes by dx/2, dy/2

File:

• codes/icosagcm/devel/Python/dynamico/meshes.py (modified) (1 diff)

codes/icosagcm/devel/Python/dynamico/meshes.py

@@ -169 +169 @@
         primal_i, primal_j = [ x.astype(np.int32) for x in lon_i, lat_i]
         edge_i, edge_j = [ x.astype(np.int32) for x in lon_e, lat_e]
-        lon_i, lon_v, lon_e = [x*dx-Lx/2 for x in lon_i, lon_v, lon_e]
-        lat_i, lat_v, lat_e = [y*dy-Ly/2 for y in lat_i, lat_v, lat_e]
+        lon_i, lon_v, lon_e = [(x+.5)*dx-.5*Lx for x in lon_i, lon_v, lon_e]
+        lat_i, lat_v, lat_e = [(y+.5)*dy-.5*Ly for y in lat_i, lat_v, lat_e]
 
         Aiv[:]=dx*dy # Ai=Av=dx*dy