Changeset 757

Timestamp:

10/10/18 12:37:11 (6 years ago)

Author:

dubos

Message:

devel/Python : added some features to Cartesian_Mesh

File:

• codes/icosagcm/devel/Python/dynamico/meshes.py (modified) (2 diffs)

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

@@ -183 +183 @@
     def field_u(self,n=1): return zeros((n,self.ny,2*self.nx,self.llm))
     def field_ps(self,n=1): return zeros((n,self.ny,self.nx))
+    def field_ucomp(self,n=1): return zeros((n,self.ny,self.nx,self.llm))
     def ucomp(self,u): 
         return u[range(0,2*self.nx,2),:] if self.ny==1 else u[:,range(0,2*self.nx,2),:]
@@ -188 +189 @@
         if self.ny==1 : uv[range(0,2*self.nx,2),:]=u
         else : uv[:,range(0,2*self.nx,2),:]=u
-    def vcomp(self,u): return u[:,range(1,2*self.nx,2),:]
+    def vcomp(self,u):
+        return u[range(1,2*self.nx,2),:] if self.ny==1 else u[:,range(1,2*self.nx,2),:]
+    def set_vcomp(self,uv,v):
+        if self.ny==1 : uv[range(1,2*self.nx,2),:]=v
+        else : uv[:,range(1,2*self.nx,2),:]=v
 
 #---------------------- DYNAMICO format for fully unstructured mesh and curvilinear meshes ----------------------