readdata.c

二维热能方程的求解

/*** ReadDATA: for reading datafile ***/
/*** Gererated by Dr. Kim's CSTART  ***/

#include <stdio.h>
#include <strings.h>

typedef  double TYPE;
typedef  double ATYPE;
typedef  double VTYPE;
typedef  float DTYPE;

void readdata(char*datafile,int n1,int n2,DTYPE DATA[][n2],int level,int*ierr)
{
   FILE *fp000;
   int i,j;
   if (level >= 1)
      printf("ReadDATA: from <%s>\n",datafile);

   if ((fp000=fopen(datafile,"r"))==NULL){
       printf("Error in ReadDATA: when try to open <%s>\n",datafile);
       *ierr=1;return;}

   for(i=0;i<n1;i++)
   for(j=0;j<n2;j++)
      fscanf(fp000,"%f",&DATA[i][j]);               
   
   if(level>=3)
   for(i=0;i<n1;i++){
   for(j=0;j<n2;j++)
         printf(" %g",DATA[i][j]);
      printf("\n");
    }

   if(level>=5)
   printf(" ReadDATA: you may have to edit it.\n");

   fclose(fp000);
   return;
}