Changeset 2433

Timestamp:

Feb 21, 2006, 5:27:25 PM (18 years ago)

Author:

ole

Message:

Coarsened okushiri mesh for use with stochastic study

Location:

development/stochastic_study

Files:

• create_mesh.py (modified) (1 diff)
• run_model.py (modified) (2 diffs)

development/stochastic_study/create_mesh.py

@@ -112 +112 @@
 
     
-#
-    
+    #res_1
+    #base_resolution = 1
+    #ocean = m.addRegionEN(xleft+.1, ybottom+.1)
+    #ocean.setMaxArea(0.1*base_resolution)
+    #mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
+    #mid.setMaxArea(0.0005*base_resolution)    
+    #inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
+    #inner.setMaxArea(0.00007*base_resolution)
+
+    #res_2 (try to improve by adding gulleys back in at a coarser res but coarsen inner and mid - this was OK)  
+    #base_resolution = 1
+    #ocean = m.addRegionEN(xleft+.1, ybottom+.1)
+    #ocean.setMaxArea(0.1*base_resolution)
+    #mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
+    #mid.setMaxArea(0.001*base_resolution)    
+    #inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
+    #inner.setMaxArea(0.0001*base_resolution)
+    #gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
+    ##gulleys.setMaxArea(0.00002*base_resolution)
+    #gulleys.setMaxArea(0.0005*base_resolution)            
+
+    #res_3 (Coarsen a bit more: mid and inner - shoulders got rounded. Lost one reflection)  
+    #base_resolution = 1
+    #ocean = m.addRegionEN(xleft+.1, ybottom+.1)
+    #ocean.setMaxArea(0.1*base_resolution)
+    #mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
+    #mid.setMaxArea(0.01*base_resolution)    
+    #inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
+    #inner.setMaxArea(0.0002*base_resolution)
+    #gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
+    #gulleys.setMaxArea(0.0005*base_resolution)            
+
+    #res_4 (Refined inner again - still rounded. Need mid to be finer)  
+    #base_resolution = 1
+    #ocean = m.addRegionEN(xleft+.1, ybottom+.1)
+    #ocean.setMaxArea(0.1*base_resolution)
+    #mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
+    #mid.setMaxArea(0.01*base_resolution)    
+    #inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
+    #inner.setMaxArea(0.0001*base_resolution)
+    #gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
+    #gulleys.setMaxArea(0.0005*base_resolution)
+
+    #res_5 (Refined mid again - coarsened inner 5 times - this almost OK)
+    #base_resolution = 1
+    #ocean = m.addRegionEN(xleft+.1, ybottom+.1)
+    #ocean.setMaxArea(0.1*base_resolution)
+    #mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
+    #mid.setMaxArea(0.001*base_resolution)    
+    #inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
+    #inner.setMaxArea(0.0005*base_resolution)
+    #gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
+    #gulleys.setMaxArea(0.0005*base_resolution)
+
+    #res_6 (Like 5 but with inner a little finer, ocean and mid a bit coarser -  no worse than 5)
+    #base_resolution = 1
+    #ocean = m.addRegionEN(xleft+.1, ybottom+.1)
+    #ocean.setMaxArea(0.2*base_resolution)
+    #mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
+    #mid.setMaxArea(0.002*base_resolution)    
+    #inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
+    #inner.setMaxArea(0.0004*base_resolution)
+    #gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
+    #gulleys.setMaxArea(0.0005*base_resolution)
+
+    #res_7 (Like 6 but with everything coarsened a tad -
+    #a bit smoothed but still OK)
     base_resolution = 1
-    base_resolution = 10    
+    ocean = m.addRegionEN(xleft+.1, ybottom+.1)
+    ocean.setMaxArea(0.3*base_resolution)
+    mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
+    mid.setMaxArea(0.003*base_resolution)    
+    inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
+    inner.setMaxArea(0.0007*base_resolution)
+    gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
+    gulleys.setMaxArea(0.0007*base_resolution)                        
 
-    ocean = m.addRegionEN(xleft+.1, ybottom+.1)
-    ocean.setMaxArea(0.1*base_resolution)
-
-    mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
-    mid.setMaxArea(0.0005*base_resolution)
-
-    inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
-    inner.setMaxArea(0.00007*base_resolution)
-
-
-    #gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
-    #gulleys.setMaxArea(0.00002*base_resolution)    
-    
     
     m.generateMesh('pzq28.0za1000000a')

development/stochastic_study/run_model.py

@@ -35 +35 @@
 
 
-#------------------------------------------------------------------------------
-# Preparation of topographic data:
-# Convert txt file to pts NetCDF
-#------------------------------------------------------------------------------
-xya2pts(project.bathymetry_filename, verbose = True,
-        z_func = lambda z: -z)
-
-
 #-----------------------------------------------------------------------------
 # Setup computational domain
@@ -67 +59 @@
 domain.set_quantity('elevation',
                     filename = project.bathymetry_filename[:-4] + '.pts',
-                    alpha = 0.0001,
+                    alpha = 0.001,
                     verbose = True,
                     use_cache = True)