Changeset 19305


Ignore:
Timestamp:
04/23/15 11:06:28 (10 years ago)
Author:
Eric.Larour
Message:

CHG: minor tweaks and improvements

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  
    244244
    245245%}}}
    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 end
    267 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 end
    276 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 end
    285 fprintf(fid,'%g];\n',contour_z(end));
    286 fprintf(fid,'<!--}}}-->\n');
    287 
    288 fclose(fid);
    289 
    290 %}}}
    291246
    292247%Come back to present directory:
  • issm/branches/trunk-jpl-ad2/src/m/plot/writejsfile.m

    r18356 r19305  
    88        nods=length(model.x);
    99        nel=size(model.index,1);
     10        nx=length(model.contourx);
    1011       
    1112        fid=fopen(filename,'w');
    1213
    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       
    1518        %write index:
    1619        fprintf(fid,'<!-- model["index"]{{{-->\n');
     
    2629        writejsfield(fid,'model["z"]',model.z,nods);
    2730        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);
    2834
    2935        results=model.results;
    30         fprintf(fid,'results=[]\n');
    31         fprintf(fid,'result=[]\n');
     36        fprintf(fid,'results={};\n');
    3237
    3338        for i=1:length(results),
     39                fprintf(fid,'result={};\n');
    3440                writejsfield(fid,'result["data"]',results(i).data,nods);
    3541                fprintf(fid,'<!--{{{-->\n');
    3642                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
    4049                fprintf(fid,'results["%i"]=result;\n',i);
    4150                fprintf(fid,'<!--}}}-->\n');
Note: See TracChangeset for help on using the changeset viewer.