# file calccr.m ; 
# purpose: to create a minimal script file
# that shows problems with plotting;

# command = diary on, echo on, more off, calccr;
   [now] = ctime(time)

# data to plot; from ccq.dat with strings "ccq" replaced by "ccr"
load  "ccr.dat"

 graphics_toolkit("gnuplot")
 
 clf
  h0 = figure(1)
  set(h0, "paperposition", [0.25, 0.25, 8.0, 10.5] )

  hax1 = axes()
  axis( "equal" )
  set(hax1, "fontweight","bold", "fontsize", 12 )
  set (hax1,"position",[.1, .5, .4, .444444 ] )
  set(hax1, "xlim", [0.94, 1.03 ] , "ylim", [ 0.96, 1.06 ] )   
  set(hax1, "xtick", [.94:.02:1.03] )
  set(hax1, "xticklabel", [])
  set(hax1, "ytick", [.96:.02:1.06] )
  set(hax1, "yticklabel", {".96",".98","1.00","1.02","1.04","1.06"})
  set(hax1, "xgrid", "on" )
  set(hax1, "ygrid", "on" )
  ylabel(ylabelstr1 , "fontsize",10,"fontweight","bold")

# doing the upper left plot
    # doing the curve for the legend, the first one;

     # duf 4 points:
  hl11=line(nomdufex(kk=[1,2,4,5]),nadufex(kk),"linewidth",3,"markersize",4,"linestyle","none","marker","+");
