| # | Line 128 | Line 128 | init_pabopto_inp(const int i, const char *fname) | |
|---|---|---|
| 128 | dv[0] = cos(M_PI/180.*inpfile[i].phi)*dv[2]; | |
| 129 | dv[1] = sin(M_PI/180.*inpfile[i].phi)*dv[2]; | |
| 130 | dv[2] = sqrt(1. - dv[2]*dv[2]); | |
| 131 | < | pos_from_vec(inpfile[i].igp, dv); |
| 131 | > | if (inpfile[i].theta <= FTINY) |
| 132 | > | inpfile[i].igp[0] = inpfile[i].igp[1] = grid_res/2 - 1; |
| 133 | > | else |
| 134 | > | pos_from_vec(inpfile[i].igp, dv); |
| 135 | return(1); | |
| 136 | } | |
| 137 | ||
| – | Removed lines |
| + | Added lines |
| < | Changed lines (old) |
| > | Changed lines (new) |