Changeset 21656

Timestamp:

04/04/17 04:53:50 (8 years ago)

Author:

bdef

Message:

NEW: adding hydrology spc to BC plot

File:

• issm/trunk-jpl/src/m/plot/plot_BC.py (modified) (1 diff)

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

@@ -44 +44 @@
                                                          y[np.where(~np.isnan(md.stressbalance.spcvz))],
                                                          marker='o',c='y',s=80,label='vz Dirichlet',linewidth=0)
-                ax.scatter(x[np.where(~np.isnan(md.hydrology.spcepl_head))],
-                                                         y[np.where(~np.isnan(md.hydrology.spcepl_head))],
-                                                         marker='v',c='r',s=240,label='EPL Head',linewidth=0)
-                ax.scatter(x[np.where(~np.isnan(md.hydrology.spcsediment_head))],
-                                                         y[np.where(~np.isnan(md.hydrology.spcsediment_head))],
-                                                         marker='^',c='b',s=240,label='IDS head',linewidth=0)
-                
+                try:
+                        ax.scatter(x[np.where(~np.isnan(md.hydrology.spcepl_head))],
+                                                                 y[np.where(~np.isnan(md.hydrology.spcepl_head))],
+                                                                 marker='v',c='r',s=240,label='EPL Head',linewidth=0)
+                        ax.scatter(x[np.where(~np.isnan(md.hydrology.spcsediment_head))],
+                                                                 y[np.where(~np.isnan(md.hydrology.spcsediment_head))],
+                                                                 marker='^',c='b',s=240,label='IDS head',linewidth=0)
+                except AttributeError:
+                        print ('Not treating Hydrologydc, skipping these boundaries')
                 ax.set_xlim(XLims)
                 ax.set_ylim(YLims)