


XYSCALE_SOUTH - calculate scaling factor
latitude is <0 in degrees.
Usage:
res=xyscale_south(latatitude)

0001 function res=xyscale_south(lat); 0002 %XYSCALE_SOUTH - calculate scaling factor 0003 % 0004 % latitude is <0 in degrees. 0005 % 0006 % Usage: 0007 % res=xyscale_south(latatitude) 0008 lat=-lat/180*pi; 0009 0010 m71_t71=1.9390295644; 0011 e2=.00669437999015; 0012 e=sqrt(e2); 0013 0014 m=cos(lat)./sqrt(1-e2*sin(lat).^2); 0015 t=tan(pi/4-lat/2)./((1-e*sin(lat))./(1+e*sin(lat))).^(e/2); 0016 k=m71_t71*t./m; 0017 res=(1./k./k);