CONVERTER PARA VB

USUARIO.EXCLUIDOS 21/11/2006 10:23:40
#185512

Bom dia a todos, estou criando um programinha em VB e tenho que incorporar a ele, um outro progarminha que está em C. Agluém poderia me ajudar, convertendo o programa de C para VB? Pois eu não entendo absolutamente nada de C. O código fonte do programinha encontra-se abaixo e todos os seus componentes em anexo.

----
/* PROGRAM MAGPOINT (GEOMAG DRIVER) */
/************************************************************************
Contact Information

Software and Model Support
National Geophysical Data Center
NOAA EGC/2
325 Broadway
Boulder, CO 80303 USA
Attn: Susan McLean or Stefan Maus
Phone: (303) 497-6478 or -6522
Email: Susan.McLean@noaa.gov or Stefan.Maus@noaa.gov
Web: http://www.ngdc.noaa.gov/seg/WMM/

Sponsoring Government Agency
National Geospatial-Intelligence Agency
PRG / CSAT, M.S. L-41
3838 Vogel Road
Arnold, MO 63010
Attn: Craig Rollins
Phone: (314) 263-4186
Email: Craig.M.Rollins@Nga.Mil

Original Program By:
Dr. John Quinn
FLEET PRODUCTS DIVISION, CODE N342
NAVAL OCEANOGRAPHIC OFFICE (NAVOCEANO)
STENNIS SPACE CENTER (SSC), MS 39522-5001

3/25/05 Stefan Maus corrected 2 bugs:
- use %c instead of %s for character read
- help text: positive inclination is downward
*/

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

/* The following include file must define a function 'isnan' */
/* This function, which returns '1' if the number is NaN and 0*/
/* otherwise, could be hand-written if not available. */
/* Comment out one of the two following lines, as applicable */
#include <math.h> /* for gcc */
#include <mathimf.h> /* for Intel icc */

#define NaN log(-1.0)

static void E0000(int IENTRY, int *maxdeg, float alt,float glat,float glon, float time, float *dec, float *dip, float *ti, float *gv);
void geomag(int *maxdeg);
void geomg1(float alt, float glat, float glon, float time, float *dec, float *dip, float *ti, float *gv);
char geomag_introduction(float epochlowlim);

