🚀 go-pugleaf

> > Hoi,
> >
> > Ik heb het volgende probleem. Als ik dit stukje programma laat lopen,
dan
> > geeft hij een waarde aan mn j van 4067520 (of toch zoiets in die grootte
> > orde...) Iemand aan idee waar hij die waarde vandaan haalt???
> > alvast bedankt
> >
> > Matthias
> >
> > -------------------------------------------
> > int j;
> >
> > FILE *ifp;
> >   FILE *ofp;
> >
> >   /*reading file*/
> >
> >   ifp = fopen("../gjelbreen","r");
> >
> >   /*reading data from file*/
> >
> >   for (j=1;j<=MAXX;j++)
> >
> >   {
> >
> >    fscanf(ifp,"%lf%lf%lf",&x[j],&hb[j],&s[j]);
> >    printf("%f",hb[j]);
> >
> >   }
> >
> >
> > exit(0);
> >
>
> 1) Bedoel je dat j die waarde heeft *na* het beeindigen van de loop en dat
> die waarde abnormaal is
> te opzichte van MAXX? Maar voor de exit(0)? Of is het na exit(0) ?
> 2) Dit stukje code is geen correct programma, geef eens jouw volledige
code.
> 3) Hoe kom je aan die waarde?
>
> Gokje:
> Na exit(0) heeft j geen gedefinieerde waarde meer, maar sommige debuggers
> zullen wel nog iets tonen bijv: 4067520
>
Idd, 't is zoiets! Ik ken er allemaal niet veel van, maar ik wou mn
programma stuk per stuk testen, dit door die exit() toe te voegen. Hieronder
staat de rest van mijn programma, alle, tot aan die exit toch....

/*standard libraries*/

#include <stdio.h>  /*standard input-output*/
#include <stdlib.h> /*standard C library*/
#include <math.h>  /*math library*/


/*definieren van de parameters*/


#define GA   9.81      /*gravity acceleration (m/s�)*/
#define K       (2.1*SEC_PER_YEAR)       /*heat conductivity (W/mK)*/
#define RHO      910.0        /*ice density (kg/m�)*/
#define G       (0.042*SEC_PER_YEAR)     /*geothermal heat flux (*/
#define CP       2009.0       /*specific heat capacity (J/kgK)*/
#define TABS     273.15       /* absolute temperature (K)*/
#define TSEA                  /*temperature at sea level (K)*/
#define BETA  0.000000098     /*constant parameter (K/Pa)*/
#define ADIAT     (-0.012)       /*Factor regarding temperature/height
(K/m)*/
#define A        1.0e-16                     /*flow parameter */

#define TIMETOT  50000.0       /*total calculation time (year)*/
#define MAXX  122       /*maximum number of elements*/
#define MAXXARRAY   (MAXX+1)
#define MAXZ  51
#define MAXZARRAY (MAXZ+1)
#define SEC_PER_YEAR 31556926     /*numbers of seconds in 1 year*/
#define DT   10.0
#define DZETA  0.02      /*interval on the vertical axes*/
#define DX    10.000      /*interval on the horizontal axes (m)*/




/*main program*/

int main(int argc, char **argv)   /*arguments to communicate with the main
system*/

 {
  int i,j,k,tlapse;
  double t0[MAXZARRAY][MAXXARRAY], t1[MAXZARRAY][MAXXARRAY],
w[MAXZARRAY][MAXXARRAY], u[MAXZARRAY][MAXXARRAY];
  double h[MAXXARRAY];      /*ice sheet thickness*/
  double s[MAXXARRAY];      /*Height of surface*/
  double hb[MAXXARRAY];      /*height of basis of glacier*/
  double x[MAXXARRAY];
  double m[MAXXARRAY];
  double zeta[MAXZARRAY];
  double tsurf[MAXXARRAY];
  double tpmp[MAXZARRAY][MAXXARRAY];
  double wsurf[MAXXARRAY];
  double f[MAXZARRAY];
  double g[MAXZARRAY];
  double alfa[MAXZARRAY];
  double beta[MAXZARRAY];
  double gamma[MAXZARRAY];
  double delta[MAXZARRAY];
  double taud[MAXXARRAY];

  FILE *ifp;
  FILE *ofp;

  /*reading file*/


  ifp = fopen("../gjelbreen","r");

  /*reading data from file*/

  for (j=1;j<=MAXX;j++)

  {

   fscanf(ifp,"%lf%lf%lf",&x[j],&hb[j],&s[j]);
   printf("%f",hb[j]);

  }


exit(0);