Changeset 666 for codes/icosagcm

Timestamp:

01/20/18 18:47:41 (6 years ago)

Author:

dubos

Message:

devel/unstructured : python fixes

Location:

codes/icosagcm/devel/Python

Files:

• dynamico/__init__.py (modified) (1 diff)
• test/py/RSW2_MPAS_W02.py (modified) (1 diff)
• test/py/slice_GW_NH.py (modified) (3 diffs)

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

@@ -5 +5 @@
 import libs
 libs.libicosa.dynamico_show_openmp()
+libs.libicosa.dynamico_init_trace()

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

@@ +1 @@
+print 'Loading DYNAMICO modules ...'
+from dynamico import unstructured as unst
+from dynamico.meshes import MPAS_Mesh as Mesh
+from dynamico import time_step
+print '...Done'
+
 print 'Loading modules ...'
 import math as math
 import matplotlib.pyplot as plt
 import numpy as np
-print '...Done'
-
-print 'Loading DYNAMICO modules ...'
-from dynamico import unstructured as unst
-from dynamico.meshes import MPAS_Mesh as Mesh
-from dynamico import time_step
 print '...Done'
 

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

@@ -10 +10 @@
 import numpy as np
 import time
-
-unst.setvar('nb_threads', 1)
 
 def run(mesh,metric,thermo,BC,caldyn_eta,u, m0ik,gas0,Phi0_il, T=300., courant=2.9):
@@ -97 +95 @@
 
 def deform(x,eta): # x and eta go from 0 to 1
-    h0, dd, xi = 2000., 5000., 4000.
+    h0, dd, xi = 500., 5000., 4000.
     cos2=lambda x: np.cos(np.pi*x)**2
     gauss=lambda x: np.exp(-x*x)
@@ -105 +103 @@
 
 #Lx, nx, llm = 3e5, 100, 1e4, 30
-Lx, nx, ztop, llm = 3e5, 300, 3e4, 80
+Lx, nx, ztop, llm = 3e5, 300, 1e4, 80
 nqdyn, ny, dx = 1, 1, Lx/nx
 mesh = Mesh(nx,ny,llm,nqdyn,nx*dx,ny*dx,0.)