- Timestamp:
- 10/08/18 18:57:20 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/test/py/RSW_2D.py
r747 r754 4 4 from dynamico import time_step 5 5 from dynamico import DCMIP 6 from dynamico.precision import zero, np_num 7 6 8 import math as math 7 9 import matplotlib.pyplot as plt … … 9 11 import time 10 12 11 nx,ny,llm,nqdyn=128,1 28,1,112 Lx,Ly,g,f = 8., 8.,1.,1.13 nx,ny,llm,nqdyn=128,192,1,1 14 Lx,Ly,g,f = 8.,12.,1.,1. 13 15 dx,dy=Lx/nx,Ly/ny 14 16 … … 30 32 dt=T/N 31 33 print N,dt,Lx/nx 32 scheme = time_step.ARK2(caldyn.bwd_fast_slow, dt )34 scheme = time_step.ARK2(caldyn.bwd_fast_slow, dt, precision=np_num) 33 35 34 36 print 'types : yy, u0, h0, u0', yy.dtype, u0.dtype, h0.dtype, u0.dtype … … 38 40 h,u=flow 39 41 print 'types : h,u', h.dtype, u.dtype 40 caldyn.bwd_fast_slow(flow, 0.)42 caldyn.bwd_fast_slow(flow, zero) 41 43 plt.figure(); plt.pcolor(mesh.x,mesh.y,caldyn.qv) 42 44 plt.colorbar(); plt.title('potential vorticity')
Note: See TracChangeset
for help on using the changeset viewer.