cd /home/disk/gamera/kwajex/dualpol/ascii/
	
test_readcfadlist("/home/disk/gamera/kwajex/dualpol/ascii/200106.output",12,3576)
test_readcfadlist("/home/disk/gamera/kwajex/dualpol/ascii/199910.output",12,3576)

#3 month data sets
test_readcfadlist("/home/disk/gamera/kwajex_ls/0207talk/1999_07_08_09.output",12,6483)
test_readcfadlist("/home/disk/gamera/kwajex_ls/0207talk/1999_10_11_12.output",12,12164)
test_readcfadlist("/home/disk/gamera/kwajex_ls/0207talk/2000_07_08_09.output",12,11423)
test_readcfadlist("/home/disk/gamera/kwajex_ls/0207talk/2000_10_11_12.output",12,9050)
test_readcfadlist("/home/disk/gamera/kwajex_ls/0207talk/2001_07_08_09.output",12,9050)
test_readcfadlist("/home/disk/gamera/kwajex_ls/0207talk/2001_10_11_12.output",12,10614)
	
test.area_computecfadlistareas(test,5)


	
showmints(test.area$kmarea/(121*121)*100,test.area$julmin,atitle="KWAJ 2001 Precip Echo Area",minint=60*1,ylab="% Echo Area > 20 dBZ within 240 km range")

kwajlowareaf_(pi*(242*242)-pi*(17*17))/100 # area of 242 km diameter circle - 17 km circle

showmints.whl(test.area$kmarea/kwajlowareaf,test.area$julmin,c(0),atitle="KWAJ 2001 Precip Echo Area",minint=60*24,ylab="% Echo Area > 20 dBZ within 240 km range",ylim=c(0,40),pretty=T)


showmints.whl(sqrt(test.area$kmarea),test.area$julmin,c(0),atitle="KWAJ 2001 Precip Echo Area",minint=60*24,ylab="Sqrt Echo Area > 20 dBZ within 240 km range",ylim=c(0,300),pretty=T)

trmm0207fig3month_function(areaval,atitle,adjdate=0){
fact_(24*60)*adjdate
showmints.whl(areaval$kmarea/kwajlowareaf,areaval$julmin-fact,c(0),atitle=atitle,minint=60*24,ylab="% Echo Area > 20 dBZ within 240 km range",ylim=c(0,40),pretty=T)
	}


trmm0207fig3month(test.area,"Jul, Aug, Sep 1999",0)
printgraph(file="/home/disk/gamera/kwajex_ls/0207talk/1999_07_08_09echoarea.ps")
trmm0207fig3month(test.area,"Oct, Nov, Dec 1999",0)
printgraph(file="/home/disk/gamera/kwajex_ls/0207talk/1999_10_11_12echoarea.ps")

 
trmm0207fig3month(test.area,"Jul, Aug, Sep 2000",1)
printgraph(file="/home/disk/gamera/kwajex_ls/0207talk/2000_07_08_09echoarea.ps")
trmm0207fig3month(test.area,"Oct, Nov, Dec 2000",1)
printgraph(file="/home/disk/gamera/kwajex_ls/0207talk/2000_10_11_12echoarea.ps")


trmm0207fig3month(test.area,"Jul, Aug, Sep 2001",1)
printgraph(file="/home/disk/gamera/kwajex_ls/0207talk/2001_07_08_09echoarea.ps")
trmm0207fig3month(test.area,"Oct, Nov, Dec 2001",0)
printgraph(file="/home/disk/gamera/kwajex_ls/0207talk/2001_10_11_12echoarea.ps")


	
calcrainareatosamplearea_function(julmin,kmarea,domainarea,begtime,endtime){
	subkmarea_subset.tepps.ts(kmarea,julmin,begtime,endtime)$data
	subdomainarea_subset.tepps.ts(domainarea,julmin,begtime,endtime)$data
	print(paste("number of samples:",length(subkmarea)))
	return(sum(subkmarea)/sum(subdomainarea)*100)}
	
calcrainareatosamplearea(kwaj.area.01.070809$julmin,kwaj.area.01.070809$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.01.070809$julmin)),
	timetojulmin(2001,7,1),timetojulmin(2001,7,31))

[1] "number of samples: 3547" [1] 5.037559 # is the same as it should be
	
# old version "number of samples: 3547" 5.037559

calcrainareatosamplearea(kwaj.area.01.070809$julmin,kwaj.area.01.070809$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.01.070809$julmin)),
	timetojulmin(2001,8,1),timetojulmin(2001,8,31)) 

[1] "number of samples: 3464" [1] 5.695382 # is the same as it should be

#old version [1] "number of samples: 3464" 5.695382

calcrainareatosamplearea(kwaj.area.01.070809$julmin,kwaj.area.01.070809$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.01.070809$julmin)),
	timetojulmin(2001,9,1),timetojulmin(2001,9,30))
[1] "number of samples: 2895" [1] 4.725553, is different since gaps filled
# old version [1] "number of samples: 1930"  4.139397

calcrainareatosamplearea(kwaj.area.01.101112$julmin,kwaj.area.01.101112$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.01.101112$julmin)),
	timetojulmin(2001,10,1),timetojulmin(2001,10,31))
[1] "number of samples: 3499" [1] 8.714166 # is the same as it should be
#old version [1] "number of samples: 3499"  8.714166

calcrainareatosamplearea(kwaj.area.01.101112$julmin,kwaj.area.01.101112$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.01.101112$julmin)),
	timetojulmin(2001,11,1),timetojulmin(2001,11,30))
