Changeset 19036
- Timestamp:
- 01/24/15 16:11:22 (10 years ago)
- Location:
- issm/trunk-jpl/src/m/plot
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/plot/export_gl.py
r19035 r19036 1 1 def export_gl(md,varargin) 2 2 3 templist=plotoptions(varargin {:});3 templist=plotoptions(varargin[:]); 4 4 optionslist=templist.list; 5 options=optionslist {1};5 options=optionslist[1]; 6 6 options=checkplotoptions(md,options); 7 7 … … 14 14 15 15 #Deal with title: 16 if exist(options,'title')16 if options['title']: 17 17 title=getfieldvalue(options,'title'); 18 else 18 else: 19 19 title=''; 20 end21 20 22 21 #initialize model: … … 26 25 #Deal with contour {{{ 27 26 28 contour_lat1=md.mesh.lat(md.mesh.segments( :,1));27 contour_lat1=md.mesh.lat(md.mesh.segments( :,1)); 29 28 contour_lat2=md.mesh.lat(md.mesh.segments(:,2)); 30 29 contour_long1=md.mesh.long(md.mesh.segments(:,1)); … … 33 32 contour_surface2=md.geometry.surface(md.mesh.segments(:,2)); 34 33 35 R1=6371000*ones(len gth(contour_surface1),1)+scaling_factor*contour_surface1;36 R2=6371000*ones(len gth(contour_surface2),1)+scaling_factor*contour_surface2;34 R1=6371000*ones(len(contour_surface1),1)+scaling_factor*contour_surface1; 35 R2=6371000*ones(len(contour_surface2),1)+scaling_factor*contour_surface2; 37 36 38 37 contourx1 = R1 .* cosd(contour_lat1) .* cosd(contour_long1); … … 79 78 #Deal with data: 80 79 results=struct([]); 81 for i=1:length(optionslist), 82 options=optionslist{i}; options=checkplotoptions(md,options); 80 for i in xrange(1,len(optionslist)), 81 options=optionslist[i]; 82 options=checkplotoptions(md,options); 83 83 data=getfieldvalue(options,'data'); 84 results (i).data=data;85 results (i).caxis=getfieldvalue(options,'caxis',[min(data(:)) max(data(:))]);84 results[i].data=data; 85 results[i].caxis=getfieldvalue(options,'caxis',[min(data) max(data)]); 86 86 87 87 label=getfieldvalue(options,'label',''); … … 89 89 #create generic label: 90 90 label=['data' num2str(i)]; 91 results (i).label=label;91 results[i].label=label; 92 92 93 93 shortlabel=getfieldvalue(options,'shortlabel',''); … … 95 95 #create generic short label: 96 96 shortlabel=['data' num2str(i)]; 97 results (i).shortlabel=shortlabel;97 results[i].shortlabel=shortlabel; 98 98 99 99 if size(data,2)>1, 100 100 time_range=getfieldvalue(options,'time_range',[0 100]); 101 results (i).time_range=time_range;101 results[i].time_range=time_range; 102 102 103 103 unit=getfieldvalue(options,'unit',''); … … 105 105 #create generic unit: 106 106 unit='SI'; 107 results (i).unit=unit;107 results[i].unit=unit; 108 108 model.results=results; 109 109 -
issm/trunk-jpl/src/m/plot/writejsfield.py
r19035 r19036 7 7 #write array: 8 8 if len(variable[2])==1: 9 fid.write('<!-- {0}{{{-->\n'.format(name)) ;10 fid.write('{0}=['.format(name)) ;9 fid.write('<!-- {0}{{{-->\n'.format(name)) 10 fid.write('{0}=['.format(name)) 11 11 for i in xrange(1, nods-1): 12 fid.write('{0},',variable[i]) ;13 fid.write('{0}];\n',variable[-1]) ;14 fid.write('<!--}}}-->\n') ;12 fid.write('{0},',variable[i]) 13 fid.write('{0}];\n',variable[-1]) 14 fid.write('<!--}}}-->\n') 15 15 else: 16 16 #multi-sized array: 17 fid.write('<!-- {0}{{{-->\n'.format(name)) ;18 fid.write('{0}=[]\n'.format(name)) ;17 fid.write('<!-- {0}{{{-->\n'.format(name)) 18 fid.write('{0}=[]\n'.format(name)) 19 19 for i in xrange(1, len(variable[2])): 20 fid.write('{0}["{1}"]=['.format(name,i)) ;20 fid.write('{0}["{1}"]=['.format(name,i)) 21 21 for j in xrange(1, nods-1): 22 fid.write('{0},'.format(variable[j,i])) ;23 fid.write('{0}];\n'.format(variable[-1,i])) ;24 fid.write('<!--}}}-->\n') ;22 fid.write('{0},'.format(variable[j,i])) 23 fid.write('{0}];\n'.format(variable[-1,i])) 24 fid.write('<!--}}}-->\n') -
issm/trunk-jpl/src/m/plot/writejsfile.py
r19035 r19036 6 6 # 7 7 8 nods=len gth(model.x);9 nel= size(model.index,1);10 nx=len gth(model.contourx1);8 nods=len(model.x) 9 nel=len(model.index[1]) 10 nx=len(model.contourx1) 11 11 12 fid=open(filename,'w') ;12 fid=open(filename,'w') 13 13 14 fid.write('model = {};\n') ;15 fid.write('model["title"]="{0}";\n'.format(model.title) ;16 fid.write('model["initialZoomFactor"]={0};\n'.format(model.initialZoomFactor) ;14 fid.write('model = {};\n') 15 fid.write('model["title"]="{0}";\n'.format(model.title) 16 fid.write('model["initialZoomFactor"]={0};\n'.format(model.initialZoomFactor) 17 17 18 18 #write index: 19 fid.write('<!-- model["index"]{{{-->\n'); 20 fid.write('model["index"]=['); 21 for i=1:nel-1, 22 fid.write('[{0}, {1}, {2}],'.format(model.index(i,1),model.index(i,2),model.index(i,3)); 23 end 24 fid.write('[{0}, {1}, {2}]];\n'.format(model.index(end,1),model.index(end,2),model.index(end,3)); 25 fid.write('<!--}}}-->\n'); 19 fid.write('<!-- model["index"]{{{-->\n') 20 fid.write('model["index"]=[') 21 for i in xrange(1,nel-1): 22 fid.write('[{0}, {1}, {2}],'.model.index[i][1],model.index[i][2],model.index[i][3]) 23 fid.write('[{0}, {1}, {2}]];\n'.format(model.index[-1][1],model.index[-1][2],model.index[-1][3]) 24 fid.write('<!--}}}-->\n') 26 25 27 writejsfield(fid,'model["x"]',model.x,nods) ;28 writejsfield(fid,'model["y"]',model.y,nods) ;29 writejsfield(fid,'model["z"]',model.z,nods) ;30 writejsfield(fid,'model["surface"]',model.surface,nods) ;31 writejsfield(fid,'model["contourx1"]',model.contourx1,nx) ;32 writejsfield(fid,'model["contoury1"]',model.contoury1,nx) ;33 writejsfield(fid,'model["contourz1"]',model.contourz1,nx) ;34 writejsfield(fid,'model["contourx2"]',model.contourx2,nx) ;35 writejsfield(fid,'model["contoury2"]',model.contoury2,nx) ;36 writejsfield(fid,'model["contourz2"]',model.contourz2,nx) ;26 writejsfield(fid,'model["x"]',model.x,nods) 27 writejsfield(fid,'model["y"]',model.y,nods) 28 writejsfield(fid,'model["z"]',model.z,nods) 29 writejsfield(fid,'model["surface"]',model.surface,nods) 30 writejsfield(fid,'model["contourx1"]',model.contourx1,nx) 31 writejsfield(fid,'model["contoury1"]',model.contoury1,nx) 32 writejsfield(fid,'model["contourz1"]',model.contourz1,nx) 33 writejsfield(fid,'model["contourx2"]',model.contourx2,nx) 34 writejsfield(fid,'model["contoury2"]',model.contoury2,nx) 35 writejsfield(fid,'model["contourz2"]',model.contourz2,nx) 37 36 38 37 39 results=model.results ;40 fid.write('results={};\n') ;38 results=model.results 39 fid.write('results={};\n') 41 40 42 for i =1:length(results):43 fid.write('result={};\n') ;44 writejsfield(fid,'result["data"]',results (i).data,nods);45 fid.write('<!--{{{-->\n') ;46 fid.write('result["caxis"]=[{0},{1}];\n'.format(results (i).caxis(1),results(i).caxis(2));47 fid.write('result["label"]="{0}";\n'.format(results (i).label);48 fid.write('result["shortlabel"]="{0}";\n'.format(results (i).shortlabel);49 fid.write('result["unit"]="{0}";\n'.format(results (i).unit);50 if size(results(i).data,2)>1:51 fid.write('result["time_range"]=[{0},{1}];\n'.format(results (i).time_range(1),results(i).time_range(2));52 fid.write('results["{0}"]=result;\n'.format (i);53 fid.write('<!--}}}-->\n') ;54 fid.write('model.results=results;\n') ;55 fid.write('models["{0}"]=model;\n'.format(keyname) ;41 for i in xrange(1,len(results)): 42 fid.write('result={};\n') 43 writejsfield(fid,'result["data"]',results[i].data,nods) 44 fid.write('<!--{{{-->\n') 45 fid.write('result["caxis"]=[{0},{1}];\n'.format(results[i].caxis[1],results[i].caxis[2]) 46 fid.write('result["label"]="{0}";\n'.format(results[i].label) 47 fid.write('result["shortlabel"]="{0}";\n'.format(results[i].shortlabel) 48 fid.write('result["unit"]="{0}";\n'.format(results[i].unit) 49 if len(results[i].data[2])>1: 50 fid.write('result["time_range"]=[{0},{1}];\n'.format(results[i].time_range[1],results[i].time_range[2]) 51 fid.write('results["{0}"]=result;\n'.format[i] 52 fid.write('<!--}}}-->\n') 53 fid.write('model.results=results;\n') 54 fid.write('models["{0}"]=model;\n'.format(keyname) 56 55 57 fid.close() ;56 fid.close()
Note:
See TracChangeset
for help on using the changeset viewer.