Changeset 19055

Timestamp:

01/30/15 18:04:17 (10 years ago)

Author:

bdef

Message:

NEW: adding netcdf export for matlab and minor change in python

Location:

issm/trunk-jpl/src/m/contrib/netCDF

Files:

• export_netCDF.m (added)
• export_netCDF.py (modified) (2 diffs)

issm/trunk-jpl/src/m/contrib/netCDF/export_netCDF.py

@@ -14 +14 @@
         def GetDim(var,shape,i,istime):
                 output=[]
-                for dim in range(0,i):
-                        if type(shape[0])==int:
+                for dim in range(0,i): #loop on the dimensions
+                        if type(shape[0])==int: 
                                 try:
-                                        output=output+[str(DimDict[shape[dim]])]
-                                except KeyError:
+                                        output=output+[str(DimDict[shape[dim]])] #test if the dimension allready exist
+                                except KeyError: #if not create it
                                         if (shape[dim])>1:
-                                                NewDim=NCData.createDimension(str(field),(shape[dim]))
-                                                DimDict[len(NewDim)]=str(field)
+                                                NewDim=NCData.createDimension(str(field)+str(dim),(shape[dim]))
+                                                DimDict[len(NewDim)]=str(field)+str(dim)
                                                 output=output+[str(DimDict[shape[dim]])]
-                                                print 'Defining dimension ' +str(shape[dim]) +' for '+str(field)
+                                                print 'Defining dimension ' +str(shape[dim]) +' for '+str(field)+str(dim)
                         elif type(shape[0])==str:#dealling with a dictionnary
                                 try:
@@ -129 +129 @@
         VertNum=NCData.createDimension('VertNum',md.mesh.numberofvertices)
         VertperElt=NCData.createDimension('VertperElt',numpy.shape(md.mesh.elements)[1])
-        if type(md.mesh) is mesh2d:
-                DimNum=NCData.createDimension('DimNum',2)
-        elif type(md.mesh) is mesh3dprisms:
-                DimNum=NCData.createDimension('DimNum',3)
-        else:
-                print 'I can not get the Dimension of the problem'
-        EdgeNum=NCData.createDimension('EdgeNum',md.mesh.numberofedges)
-        Time=NCData.createDimension('Time',StepNum)
-        SegNum=NCData.createDimension('SegNum',numpy.shape(md.mesh.segmentmarkers)[0])
-        InvSteps=NCData.createDimension('InvSteps',md.inversion.nsteps)
-        DictDim=NCData.createDimension('DictDim',2)
 
         DimDict = {len(EltNum):'EltNum',
                                                  len(VertNum):'VertNum',
-                                                 len(VertperElt):'VertperElt',
-                                                 len(DimNum):'DimNum',
-                                                 len(EdgeNum):'EdgeNum',
-                                                 len(Time):'Time',
-                                                 len(SegNum):'SegNum',
-                                                 len(InvSteps):'InvSteps'}
+                                                 len(VertperElt):'VertperElt'}
 
         TypeDict = {float:'f8',