- Timestamp:
- 11/14/18 15:07:09 (6 years ago)
- Location:
- codes/icosagcm/devel/Python
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/dynamico/meshes.py
r770 r773 178 178 dual_deg,dual_edge,dual_ne,dual_vertex, 179 179 left,right,down,up,trisk_deg,trisk, 180 Aiv,Aiv,f+0.*Aiv,le_de,Riv2, -wee)180 Aiv,Aiv,f+0.*Aiv,le_de,Riv2,wee) 181 181 self.set_members(locals(), 'primal_deg', 'primal_edge', 'primal_ne', 'primal_vertex', 182 182 'dual_deg', 'dual_edge', 'dual_ne', 'dual_vertex', -
codes/icosagcm/devel/Python/test/py/RSW_2D.py
r754 r773 23 23 h0[:] = 1+0.1*(np.exp(-2.*(x1*x1+yy*yy))+ 24 24 np.exp(-2.*(x2*x2+yy*yy))) 25 #h0[:] = 1+0.1*(np.exp(-2.*yy*yy)) 26 25 27 flow0=(h0,u0) 26 28 -
codes/icosagcm/devel/Python/test/py/RSW_2D_mesh.py
r759 r773 35 35 h0 = 1+0.1*(np.exp(-2.*(x1*x1+yy*yy))+ 36 36 np.exp(-2.*(x2*x2+yy*yy))) 37 #h0 = 1+0.1*(np.exp(-5.*yy*yy)) 38 37 39 flow0=prec.asnum([h0,u0]) 38 40 … … 40 42 dt = cfl/math.sqrt((nx/Lx)**2+(ny/Ly)**2) 41 43 42 T=1 .44 T=10. 43 45 N=int(T/dt)+1 44 46 dt=T/N -
codes/icosagcm/devel/Python/test/py/write_Cartesian_mesh.py
r758 r773 19 19 default=8., choices=None, 20 20 help="Ly") 21 parser.add_argument("-llm", type=int,22 default=1, choices=[1],23 help="number of vertical levels")24 21 args = parser.parse_args() 25 nx, ny, Lx, Ly, llm, nqdyn = args.nx, args.ny,args.Lx, args.Ly, args.llm, 1 22 23 nx, ny, Lx, Ly, llm, nqdyn = args.nx, args.ny,args.Lx, args.Ly, 1, 1 26 24 27 25 dx,dy=Lx/nx,Ly/ny
Note: See TracChangeset
for help on using the changeset viewer.