fixed a bug in plm2rnd demo2

csdms-contrib · Sep 29, 2020 · 1a003ad · 1a003ad
1 parent 3abd8b8
commit 1a003ad
Showing 1 changed file with 7 additions and 4 deletions.
diff --git a/plm2rnd.m b/plm2rnd.m
@@ -32,6 +32,7 @@
 % plm2rnd('demo1') % Verifies input and output
 % plm2rnd('demo2') % Compares EGM96 to random structures
 %
+% Last modified by charig-at-arizona.edu, 09/29/2020
 % Last modified by fjsimons-at-alum.mit.edu, 03/13/2013
 
 defval('bta',-4.0361)
@@ -91,21 +92,23 @@
     axes(ah(1))
     plotplm(egm(4:end,:),[],[],1)
     axes(ah(2))
-    dbt=plotplm(egm,[],[],3);
-
+    %dbt=plotplm(egm,[],[],3);
+    egm = [0 0 0 0; 1 0 0 0; 1 1 0 0; egm(:,1:4)];
+    plotplm(egm,[],[],3);
+
     syn1=plm2rnd(max(egm(:,1)),bta);
     syn2=plm2rnd(max(egm(:,1)),bta);
 
     axes(ah(3))
     plotplm(syn1(7:end,:),[],[],1)
 
     axes(ah(4))
-    dbt1=plotplm(syn1(4:end,:),[],[],3);
+    dbt1=plotplm(syn1,[],[],3);
 
     axes(ah(5))
     plotplm(syn2(7:end,:),[],[],1)
     axes(ah(6))
-    dbt2=plotplm(syn1(4:end,:),[],[],3);
+    dbt2=plotplm(syn2,[],[],3);
     figdisp([],'demo2')
   end
 end