# only 1 "line" above appears in legend;
  hleg1= legend({"duf"})
  legend location northwest
   set(hleg1,"fontsize",12,"fontweight","bold")
   legend boxoff   
 # from here down plotting doesn't go into legend;
     # duf base point:
  hl12=line(nomdufex(3),nadufex(3),"linewidth",4,"markersize",7,"marker","+");
     # duf line:
  hl13=line(om0/omref, a0/aref, "linewidth",1,"linestyle",":"    );
     # vw axes origin:
  hl14=line(vwh1om,vwh1a,"linewidth",3,"marker","x","markersize",6);
     # om* axis:
  hl123=line([xtail1,xhead1],[ytail1,yhead1],"linewidth",2,"linestyle","-"    );
  hl124=line(xarrow1,yarrow1,"linewidth",4,"linestyle","-"    );
     # a* axis:
  hl125=line([xtail2,xhead2],[ytail2,yhead2],"linewidth",2,"linestyle","-"    );
  hl126=line(xarrow2,yarrow2,"linewidth",4,"linestyle","-"    );
     # v axis:
  hl127=line([xtail3,xhead3],[ytail3,yhead3],"linewidth",2,"linestyle","-"    );
  hl128=line(xarrow3,yarrow3,"linewidth",4,"linestyle","-"    );
     # w axis:
  hl129=line([xtail4,xhead4],[ytail4,yhead4],"linewidth",2,"linestyle","-"    );
  hl130=line(xarrow4,yarrow4,"linewidth",4,"linestyle","-"    );
  hl132=line(xarrow15,yarrow15,"linewidth",3,"linestyle","-" );
    # dimension arc/arrow for al:
  hl137=line(xarc,yarc,"linewidth",1,"linestyle","-" );
  hl138=line(xarrow20,yarrow20,"linewidth",4,"linestyle","-" );

  stromstar = "{/Symbol w}*"
  text ( 0.640, 0.275, stromstar, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  strastar = "a*"
  text ( 0.450, 0.422, strastar, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  strv = "v"
  text ( 0.644, 0.350, strv, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  strw = "w"
  text ( 0.311, 0.396, strw, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  strdelstar = "{/Symbol d}*"
  text ( 0.440, 0.223, strdelstar, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  stralpha = "{/Symbol a}"
  text ( 0.556, 0.276, stralpha, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

# doing the upper right plot:
  hax2 = axes()
  axis( "equal" )
  set(hax2, "fontweight","bold", "fontsize", 10 )
  set (hax2,"position",[.5, .5, .4, .444444 ] )
  set (hax2, "xlim",[.53, .62],"ylim", [ 0.96, 1.06 ] ) 
  set(hax2, "ytick", [.96:.02:1.06])
  set (hax2, "yticklabel", [] )
  set( hax2, "xtick", [.54:.02:.60] )
  set (hax2, "xticklabel", {".54",".56",".58",".60"} )
  grid on ;
  xlabelstr2 = "Phase Lag, {/Symbol d}*"
  xlabel(xlabelstr2 , "fontsize",10,"fontweight","bold")

     # duf 4 points:
  hl21=line(ndeldufex(kk),nadufex(kk),"linewidth",2,"markersize",3,"linestyle","none","marker","+")
     # duf base point:
  hl23=line(ndeldufex(3),nadufex(3),"linewidth",4,"markersize",7,"linestyle","none","marker","+")
     # duf line:
  hl24=line(del0/delref, a0/aref, "linewidth",1,"linestyle",":" );
     # vw axes origin:
  hl25=line(vwh2del,vwh2a,"linewidth",3,"marker",".","markersize",14);
     # delta* axis:
  hl232=line([xtail5,xhead5],[ytail5,yhead5],"linewidth",2,"linestyle","-"    );
  hl233=line(xarrow5,yarrow5,"linewidth",4,"linestyle","-" );
     # a* axis:
  hl234=line([xtail6,xhead6],[ytail6,yhead6],"linewidth",2,"linestyle","-"    );
  hl235=line(xarrow6,yarrow6,"linewidth",4,"linestyle","-" );

  stromstar = "{/Symbol w}*"
  text ( 0.560, 0.223, stromstar, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  strdelstar = "{/Symbol d}*"
  text ( 0.760, 0.275, strdelstar, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  text ( 0.558, 0.422, strastar, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

 # doing the lower left plot:
  hax3 = axes()
  axis( "equal" )
  set(hax3, "ydir", "reverse" )
  set(hax3, "fontweight","bold", "fontsize", 12 )
  set (hax3,"position",[.1, .20  , .4, .4  ] )
  set (hax3, "xlim",[0.94,1.03] )
  set (hax3, "ylim",[.53,.62])
  set( hax3, "ytick", [.54:.02:.62] )
  set (hax3, "yticklabel", {".54",".56",".58",".60",".62"} )
  set( hax3, "xtick", [.94:.02:1.02] )
  set (hax3, "xticklabel", [{".94",".96",".98","1.00","1.02"}])
  ylabelstr3 = "Phase Lag, {/Symbol d}* \\equiv  {/Symbol d}/{/Symbol d}_{ref} "
  ylabel(ylabelstr3 , "fontsize",12,"fontweight","bold")
  xlabelstr3 = "Frequency"
  xlabel(xlabelstr3 , "fontsize",12,"fontweight","bold")
  grid on ;

     # duf 4 points:
  hl31=line(nomdufex(kk),ndeldufex(kk),"linewidth",2,"markersize",4,"linestyle","none","marker","+");
     # duf base point:
  hl33=line(nomdufex(3),ndeldufex(3),"linewidth",4,"markersize",7,"linestyle","none","marker","+")
     # duf line:
  hl34=line(om0/omref, del0/delref, "linewidth",1,"linestyle",":" );
     # vw axes origin:
  hl35=line(vwh3om,vwh3del,"linewidth",3,"marker","x","markersize",6);
     # om* axis:
  hl323=line([xtail7,xhead7],[ytail7,yhead7],"linewidth",2,"linestyle","-");
  hl324=line(xarrow7,yarrow7,"linewidth",4,"linestyle","-" );
     # delta* axis:
  hl325=line([xtail8,xhead8],[ytail8,yhead8],"linewidth",2,"linestyle","-");
  hl326=line(xarrow8,yarrow8,"linewidth",4,"linestyle","-" );

    xlabelstr3 = strcat("Dimensionless Frequency,"\
            , " {{/Symbol w}* {\\equiv} {/Symbol w}/{/Symbol w}_{ref} }")
  xlabel(xlabelstr3 , "fontsize",12,"fontweight","bold")

  strdelstar = "{/Symbol d}*"
  text ( 0.442, 0.246, strdelstar, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  strastar = "a*"
  text ( 0.406, 0.440, strastar, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  stromstar = "{/Symbol w}*"
  text ( 0.640, 0.468, stromstar, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  # doing the leftmost part of the lower right plot:

  hax4 = axes()
  axis( "equal" )
  set(hax4, "fontweight","bold", "fontsize", 12 )
  set (hax4,"position",[.1  , .149, .80, .84 ],"color","none" )
  set(hax4, "xlim", [0.94, 1.12 ] )
  set(hax4, "ylim", [0.87, 1.06 ] )
  set( hax4, "ytick", [] )
  set( hax4, "xtick", [] )
  set (hax4, "xticklabel", [] )
  get(hax4,"position")
  xlabelstr4 = strcat("plotccr.ps,","\t",  now )
  text ( 0.75, 0.06, xlabelstr4, "units","normalized"...
        , "horizontalalignment","center"...
        , "verticalalignment", "baseline"...
        , "fontweight",  "bold"...
        , "interpreter", "tex"...
        , "fontsize", 12 ...
       ) 

f_amplstr = "1" 
m_massstr = "1"
c_viscstr = "0.2"
k_sprstr  = "0"
h_sprstr  = "1"
omrefstr  = "1.5"

    titlestr = strcat("Equally Scaled Resonance Curve Projections for:"...
            , "\n"
            , "f = "...
            , "\t"...
            , f_amplstr...
            , ", \t"...
            , " m = \t"...
            , m_massstr...
            , ", \t"...
            , " c = \t"...
            , c_viscstr...
            , ", \t"...
            , " k = \t"...
            , k_sprstr...
            , ", \t"...
            , " and"...
            , "\t"...
            , " h = \t"...
            , h_sprstr...
            , " on RES Branch;"...
            , "\n"
            , "{/Symbol w}_{ref} = "...
            , "\t"...
            , omrefstr...
            , ", \t"...
            , " a_{ref} \\approx "...
            , "\t"...
            , arefstr...
            , ", \t"...
            , "( {/Symbol d}_{a}  \\approx "...
            , "\t"...
            , ndelduf3str...
            , ") \t"...
            )
  text ( 0.5, 0.94 ...
        , titlestr ...
        , "units","normalized"...
        , "horizontalalignment","center"...
        , "verticalalignment", "baseline"...
        , "fontweight",  "bold"...
        , "interpreter", "tex"...
        , "fontsize", 12 ...
       ) 
     # vw orig, connecting grid line:
  hl42=line(xorigin, yorigin, "linewidth", .5,"linestyle",":","markersize",3);
     # vw orig:
  hl43=line(vwh4om, vwh4a, "linewidth",3,"marker","x","linestyle","none","markersize",6);
     # vw orig2, connecting grid line:
  hl44=line(xorigin2, yorigin2, "linewidth",.5,"marker","x","linestyle",":","markersize",3);
     # vw orig2:
  hl45=line(vwh4om2, vwh4a2, "linewidth",3,"marker",".","linestyle","none","markersize",14);
     # bp (duf):
  hl46=line(bph4om,bph4a,"linewidth",4,"marker","+","linestyle","none","markersize",9);
     # bp2 (duf):
  hl47=line(bph4om2,bph4a2,"linewidth",4,"marker","+","linestyle","none","markersize",9);
     # bp-bp2 vertical boundary line:
  hl48=line(xblside,yblside,"linewidth",3,"linestyle","-");
     # bp-bp2 top boundary line:
  hl49=line(xbltop,ybltop,"linewidth",3,"linestyle","-");
     # bp-bp2 diag boundary line:
  hl410=line(xbldiag,ybldiag,"linewidth",3,"linestyle","-");
     # bp4ex(duf):
  hl411=line(bp4exh4om,bp4exh4a,"linewidth",2,"marker","+","linestyle","none","markersize",4);
     # duf line:
  hl412=line(dufh4om,dufh4a,"linewidth",1,"linestyle",":");
     # grid lines:
  hl427=line(xgridh4(:,1),ygridh4(:,1),"linewidth",0.5,"linestyle",":");
  hl428=line(xgridh4(:,2),ygridh4(:,2),"linewidth",0.5,"linestyle",":");
  hl429=line(xgridh4(:,3),ygridh4(:,3),"linewidth",0.5,"linestyle",":");
  hl430=line(xgridh4(:,4),ygridh4(:,4),"linewidth",0.5,"linestyle",":");
  hl431=line(xgridh4(:,5),ygridh4(:,5),"linewidth",0.5,"linestyle",":");
     # v axis:
  hl441=line([xtail9,xhead9],[ytail9,yhead9],"linewidth",2,"linestyle","-");
  hl442=line(xarrow9,yarrow9,"linewidth",4,"linestyle","-" );
     # del axis:
  hl443=line([xtail10,xhead10],[ytail10,yhead10],"linewidth",2,"linestyle","-");
  hl444=line(xarrow10,yarrow10,"linewidth",4,"linestyle","-" );
     # 't' axis:
  hl445=line([xtail11,xhead11],[ytail11,yhead11],"linewidth",2,"linestyle","-");
  hl446=line(xarrow11,yarrow11,"linewidth",4,"linestyle","-" );
     # 'u' axis:
  hl447=line([xtail12,xhead12],[ytail12,yhead12],"linewidth",2,"linestyle","-");
  hl448=line(xarrow12,yarrow12,"linewidth",4,"linestyle","-" );
     # 'u2' axis:
  hl449=line([xtail13,xhead13],[ytail13,yhead13],"linewidth",2,"linestyle","-");
  hl450=line(xarrow13,yarrow13,"linewidth",4,"linestyle","-" );
     # 'w' axis:
  hl451=line([xtail14,xhead14],[ytail14,yhead14],"linewidth",2,"linestyle","-");
  hl452=line(xarrow14,yarrow14,"linewidth",4,"linestyle","-" );
    # dimension arc/arrow for bet:
  hl459=line(xarc21,yarc21,"linewidth",1,"linestyle","-" );
  hl460=line(xarrow21,yarrow21,"linewidth",4,"linestyle","-" );

  str56 = ".56"
  text ( 0.612, 0.404, str56, "units","normalized"...
        , "horizontalalignment","center"...
        , "verticalalignment", "baseline"...
        , "fontweight",  "bold"...
        , "interpreter", "tex"...
        , "fontsize", 10 ...
       ) 

  str58 = ".58"
  text ( 0.800, 0.404, str58, "units","normalized"...
        , "horizontalalignment","center"...
        , "verticalalignment", "baseline"...
        , "fontweight",  "bold"...
        , "interpreter", "tex"...
        , "fontsize", 10 ...
       ) 

  str60 = ".60"
  text ( 0.982, 0.404, str60, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       ) 

  str62 = ".62"
  text ( 0.782, 0.200, str62, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       ) 

  strn = "u"
  text ( 0.526, 0.344, strn, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       ) 

  strv = "v"
  text ( 0.598, 0.312, strv, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       ) 

  strt = "t"
  text ( 0.620, 0.239, strt, "units","normalized"\
        , "horizontalalignment","center"\
        , "verticalalignment", "baseline"\
        , "fontweight",  "bold"\
        , "interpreter", "tex"\
        , "fontsize", 10\
       ) 

  strdelstar = "{/Symbol d}*"
  text ( 0.586, 0.189, strdelstar, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       )

  strw = "w"
  text ( 0.512, 0.248, strw, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       ) 

  strt = "t"
  text ( 0.809, 0.214, strt, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       ) 

  strn = "u"
  text ( 0.834, 0.296, strn, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       ) 

  strw = "w"
  text ( 0.918, 0.208, strw, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       ) 

  strbeta = "{/Symbol b}"
  text ( 0.560, 0.240, strbeta, "units","normalized" ...
        , "horizontalalignment","center" ...
        , "verticalalignment", "baseline" ...
        , "fontweight",  "bold" ...
        , "interpreter", "tex" ...
        , "fontsize", 10 ...
       )

  print(h0, "plotccr.ps")

  [now] = ctime(time)

 uname ()
 version

  diary off;