#
# $Id: poldat.dem,v 1.7 2003/10/28 05:35:54 sfeam Exp $
#

print = 0
#print = 1
table = 0
#table = 1

set key left bottom
if(print) set term post landscape enhanced color solid defaultplex "HelveticaBold" 12; set out'CP_CFD_Polar.ps'



#******Do not chaqnge ****
set parametric
set size square
set size 1,1
set size ratio -1
set clip
set xtics axis mirror scale 0.5 offset 0.5,0.5
set ytics axis mirror scale 0.5 offset 0.5,0.5
unset grid
set angles degrees
set polar
set grid polar 10.


unset border
unset param
set style data linespoints 

set tmargin 3.5; set bmargin 1.5


set key outside left box

#******END of Do not chaqnge  settings****

set title "C_P along circumference" offset 1.5,1.5
set key at -1.75,2.25 box

#####**** Change settings here afterwards*******####
#vinodkumar_g@vssc.gov.in,, vinodkgul@gmail.com
#The data is within -25 and 10, so  decide to have origin at  -30
#So total distance along radial line is,  R = rmax - rmin = 10 - (-30) = 40
# so set xrange and yrange to -40 to 40
# del_r = xtics=ytics = (R/No of Div)
#*****
# No of sectors in radial dira, Nsector
Nsector	= 5 # 10 tics along the radial line


origin	= -30.0 # Decided based on the min value of data set
Max	= 20.0 # Decided based on the miax value of data set
xy_max 	= (Max - origin) # Max value - Min value(origin)
xy_tics	= xy_max/Nsector
print "Origin is at = ",origin
print "RMax value is = ",Max
print "Total Radial length is = ",xy_max
print "No of sectors is = ",Nsector
print "Radial tics = ",xy_tics

set yrange [-xy_max:xy_max]  
set xrange [-xy_max:xy_max] 
set ytics abs(xy_tics)
set xtics abs(xy_tics)
set rrange [origin:Max]    # rmin = origin = -30, rmax = 0 = rmin + max_xrange*del_ro
#set rrange [Max:origin]    # rmin = origin = -30, rmax = 0 = rmin + max_xrange*del_ro

set title "C_P along circumference" offset 1.5,1.5
set key outside left box
set key at -xy_max,xy_max box # top left corner


#****
set label 100 "Path `pwd`" at graph 1.0, -0.18 right font "Arial,8" back tc lt 1
set angle radians
	factor = 1.0 #(xy_max/1.1)
set label 50 " Y-" at factor*xy_max*cos(0*pi),factor*xy_max*sin(0*pi) font "ArialBold,10" offset 3.5, -1.2
set label 51 " P-" at factor*xy_max*cos(0.5*pi)-0.025,factor*xy_max*sin(0.5*pi) font "ArialBold,10" offset 2.0,1.2
set label 52 " Y+" at factor*xy_max*cos(1*pi)-0.0,factor*xy_max*sin(1*pi) font "ArialBold,10" offset -3.5,-1.2 
set label 53 " P+" at factor*xy_max*cos(1.5*pi)-0.025,factor*xy_max*sin(1.5*pi) font "ArialBold,10" offset 2.5,-1.0


   set label 1 "0^o" at factor*xy_max*cos(0*pi)-0.025,factor*xy_max*sin(0*pi) font "Arial,10" offset 3.5,0
    set label 2 "45^o" at factor*xy_max*cos(0.25*pi)-0.025,factor*xy_max*sin(0.25*pi) rotate by -45 font "Arial,10" offset 1.2,1.2
    set label 3 "90^o" at factor*xy_max*cos(0.5*pi)-0.025,factor*xy_max*sin(0.5*pi) font "Arial,10" offset 0,1.3 
    set label 4 "135^o" at factor*xy_max*cos(0.75*pi)-0.025,factor*xy_max*sin(0.75*pi) rotate by 45 font "Arial,10" offset 2.2,2.2 
#    set label 5 "0^o" at factor*1.05*cos(1*pi)-0.025,factor*1.05*sin(1*pi) font "Arial,10"
    set label 5 "180^o" at factor*xy_max*cos(1*pi)-0.0,factor*xy_max*sin(1*pi) font "Arial,10" offset -3.5,0
    set label 6 "225^o" at factor*xy_max*cos(1.25*pi)-0.025,factor*xy_max*sin(1.25*pi) rotate by 135 font "Arial,10" offset -1.2,-1.2
    set label 7 "270^o" at factor*xy_max*cos(1.5*pi)-0.025,factor*xy_max*sin(1.5*pi) font "Arial,10" offset 0,-1.0
    

set angle degrees

maxxloc = 1
ind = 0

call 'plot_polar.gnu' "pl.dat"