EPANET中读取INPUT文件的函数文件——INPUT1.C/INPUT2.C/INPUT3.C
首先介绍下这3个文件的关系：可以说INPUT1.C的函数粒度最大，它的函数getdata()就完成了整个INPUT文件数据的读入，该函数又调用了INPUT2.C中的部分函数，INPUT2.C文件中的函数分为两类，一类是解析INPUT文件的函数；另一类是工具函数，包括将节点、管段数据存入哈希表中，构建哈希索引、字符串的匹配，从哈希表中查找数据等，该文件的粒度中等。而INPUT3.C的粒度是最细的，他是用于解析单行INPUT文件的对象内容并保存到内存中。INPUT3.C的组织结构也很简单,各个函数具有较大的相似性，分别是对INPUT文件中的不同段落数据的处理。
其中INPUT2.C文件中的unlinked()函数是检查孤立管点，这个数据检查的算法比较高效值得借鉴，该算法的巧妙之处可以参考代码中的注释：
1 int unlinked() 2 /* 3 **-------------------------------------------------------------- 4 ** Input: none ；输入：无 5 ** Output: returns error code if any unlinked junctions found ；输出：如果存在孤立点则返回错误代码 6 ** Purpose: checks for unlinked junctions in network ；目的：检查管网中的孤立点 7 ** 8 ** NOTE: unlinked tanks have no effect on computations. ；注意：没有连接的蓄水池对计算无效 9 **-------------------------------------------------------------- 10 */ 11 { 12 char *marked; 13 int i,err, errcode; 14 errcode = 0; 15 err = 0; 16 marked = (char *) calloc(Nnodes+1,sizeof(char)); //开辟连续的所有节点内存空间 17 ERRCODE(MEMCHECK(marked)); 18 if (!errcode) 19 { 20 memset(marked,0,(Nnodes+1)*sizeof(char));//给这块内存空间置0 21 for (i=1; i<=Nlinks; i++) /* Mark end nodes of each link */ 22 { 23 marked[Link[i].N1]++;//以marked的内存的逻辑地址从1、2、3...开始去一一对应管段节点在Node中的索引值，从而标记出所有被管段引用作为始末节点的索引值，而所有未标记成功的就是孤立的点。 24 marked[Link[i].N2]++; 25 } 26 for (i=1; i<=Njuncs; i++) /* Check each junction */ 27 { 28 if (marked[i] == 0) /* If not marked then error */ 29 { 30 err++; 31 sprintf(Msg,ERR233,Node[i].ID); 32 writeline(Msg); 33 } 34 if (err >= MAXERRS) break; 35 } 36 if (err > 0) errcode = 200; 37 } 38 free(marked); 39 return(errcode); 40 } /* End of unlinked */
在INPUT1.C这个文件中，我们重点看下getdata()，就知道了它的功能：
1 int getdata() 2 /* 3 **---------------------------------------------------------------- 4 ** Input: none ；输入：无 5 ** Output: returns error code ；输出：错误编号（大于100表示严重错误） 6 ** Purpose: reads in network data from disk file；作用：从磁盘文件（*.inp）中读入管网模型数据 7 **---------------------------------------------------------------- 8 */ 9 { 10 int errcode = 0; 11 setdefaults(); /* Assign default data values ；给vars.h中的部分全局变量赋值*/ 12 initreport(); /* Initialize reporting options ；初始化报表选项*/ 13 rewind(InFile); /* Rewind input file ；将INPUT文件的指针指向文件的开始*/ 14 ERRCODE(readdata()); /* Read in network data ；读取整个管网数据，该部分的具体实现是依赖INPUT2.C与INPUT3.C*/ 15 if (!errcode) adjustdata(); /* Adjust data for default values ；调整数据的默认值，当INPUT文件对部分关键参数缺少描述时，这里将重新赋值*/ 16 if (!errcode) initunits(); /* Initialize units on input data ；确定单位转换因子，EPANET中的单位体系就2种，一种是英美制US，一种是米制，根据确定的单位体系，来设置各节点、管段等的字段的转换参数*/ 17 ERRCODE(inittanks()); /* Initialize tank volumes ；初始化水库的容积*/ 18 if (!errcode) convertunits(); /* Convert units on input data ；在进行水力解算前对管段、节点、水泵、阀门的相关字段的单位数据进行转换*/ 19 return(errcode); 20 } /* End of getdata */
