Alexandre Schimel

FYI. I just submitted a new version of the wgs2utm allowing to force UTM zone and hemisphere in input, as per requests.

You can still use it like the old one though, ie with no forcing and let the function choose the appropriate zone. The new function should be compatible with any code using the old version. Check function's help for more info on use.

Comment about forcing output UTM zone.

Two parameters change for each UTM zone: the reference longitude and the false northing. In the current code, these values are caculated from the input data at lines 63 and 68:
Lon0 = floor(Lon./6).*6+3; % reference longitude in degrees
FN = (Lat < 0).*10000000; % false northing

In order to force the output UTM zone, one can force these two parameters to their desired value. If you know your reference longitude, just replace Lon0 by this value. If you know the zone number ii, replace Lon0 by 6*ii-183

Set FN to 0 for northern hemisphere, or 10000000 for southern hemisphere.

Hi and thanks for your code,
I just did a quick comparison between your code and Matlab's built-in functions to calculate UTM values based on LatLon, here are the results:
Matlab:
mstruct = defaultm('utm');
mstruct.zone=utmzone(-43.5,172.5);
mstruct = defaultm(mstruct);
[x,y] = mfwdtran(mstruct,-43.5,172.5)
Result: 621272.828768645 5182478.29916538
Your code:
[x,y,zone1,utmhemi1] = wgs2utm(-43.5, 172.5)
Result: 621267.241003134 5182565.80639284
(there is almost 100 meter diff in y).
By the way, you have used utmzone as variable name in the output which conflicts with Matlab's built-in function name, it took me a while why I am getting weird errors after trying your function (output from your function stops utmzone function from working). Good luck

FYI. I just submitted a new version of the wgs2utm allowing to force UTM zone and hemisphere in input, as per requests.

You can still use it like the old one though, ie with no forcing and let the function choose the appropriate zone. The new function should be compatible with any code using the old version. Check function's help for more info on use.