Changeset 753

Timestamp:

10/08/18 18:46:28 (6 years ago)

Author:

dubos

Message:

devel/unstructured : fixed indexing bug in Cartesian mesh

Location:

codes/icosagcm/devel/Python

Files:

• dynamico/meshes.py (modified) (1 diff)
• test/py/NH_3D_bubble.py (modified) (1 diff)

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

@@ -94 +94 @@
         self.x, self.y, self.lev, self.levp1 = x,y,lev,levp1
         # 3D coordinate arrays
-        self.xx,self.yy,self.ll = np.meshgrid(x,y,lev, indexing='ij')       
-        self.xxp1,self.yyp1,self.llp1 = np.meshgrid(x,y,levp1, indexing='ij')
+        self.yy,self.xx,self.ll = np.meshgrid(y,x,lev, indexing='ij')       
+        self.yyp1,self.xxp1,self.llp1 = np.meshgrid(y,x,levp1, indexing='ij')
         # beware conventions for indexing
         # Fortran order : llm,nx*ny,nqdyn  / indices start at 1

codes/icosagcm/devel/Python/test/py/NH_3D_bubble.py

@@ -117 +117 @@
     print 'ptop, model top (m) :', unst.getvar('ptop'), Phi.max()/unst.getvar('g')
     jj=ny/2
-    xx,zz,ss,ww = mesh.xx[:,jj,:]/1000., z[:,jj,:]/1000., s[:,jj,:], w[:,jj,:]
+    xx,zz,ss,ww = mesh.xx[jj,:,:]/1000., z[jj,:,:]/1000., s[jj,:,:], w[jj,:,:]
 
     f, (ax1, ax2) = plt.subplots(1, 2, sharey=True, figsize=(12,4))