- Timestamp:
- 11/19/14 13:23:30 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/STGrid.cpp
r3011 r3015 1136 1136 Double wx = max( abs(xmax-center_(0)), abs(xmin-center_(0)) ) * 2 ; 1137 1137 Double wy = max( abs(ymax-center_(1)), abs(ymin-center_(1)) ) * 2 ; 1138 // take 10% margin1139 wx *= 1.10 ;1140 wy *= 1.10 ;1141 1138 1142 1139 Quantum<Double> qcellx ; … … 1172 1169 wy = 0.00290888 ; 1173 1170 } 1174 qcellx = Quantum<Double>( wx/nx_*cos(center_[1]), "rad" ) ; 1175 qcelly = Quantum<Double>( wy/ny_, "rad" ) ; 1171 if (nx_ > 1) { 1172 qcellx = Quantum<Double>(wx / (nx_ - 1) * cos(center_[1]), "rad"); 1173 } 1174 else { 1175 qcellx = Quantum<Double>( 1.1f * wx / nx_ * cos(center_[1]), "rad" ); 1176 } 1177 if (ny_ > 1) { 1178 qcelly = Quantum<Double>(wy / (ny_ - 1), "rad"); 1179 } 1180 else { 1181 qcelly = Quantum<Double>( 1.1f * wy / ny_, "rad" ); 1182 } 1176 1183 } 1177 1184 }
Note:
See TracChangeset
for help on using the changeset viewer.