Index: /trunk/DataCheck/Tools/fill_r750.sh
===================================================================
--- /trunk/DataCheck/Tools/fill_r750.sh	(revision 18938)
+++ /trunk/DataCheck/Tools/fill_r750.sh	(revision 18939)
@@ -12,75 +12,75 @@
          0.0 # 5
          0.0 # 6
-       267.1 # 7
-       267.1 # 8
-       271.8 # 9
-       295.0 #10
-       257.3 #11
-       248.6 #12
-       237.1 #13
-       244.0 #14
-       241.0 #15
-       263.8 #16
-       245.4 #17
-       254.3 #18
-       263.3 #19
-       278.1 #20
-       267.1 #21
-       234.9 #22
-       259.9 #23
-       268.8 #24
-       281.6 #25
-       272.1 #26
-       269.5 #27
-       266.6 #28
+       260.7 # 7
+       261.2 # 8
+       267.9 # 9
+       287.0 #10
+       251.3 #11
+       251.3 #12
+       237.7 #13
+       246.6 #14
+       243.4 #15
+       267.9 #16
+       249.4 #17
+       254.8 #18
+       261.5 #19
+       274.2 #20
+       264.5 #21
+       228.9 #22
+       255.1 #23
+       269.0 #24
+       281.9 #25
+       271.4 #26
+       270.5 #27
+       267.9 #28
          0.0 #29
-       246.5 #30
-       246.5 #31
-       286.9 #32
-       289.4 #33
-       278.7 #34
-       266.3 #35
-       258.4 #36
-       273.8 #37
-       276.3 #38
-       275.4 #39
-       279.9 #40
-       286.2 #41
-       289.1 #42
-       276.4 #43
-       279.3 #44
-       271.8 #45
-       252.9 #46b
-       254.3 #47
-       252.0 #48
-       255.4 #49
-       265.3 #50
-       263.8 #51
-       261.3 #52
-       253.9 #53
-       242.5 #54
-       242.6 #55
-       241.5 #56
-       248.7 #57
-       248.7 #58
-       243.4 #59
-       249.3 #60b
-       248.4 #61
-       261.3 #62
-       257.6 #63
-       244.3 #64
+       246.9 #30
+       246.9 #31
+       282.8 #32
+       282.7 #33
+       272.0 #34
+       261.1 #35
+       256.0 #36
+       273.6 #37
+       273.2 #38
+       280.8 #39
+       282.9 #40
+       286.7 #41
+       290.7 #42
+       275.5 #43
+       276.1 #44
+       267.2 #45
+       246.9 #46b
+       248.1 #47
+       246.2 #48
+       253.3 #49
+       265.1 #50
+       266.2 #51
+       264.7 #52
+       257.6 #53
+       244.9 #54
+       244.7 #55
+       239.6 #56
+       244.6 #57
+       242.6 #58
+       237.1 #59
+       247.3 #60b
+       245.5 #61
+       260.3 #62
+       258.7 #63
+       248.7 #64
        248.2 #65
-       241.6 #66
-       234.8 #67
-       234.1 #68
-       291.2 #69
-       248.8 #70
-       233.1 #71
-       242.0 #72
-       226.2 #73
-       215.1 #74
-       244.6 #75
-       246.7 #76
-       249.9 #77
+       245.9 #66
+       235.9 #67
+       233.2 #68
+       287.2 #69
+       243.1 #70
+       231.4 #71
+       237.8 #72
+       222.3 #73
+       214.0 #74
+       244.0 #75
+       249.6 #76
+       253.4 #77
          0.0 #78
        )
@@ -101,15 +101,18 @@
 # periods 42 and 60 are split: 
 # 60a
-query="UPDATE RunInfo SET fR750Ref=243.4 WHERE fPeriod=60 AND fNight<20170901"
+query="UPDATE RunInfo SET fR750Ref=237.1 WHERE fPeriod=60 AND fNight<20170901"
 mysql --defaults-file=$sqlpw -u operator -h 10.0.100.21 $dbname -s -e "$query" 
 # 46a
-query="UPDATE RunInfo SET fR750Ref=275 WHERE fPeriod=46 AND fNight<20150716"
+query="UPDATE RunInfo SET fR750Ref=271 WHERE fPeriod=46 AND fNight<20150716"
 mysql --defaults-file=$sqlpw -u operator -h 10.0.100.21 $dbname -s -e "$query" 
 
+zd="cos(fZenithDistanceMean/180*Pi())"
+min="TIME_TO_SEC(TIMEDIFF(fRunStop,fRunStart))*fEffectiveOn/60"
+cor="pow("$zd", 1.091)*exp((1-1/"$zd")*0.593)*exp(1.829*(1-"$zd"))"
 
 # fill corrected R750
 query="UPDATE RunInfo SET "
-query=$query" fR750Cor=round(fNumThreshold750/(TIME_TO_SEC(TIMEDIFF(fRunStop,fRunStart))*fEffectiveOn/60)/pow(1/cos(fZenithDistanceMean/180*Pi())*exp(1-1/cos(fZenithDistanceMean/180*Pi())), 1.22),2) "
-query=$query" WHERE fRunTypeKey=1 AND NOT ((TIME_TO_SEC(TIMEDIFF(fRunStop,fRunStart))*fEffectiveOn/60)/pow(1/cos(fZenithDistanceMean/180*Pi())*exp(1-1/cos(fZenithDistanceMean/180*Pi())), 1.22))=0"
+query=$query" fR750Cor=round(fNumThreshold750/("$min")/("$cor"),2) "
+query=$query" WHERE fRunTypeKey=1 AND NOT (("$min")/("$cor"))=0"
 echo $query
 mysql --defaults-file=$sqlpw -u operator -h 10.0.100.21 $dbname -s -e "$query" 