int main()
{
int warn_H, warn_H_strong, warn_P;
static int maxdeg;
static float altm, dlat, dlon;
static float ati, adec, adip;
static float alt, time, dec, dip, ti, gv;
static float time1, dec1, dip1, ti1;
static float time2, dec2, dip2, ti2;
char answer, ans;
char decd[5], dipd[5],d_str[81],modl[20];
char goodbye[81];
float x1,x2,y1,y2,z1,z2,h1,h2;
float ax,ay,az,ah;
float rTd=0.017453292;
float epochlowlim,epochuplim;
float epochrange = 5.0;
float warn_H_val, warn_H_strong_val;
float dmin, imin, ddeg, ideg;
FILE *wmmtemp;

wmmtemp = fopen("WMM.COF","r");
fgets(d_str, 80, wmmtemp);
sscanf(d_str,"%f%s",&epochlowlim,modl);
fclose(wmmtemp);
strcpy (goodbye,"
-- End of WMM Point Calculation Program --

");
ans = geomag_introduction(epochlowlim);

if ((ans == 'y') || (ans == 'Y'))
{

/* INITIALIZE GEOMAG ROUTINE */

S1:

maxdeg = 12;
warn_H = 0;
warn_H_val = 99999.0;
warn_H_strong = 0;
warn_H_strong_val = 99999.0;
warn_P = 0;

geomag(&maxdeg);


printf("


ENTER LATITUDE IN DECIMAL DEGREES ");
printf("
(North latitude positive, South latitude negative
");
printf("i.e. 25.5 for 25 degrees 30 minutes north.)
");
scanf("%f%*[^
]", &dlat);
getchar();

printf("ENTER LONGITUDE IN DECIMAL DEGREES");
printf("(East longitude positive, West negative
");
printf("i.e.- 100.0 for 100.0 degrees west.)
");
scanf("%f%*[^
]", &dlon);
getchar();

printf("ENTER ALTITUDE IN METERS ABOVE MEAN SEA LEVEL (WGS84)
");
scanf("%f%*[^
]", &altm);
getchar();
alt = altm/1000;

epochuplim = epochlowlim + epochrange;
printf("ENTER TIME IN DECIMAL YEAR (%-7.2f - %-7.2f)
",epochlowlim,epochuplim);
scanf("%f%*[^
]",&time);
getchar();

geomg1(alt,dlat,dlon,time,&dec,&dip,&ti,&gv);
time1 = time;
dec1 = dec;
dip1 = dip;
ti1 = ti;
time = time1 + 1.0;

geomg1(alt,dlat,dlon,time,&dec,&dip,&ti,&gv);
time2 = time;
dec2 = dec;
dip2 = dip;
ti2 = ti;

/*COMPUTE X, Y, Z, AND H COMPONENTS OF THE MAGNETIC FIELD*/

x1=ti1*(cos((dec1*rTd))*cos((dip1*rTd)));
x2=ti2*(cos((dec2*rTd))*cos((dip2*rTd)));
y1=ti1*(cos((dip1*rTd))*sin((dec1*rTd)));
y2=ti2*(cos((dip2*rTd))*sin((dec2*rTd)));
z1=ti1*(sin((dip1*rTd)));
z2=ti2*(sin((dip2*rTd)));
h1=ti1*(cos((dip1*rTd)));
h2=ti2*(cos((dip2*rTd)));

/* COMPUTE ANNUAL CHANGE FOR TOTAL INTENSITY */
ati = ti2 - ti1;

/* COMPUTE ANNUAL CHANGE FOR DIP & DEC */
adip = (dip2 - dip1) * 60.;
adec = (dec2 - dec1) * 60.;


/* COMPUTE ANNUAL CHANGE FOR X, Y, Z, AND H */
ax = x2-x1;
ay = y2-y1;
az = z2-z1;
ah = h2-h1;


if (dec1 < 0.0) {
strcpy (decd,"(WEST)");
}
else
{
strcpy(decd,"(EAST)");
}

if (dip1 < 0.0)
{
strcpy(dipd,"(UP) ");
}
else
{
strcpy(dipd,"(DOWN)");
}

/* deal with geographic and magnetic poles */

if (h1 < 100.0) /* at magnetic poles */
{
dec1 = NaN;
adec = NaN;
strcpy(decd,"(VOID)");
/* while rest is ok */
}

if (h1 < 1000.0)
{
warn_H = 0;
warn_H_strong = 1;
warn_H_strong_val = h1;
}
else if (h1 < 5000.0 && !warn_H_strong)
{
warn_H = 1;
warn_H_val = h1;
}

if (90.0-fabs(dlat) <= 0.001) /* at geographic poles */
{
x1 = NaN;
y1 = NaN;
dec1 = NaN;
ax = NaN;
ay = NaN;
adec = NaN;
strcpy(decd,"(VOID)");
warn_P = 1;
warn_H = 0;
warn_H_strong = 0;
/* while rest is ok */
}

/* convert D and I to deg and min */
if (isnan(dec1)) ddeg = dec1; else ddeg=(int)dec1;
dmin=(dec1-(float)ddeg)*60;
if(ddeg!=0) dmin=fabs(dmin);
if (isnan(dip1)) ideg = dip1; else ideg=(int)dip1;
imin=(dip1-(float)ideg)*60;
if(ideg!=0) imin=fabs(imin);

printf("
Results For
");
if (dlat < 0)
printf("
LATITUDE: %7.2fS",-dlat);
else
printf("
LATITUDE: %7.2fN",dlat);
if (dlon < 0)
printf("
LONGITUDE: %7.2fW",-dlon);
else
printf("
LONGITUDE: %7.2fE",dlon);

printf("
ALTITUDE: %8.2f METERS AMSL (WGS84)",altm);
printf("
DATE: %6.1f
",time1);


printf("
Main Field              Secular Change");

printf("
F = %-9.1f nT         dF = %-8.1f nT/yr",ti1,ati);
if (isnan(h1))
printf("
H = NaN          dH = NaN");
else
printf("
H = %-9.1f nT         dH = %-8.1f nT/yr",h1,ah);
if (isnan(x1))
printf("
X = NaN          dX = NaN");
else
printf("
X = %-9.1f nT         dX = %-8.1f nT/yr ",x1,ax);
if (isnan(y1))
printf("
Y = NaN          dY = NaN");
else
printf("
Y = %-9.1f nT         dY = %-8.1f nT/yr ",y1,ay);
printf("
Z = %-9.1f nT         dZ = %-8.1f nT/yr ",z1,az);
if (isnan(dec1))
printf("
D = NaN          dD = NaN");
else
printf("
D = %4.0f Deg %3.0f Min %s     dD = %-8.1f Min/yr",ddeg,dmin,decd,adec);
printf("
I = %4.0f Deg %3.0f Min %s     dI = %-8.1f Min/yr",ideg,imin,dipd,adip);

if (warn_H)
{
printf("

Warning: The horizontal field strength at this location is only %6.1f nT
",warn_H_val);
printf(" Compass readings have large uncertainties in areas where H is
");
printf(" smaller than 5000 nT
");
}
if (warn_H_strong)
{
printf("

Warning: The horizontal field strength at this location is only %6.1f nT
",warn_H_strong_val);
printf(" Compass readings have VERY LARGE uncertainties in areas where H is
");
printf(" smaller than 1000 nT
");
}
if (warn_P)
{
printf("

Warning: Location is at geographic pole where X, Y, and Decl are undefined
");
}

printf("

DO YOU NEED MORE POINT DATA? (y or n) ");
scanf("%c%*[^
]", &answer);
getchar();

if ((answer =='y')||(answer == 'Y')) goto S1;
else
{
printf("%s",goodbye);
}

}

else
{
printf("%s",goodbye);
}

exit(0);

}
/*************************************************************************/

static void E0000(int IENTRY, int *maxdeg, float alt, float glat, float glon, float time, float *dec, float *dip, float *ti, float *gv)
{
static int maxord,i,icomp,n,m,j,D1,D2,D3,D4;
static float c[13][13],cd[13][13],tc[13][13],dp[13][13],snorm[169],
sp[13],cp[13],fn[13],fm[13],pp[13],k[13][13],pi,dtr,a,b,re,
a2,b2,c2,a4,b4,c4,epoch,gnm,hnm,dgnm,dhnm,flnmj,otime,oalt,
olat,olon,dt,rlon,rlat,srlon,srlat,crlon,crlat,srlat2,
crlat2,q,q1,q2,ct,st,r2,r,d,ca,sa,aor,ar,br,bt,bp,bpp,
par,temp1,temp2,parp,bx,by,bz,bh;
static char model[20], c_str[81], c_new[5];
static float *p = snorm;
char answer;

FILE *wmmdat;

switch(IENTRY){case 0: goto GEOMAG; case 1: goto GEOMG1;}

GEOMAG:
wmmdat = fopen("WMM.COF","r");

/* INITIALIZE CONSTANTS */
maxord = *maxdeg;
sp[0] = 0.0;
cp[0] = *p = pp[0] = 1.0;
dp[0][0] = 0.0;
a = 6378.137;
b = 6356.7523142;
re = 6371.2;
a2 = a*a;
b2 = b*b;
c2 = a2-b2;
a4 = a2*a2;
b4 = b2*b2;
c4 = a4 - b4;

/* READ WORLD MAGNETIC MODEL SPHERICAL HARMONIC COEFFICIENTS */
c[0][0] = 0.0;
cd[0][0] = 0.0;

fgets(c_str, 80, wmmdat);
sscanf(c_str,"%f%s",&epoch,model);
S3:
fgets(c_str, 80, wmmdat);
/* CHECK FOR LAST LINE IN FILE */
for (i=0; i<4 && (c_str[i] != '\0'); i++)
{
c_new[i] = c_str[i];
c_new[i+1] = '\0';
}
icomp = strcmp("9999", c_new);
if (icomp == 0) goto S4;
/* END OF FILE NOT ENCOUNTERED, GET VALUES */
sscanf(c_str,"%d%d%f%f%f%f",&n,&m,&gnm,&hnm,&dgnm,&dhnm);
if (m <= n)
{
c[m][n] = gnm;
cd[m][n] = dgnm;
if (m != 0)
{
c[n][m-1] = hnm;
cd[n][m-1] = dhnm;
}
}
goto S3;

/* CONVERT SCHMIDT NORMALIZED GAUSS COEFFICIENTS TO UNNORMALIZED */
S4:
*snorm = 1.0;
for (n=1; n<=maxord; n++)
{
*(snorm+n) = *(snorm+n-1)*(float)(2*n-1)/(float)n;
j = 2;
for (m=0,D1=1,D2=(n-m+D1)/D1; D2>0; D2--,m+=D1)
{
k[m][n] = (float)(((n-1)*(n-1))-(m*m))/(float)((2*n-1)*(2*n-3));
if (m > 0)
{
flnmj = (float)((n-m+1)*j)/(float)(n+m);
*(snorm+n+m*13) = *(snorm+n+(m-1)*13)*sqrt(flnmj);
j = 1;
c[n][m-1] = *(snorm+n+m*13)*c[n][m-1];
cd[n][m-1] = *(snorm+n+m*13)*cd[n][m-1];
}
c[m][n] = *(snorm+n+m*13)*c[m][n];
cd[m][n] = *(snorm+n+m*13)*cd[m][n];
}
fn[n] = (float)(n+1);
fm[n] = (float)n;
}
k[1][1] = 0.0;

otime = oalt = olat = olon = -1000.0;
fclose(wmmdat);
return;

/*************************************************************************/

GEOMG1:

dt = time - epoch;
if (otime < 0.0 && (dt < 0.0 || dt > 5.0))
{
printf("

WARNING - TIME EXTENDS BEYOND MODEL 5-YEAR LIFE SPAN");
printf("
CONTACT NGDC FOR PRODUCT UPDATES:");
printf("
National Geophysical Data Center");
printf("
NOAA EGC/2");
printf("
325 Broadway");
printf("
Boulder, CO 80303 USA");
printf("
Attn: Susan McLean or Stefan Maus");
printf("
Phone: (303) 497-6478 or -6522");
printf("
Email: Susan.McLean@noaa.gov");
printf("
or");
printf("
Stefan.Maus@noaa.gov");
printf("
Web: http://www.ngdc.noaa.gov/seg/WMM/");
printf("

EPOCH = %.3lf",epoch);
printf("
TIME = %.3lf",time);
printf("
Do you wish to continue? (y or n) ");
scanf("%c%*[^
]",&answer);
getchar();
if ((answer == 'n') || (answer == 'N'))
{
printf("
Do you wish to enter more point data? (y or n) ");
scanf("%c%*[^
]",&answer);
getchar();
if ((answer == 'y')||(answer == 'Y')) goto GEOMG1;
else exit (0);
}
}

pi = 3.14159265359;
dtr = pi/180.0;
rlon = glon*dtr;
rlat = glat*dtr;
srlon = sin(rlon);
srlat = sin(rlat);
crlon = cos(rlon);
crlat = cos(rlat);
srlat2 = srlat*srlat;
crlat2 = crlat*crlat;
sp[1] = srlon;
cp[1] = crlon;

/* CONVERT FROM GEODETIC COORDS. TO SPHERICAL COORDS. */
if (alt != oalt || glat != olat)
{
q = sqrt(a2-c2*srlat2);
q1 = alt*q;
q2 = ((q1+a2)/(q1+b2))*((q1+a2)/(q1+b2));
ct = srlat/sqrt(q2*crlat2+srlat2);
st = sqrt(1.0-(ct*ct));
r2 = (alt*alt)+2.0*q1+(a4-c4*srlat2)/(q*q);
r = sqrt(r2);
d = sqrt(a2*crlat2+b2*srlat2);
ca = (alt+d)/r;
sa = c2*crlat*srlat/(r*d);
}
if (glon != olon)
{
for (m=2; m<=maxord; m++)
{
sp[m] = sp[1]*cp[m-1]+cp[1]*sp[m-1];
cp[m] = cp[1]*cp[m-1]-sp[1]*sp[m-1];
}
}
aor = re/r;
ar = aor*aor;
br = bt = bp = bpp = 0.0;
for (n=1; n<=maxord; n++)
{
ar = ar*aor;
for (m=0,D3=1,D4=(n+m+D3)/D3; D4>0; D4--,m+=D3)
{
/*
COMPUTE UNNORMALIZED ASSOCIATED LEGENDRE POLYNOMIALS
AND DERIVATIVES VIA RECURSION RELATIONS
*/
if (alt != oalt || glat != olat)
{
if (n == m)
{
*(p+n+m*13) = st**(p+n-1+(m-1)*13);
dp[m][n] = st*dp[m-1][n-1]+ct**(p+n-1+(m-1)*13);
goto S50;
}
if (n == 1 && m == 0)
{
*(p+n+m*13) = ct**(p+n-1+m*13);
dp[m][n] = ct*dp[m][n-1]-st**(p+n-1+m*13);
goto S50;
}
if (n > 1 && n != m)
{
if (m > n-2) *(p+n-2+m*13) = 0.0;
if (m > n-2) dp[m][n-2] = 0.0;
*(p+n+m*13) = ct**(p+n-1+m*13)-k[m][n]**(p+n-2+m*13);
dp[m][n] = ct*dp[m][n-1] - st**(p+n-1+m*13)-k[m][n]*dp[m][n-2];
}
}
S50:
/*
TIME ADJUST THE GAUSS COEFFICIENTS
*/
if (time != otime)
{
tc[m][n] = c[m][n]+dt*cd[m][n];
if (m != 0) tc[n][m-1] = c[n][m-1]+dt*cd[n][m-1];
}
/*
ACCUMULATE TERMS OF THE SPHERICAL HARMONIC EXPANSIONS
*/
par = ar**(p+n+m*13);
if (m == 0)
{
temp1 = tc[m][n]*cp[m];
temp2 = tc[m][n]*sp[m];
}
else
{
temp1 = tc[m][n]*cp[m]+tc[n][m-1]*sp[m];
temp2 = tc[m][n]*sp[m]-tc[n][m-1]*cp[m];
}
bt = bt-ar*temp1*dp[m][n];
bp += (fm[m]*temp2*par);
br += (fn[n]*temp1*par);
/*
SPECIAL CASE: NORTH/SOUTH GEOGRAPHIC POLES
*/
if (st == 0.0 && m == 1)
{
if (n == 1) pp[n] = pp[n-1];
else pp[n] = ct*pp[n-1]-k[m][n]*pp[n-2];
parp = ar*pp[n];
bpp += (fm[m]*temp2*parp);
}
}
}
if (st == 0.0) bp = bpp;
else bp /= st;
/*
ROTATE MAGNETIC VECTOR COMPONENTS FROM SPHERICAL TO
GEODETIC COORDINATES
*/
bx = -bt*ca-br*sa;
by = bp;
bz = bt*sa-br*ca;
/*
COMPUTE DECLINATION (DEC), INCLINATION (DIP) AND
TOTAL INTENSITY (TI)
*/
bh = sqrt((bx*bx)+(by*by));
*ti = sqrt((bh*bh)+(bz*bz));
*dec = atan2(by,bx)/dtr;
*dip = atan2(bz,bh)/dtr;
/*
COMPUTE MAGNETIC GRID VARIATION IF THE CURRENT
GEODETIC POSITION IS IN THE ARCTIC OR ANTARCTIC
(I.E. GLAT > +55 DEGREES OR GLAT < -55 DEGREES)

OTHERWISE, SET MAGNETIC GRID VARIATION TO -999.0
*/
*gv = -999.0;
if (fabs(glat) >= 55.)
{
if (glat > 0.0 && glon >= 0.0) *gv = *dec-glon;
if (glat > 0.0 && glon < 0.0) *gv = *dec+fabs(glon);
if (glat < 0.0 && glon >= 0.0) *gv = *dec+glon;
if (glat < 0.0 && glon < 0.0) *gv = *dec-fabs(glon);
if (*gv > +180.0) *gv -= 360.0;
if (*gv < -180.0) *gv += 360.0;
}
otime = time;
oalt = alt;
olat = glat;
olon = glon;
return;
}

/*************************************************************************/

void geomag(int *maxdeg)
{
E0000(0,maxdeg,0.0,0.0,0.0,0.0,NULL,NULL,NULL,NULL);
}

/*************************************************************************/

void geomg1(float alt, float glat, float glon, float time, float *dec, float *dip, float *ti, float *gv)
{
E0000(1,NULL,alt,glat,glon,time,dec,dip,ti,gv);
}

/*************************************************************************/

char geomag_introduction(float epochlowlim)
{
char help;
static char ans;

printf("

Welcome to the World Magnetic Model (WMM) %4.0f C-Program

", epochlowlim);
printf(" --- Version 2.0, September 2005 ---

");
printf("
This program estimates the strength and direction of ");
printf("
Earth's main magnetic field for a given point/area.");
printf("
Enter h for help and contact information or c to continue.");
printf ("
>");
scanf("%c%*[^
]",&help);
getchar();

if ((help == 'h') || (help == 'H'))
{
printf("
Help information ");

printf("
The World Magnetic Model (WMM) for %7.2f", epochlowlim);
printf("
is a model of Earth's main magnetic field. The WMM");
printf("
is recomputed every five (5) years, in years divisible by ");
printf("
five (i.e. 2005, 2010). See the contact information below");
printf("
to obtain more information on the WMM and associated software.");
printf("
");
printf("
Input required is the location in geodetic latitude and");
printf("
longitude (positive for northern latitudes and eastern ");
printf("
longitudes), geodetic altitude in meters, and the date of ");
printf("
interest in years.");

printf("


The program computes the estimated magnetic Declination");
printf("
(D) which is sometimes called MAGVAR, Inclination (I), Total");
printf("
Intensity (F or TI), Horizontal Intensity (H or HI), Vertical");
printf("
Intensity (Z), and Grid Variation (GV). Declination and Grid");
printf("
Variation are measured in units of degrees and are considered");
printf("
positive when east or north. Inclination is measured in units");
printf("
of degrees and is considered positive when pointing down (into");
printf("
the Earth). The WMM is reference to the WGS-84 ellipsoid and");
printf("
is valid for 5 years after the base epoch.");

printf("


It is very important to note that a degree and order 12 model,");
printf("
such as WMM, describes only the long wavelength spatial magnetic ");
printf("
fluctuations due to Earth's core. Not included in the WMM series");
printf("
models are intermediate and short wavelength spatial fluctuations ");
printf("
that originate in Earth's mantle and crust. Consequently, isolated");
printf("
angular errors at various positions on the surface (primarily over");
printf("
land, incontinental margins and over oceanic seamounts, ridges and");
printf("
trenches) of several degrees may be expected. Also not included in");
printf("
the model are temporal fluctuations of magnetospheric and ionospheric");
printf("
origin. On the days during and immediately following magnetic storms,");
printf("
temporal fluctuations can cause substantial deviations of the geomagnetic");
printf("
field from model values. If the required declination accuracy is");
printf("
more stringent than the WMM series of models provide, the user is");
printf("
advised to request special (regional or local) surveys be performed");
printf("
and models prepared. Please make requests of this nature to the");
printf("
National Geospatial-Intelligence Agency (NGA) at the address below.");

printf("


Contact Information");

printf("
Software and Model Support");
printf("
National Geophysical Data Center");
printf("
NOAA EGC/2");
printf("
325 Broadway");
printf("
Boulder, CO 80303 USA");
printf("
Attn: Susan McLean or Stefan Maus");
printf("
Phone: (303) 497-6478 or -6522");
printf("
Email: Susan.McLean@noaa.gov or Stefan.Maus@noaa.gov ");

printf("


Continue with program? (y or n) ");
scanf("%c%*[^
]", &ans);
getchar();
}
else
{
ans = 'y';
}

return(ans);
}
----

Grato pela atenção;
Bruno Lelli