[cig-commits] [commit] devel: update compatibility with cubit13+ (0f3b542)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d/compare/3b04b68452f55959366a6a10017252704d6e711d...0f3b542a91577fd3452edc0a5d46c4d8e62bc914

>---------------------------------------------------------------

commit 0f3b542a91577fd3452edc0a5d46c4d8e62bc914
Author: emanuele casarotti <casarotti at tiscali.it>
Date:   Thu Jan 29 16:24:19 2015 +0100

    update compatibility with cubit13+
    
    fix compatibility issue with the creation of volumes
    
    it solves this error:
    
    File "specfem3d/CUBIT_GEOCUBIT/GEOCUBIT.py", line 56, in
    mesh_volume.mesh()
    
    File "specfem3d/CUBIT_GEOCUBIT/geocubitlib/mesh_volume.py", line 44, in mesh
    mesh_layercake_regularmap(filename=filename)
    File "specfem3d/CUBIT_GEOCUBIT/geocubitlib/mesh_volume.py", line 100, in mesh_layercake_regularmap
    vol[id_vol].intervalv=cfg.iv_interval[id_vol]
    IndexError: list index out of range


>---------------------------------------------------------------

0f3b542a91577fd3452edc0a5d46c4d8e62bc914
 CUBIT_GEOCUBIT/geocubitlib/volumes.py | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/CUBIT_GEOCUBIT/geocubitlib/volumes.py b/CUBIT_GEOCUBIT/geocubitlib/volumes.py
index 739d998..3967e43 100644
--- a/CUBIT_GEOCUBIT/geocubitlib/volumes.py
+++ b/CUBIT_GEOCUBIT/geocubitlib/volumes.py
@@ -433,6 +433,10 @@ def layercake_volume_ascii_regulargrid_mpiregularmap(filename=None,verticalsandw
             cubitcommand= 'del surface all'
             cubit.cmd(cubitcommand)
             list_vol=cubit.parse_cubit_list("volume","all")
+            for v in list_vol:
+                sv=cubit.get_relatives("volume",v,"surface")
+                if len(sv) <=1: cubit.cmd("del vol "+str(v))
+            list_vol=cubit.parse_cubit_list("volume","all")
             if len(list_vol) > 1:     
                 cubitcommand= 'imprint volume all'
                 cubit.cmd(cubitcommand)
@@ -637,7 +641,7 @@ def create_volume(surf1,surf2,method='loft'):
             cmd='create surface skin curve '+str(cs[1])+' '+str(cs[0])
             cubit.cmd(cmd)
         is_stop=cubit.get_last_id('surface')
-        cmd="create volume surface "+str(surf1)+' '+str(surf2)+' '+str(is_start)+' to '+str(is_stop)+"  heal keep"
+        cmd="create volume surface "+str(surf1)+' '+str(surf2)+' '+str(is_start)+' to '+str(is_stop)+"  heal"
         cubit.cmd(cmd)