Changeset 35204 for trunk/ippToPsps/jython/queue.py
- Timestamp:
- Feb 26, 2013, 8:50:05 AM (13 years ago)
- File:
-
- 1 edited
-
trunk/ippToPsps/jython/queue.py (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ippToPsps/jython/queue.py
r35203 r35204 8 8 import sys 9 9 import os 10 import math 10 11 import logging.config 11 12 … … 117 118 # starting positions 118 119 119 dD = 0.5*s kychunk.boxSize120 dD = 0.5*self.skychunk.boxSize 120 121 dec = self.skychunk.minDec + dD 122 if (dec > 85): raise 123 124 # print "dec, dD: ", dec, dD 121 125 122 126 # a given box centered at RA,DEC with (linear) size of SIDE has the following bounds 123 dR = 0.5*s kychunk.boxSize / math.cos(math.radians(dec))127 dR = 0.5*self.skychunk.boxSize / math.cos(math.radians(dec)) 124 128 ra = self.skychunk.minRa + dR 125 129 130 # print "ra, dR: ", ra, dR 131 132 # XXX not sure why this is not done with a for-loop... 126 133 while dec <= self.skychunk.maxDec: 127 128 # XXX EAM : need to do something special at the poles...129 if (dec > 85): raise130 131 134 while ra <= self.skychunk.maxRa: 132 135 … … 137 140 maxDEC = dec + dD 138 141 139 XXX EAM:142 # this function inserts the box with given center and (linear) skychunk.boxSize 140 143 box_id = self.ippToPspsDb.insertBox(ra, dec) 141 144 142 XXX EAM: 143 ids = self.ippToPspsDb.getItemsInThisThisBox(ra, dec) 145 ids = self.ippToPspsDb.getItemsInThisBox(minRA, maxRA, minDEC, maxDEC) 144 146 145 147 self.logger.info("| %5.1f | %5.1f | %9d | %9d |" % ( … … 149 151 len(ids))) 150 152 153 # print "box_id: ", box_id 154 151 155 if len(ids) > 0: self.ippToPspsDb.insertPending(box_id, batchType, ids) 156 # print "inserted pending" 152 157 ra = ra + 2*dR 153 154 dec = dec + 2*dD 155 156 dR = 0.5*skychunk.boxSize / math.cos(math.radians(dec)) 157 ra = self.skychunk.minRa + dR 158 # print "new ra: ", ra 159 160 dec = dec + 2*dD 161 if (dec > 85): raise 162 163 dR = 0.5*self.skychunk.boxSize / math.cos(math.radians(dec)) 164 ra = self.skychunk.minRa + dR 158 165 159 166 self.logger.info("+-------------+-------------+-------------+-------------+")
Note:
See TracChangeset
for help on using the changeset viewer.
