- Timestamp:
- 10/10/18 12:37:11 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/dynamico/meshes.py
r756 r757 183 183 def field_u(self,n=1): return zeros((n,self.ny,2*self.nx,self.llm)) 184 184 def field_ps(self,n=1): return zeros((n,self.ny,self.nx)) 185 def field_ucomp(self,n=1): return zeros((n,self.ny,self.nx,self.llm)) 185 186 def ucomp(self,u): 186 187 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 189 190 else : uv[:,range(0,2*self.nx,2),:]=u 190 def vcomp(self,u): return u[:,range(1,2*self.nx,2),:] 191 def vcomp(self,u): 192 return u[range(1,2*self.nx,2),:] if self.ny==1 else u[:,range(1,2*self.nx,2),:] 193 def set_vcomp(self,uv,v): 194 if self.ny==1 : uv[range(1,2*self.nx,2),:]=v 195 else : uv[:,range(1,2*self.nx,2),:]=v 191 196 192 197 #---------------------- DYNAMICO format for fully unstructured mesh and curvilinear meshes ----------------------
Note: See TracChangeset
for help on using the changeset viewer.