[1] "number of samples: 3426" [1] 7.733 # is the same as it should be
# old version [1] "number of samples: 3426"  7.733

calcrainareatosamplearea(kwaj.area.01.101112$julmin,kwaj.area.01.101112$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.01.101112$julmin)),
	timetojulmin(2001,12,1),timetojulmin(2001,12,31)) 
# [1] "number of samples: 3330" [1] 3.006504 is the same as it should be
# old version [1] "number of samples: 3330"3.006504

calcrainareatosamplearea(kwaj.area.00070809$julmin,kwaj.area.00070809$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.00070809$julmin)),
	timetojulmin(2000,7,1),timetojulmin(2000,7,31)) 
[1] "number of samples: 3846" [1] 6.65887
# old version [1] "number of samples:\t 3846" [1] 6.65887

calcrainareatosamplearea(kwaj.area.00070809$julmin,kwaj.area.00070809$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.00070809$julmin)),
	timetojulmin(2000,8,1),timetojulmin(2000,8,31))

#[1] "number of samples: 3703" [1] 9.070178

# old version [1] "number of samples: 3703" [1] 9.702242

calcrainareatosamplearea(kwaj.area.00070809$julmin,kwaj.area.00070809$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.00070809$julmin)),
	timetojulmin(2000,9,1),timetojulmin(2000,9,30))
[1] "number of samples: 3476" [1] 5.583174
# old version [1] "number of samples: 3476" [1] 7.487812


calcrainareatosamplearea(kwaj.area.00101112$julmin,kwaj.area.00101112$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.00101112$julmin)),
	timetojulmin(2000,10,1),timetojulmin(2000,10,31))

[1] "number of samples: 2322" [1] 6.843451
# old version [1] "number of samples: 2635" [1] 8.819308

calcrainareatosamplearea(kwaj.area.00101112$julmin,kwaj.area.00101112$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.00101112$julmin)),
	timetojulmin(2000,11,1),timetojulmin(2000,11,30))  

[1] "number of samples: 2714" [1] 5.654009

# oldversion [1] "number of samples: 2714" [1] 7.472861

calcrainareatosamplearea(kwaj.area.00101112$julmin,kwaj.area.00101112$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.00101112$julmin)),
	timetojulmin(2000,12,1),timetojulmin(2000,12,31))
[1] "number of samples: 3392" [1] 6.15623 
# old version [1] "number of samples: 3392" [1] 5.896987

calcrainareatosamplearea(kwaj.area.99070809$julmin,kwaj.area.99070809$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.99070809$julmin)),
	timetojulmin(1999,7,1),timetojulmin(1999,7,31)) 
[1] "number of samples: 519" [1] 3.534944 
# old version [1] "number of samples: 856" [1] 7.467332


calcrainareatosamplearea(kwaj.area.99070809$julmin,kwaj.area.99070809$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.99070809$julmin)),
	timetojulmin(1999,8,1),timetojulmin(1999,8,31))
	
[1] "number of samples: 3581" [1] 4.327397 # same as before as it should be
# old version [1] "number of samples: 3581" [1] 4.327397

calcrainareatosamplearea(kwaj.area.99070809$julmin,kwaj.area.99070809$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.99070809$julmin)),
	timetojulmin(1999,9,1),timetojulmin(1999,9,30))
	
[1] "number of samples: 1763" [1] 5.662974
# old version [1] "number of samples: 1763" [1] 3.113239

calcrainareatosamplearea(kwaj.area.99.101112$julmin,kwaj.area.99.101112$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.99.101112$julmin)),
	timetojulmin(1999,10,1),timetojulmin(1999,10,31))
	
[1] "number of samples: 3801" [1] 5.668643
# old version [1] "number of samples: 3845" [1] 3.883527

calcrainareatosamplearea(kwaj.area.99.101112$julmin,kwaj.area.99.101112$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.99.101112$julmin)),
	timetojulmin(1999,11,1),timetojulmin(1999,11,30)) 

[1] "number of samples: 3946" [1] 6.473102
# old version [1] "number of samples: 3946" [1] 4.497403


calcrainareatosamplearea(kwaj.area.99.101112$julmin,kwaj.area.99.101112$kmarea,
	rep(kwajlowareaf*100,length(kwaj.area.99.101112$julmin)),
	timetojulmin(1999,12,1),timetojulmin(1999,12,31))
[1] "number of samples: 3957" [1] 2.73825
# old version [1] "number of samples: 3957" [1] 1.740464

#redone 0716# kwajpercentrainarea99_c(0, 4.327397, 5.662974, 5.668643, 6.473102, 2.73825)
#redone 0716# kwajpercentrainarea00_c(6.65887,9.070178, 5.583174,  6.843451, 5.654009, 6.15623 )
#redone 0716# kwajpercentrainarea01_c(5.037559, 5.695382, 4.725553, 8.714166, 7.733, 3.006504)
	

par(mfrow=c(1,3))
barplot(kwajpercentrainarea99,names=c("Jul","Aug","Sep","Oct","Nov","Dec"),
		ylim=c(0,10),cex=1.2)
barplot(kwajpercentrainarea00,names=c("Jul","Aug","Sep","Oct","Nov","Dec"),
		ylim=c(0,10),cex=1.2)
barplot(kwajpercentrainarea01,names=c("Jul","Aug","Sep","Oct","Nov","Dec"),
		ylim=c(0,10),cex=1.2)
par(mfrow=c(1,1))
# old version, see useful.s in gamera/kwajex_ls/0207talk printgraph(file="/home/disk/gamera/kwajex_ls/0207talk/monthlyrainareabar.ps")