BookGPU/Chapters/chapter7/figures/.svn/text-base/Series60CB06HalfModelforcex02.m.svn-base

   1 clear all; close all
   2
   3 path = '/Users/apek/SVN/bookgpu/Chapters/chapter7/figures/';
   4 printONOFF = 1;
   5
   6 %set grid
   7 %set xlabel "Fr"
   8 %set ylabel "Cw=Rw/(0.5*rho*S*U^2)"
   9 %L = 3.048
  10 %Fsc(U,S) = 0.5*1000*S*U*U
  11 data0   = load('Series60CB06HalfModelforcex.dat');
  12 froude  = data0(:,1);
  13 time    = data0(:,2);
  14 forcex  = data0(:,3);
  15
  16 U = froude*sqrt(9.81);
  17 S = 0.0797;
  18
  19 force_scale = 0.5*1000*S*U.*U;
  20
  21 data1       =   load('S60b_resistance.dat');
  22 froude1     =   data1(:,3);
  23 forcex1     =   data1(:,7);
  24 plot(froude,forcex./force_scale,'ks',froude1,forcex1,'kx')
  25 grid on
  26
  27 xlabel('Fr')
  28 ylabel('Cw=Rw/(0.5*rho*S*U^2)')
  29 legend('N=513x257x9','Cr, Toda et al 1992')
  30
  31 %plot "Series60CB06HalfModelforcex.dat" using 1:(2.0*abs($3)/Fsc(($1)*sqrt(9.81),2.0*0.0797)) title "N=513x257" pointsize 2 linewidth 2 lc rgb "black",\
  32 %     "/home/ole/postdoc/Series60/S60b_resistance.dat" using 3:7 title "Cr, Toda et al 1992" linewidth 2 lc rgb "gray"
  33
  34 %set term postscript eps enhanced color
  35 %set output "Series60CB06Cw.eps"
  36 %replot
  37
  38 %set terminal png enhanced large
  39 %set output "Series60CB06Cw.png"
  40 %replot
  41
  42
  43 if printONOFF
  44     filename = 'figSeries60CB06Type7Resistance.eps';
  45     title('')
  46     disp(filename)
  47     print('-depsc',[path filename])
  48     disp('saved!')
  49 end
  50