計程作業第五章-2

[b7505002 (熱~~~~~~)]

#include <stdio.h>
#include <process.h>
double read_data(double *, double *);
int main()
{
 double i=0.,A=0.,C,Q;
 C=read_data(&i,&A);
 Q=C*i*A;
 printf("Total watershed area A=%lf\n",A);
 printf("Weighted average runoff coeff. C=%lf\n",C);
 printf("Rate of runoff Q=%lf\n",Q);
 return 0;
 }

 double read_data(double *i,double *A)
 {
  FILE *pf;
  double ai,ci,sumc=0.0;
  int j;
  pf=fopen("runoff.dat","r");
  if(pf==NULL){
         printf("Can not open file \n");
         exit(1);
  }
  fscanf(pf,"%lf",i);
  printf("Rainfall intensity i=%lf\n",*i);
  *A=0.0;
  for(j=0; j<4; ++j){
         fscanf(pf,"%lf %lf",&ai,&ci);
         printf(" ai= %lf  ci=%lf\n",ai,ci);
         sumc=sumc+ai*ci;
         *A=*A+ai;
  }
  return sumc/(*A);
  fclose(pf);
 }


--
※ 發信站: 批踢踢實業坊(ptt.twbbs.org) 
◆ From: pc86.na.ntu.edu