Changeset 754

Timestamp:

10/08/18 18:57:20 (6 years ago)

Author:

dubos

Message:

devel/unstructured : fix RSW_2D test case for mixed precision

File:

• codes/icosagcm/devel/Python/test/py/RSW_2D.py (modified) (4 diffs)

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

@@ -4 +4 @@
 from dynamico import time_step
 from dynamico import DCMIP
+from dynamico.precision import zero, np_num
+
 import math as math
 import matplotlib.pyplot as plt
@@ -9 +11 @@
 import time
 
-nx,ny,llm,nqdyn=128,128,1,1
-Lx,Ly,g,f = 8.,8.,1.,1.
+nx,ny,llm,nqdyn=128,192,1,1
+Lx,Ly,g,f = 8.,12.,1.,1.
 dx,dy=Lx/nx,Ly/ny
 
@@ -30 +32 @@
 dt=T/N
 print N,dt,Lx/nx
-scheme = time_step.ARK2(caldyn.bwd_fast_slow, dt)
+scheme = time_step.ARK2(caldyn.bwd_fast_slow, dt, precision=np_num)
 
 print 'types : yy, u0, h0, u0', yy.dtype, u0.dtype, h0.dtype, u0.dtype
@@ -38 +40 @@
     h,u=flow
     print 'types : h,u', h.dtype, u.dtype
-    caldyn.bwd_fast_slow(flow, 0.)
+    caldyn.bwd_fast_slow(flow, zero)
     plt.figure(); plt.pcolor(mesh.x,mesh.y,caldyn.qv)
     plt.colorbar(); plt.title('potential vorticity')