makegrid_tilteddipole_varx2_oneside_3D for Southern Hemisphere

[paulinchin]

While using
makegrid_tilteddipole_varx2_oneside_3D

have an error for preparing grid for southern hemisphere:
% Iowa
p.dtheta=34.5; % 30.5
p.dphi=70; % 70
p.lp=432; % USE FOR VEGA (36)
p.lq=360;
p.lphi=360; % USE FOR VEGA (10) 360
p.altmin=80e3;
p.glat=-40;
p.glon=360-93;
p.gridflag=0; % 1 - works
p.flagsource=0;
p.iscurv=true;

Index exceeds the number of array elements (1).

Error in gemini3d.grid.makegrid_tilteddipole_varx2_oneside_3D (line 177)
pstride=p(2)-p(1);

Error in iowaeqmy (line 19)
xg=gemini3d.grid.makegrid_tilteddipole_varx2_oneside_3D(p.dtheta,p.dphi,p.lp,p.lq,p.lphi,p.altmin,p.glat,p.glon,p.gridflag);

[scivision]

I get the same error. I updated numerous mat_gemini-scripts to have the correct references to mat_gemini package and I get this same error message

[scivision]

p(1) = pmin  % 4.005

pmax is 1.06 so the while p(ip)<pmax never runs. I just added an assert() check for this so it errors early.

[scivision]

for the many similar cases such as tilted_dipole.m in the same directory, the linspace(pmin,pmax) works to make a decreasing vector. For this specific makegrid function, a distinct while loop is used that expected an increasing p vector.
@mattzett

[mattzett]

It's very likely this is a use that has never been tested before, @paulinchin did we fix this on the Dropbox version or do I still need to look over it?

[mattzett]

There are numerous things that need to be fixed here; I'm working on it :/