![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 21427
- Timestamp:
- 12/02/16 02:14:12 (8 years ago)
- Location:
- issm/trunk-jpl/src/m/plot
- Files:
-
- 3 edited
-
plot_unit.py (modified) (4 diffs)
-
plotdoc.py (modified) (2 diffs)
-
processmesh.py (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/plot/plot_unit.py
r21426 r21427 42 42 options.addfield('colormap',cmap) 43 43 # }}} 44 # {{{ if plotting only one of several layers reduce dataset 45 if options.getfieldvalue('layer',0)>=1: 46 plotlayer=options.getfieldvalue('layer',1) 47 if datatype==1: 48 slicesize=np.shape(elements)[0] 49 elif datatype in [2,3]: 50 slicesize=len(x) 51 data=data[(plotlayer-1)*slicesize:plotlayer*slicesize] 52 # }}} 44 53 # {{{ Get the colormap limits 45 54 if options.exist('clim'): … … 86 95 else: 87 96 triangles=mpl.tri.Triangulation(x,y,elements) 88 ax.tricontourf(triangles,data,colorlevels,cmap=cmap,norm=norm,alpha=alpha)97 tri=ax.tricontourf(triangles,data,colorlevels,cmap=cmap,norm=norm,alpha=alpha) 89 98 if edgecolor != 'None': 90 99 ax.triplot(x,y,elements,color=edgecolor) … … 93 102 return 94 103 # }}} 95 # {{{ plotting quiver (TODO)104 # {{{ plotting quiver 96 105 elif datatype==3: 97 106 if is2d: … … 103 112 # }}} 104 113 # {{{ plotting P1 Patch (TODO) 114 105 115 elif datatype==4: 106 116 print 'plot_unit message: P1 patch plot not implemented yet' 107 117 return 118 108 119 # }}} 109 120 # {{{ plotting P0 Patch (TODO) 121 110 122 elif datatype==5: 111 123 print 'plot_unit message: P0 patch plot not implemented yet' 112 124 return 125 113 126 # }}} 114 127 else: -
issm/trunk-jpl/src/m/plot/plotdoc.py
r21355 r21427 10 10 'elementnumbering':' numbering of elements (matlab indices)'} 11 11 TODOdata={'basal_drag':' plot the basal drag on the bed (in kPa) based on the velocity in md.initialization', 12 'basal_dragx or basal_dragy' :' plot a component of the basal drag on the bed (in kPa)',13 'boundaries':' this will draw all the segment boundaries to the model, including rifts.',14 'icefront':' this will show segments that are used to define the icefront of the model (Neumann boundary conditions).',15 'deviatoricstress_tensor':' plot the components of the deviatoric stress tensor (tauxx,tauyy,tauzz,tauxy,tauxz,tauyz, if computed',16 'deviatoricstress_principal':' plot the deviatoricstress tensor principal axis and principal values',17 'deviatoricstress_principalaxis1':' arrow plot the first principal axis of the deviatoricstress tensor(replace 1 by 2 or 3 if needed)',18 'driving_stress':' plot the driving stress (in kPa)',19 'elements_type':' model used for each element',20 'vertexnumbering':' numbering of vertices',21 'highlightelements':' to highlight elements to highlight the element list',22 'highlightvertices':' to highlight vertices (use highlight option to enter the vertex list',23 'referential':' stressbalance referential',24 'riftvel':' velocities along rifts',25 'riftrelvel':' relative velocities along rifts',26 'riftpenetration':' penetration levels for a fault',27 'riftfraction':' fill fractions for every node of the rifts',28 'rifts':' plot mesh with an offset so that rifts are visible',29 'strainrate_tensor':' plot the components of the strainrate tensor (exx,eyy,ezz,exy,exz,eyz) if computed',30 'strainrate_principal':' plot the strainrate tensor principal axis and principal values)',31 'strainrate_principalaxis1':' arrow plot the first principal axis of the strainrate tensor(replace 1 by 2 or 3 if needed)',32 'stress_tensor':' plot the components of stress tensor (sxx,syy,szz,sxy,sxz,syz) if computed',33 'stress_principal':' plot the stress tensor principal axis and principal values',34 'stress_principalaxis1':' arrow plot the first principal axis of the stress tensor(replace 1 by 2 or 3 if needed)',35 'transient_results':' this will printlay all the time steps of a transient run (use steps to specify the steps requested)',36 'transient_vel':' this will printlay the velocity for the time steps requested in ''steps'' of a transient run',37 'transient_vel':' vel can be by any field of the transient results (vx, vy, vz, vel, temperature, melting, pressure, bed, thickness, surface)',38 'transient_field':' dynamic plot of results. specify ''steps'' option, as fell as ''field'' (defaults are all steps, for ''Vel'' field)',39 'transient_movie':' this will printlay the time steps of a given field of a transient run',40 'transient_movie_field':' field to be printlayed when doing transient_movie data printlay',41 'transient_movie_output':' filename if output is desired for movie',42 'transient_movie_time':' time for each image (default 2 seconds)',43 'thermaltransient_results':' this will printlay all the time steps of a thermal transient run',44 'qmuhistnorm':' histogram normal distribution. needs option qmudata',45 'qmumean':' plot of mean distribution in sampling analysis with scaled response. needs option qmudata for descriptor',46 'qmustddev':' plot of stddev distribution in sampling analysis with scaled response. needs option qmudata for descriptor',47 'part_hist':' partitioning node and area histogram'}12 'basal_dragx or basal_dragy' :' plot a component of the basal drag on the bed (in kPa)', 13 'boundaries':' this will draw all the segment boundaries to the model, including rifts.', 14 'icefront':' this will show segments that are used to define the icefront of the model (Neumann boundary conditions).', 15 'deviatoricstress_tensor':' plot the components of the deviatoric stress tensor (tauxx,tauyy,tauzz,tauxy,tauxz,tauyz, if computed', 16 'deviatoricstress_principal':' plot the deviatoricstress tensor principal axis and principal values', 17 'deviatoricstress_principalaxis1':' arrow plot the first principal axis of the deviatoricstress tensor(replace 1 by 2 or 3 if needed)', 18 'driving_stress':' plot the driving stress (in kPa)', 19 'elements_type':' model used for each element', 20 'vertexnumbering':' numbering of vertices', 21 'highlightelements':' to highlight elements to highlight the element list', 22 'highlightvertices':' to highlight vertices (use highlight option to enter the vertex list', 23 'referential':' stressbalance referential', 24 'riftvel':' velocities along rifts', 25 'riftrelvel':' relative velocities along rifts', 26 'riftpenetration':' penetration levels for a fault', 27 'riftfraction':' fill fractions for every node of the rifts', 28 'rifts':' plot mesh with an offset so that rifts are visible', 29 'strainrate_tensor':' plot the components of the strainrate tensor (exx,eyy,ezz,exy,exz,eyz) if computed', 30 'strainrate_principal':' plot the strainrate tensor principal axis and principal values)', 31 'strainrate_principalaxis1':' arrow plot the first principal axis of the strainrate tensor(replace 1 by 2 or 3 if needed)', 32 'stress_tensor':' plot the components of stress tensor (sxx,syy,szz,sxy,sxz,syz) if computed', 33 'stress_principal':' plot the stress tensor principal axis and principal values', 34 'stress_principalaxis1':' arrow plot the first principal axis of the stress tensor(replace 1 by 2 or 3 if needed)', 35 'transient_results':' this will printlay all the time steps of a transient run (use steps to specify the steps requested)', 36 'transient_vel':' this will printlay the velocity for the time steps requested in ''steps'' of a transient run', 37 'transient_vel':' vel can be by any field of the transient results (vx, vy, vz, vel, temperature, melting, pressure, bed, thickness, surface)', 38 'transient_field':' dynamic plot of results. specify ''steps'' option, as fell as ''field'' (defaults are all steps, for ''Vel'' field)', 39 'transient_movie':' this will printlay the time steps of a given field of a transient run', 40 'transient_movie_field':' field to be printlayed when doing transient_movie data printlay', 41 'transient_movie_output':' filename if output is desired for movie', 42 'transient_movie_time':' time for each image (default 2 seconds)', 43 'thermaltransient_results':' this will printlay all the time steps of a thermal transient run', 44 'qmuhistnorm':' histogram normal distribution. needs option qmudata', 45 'qmumean':' plot of mean distribution in sampling analysis with scaled response. needs option qmudata for descriptor', 46 'qmustddev':' plot of stddev distribution in sampling analysis with scaled response. needs option qmudata for descriptor', 47 'part_hist':' partitioning node and area histogram'} 48 48 49 49 pyoptions={'axis':" show ('on') or hide ('off') axes", … … 77 77 'log':" cutoff value for log", 78 78 'backgroundcolor':" plot background color. RGB tuple or standard string", 79 80 81 82 83 84 85 86 87 88 89 90 'alpha':" set transparency of plotted data (in [0,1]) " 91 79 'expdisp':" path (or list of paths) to the exp file to be plotted ", 80 'explinewidth':" linewidth ", 81 'explinestyle':" matplotlib linestyle string ", 82 'explinecolor':" matplotlib color string ", 83 'expfill':" (True/False) fill a closed contour ", 84 'expfillcolor':" Color for a filled contour, only used if expfill is True ", 85 'expfillalpha':" alpha transparency for filled contour ", 86 'overlay':" True/False. Overlay a georeferenced image (radar/visible) ", 87 'overlay_image':" path to overlay image ", 88 'overlayhist':" plot a histogram of overlay image, used for setting overlaylims ", 89 'overlaylims':" normalized limits to clip and stretch contrast of overlay image (in [0,1], ex. [0.25,0.75]) ", 90 'alpha':" set transparency of plotted data (in [0,1]) ", 91 'layer':"number of the layer to display for 3D runs"} 92 92 93 93 TODOoptions={'basin':" zoom on a given basin ('pineislandglacier','ronneiceshelf', use isbasin to identify a basin", 94 'figurebackgroundcolor':" figure background color. (default is 'none',", 95 'coord':" 'xy' (default) or 'latlon'", 96 'colorbarpos':" [x,y,dx,dy] where x,y,dx and dy are within [0 1]", 97 'colorbarcornerposition':" 'West','North',etc ...", 98 'colorbartitlerotation':" -90, etc ...", 99 'colorbarfontsize':" specify colorbar fontsize", 100 'colorbarwidth':" multiplier (default 1) to the default width colorbar", 101 'colorbarheight':" multiplier (default 1) to the default height colorbar", 102 'contourticks':" 'on' or 'off' to printlay the ticks of the contours", 103 'contouronly':" 'on' or 'off' to printlay the contours on a white background", 104 'contourcolor':" ticks and contour color", 105 'density':" density of quivers (one arrow every N nodes, N integer)", 106 'inset':" add an inset (zoom) of the current figure if 1 (use 'insetx', 'insety' and 'insetpos' to determine the inset position and content)", 107 'insetx':" [min(x) max(x)] where min(x) and max(x) are values determining the inset content", 108 'insety':" [min(y) max(y)] where min(y) and max(y) are values determining the inset content", 109 'insetpos':" [x,y,dx,dy] where x,y,dx and dy are within [0 1]", 110 'resolution':" resolution used by section value (array of type [horizontal_resolution vertical_resolution])", 111 'showsection':" show section used by 'sectionvalue' (string 'on' or a number of labels)", 112 'sectionvalue':" give the value of data on a profile given by an Argus file (string 'Argusfile_name.exp',", 113 'profile':" give the value of data along a vertical profile ([xlocation ylocation])", 114 'smooth':" smooth element data (string 'yes' or integer)", 115 116 'view':" same as standard matlab option (ex: 2, 3 or [90 180]", 117 'zlim':" same as standard matlab option", 118 'xticklabel':" specifiy xticklabel", 119 'yticklabel':" specifiy yticklabel", 120 'contrast':" (default 1) coefficient to add contrast to the radar amplitude image used in overlays", 121 'highres':" resolution of overlayed radar amplitude image (default is 0, high resolution is 1).", 122 'alpha':" transparency coefficient (the higher, the more transparent). Default is 1.5", 123 'scaling':" scaling factor used by quiver plots. Default is 0.4", 124 'autoscale':" set to 'off' to have all the quivers with the same size. Default is 'on'", 125 'linewidth':" line width for expprint plot (use a cell of strings if more than one)", 126 'border':" size of printlay border (in pixels). active only for overlay plots", 127 128 'nan':" value assigned to NaNs (convenient when plotting BC)", 129 'partitionedges':" 'off' by default. overlay plot of partition edges", 130 131 'latlon':" 'on' or {latstep lonstep [resolution [color]]} where latstep,longstep and resolution are in degrees, color is a [r g b] array", 132 'latlonnumbering':" 'on' or {latgap longap colornumber latangle lonangle} where latgap and longap are pixel gaps for the numbers", 133 'latlonclick':" 'on' to click on latlon ticks positions colornumber is a [r g b] array and latangle and lonangle are angles to flip the numbers", 134 'northarrow':" add an arrow pointing north, 'on' for default value or [x0 y0 length [ratio width fontsize]] where (x0,y0) are the coordinates of the base, ratio=headlength/length", 135 'offset':" mesh offset used by 'rifts', default is 500", 136 'scaleruler':" add a scale ruler, 'on' for default value or [x0 y0 length width numberofticks] where (x0,y0) are the coordinates of the lower left corner", 137 'showregion':" show domain in Antarctica on an inset, use 'insetpos' properties", 138 'visible':" 'off' to make figure unvisible, default is 'on'", 139 'wrapping':" repeat 'n' times the colormap ('n' must be an integer)", 140 'unit':" by default, in m, otherwise, 'km' is available", 141 'legend_position':" by default, 'NorthEasth'", 142 'qmudata':" ", 143 'figposition':" position of figure: 'fullscreen', 'halfright', 'halfleft', 'portrait', 'landscape',... (hardcoded in applyoptions.m)", 144 'offsetaxispos':" offset of current axis position to get more space (ex: [-0.02 0 0.04 0])", 145 'axispos':" axis position to get more space", 146 'hmin':" (numeric, minimum for histogram)", 147 'hmax':" (numeric, maximum for histogram)", 148 'hnint':" (numeric, number of intervals for histogram)", 149 'ymin1':" (numeric, minimum of histogram y-axis)", 150 'ymax1':" (numeric, maximum of histogram y-axis)", 151 'ymin2':" (numeric, minimum of cdf y-axis)", 152 'ymax2':" (numeric, maximum of cdf y-axis)", 153 'cdfplt':" (char, 'off' to turn off cdf line plots)", 154 'cdfleg':" (char, 'off' to turn off cdf legends)", 155 'segmentnumbering':" ('off' by default)", 156 'kmlgroundoverlay':" ('off' by default)", 157 'kmlfilename':" ('tempfile.kml' by default)", 158 'kmlroot':" ('./' by default)", 159 'kmlimagename':" ('tempimage' by default)", 160 'kmlimagetype':" ('png' by default)", 161 'kmlresolution':" (1 by default)", 162 'kmlfolder':" ('Ground Overlay' by default)", 163 'kmlfolderdescription':" ('' by default)", 164 'kmlgroundoverlayname':" ('' by default)", 165 'kmlgroundoverlaydescription':"N/A by default')"} 94 'figurebackgroundcolor':" figure background color. (default is 'none',", 95 'coord':" 'xy' (default) or 'latlon'", 96 'colorbarpos':" [x,y,dx,dy] where x,y,dx and dy are within [0 1]", 97 'colorbarcornerposition':" 'West','North',etc ...", 98 'colorbartitlerotation':" -90, etc ...", 99 'colorbarfontsize':" specify colorbar fontsize", 100 'colorbarwidth':" multiplier (default 1) to the default width colorbar", 101 'colorbarheight':" multiplier (default 1) to the default height colorbar", 102 'contourticks':" 'on' or 'off' to printlay the ticks of the contours", 103 'contouronly':" 'on' or 'off' to printlay the contours on a white background", 104 'contourcolor':" ticks and contour color", 105 'density':" density of quivers (one arrow every N nodes, N integer)", 106 'inset':" add an inset (zoom) of the current figure if 1 (use 'insetx', 'insety' and 'insetpos' to determine the inset position and content)", 107 'insetx':" [min(x) max(x)] where min(x) and max(x) are values determining the inset content", 108 'insety':" [min(y) max(y)] where min(y) and max(y) are values determining the inset content", 109 'insetpos':" [x,y,dx,dy] where x,y,dx and dy are within [0 1]", 110 'resolution':" resolution used by section value (array of type [horizontal_resolution vertical_resolution])", 111 'showsection':" show section used by 'sectionvalue' (string 'on' or a number of labels)", 112 'sectionvalue':" give the value of data on a profile given by an Argus file (string 'Argusfile_name.exp',", 113 'profile':" give the value of data along a vertical profile ([xlocation ylocation])", 114 'smooth':" smooth element data (string 'yes' or integer)", 115 'view':" same as standard matlab option (ex: 2, 3 or [90 180]", 116 'zlim':" same as standard matlab option", 117 'xticklabel':" specifiy xticklabel", 118 'yticklabel':" specifiy yticklabel", 119 'contrast':" (default 1) coefficient to add contrast to the radar amplitude image used in overlays", 120 'highres':" resolution of overlayed radar amplitude image (default is 0, high resolution is 1).", 121 'alpha':" transparency coefficient (the higher, the more transparent). Default is 1.5", 122 'scaling':" scaling factor used by quiver plots. Default is 0.4", 123 'autoscale':" set to 'off' to have all the quivers with the same size. Default is 'on'", 124 'linewidth':" line width for expprint plot (use a cell of strings if more than one)", 125 'border':" size of printlay border (in pixels). active only for overlay plots", 126 'nan':" value assigned to NaNs (convenient when plotting BC)", 127 'partitionedges':" 'off' by default. overlay plot of partition edges", 128 'latlon':" 'on' or {latstep lonstep [resolution [color]]} where latstep,longstep and resolution are in degrees, color is a [r g b] array", 129 'latlonnumbering':" 'on' or {latgap longap colornumber latangle lonangle} where latgap and longap are pixel gaps for the numbers", 130 'latlonclick':" 'on' to click on latlon ticks positions colornumber is a [r g b] array and latangle and lonangle are angles to flip the numbers", 131 'northarrow':" add an arrow pointing north, 'on' for default value or [x0 y0 length [ratio width fontsize]] where (x0,y0) are the coordinates of the base, ratio=headlength/length", 132 'offset':" mesh offset used by 'rifts', default is 500", 133 'scaleruler':" add a scale ruler, 'on' for default value or [x0 y0 length width numberofticks] where (x0,y0) are the coordinates of the lower left corner", 134 'showregion':" show domain in Antarctica on an inset, use 'insetpos' properties", 135 'visible':" 'off' to make figure unvisible, default is 'on'", 136 'wrapping':" repeat 'n' times the colormap ('n' must be an integer)", 137 'unit':" by default, in m, otherwise, 'km' is available", 138 'legend_position':" by default, 'NorthEasth'", 139 'qmudata':" ", 140 'figposition':" position of figure: 'fullscreen', 'halfright', 'halfleft', 'portrait', 'landscape',... (hardcoded in applyoptions.m)", 141 'offsetaxispos':" offset of current axis position to get more space (ex: [-0.02 0 0.04 0])", 142 'axispos':" axis position to get more space", 143 'hmin':" (numeric, minimum for histogram)", 144 'hmax':" (numeric, maximum for histogram)", 145 'hnint':" (numeric, number of intervals for histogram)", 146 'ymin1':" (numeric, minimum of histogram y-axis)", 147 'ymax1':" (numeric, maximum of histogram y-axis)", 148 'ymin2':" (numeric, minimum of cdf y-axis)", 149 'ymax2':" (numeric, maximum of cdf y-axis)", 150 'cdfplt':" (char, 'off' to turn off cdf line plots)", 151 'cdfleg':" (char, 'off' to turn off cdf legends)", 152 'segmentnumbering':" ('off' by default)", 153 'kmlgroundoverlay':" ('off' by default)", 154 'kmlfilename':" ('tempfile.kml' by default)", 155 'kmlroot':" ('./' by default)", 156 'kmlimagename':" ('tempimage' by default)", 157 'kmlimagetype':" ('png' by default)", 158 'kmlresolution':" (1 by default)", 159 'kmlfolder':" ('Ground Overlay' by default)", 160 'kmlfolderdescription':" ('' by default)", 161 'kmlgroundoverlayname':" ('' by default)", 162 'kmlgroundoverlaydescription':"N/A by default')"} 166 163 167 164 -
issm/trunk-jpl/src/m/plot/processmesh.py
r21410 r21427 5 5 """ 6 6 PROCESSMESH - process the mesh for plotting 7 7 8 8 Usage: 9 10 9 x,y,z,elements,is2d=processmech(md,data,options) 10 11 11 See also: PLOTMODEL, PROCESSDATA 12 12 """ 13 13 14 14 # {{{ check mesh size parameters 15 15 if md.mesh.numberofvertices==0: 16 16 raise ValueError('processmesh error: mesh is empty') 17 if md.mesh.numberofvertices==md.mesh.numberofelements:18 raise ValueError('processmesh error: the number of elements is the same as the number of nodes')17 if md.mesh.numberofvertices==md.mesh.numberofelements: 18 raise ValueError('processmesh error: the number of elements is the same as the number of nodes') 19 19 # }}} 20 20 # {{{ treating non data plots mesh … … 26 26 x=md.mesh.x 27 27 try: 28 y=md.mesh. x2d28 y=md.mesh.y2d 29 29 except AttributeError: 30 30 y=md.mesh.y … … 53 53 raise ValueError('processmesh error: cannot work with 3D mesh in lat-lon coords') 54 54 #we modify the mesh temporarily to a 2D mesh from which the 3D mesh was extruded 55 #Basile: does not seem necessary as we already picked x as x2d 56 # x=x2d 57 # y=y2d 58 # z=zeros(size(x2d)) 59 # elements=elements2d 55 z=np.zeros(np.size(x)) 56 elements=elements 60 57 else: 61 58 #Process mesh for plotting … … 65 62 # process polycollection here for 3D plot 66 63 is2d=0 67 68 64 #units 69 65 if options.exist('unit'):
Note:
See TracChangeset
for help on using the changeset viewer.
