Skip to content

Commit 224a5d4

Browse files
committed
BUG: um, dk needs to be POSITIVE for dipole. nan_checks pass now
1 parent 2274a46 commit 224a5d4

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

src/grid_spacing.cpp

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -148,11 +148,11 @@ void Grid::calc_k_grid_spacing() {
148148
}
149149
// This needs to be turned into a distance for the dipole:
150150
if (iGridShape_ == iDipole_){
151-
// the dk's may not be negative (like the sphere). make sure they are
151+
// the dk's may be negative (not allowed). make sure they are positive
152152
// this gets rid of SO many errors...
153-
dk_center_scgc = -1.0 * abs(dk_center_scgc);
154-
dk_edge = -1.0 * abs(dk_edge);
155-
dr_edge = -1.0 * abs(dr_edge);
153+
dk_center_scgc = abs(dk_center_scgc);
154+
dk_edge = abs(dk_edge);
155+
dr_edge = abs(dr_edge);
156156
dk_center_m_scgc = pow(magAlt_scgc, 3) % dk_center_scgc / delTheta(magLat_scgc);
157157
dk_edge_m = pow(magAlt_scgc, 3) % dk_edge / delTheta(magLat_scgc);
158158
}

0 commit comments

Comments
 (0)