Changeset 19038

Timestamp:

01/26/15 11:32:43 (10 years ago)

Author:

dlcheng

Message:

export_gl

File:

• issm/trunk-jpl/src/m/plot/export_gl.py (modified) (7 diffs)

issm/trunk-jpl/src/m/plot/export_gl.py

@@ -1 @@
-def export_gl(md,varargin)
+def export_gl(md,varargin):
 
         templist=plotoptions(varargin[:]); 
@@ -25 +25 @@
         #Deal with contour {{{
 
-        contour_lat1=md.mesh.lat(md.mesh.segments(      :,1));
+        contour_lat1=md.mesh.lat(md.mesh.segments(:,1));
         contour_lat2=md.mesh.lat(md.mesh.segments(:,2));
         contour_long1=md.mesh.long(md.mesh.segments(:,1));
@@ -35 +35 @@
         R2=6371000*ones(len(contour_surface2),1)+scaling_factor*contour_surface2;
 
-        contourx1 = R1 .* cosd(contour_lat1) .* cosd(contour_long1);
-        contoury1 = R1 .* cosd(contour_lat1) .* sind(contour_long1);
-        contourz1 = R1 .* sind(contour_lat1);
+        contourx1 = R1 * cosd(contour_lat1) * cosd(contour_long1;
+        contoury1 = R1 * cosd(contour_lat1) * sind(contour_long1);
+        contourz1 = R1 * sind(contour_lat1);
         
-        contourx2 = R2 .* cosd(contour_lat2) .* cosd(contour_long2);
-        contoury2 = R2 .* cosd(contour_lat2) .* sind(contour_long2);
-        contourz2 = R2 .* sind(contour_lat2);
+        contourx2 = R2 * cosd(contour_lat2) * cosd(contour_long2);
+        contoury2 = R2 * cosd(contour_lat2) * sind(contour_long2);
+        contourz2 = R2 * sind(contour_lat2);
 
 
@@ -64 +64 @@
         R=6371000*ones(numberofvertices,1)+scaling_factor*surface;
 
-        x = R .* cosd(lat) .* cosd(long);
-        y = R .* cosd(lat) .* sind(long);
-        z = R .* sind(lat);
+        x = R * cosd(lat) * cosd(long);
+        y = R * cosd(lat) * sind(long);
+        z = R * sind(lat);
 
 
@@ -78 +78 @@
         #Deal with data: 
         results=struct([]);
-        for i in xrange(1,len(optionslist)),
+        for i in xrange(1,len(optionslist)):
                 options=optionslist[i]; 
                 options=checkplotoptions(md,options);
@@ -86 +86 @@
 
                 label=getfieldvalue(options,'label','');
-                if strcmpi(label,''),
+                if strcmpi(label,''):
                         #create generic label: 
                         label=['data' num2str(i)];
@@ -92 +92 @@
 
                 shortlabel=getfieldvalue(options,'shortlabel','');
-                if strcmpi(shortlabel,''),
+                if strcmpi(shortlabel,''):
                         #create generic short label: 
                         shortlabel=['data' num2str(i)];
                 results[i].shortlabel=shortlabel;
                 
-                if size(data,2)>1,
+                if len(data[2])>1:
                         time_range=getfieldvalue(options,'time_range',[0 100]);
                         results[i].time_range=time_range;
 
                 unit=getfieldvalue(options,'unit','');
-                if strcmpi(unit,''),
+                if strcmpi(unit,''):
                         #create generic unit: 
                         unit='SI';