Changeset 19305
- Timestamp:
- 04/23/15 11:06:28 (10 years ago)
- Location:
- issm/branches/trunk-jpl-ad2/src/m/plot
- Files:
-
- 1 added
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/branches/trunk-jpl-ad2/src/m/plot/prepare_gl.m
r18352 r19305 244 244 245 245 %}}} 246 %handle rendering of contour {{{247 248 249 fid=fopen([database '/contour_coords.js'],'w');250 251 contour_lat=md.mesh.lat(md.mesh.segments(:,1));252 contour_long=md.mesh.long(md.mesh.segments(:,1));253 contour_surface=md.geometry.surface(md.mesh.segments(:,1));254 255 R=6371000*ones(length(contour_surface),1)+scaling_factor*contour_surface;256 257 contour_x = R .* cosd(contour_lat) .* cosd(contour_long);258 contour_y = R .* cosd(contour_lat) .* sind(contour_long);259 contour_z = R .* sind(contour_lat);260 261 %write contour_x:262 fprintf(fid,'<!-- contour_x{{{-->\n');263 fprintf(fid,'var contour_x=[');264 for i=1:length(contour_x)-1,265 fprintf(fid,'%g,',contour_x(i));266 end267 fprintf(fid,'%g];\n',contour_x(end));268 fprintf(fid,'<!--}}}-->\n');269 270 %write contour_y:271 fprintf(fid,'<!-- contour_y{{{-->\n');272 fprintf(fid,'var contour_y=[');273 for i=1:length(contour_y)-1,274 fprintf(fid,'%g,',contour_y(i));275 end276 fprintf(fid,'%g];\n',contour_y(end));277 fprintf(fid,'<!--}}}-->\n');278 279 %write contour_surface:280 fprintf(fid,'<!-- contour_z{{{-->\n');281 fprintf(fid,'var contour_z=[');282 for i=1:length(contour_z)-1,283 fprintf(fid,'%g,',contour_z(i));284 end285 fprintf(fid,'%g];\n',contour_z(end));286 fprintf(fid,'<!--}}}-->\n');287 288 fclose(fid);289 290 %}}}291 246 292 247 %Come back to present directory: -
issm/branches/trunk-jpl-ad2/src/m/plot/writejsfile.m
r18356 r19305 8 8 nods=length(model.x); 9 9 nel=size(model.index,1); 10 nx=length(model.contourx); 10 11 11 12 fid=fopen(filename,'w'); 12 13 13 fprintf(fid,'model = []\n'); 14 14 fprintf(fid,'model = {};\n'); 15 fprintf(fid,'model["title"]="%s";\n',model.title); 16 fprintf(fid,'model["initialZoomFactor"]=%s;\n',model.initialZoomFactor); 17 15 18 %write index: 16 19 fprintf(fid,'<!-- model["index"]{{{-->\n'); … … 26 29 writejsfield(fid,'model["z"]',model.z,nods); 27 30 writejsfield(fid,'model["surface"]',model.surface,nods); 31 writejsfield(fid,'model["contourx"]',model.contourx,nx); 32 writejsfield(fid,'model["contoury"]',model.contoury,nx); 33 writejsfield(fid,'model["contourz"]',model.contourz,nx); 28 34 29 35 results=model.results; 30 fprintf(fid,'results=[]\n'); 31 fprintf(fid,'result=[]\n'); 36 fprintf(fid,'results={};\n'); 32 37 33 38 for i=1:length(results), 39 fprintf(fid,'result={};\n'); 34 40 writejsfield(fid,'result["data"]',results(i).data,nods); 35 41 fprintf(fid,'<!--{{{-->\n'); 36 42 fprintf(fid,'result["caxis"]=[%g,%g];\n',results(i).caxis(1),results(i).caxis(2)); 37 fprintf(fid,'result["label"]="%s"\n',results(i).label); 38 fprintf(fid,'result["shortlabel"]="%s"\n',results(i).shortlabel); 39 fprintf(fid,'result["unit"]="%s"\n',results(i).unit); 43 fprintf(fid,'result["label"]="%s";\n',results(i).label); 44 fprintf(fid,'result["shortlabel"]="%s";\n',results(i).shortlabel); 45 fprintf(fid,'result["unit"]="%s";\n',results(i).unit); 46 if size(results(i).data,2)>1, 47 fprintf(fid,'result["time_range"]=[%g,%g];\n',results(i).time_range(1),results(i).time_range(2)); 48 end 40 49 fprintf(fid,'results["%i"]=result;\n',i); 41 50 fprintf(fid,'<!--}}}-->\n');
Note:
See TracChangeset
for help on using the changeset viewer.