Changeset 24214 for issm/trunk-jpl/test/NightlyRun/test1205.py

Timestamp:

10/11/19 00:27:00 (5 years ago)

Author:

bdef

Message:

CHG: syntax cahnge to meet most of Pep8 requirement

File:

• issm/trunk-jpl/test/NightlyRun/test1205.py (modified) (5 diffs)

issm/trunk-jpl/test/NightlyRun/test1205.py

@@ -22 +22 @@
 md = model()
 md = roundmesh(md, 750000., resolution)
-md = setmask(md, '', '')    #We can not test iceshelves nor ice rises with this analytical solution
+md = setmask(md, '', '')  #We can not test iceshelves nor ice rises with this analytical solution
 md = parameterize(md, '../Par/RoundSheetStaticEISMINT.py')
 
 #Calculation of the analytical 2d velocity field
 constant = 0.3
-vx_obs = constant / 2. * md.mesh.x * (md.geometry.thickness)**-1
-vy_obs = constant / 2. * md.mesh.y * (md.geometry.thickness)**-1
+vx_obs = constant / 2. * md.mesh.x * (md.geometry.thickness)**- 1
+vy_obs = constant / 2. * md.mesh.y * (md.geometry.thickness)**- 1
 vel_obs = np.sqrt((md.inversion.vx_obs)**2 + (md.inversion.vy_obs)**2)
 
@@ -58 +58 @@
 #Plot of the velocity from the exact and calculated solutions
 #figure(1)
-#set(gcf, 'Position',[1 1 1580 1150])
+#set(gcf, 'Position', [1 1 1580 1150])
 #subplot(2, 2, 1)
-#p = patch('Faces', md.mesh.elements2d, 'Vertices',[md.mesh.x2d md.mesh.y2d], 'FaceVertexCData',...
+#p = patch('Faces', md.mesh.elements2d, 'Vertices', [md.mesh.x2d md.mesh.y2d], 'FaceVertexCData', ...
 #vel, 'FaceColor', 'interp', 'EdgeColor', 'none')
 #title('Modelled velocity', 'FontSize', 14, 'FontWeight', 'bold')
@@ -67 +67 @@
 
 #subplot(2, 2, 2)
-#p = patch('Faces', md.mesh.elements2d, 'Vertices',[md.mesh.x2d md.mesh.y2d], 'FaceVertexCData',...
+#p = patch('Faces', md.mesh.elements2d, 'Vertices', [md.mesh.x2d md.mesh.y2d], 'FaceVertexCData', ...
 #vel_obs, 'FaceColor', 'interp', 'EdgeColor', 'none')
 #title('Analytical velocity', 'FontSize', 14, 'FontWeight', 'bold')
@@ -75 +75 @@
 #subplot(2, 2, 3)
 #hold on
-#plot(sqrt((md.mesh.x(1:md.mesh.numberofvertices2d)).^2+(md.mesh.y(1:md.mesh.numberofvertices2d)).^2), vel, 'r.')
-#plot(sqrt((md.mesh.x2d).^2+(md.mesh.y2d).^2), vel_obs, 'b.')
+#plot(sqrt((md.mesh.x(1:md.mesh.numberofvertices2d)).^2 + (md.mesh.y(1:md.mesh.numberofvertices2d)).^2), vel, 'r.')
+#plot(sqrt((md.mesh.x2d).^2 + (md.mesh.y2d).^2), vel_obs, 'b.')
 #title('Analytical vs calculated velocity', 'FontSize', 14, 'FontWeight', 'bold')
 #xlabel('distance to the center of the icesheet [m]', 'FontSize', 14, 'FontWeight', 'bold')
-#ylabel('velocity [m/yr]', 'FontSize', 14, 'FontWeight', 'bold')
+#ylabel('velocity [m / yr]', 'FontSize', 14, 'FontWeight', 'bold')
 #legend('calculated velocity', 'exact velocity')
 #axis([0 750000 0 200])
@@ -85 +85 @@
 
 #subplot(2, 2, 4)
-#p = patch('Faces', md.mesh.elements2d, 'Vertices',[md.mesh.x2d md.mesh.y2d], 'FaceVertexCData',...
-#abs(vel-vel_obs)./vel_obs*100, 'FaceColor', 'interp', 'EdgeColor', 'none')
+#p = patch('Faces', md.mesh.elements2d, 'Vertices', [md.mesh.x2d md.mesh.y2d], 'FaceVertexCData', ...
+#abs(vel - vel_obs). / vel_obs * 100, 'FaceColor', 'interp', 'EdgeColor', 'none')
 #title('Relative misfit [%]', 'FontSize', 14, 'FontWeight', 'bold')
 #colorbar