炒股，大多数人看见有人炒股票挣钱了都非常想试试，但是要做好股票交易，真的是很困难的一件事情，所以有朋友在问如何做股票这个问题，小编真的不知道如何说起，因此写了《sss988》这篇文章，从基本的知识开始学习吧！
现在，年轻女性喜欢喝酒的越来越多，这是正常的。不过，女性在月经来临时或月经期间，最好别喝酒，否则更容易造成身体伤害。但是很多人还是不注意，下面还是说说月经喝酒的危害吧，来看看喝酒对月经有什么影响吗？
喝酒对月经有什么影响吗
月经来临前及期间，女性受荷尔蒙分泌影响，体内分解酶的活动能力低下，酒精代谢能力下降，结果使得酒精不易迅速从血液中排泄出去，而是变成了对身体有害的“酸性物质”。为清除这些酸性物质，肝脏就要不断地制造出酶。其结果就会加重肝脏的负担，使引发肝脏机能障碍的可能性也加大，。
有报道称，同样是喝酒，女性经期饮酒引发肝损害或酒精中毒的几率将比男性多一半。因为是女性在月经期间，体内缺乏分解酶，如果一时喝得过多，将使处于醉酒状态的时间延长、酒醉感觉或症状也会更严重。这就是月经时饮酒容易上瘾、容易引发酒精中毒原因。
另外，经期由于不断流血，身体虚弱，抵抗力较差，喝酒会加快血液循环，此时有可能导致月经量增多，如饮凉啤酒，还可能引起痛经等。所以，月经临近或月经期间，原则上应当禁饮白酒。当然，可以少量喝点葡萄酒(50ml左右为宜)，但也不宜过量。
女性喝酒对月经有什么影响
1、女性在月经期间喝酒的话，就会增加肝脏的负担，因为女性受到荷尔蒙分泌的影响，体内的一些活动能力低下，所以酒精代谢能力也会下降的结果，导致酒精不能从血液里排泄出现，而就会对身体产生一种有害的酸性物质，结果就会加重了肝脏的负担。
2、在月经期间受到内分泌功能的影响，肝脏代谢酒精能力就会降低的，所以导致酒精残留在体内的时间非常的长，不断的喝酒就会造成酒精中毒的几率是非常大的，所以很容易导致出现酒精中毒的情况，月经期间身体是比较虚弱的，抵抗力也是非常的差，月经期间喝酒就会加快血液循环，会导致月经量增多。
3、在月经期间喝酒的话就会导致月经血液排出，去的时候受到一定的阻碍，所以会导致出现痛经的情况，所以女性朋友在月经期间不要喝酒，不然会给身体带来很多的危害，严重的话会导致女性出现不孕的情况。
注意事项：
在月经期间不仅不能喝酒，也不能喝浓茶，不能喝咖啡，不能碰一些冰凉的东西，禁止吃一些辛辣刺激性的食物，一定要保持乐观的心态，保障充足的睡眠，避免过于劳累，注意劳逸结合，月经期间避免剧烈的运动。
1 过两点有且只有一条直线
2 两点之间线段最短
3 同角或等角的补角相等
4 同角或等角的余角相等
5 过一点有且只有一条直线和已知直线垂直
6 直线外一点与直线上各点连接的所有线段中，垂线段最短
7 平行公理 经过直线外一点，有且只有一条直线与这条直线平行
8 如果两条直线都和第三条直线平行，这两条直线也互相平行
9 同位角相等，两直线平行
10 内错角相等，两直线平行
11 同旁内角互补，两直线平行
12两直线平行，同位角相等
13 两直线平行，内错角相等
14 两直线平行，同旁内角互补
15 定理 三角形两边的和大于第三边
16 推论 三角形两边的差小于第三边
17 三角形内角和定理 三角形三个内角的和等于180°
18 推论1 直角三角形的两个锐角互余
19 推论2 三角形的一个外角等于和它不相邻的两个内角的和
20 推论3 三角形的一个外角大于任何一个和它不相邻的内角
21 全等三角形的对应边、对应角相等
22边角边公理(SAS) 有两边和它们的夹角对应相等的两个三角形全等
23 角边角公理( ASA)有两角和它们的夹边对应相等的两个三角形全等
24 推论(AAS) 有两角和其中一角的对边对应相等的两个三角形全等
25 边边边公理(SSS) 有三边对应相等的两个三角形全等
26 斜边、直角边公理(HL) 有斜边和一条直角边对应相等的两个直角三角形全等
27 定理1 在角的平分线上的点到这个角的两边的距离相等
28 定理2 到一个角的两边的距离相同的点，在这个角的平分线上
29 角的平分线是到角的两边距离相等的所有点的集合
30 等腰三角形的性质定理 等腰三角形的两个底角相等 (即等边对等角）
31 推论1 等腰三角形顶角的平分线平分底边并且垂直于底边
32 等腰三角形的顶角平分线、底边上的中线和底边上的高互相重合
33 推论3 等边三角形的各角都相等，并且每一个角都等于60°
34 等腰三角形的判定定理 如果一个三角形有两个角相等，那么这两个角所对的边也相等（等角对等边）
35 推论1 三个角都相等的三角形是等边三角形
36 推论 2 有一个角等于60°的等腰三角形是等边三角形
37 在直角三角形中，如果一个锐角等于30°那么它所对的直角边等于斜边的一半
38 直角三角形斜边上的中线等于斜边上的一半
39 定理 线段垂直平分线上的点和这条线段两个端点的距离相等
40 逆定理 和一条线段两个端点距离相等的点，在这条线段的垂直平分线上
41 线段的垂直平分线可看作和线段两端点距离相等的所有点的集合
42 定理1 关于某条直线对称的两个图形是全等形
43 定理 2 如果两个图形关于某直线对称，那么对称轴是对应点连线的垂直平分线
44定理3 两个图形关于某直线对称，如果它们的对应线段或延长线相交，那么交点在对称轴上
45逆定理 如果两个图形的对应点连线被同一条直线垂直平分，那么这两个图形关于这条直线对称
46勾股定理 直角三角形两直角边a、b的平方和、等于斜边c的平方，即a^2+b^2=c^2
47勾股定理的逆定理 如果三角形的三边长a、b、c有关系a^2+b^2=c^2 ，那么这个三角形是直角三角形
48定理 四边形的内角和等于360°
49四边形的外角和等于360°
50多边形内角和定理 n边形的内角的和等于（n-2）×180°
51推论 任意多边的外角和等于360°
52平行四边形性质定理1 平行四边形的对角相等
53平行四边形性质定理2 平行四边形的对边相等
54推论 夹在两条平行线间的平行线段相等
55平行四边形性质定理3 平行四边形的对角线互相平分
56平行四边形判定定理1 两组对角分别相等的四边形是平行四边形
57平行四边形判定定理2 两组对边分别相等的四边形是平行四边形
58平行四边形判定定理3 对角线互相平分的四边形是平行四边形
59平行四边形判定定理4 一组对边平行相等的四边形是平行四边形
60矩形性质定理1 矩形的四个角都是直角 61矩形性质定理2 矩形的对角线相等
62矩形判定定理1 有三个角是直角的四边形是矩形
63矩形判定定理2 对角线相等的平行四边形是矩形
64菱形性质定理1 菱形的四条边都相等
65菱形性质定理2 菱形的对角线互相垂直，并且每一条对角线平分一组对角
66菱形面积=对角线乘积的一半，即S=（a×b）÷2
67菱形判定定理1 四边都相等的四边形是菱形
68菱形判定定理2 对角线互相垂直的平行四边形是菱形
69正方形性质定理1 正方形的四个角都是直角，四条边都相等
70正方形性质定理2正方形的两条对角线相等，并且互相垂直平分，每条对角线平分一组对角
71定理1 关于中心对称的两个图形是全等的
72定理2 关于中心对称的两个图形，对称点连线都经过对称中心，并且被对称中心平分
73逆定理 如果两个图形的对应点连线都经过某一点，并且被这一
点平分，那么这两个图形关于这一点对称
74等腰梯形性质定理 等腰梯形在同一底上的两个角相等
75等腰梯形的两条对角线相等
76等腰梯形判定定理 在同一底上的两个角相等的梯形是等腰梯形
77对角线相等的梯形是等腰梯形
78平行线等分线段定理 如果一组平行线在一条直线上截得的线段
相等，那么在其他直线上截得的线段也相等
79 推论1 经过梯形一腰的中点与底平行的直线，必平分另一腰
80 推论2 经过三角形一边的中点与另一边平行的直线，必平分第
三边
81 三角形中位线定理 三角形的中位线平行于第三边，并且等于它
的一半
82 梯形中位线定理 梯形的中位线平行于两底，并且等于两底和的
一半 L=（a+b）÷2 S=L×h
83 (1)比例的基本性质 如果a:b=c:d,那么ad=bc
如果ad=bc,那么a:b=c:d
84 (2)合比性质 如果a／b=c／d,那么(a±b)／b=(c±d)／d
85 (3)等比性质 如果a／b=c／d=…=m／n(b+d+…+n≠0),那么
(a+c+…+m)／(b+d+…+n)=a／b
86 平行线分线段成比例定理 三条平行线截两条直线，所得的对应
线段成比例
87 推论 平行于三角形一边的直线截其他两边（或两边的延长线），所得的对应线段成比例
88 定理 如果一条直线截三角形的两边（或两边的延长线）所得的对应线段成比例，那么这条直线平行于三角形的第三边
89 平行于三角形的一边，并且和其他两边相交的直线，所截得的三角形的三边与原三角形三边对应成比例
90 定理 平行于三角形一边的直线和其他两边（或两边的延长线）相交，所构成的三角形与原三角形相似
91 相似三角形判定定理1 两角对应相等，两三角形相似（ASA）
92 直角三角形被斜边上的高分成的两个直角三角形和原三角形相似
93 判定定理2 两边对应成比例且夹角相等，两三角形相似（SAS）
94 判定定理3 三边对应成比例，两三角形相似（SSS）
95 定理 如果一个直角三角形的斜边和一条直角边与另一个直角三
角形的斜边和一条直角边对应成比例，那么这两个直角三角形相似
96 性质定理1 相似三角形对应高的比，对应中线的比与对应角平
分线的比都等于相似比
97 性质定理2 相似三角形周长的比等于相似比
98 性质定理3 相似三角形面积的比等于相似比的平方
99 任意锐角的正弦值等于它的余角的余弦值，任意锐角的余弦值等
于它的余角的正弦值
100任意锐角的正切值等于它的余角的余切值，任意锐角的余切值等
于它的余角的正切值
101圆是定点的距离等于定长的点的集合
102圆的内部可以看作是圆心的距离小于半径的点的集合
103圆的外部可以看作是圆心的距离大于半径的点的集合
104同圆或等圆的半径相等
105到定点的距离等于定长的点的轨迹，是以定点为圆心，定长为半
径的圆
106和已知线段两个端点的距离相等的点的轨迹，是着条线段的垂直
平分线
107到已知角的两边距离相等的点的轨迹，是这个角的平分线
108到两条平行线距离相等的点的轨迹，是和这两条平行线平行且距离相等的一条直线
109定理 不在同一直线上的三点确定一个圆。
110垂径定理 垂直于弦的直径平分这条弦并且平分弦所对的两条弧
111推论1 ①平分弦（不是直径）的直径垂直于弦，并且平分弦所对的两条弧
②弦的垂直平分线经过圆心，并且平分弦所对的两条弧
③平分弦所对的一条弧的直径，垂直平分弦，并且平分弦所对的另一条弧
112推论2 圆的两条平行弦所夹的弧相等
113圆是以圆心为对称中心的中心对称图形
114定理 在同圆或等圆中，相等的圆心角所对的弧相等，所对的弦
相等，所对的弦的弦心距相等
115推论 在同圆或等圆中，如果两个圆心角、两条弧、两条弦或两
弦的弦心距中有一组量相等那么它们所对应的其余各组量都相等
116定理 一条弧所对的圆周角等于它所对的圆心角的一半
117推论1 同弧或等弧所对的圆周角相等；同圆或等圆中，相等的圆周角所对的弧也相等
118推论2 半圆（或直径）所对的圆周角是直角；90°的圆周角所
对的弦是直径
119推论3 如果三角形一边上的中线等于这边的一半，那么这个三角形是直角三角形
120定理 圆的内接四边形的对角互补，并且任何一个外角都等于它
的内对角
121①直线L和⊙O相交 d＜r
②直线L和⊙O相切 d=r
③直线L和⊙O相离 d＞r
122切线的判定定理 经过半径的外端并且垂直于这条半径的直线是圆的切线
123切线的性质定理 圆的切线垂直于经过切点的半径
124推论1 经过圆心且垂直于切线的直线必经过切点
125推论2 经过切点且垂直于切线的直线必经过圆心
126切线长定理 从圆外一点引圆的两条切线，它们的切线长相等，
圆心和这一点的连线平分两条切线的夹角
127圆的外切四边形的两组对边的和相等
128弦切角定理 弦切角等于它所夹的弧对的圆周角
129推论 如果两个弦切角所夹的弧相等，那么这两个弦切角也相等
130相交弦定理 圆内的两条相交弦，被交点分成的两条线段长的积
相等
131推论 如果弦与直径垂直相交，那么弦的一半是它分直径所成的
两条线段的比例中项
132切割线定理 从圆外一点引圆的切线和割线，切线长是这点到割
线与圆交点的两条线段长的比例中项
133推论 从圆外一点引圆的两条割线，这一点到每条割线与圆的交点的两条线段长的积相等
134如果两个圆相切，那么切点一定在连心线上
135①两圆外离 d＞R+r ②两圆外切 d=R+r
③两圆相交 R-r＜d＜R+r(R＞r)
④两圆内切 d=R-r(R＞r) ⑤两圆内含d＜R-r(R＞r)
136定理 相交两圆的连心线垂直平分两圆的公共弦
137定理 把圆分成n(n≥3):
⑴依次连结各分点所得的多边形是这个圆的内接正n边形
⑵经过各分点作圆的切线，以相邻切线的交点为顶点的多边形是这个圆的外切正n边形
138定理 任何正多边形都有一个外接圆和一个内切圆，这两个圆是同心圆
139正n边形的每个内角都等于（n-2）×180°／n
140定理 正n边形的半径和边心距把正n边形分成2n个全等的直角三角形
141正n边形的面积Sn=pnrn／2 p表示正n边形的周长
142正三角形面积√3a／4 a表示边长
143如果在一个顶点周围有k个正n边形的角，由于这些角的和应为
360°，因此k×(n-2)180°／n=360°化为（n-2）(k-2)=4
144弧长计算公式：L=n兀R／180
145扇形面积公式：S扇形=n兀R^2／360=LR／2
146内公切线长= d-(R-r) 外公切线长= d-(R+r)
实用工具:常用数学公式
公式分类 公式表达式
乘法与因式分 a2-b2=(a+b)(a-b) a3+b3=(a+b)(a2-ab+b2) a3-b3=(a-b(a2+ab+b2)
三角不等式 |a+b|≤|a|+|b| |a-b|≤|a|+|b| |a|≤b=-b≤a≤b
|a-b|≥|a|-|b| -|a|≤a≤|a|
一元二次方程的解 -b+√(b2-4ac)/2a -b-√(b2-4ac)/2a
根与系数的关系 X1+X2=-b/a X1*X2=c/a 注：韦达定理
判别式
b2-4ac=0 注：方程有两个相等的实根
b2-4ac0 注：方程有两个不等的实根
b2-4ac0 注：方程没有实根，有共轭复数根
三角函数公式
两角和公式
sin(A+B)=sinAcosB+cosAsinB sin(A-B)=sinAcosB-sinBcosA
cos(A+B)=cosAcosB-sinAsinB cos(A-B)=cosAcosB+sinAsinB
tan(A+B)=(tanA+tanB)/(1-tanAtanB) tan(A-B)=(tanA-tanB)/(1+tanAtanB)
ctg(A+B)=(ctgActgB-1)/(ctgB+ctgA) ctg(A-B)=(ctgActgB+1)/(ctgB-ctgA)
倍角公式
tan2A=2tanA/(1-tan2A) ctg2A=(ctg2A-1)/2ctga
cos2a=cos2a-sin2a=2cos2a-1=1-2sin2a
半角公式
sin(A/2)=√((1-cosA)/2) sin(A/2)=-√((1-cosA)/2)
cos(A/2)=√((1+cosA)/2) cos(A/2)=-√((1+cosA)/2)
tan(A/2)=√((1-cosA)/((1+cosA)) tan(A/2)=-√((1-cosA)/((1+cosA))
ctg(A/2)=√((1+cosA)/((1-cosA)) ctg(A/2)=-√((1+cosA)/((1-cosA))
和差化积
2sinAcosB=sin(A+B)+sin(A-B) 2cosAsinB=sin(A+B)-sin(A-B)
2cosAcosB=cos(A+B)-sin(A-B) -2sinAsinB=cos(A+B)-cos(A-B)
sinA+sinB=2sin((A+B)/2)cos((A-B)/2 cosA+cosB=2cos((A+B)/2)sin((A-B)/2)
tanA+tanB=sin(A+B)/cosAcosB tanA-tanB=sin(A-B)/cosAcosB
ctgA+ctgBsin(A+B)/sinAsinB -ctgA+ctgBsin(A+B)/sinAsinB
某些数列前n项和
1+2+3+4+5+6+7+8+9+…+n=n(n+1)/2 1+3+5+7+9+11+13+15+…+(2n-1)=n2
2+4+6+8+10+12+14+…+(2n)=n(n+1) 12+22+32+42+52+62+72+82+…+n2=n(n+1)(2n+1)/6
13+23+33+43+53+63+…n3=n2(n+1)2/4 1*2+2*3+3*4+4*5+5*6+6*7+…+n(n+1)=n(n+1)(n+2)/3
正弦定理 a/sinA=b/sinB=c/sinC=2R 注： 其中 R 表示三角形的外接圆半径
余弦定理 b2=a2+c2-2accosB 注：角B是边a和边c的夹角
圆的标准方程 (x-a)2+(y-b)2=r2 注：（a,b）是圆心坐标
圆的一般方程 x2+y2+Dx+Ey+F=0 注：D2+E2-4F0
抛物线标准方程 y2=2px y2=-2px x2=2py x2=-2py
直棱柱侧面积 S=c*h 斜棱柱侧面积 S=c'*h
正棱锥侧面积 S=1/2c*h' 正棱台侧面积 S=1/2(c+c')h'
圆台侧面积 S=1/2(c+c')l=pi(R+r)l 球的表面积 S=4pi*r2
圆柱侧面积 S=c*h=2pi*h 圆锥侧面积 S=1/2*c*l=pi*r*l
弧长公式 l=a*r a是圆心角的弧度数r 0 扇形面积公式 s=1/2*l*r
锥体体积公式 V=1/3*S*H 圆锥体体积公式 V=1/3*pi*r2h
斜棱柱体积 V=S'L 注：其中,S'是直截面面积， L是侧棱长
柱体体积公式 V=s*h 圆柱体 V=pi*r2h
圆锥曲线包括椭圆，双曲线，抛物线
1. 椭圆：到两个定点的距离之和等于定长（定长大于两个定点间的距离）的动点的轨迹叫做椭圆。即：{P| |PF1|+|PF2|=2a, (2a|F1F2|)}。
2. 双曲线：到两个定点的距离的差的绝对值为定值（定值小于两个定点的距离）的动点轨迹叫做双曲线。即{P|||PF1|-|PF2||=2a, (2a|F1F2|)}。
3. 抛物线：到一个定点和一条定直线的距离相等的动点轨迹叫做抛物线。
4. 圆锥曲线的统一定义：到定点的距离与到定直线的距离的比e是常数的点的轨迹叫做圆锥曲线。当0e1时为椭圆：当e=1时为抛物线；当e1时为双曲线。
·圆锥曲线由来：圆，椭圆，双曲线，抛物线同属于圆锥曲线。早在两千多年前，古希腊数学家对它们已经很熟悉了。古希腊数学家阿波罗尼采用平面切割圆锥的方法来研究这几种曲线。用垂直与锥轴的平面去截圆锥，得到的是圆；把平面渐渐倾斜，得到椭圆；当平面和圆锥的一条母线平行时，得到抛物线；当平面再倾斜一些就可以得到双曲线。阿波罗尼曾把椭圆叫“亏曲线”，把双曲线叫做“超曲线”，把抛物线叫做“齐曲线”。
·圆锥曲线的参数方程和直角坐标方程：
1）直线
参数方程：x=X+tcosθ y=Y+tsinθ (t为参数）
直角坐标：y=ax+b
2）圆
参数方程：x=X+rcosθ y=Y+rsinθ (θ为参数 )
直角坐标：x^2+y^2=r^2 (r 为半径）
3）椭圆
参数方程：x=X+acosθ y=Y+bsinθ (θ为参数 )
直角坐标（中心为原点）：x^2/a^2 + y^2/b^2 = 1
4）双曲线
参数方程：x=X+asecθ y=Y+btanθ (θ为参数 )
直角坐标（中心为原点）：x^2/a^2 - y^2/b^2 = 1 (开口方向为x轴） y^2/a^2 - x^2/b^2 = 1 (开口方向为y轴）
5）抛物线
参数方程：x=2pt^2 y=2pt (t为参数)
直角坐标：y=ax^2+bx+c (开口方向为y轴, a0 ） x=ay^2+by+c （开口方向为x轴, a0 )
圆锥曲线（二次非圆曲线）的统一极坐标方程为
ρ=ep/(1-e·cosθ)
其中e表示离心率，p为焦点到准线的距离。
我是高考过来的，一般我们省是自主命题，最后一道大题通常就是圆锥曲线的综合型题目，这种题目的分值大约18分左右但是计算量相当的巨大，一般会设几个小问题，建议楼主视自己的情况而定，有取舍的做这些题目，而所谓的重点就是平常练习中的熟练程度了，高考的数学还是考察个人的解题熟练程度，所以想要取得高分还是要做一些有代表性的题目在注意总结考120以上应该没有问题，最后祝你金榜题名哈！ 1、空间一点P位于平面MAB的充要条件是存在唯一的有序实数对x、y，使得PM=xPA+yPB(其中PM等为向量，由于图不方便做就如此代替，下同）
2、对空间任一点O和不共线的三点A，B，C，若：OP=xOA+yOB+zOC （其中x＋y＋z=1），则四点P、A、B、C共面．
3、利用向量证a‖b，就是分别在a，b上取向量 （k∈R）．
4、利用向量证在线a⊥b，就是分别在a，b上取向量 ．
5、利用向量求两直线a与b的夹角，就是分别在a，b上取 ，求： 的问题．
6、利用向量求距离就是转化成求向量的模问题： ．
7、利用坐标法研究线面关系或求角和距离，关键是建立正确的空间直角坐标系，正确表达已知点的坐标．
首先该图形能建坐标系
如果能建
则先要会求面的法向量
求面的法向量的方法是
1。尽量在空中找到与面垂直的向量
2。如果找不到，那么就设n=（x,y,z）
然后因为法向量垂直于面
所以n垂直于面内两相交直线
可列出两个方程
两个方程，三个未知数
然后根据计算方便
取z（或x或y）等于一个数
然后就求出面的一个法向量了
会求法向量后
1。二面角的求法就是求出两个平面的法向量
可以求出两个法向量的夹角为两向量的数量积除以两向量模的乘积 ：cosa,b=|n·n1|/|n|
如过在两面的同一边可以看到两向量的箭头或箭尾相交
那么二面角就是上面求的两法向量的夹角的补角
2。点到平面的距离就是求出该面的法向量 在平面上任取(除被求点在该平面的射影外）一点，
求出平面外那点和你所取的那点所构成的向量记为n1
点到平面的距离就是法向量与n1的数量积的绝对值除以法向量的模即得所求
设直线l,m的方向向量分别为a,b，平面α,β的法向量分别为μ,ν 则
线线平行 l‖m = a‖b = a=kb；
线面平行 l‖α = a⊥μ = a·μ=0；
面面平行 α‖β = μ‖ν = μ=kν
线线垂直 l⊥m = a⊥b =a·b=0；
线面垂直 l⊥α = a‖μ = a=kμ；
面面垂直 α⊥β = μ⊥ν = μ·ν=0
很高兴认识你，祝你天天开心!
c feng.for五层准三维流（潜/弱/承/弱/承，对弱透水层只考虑越流不考虑弹性释放），调参用
c 使潜、承、承上下节点、上下单元完全重合，且总单元数ms＝3×ms1，总节点数ids＝3×ids1
c 潜单元编号由1-ms1，中承单元编号由ms1+1-ms2，深承单元编号由ms2+1-ms
c 10-13行：可调数组的各个参数预先赋值，不同计算区只要改变这些参数即可用此程序，其中潜单元个数＝ms1个，中承单元个数＝ms2-ms1个，深承单元个数＝ms-ms2个；潜节点编号由1-ids1，中承节点编号由ids1+1ids2，深承节点编号由ids2+1-ids，n：未知水头节点数；mi：开采节点数；igs：拟合点总数；igs1：潜拟合点数，igs2igs1：中承拟合点数，icqs：参数区个数：ihv：计算时段数，jbn：纯外引地表水灌溉节点个数；nn1：潜已知水头节点个数；nn：潜+中+深已知水头节点总数
parameter（ms1＝108，ids1＝67，n1＝59，mi1＝17，igs1＝11，icqs1＝8）
parameter（ms2＝216，ids2＝134，n2＝126，mi2＝47，igs2＝25，icqs2＝16）
parameter（ms＝324，ids＝201，n＝193，mi＝69，igs＝35，icqs＝24）
parameter（ihv＝24，jbn＝9，nn1＝8，nn＝24）
c h0：时段初水头（m）；hed：时段末水头（m）；in：各单元三节点编号（必须按小中大顺序排）
c fh：拟合点各时段计算水位（m）；sh：拟合点各时段实测水位（m）
dimension h0（ids），hed（ids），in（ms，3），fh（igs，ihv），sh（igs，ihv）
c idyh：导水矩阵工作单元；zb：各节点x，y坐标（输入坐标为图面mm数）；igdh：拟合点节点号；q：开采井各时段水量（抽为正，注为负，流出边界为正，流入为负）（m3/d）
dimension idyh（ids，9），zb（ids，2），igdh（igs），q（mi，ihv）
c d：各节点导水矩阵；ifdh：开采井节点号；s：拟合点拟合误差（m）；mqh：各单元的参数区号；e：各节点释（储）水矩阵dimension ifdh（mi），d（ids，9），s（igs，ihv），mqh（ms），e（ids）
c fqc：各参数区最多4个参数之值；sss：各节点水头总变幅，由t0时刻起算，（m）；h：初始流场，gh潜水层独有的灌溉入渗系数
dimension fqc（icqs，4），sss（ids），h（ids），gh（icqs1）
c yo：各节点越流矩阵工作单元；jbiao：纯引地表水灌溉的节点编号，共有jbn个节点；qxia：灌溉抽取地下水均铺在灌溉区的水层厚度；gq：灌溉抽出的地下水水层厚度；gq＝qxia；w：各时段潜水各节点的降水灌溉回归水量
dimension yo（ids），jbiao（jbn），qxia（ihv），gq（ihv），w（ids1，ihv）
c zzz：潜水各节点底板标高（m）；ep0：潜水各节点含水层厚（m）；ep1：潜水各单元含水层厚均值（m）；x：各时段降水量（m/时段）
dimension zzz（ids1），ep0（ids1），ep1（ms1），x（ihv）
c tl：各时段步长值（天）；bc：三角单元几何量（bi，ci，bj，cj，bk，ck）
dimension tl（ihv），bc（3，2）
c 分别是潜（h0nn1），中（h0nn2），深（h0nn3）层已知变水头节点各时段的水头
dimension h0nn1（nn1，ihv），h0nn2（nn1，ihv），h0nn3（nn1，ihv）
c file1（igs）*8：定义8字符的字符串igs个；filee*8：定义8字符的文件名
character file1（igs）*8，file2（igs）*8，filee*8
write（*，23）
23 format（1x，’zz1＝？zz2＝？迭代因子：zz1潜水0—1，zz2承压水1—2’）
c zz1：潜水为亚松弛系数，一般取0.85为好；zz2：承压水为超松弛系数，当ids＜300时取1.2-1.3，当ids＞300时取1.3-1.5为好
zz1＝0.85
zz2＝1.25
write（*，24）ms1，ids1，n1，mi1，igs1，icqs1
write（*，24）ms2，ids2，n2，mi2，igs2，icqs2，ihv
write（*，24）ms，ids，n，mi，igs，icqs
24 format（5x，12i5）
open（1，file＝’fqc’，status＝’old’）
c 首先从’fqc’中读取潜水的区号m10（空读），参数1（K），参数2（μ），参数3潜-中的越流系数（k’/m’），参数4降水入渗系数（α），gh灌溉入渗系数
read（1，*）（m10，fqc（i，1），fqc（i，2），fqc（i，3），fqc（i，4），gh（i），
* i＝1，icqs1）
c 接着从’fqc’中读取中层水的区号m10（空读），参数1（T），参数2（μ*）
read（1，*）（m10，fqc（i，1），fqc（i，2），i＝icqs1+1，icqs2）
c 接着从’fqc’中读取深层水的区号m10（空读），参数1（T），参数2（μ*），参数3：中-深的越流系数（k’/m’）
read（1，*）（m10，fqc（i，1），fqc（i，2），fqc（i，3），i＝icqs2+1，icqs）
close（1）
c 参数3：中-深的越流系数（k’/m’）同样用于中层水
do 54 i＝1，icqs1
fqc（icqs1+i，3）＝fqc（i，3）
54 continue
c qxia：灌溉抽取地下水平均铺在灌溉区的水层厚度
open（1，file＝’gq’，status＝’old’）
read（1，*）（qxia（i），i＝1，ihv）
close（1）
c gq：灌溉抽出的地下水水层厚度
open（1，file＝’gq’，status＝’old’）
read（1，*）（gq（i），i＝1，ihv）
close（1）
c jbiao：纯引地表水灌溉的节点编号，共有jbn个节点
open（1，file＝’jbiao’，status＝’old’）
read（1，*）（jbiao（i），i＝1，jbn）
close（1）
c’in’：各单元三节点编号文件（必须按小中大顺序排），ia为单元号空读
open（1，file＝’in’，status＝’old’）
read（1，*）（ia，（in（i，j），j＝1，3），i＝1，ms）
close（1）
c igdh：拟合点节点编号
open（1，file＝’igdh’，status＝’old’）
read（1，*）（igdh（i），i＝1，igs）
close（1）
c mqh：各单元的参数区号，ia为单元号空读
open（1，file＝’mqh’，status＝’old’）
read（1，*）（ia，mqh（i），i＝1，ms）
close（1）
c ifdh：开采井节点号
open（1，file＝’ifdh’，status＝’old’）
read（1，*）（ifdh（i），i＝1，mi）
close（1）
c tl（i）：各时段值（天/时段，每月为1个时段）
open（1，file＝’tl1’，status＝’old’）
read（1，*）（tl（i），i＝1，ihv）
close（1）
c 每个时段的大气降水量（m），x（i）/dt则换算为（m/d）
open（1，file＝’x’，status＝’old’）
read（1，*）（x（i），i＝1，ihv）
close（1）
c’zb’：各节点x，y坐标文件（此处为1：2.5万图面的mm数），ia为节点号空读
open（1，file＝’zb’，status＝’old’）
read（1，*）（ia，（zb（i，j），j＝1，2），i＝1，ids）
close（1）
c 把各节点x，y坐标1：2.5万图面的mm数换算为实地的m数
do 888 i＝1，ids
do 888 j＝1，2
zb（i，j）＝zb（i，j）*25
888 continue
c 潜水开采井各时段的开采量（m3/d），空读ia开采井个数，空读ib开采井所在节点编号
open（1，file＝’qq’，status＝’old’）
read（1，*）ia，（（ib，q（i，j），i＝1，mi1），j＝1，ihv）
close（1）
c 中层水开采井第一时段的开采量（m3/d），空读ia开采井个数，空读ib开采井所在节点编号
open（1，file＝’qz1’，status＝’old’）
read（1，*）ia，（ib，q（i，1），i＝mi1+1，mi2）
close（1）
c 中层水开采井第13时段的开采量（m3/d），空读ia开采井个数，空读ib开采井所在节点编号
open（1，file＝’qz2’，status＝’old’）
read（1，*）ia，（ib，q（i，13），i＝mi1+1，mi2）
close（1）
c 深层水开采井第一时段的开采量（m3/d），空读ia开采井个数，空读ib开采井所在节点编号
open（1，file＝’qs1’，status＝’old’）
read（1，*）ia，（ib，q（i，1），i＝mi2+1，mi）
close（1）
c 深层水开采井第13时段的开采量（m3/d），空读ia开采井个数，空读ib开采井所在节点编号
open（1，file＝’qs2’，status＝’old’）
read（1，*）ia，（ib，q（i，13），i＝mi2+1，mi）
close（1）
c 中、深层水开采井第2-12时段的开采量等于第一时段的开采量
do 887 j＝2，12
do 886 i＝mi1+1，mi
q（i，j）＝q（i，1）
886 continue
887 continue
c 中、深层水开采井第14-ihv时段的开采量等于第13时段的开采量
do 884 j＝14，ihv
do 883 i＝mi1+1，mi
q（i，j）＝q（i，13）
883 continue
884 continue
c sh：拟合点各时段实测水位（m），空读ia拟合点所在节点编号
open（3，file＝’sh’，status＝’old’）
read（3，*）（ia，（sh（i，j），j＝1，ihv），i＝1，igs）
close（3）
c qc：潜水初始流场，空读is节点编号
open（1，file＝’qc’，status＝’old’）
read（1，*）（is，h（i），i＝1，ids1）
close（1）
c zc：中层水初始流场，空读is节点编号
open（1，file＝’zc’，status＝’old’）
read（1，*）（is，h（i），i＝ids1+1，ids2）
close（1）
c sc：深层水初始流场，空读is节点编号
open（1，file＝’sc’，status＝’old’）
read（1，*）（is，h（i），i＝ids2+1，ids）
close（1）
c h0nn：分别是潜（h0nn1），中（h0nn2），深（h0nn3）层已知变水头节点各时段的水头，空读ii节点编号
open（1，file＝’h0nn’）
read（1，*）（ii，（h0nn1（i，ikv），ikv＝1，ihv），i＝1，nn1）
read（1，*）（ii，（h0nn2（i，ikv），ikv＝1，ihv），i＝1，nn1）
read（1，*）（ii，（h0nn3（i，ikv），ikv＝1，ihv），i＝1，nn1）
close（1）
c zzz：潜水各节点底板标高（m），空读ii节点编号
open（1，file＝’zzz’）
read（1，*）（ii，zzz（i），i＝1，ids1）
close（1）
c’file1’存放igs个字符串（拟合节点名称占4字符，及观测井的原编号再占4字符）
open（1，file＝’file1’，status＝’old’）
read（1，110）（file1（i），i＝1，igs）
close（1）
c’file2’存放igs个字符串（拟合节点名称占4字符，加.dat后缀再占4字符）
open（1，file＝’file2’，status＝’old’）
read（1，110）（file2（i），i＝1，igs）
close（1）
110 format（10a8）
c 计算开始前先把全部节点的时段末刻水头hed（i），时段初刻水头h0（i）赋初值h（i），以使开始迭代时hed（i），h0（i）不为零
do 1993 i＝1，ids
hed（i）＝h（i）
1993 h0（i）＝h（i）
c 对中、深层承压水导水矩阵工作单元idyh，释水矩阵e，越流矩阵yo，导水矩阵d，赋0
do 25 i＝1+ids1，ids
do 25 j＝1，9
idyh（i，j）＝0
e（i）＝0.0
yo（i）＝0.0
25 d（i，j）＝0.0
c sum2用来累计计算区总面积（用承压水总面积代替），先赋零
sum2＝0
c 对承压水逐个单元计算几何量及导水、释水、越流矩阵
c 对承压水第ip单元的三个节点号依次赋给i，j，k及i1，j1，k1
do 80 ip＝ms1+1，ms
i＝in（ip，1）
j＝in（ip，2）
k＝in（ip，3）
i1＝i
j1＝j
k1＝k
c idyh（i1，9）存放与i1点同单元的所有节点号，最多9个，可以用不完，即i1点的idyh可以有几个idyh（i1，i2）＝0；当计算第ip单元时，i1点的idyh由i1占第1个（i2＝1）位置，j1，k1只能占i2＝2，3，…，9 中的位置；且先占者排前，193行：使j1，k1分两轮到idyh中找位置；当j1找时，195行发现i1的idyh中已有j1，则跳到j1对i1的96循环体头，换为k1对i1循环；当196行没发现已有j1，且196行判断此位不空时返回96循环体头找下一个位置，当碰到第1个空位时，由j1占据，然后跳到j1对i1循环体头，换为k1对i1循环
do 90 iv＝1，3
idyh（i1，1）＝i1
do 94 iu＝j1，k1，k1-j1
do 96 i2＝2，9
if（idyh（i1，i2）.eq.iu）goto 94
if（idyh（i1，i2）.ne.0）goto 96
if（idyh（i1，i2）.eq.0）idyh（i1，i2）＝iu
goto 94
96 continue
94 continue
c 把ip单元的下一个节点j提为i1，k提为j1，i降为k1，然后返回90循环头处理ip单元的下一个节点，循环3次则ip单元中i，j，k 3个点的idyh全部找完
iuu＝i1
i1＝j1
j1＝k1
k1＝iuu
90 continue
c 第ip单元中i，j，k三节点的x坐标赋给xi，xj，xk，y坐标赋给yi，yj，yk
xi＝zb（i，1）
xj＝zb（j，1）
xk＝zb（k，1）
yi＝zb（i，2）
yj＝zb（j，2）
yk＝zb（k，2）
c 第ip单元三角形面积ss
ss＝abs（（xj-xi）*（yk-yi）-（xk-xi）*（yj-yi））*0.5
c 累加各单元面积（这里的sum2累加之后是中、深层两层面积之和）
sum2＝sum2+ss
c 第ip单元的bi～bc（1，1），ci～bc（1，2），bj～b（2，1），cj～b（2，2），bk～b（3，1），ck～b（3，2）
bc（1，1）＝yj-yk
bc（1，2）＝xj-xk
bc（2，1）＝yk-yi
bc（2，2）＝xk-xi
bc（3，1）＝yi-yj
bc（3，2）＝xi-xj
c 第ip单元所在参数区号赋给jv，参数1（T）赋给txy，参数2（μ）赋给ts，参数3（k’/m’）赋给rmk
jv＝mqh（ip）
txy＝fqc（jv，1）
ts＝fqc（jv，2）
rmk＝fqc（jv，3）
c 第ip单元三节点i，j，k的释水矩阵元素c（ii，ii）及c（ii，jj）……式中没包括1/（2Δt）时，随着ip＝ms1+1，ms的单元循环而对有关单元求其和
e（i）＝-ts*ss/3.0+e（i）
e（j）＝-ts*ss/3.0+e（j）
e（k）＝-ts*ss/3.0+e（k）
c越流矩阵元素，第ip单元1m水头差时的越流量（m3/d/m）均分给三节点i，j，k，随着ip＝ms1+1，ms的单元循环而对有关单元求其和
yo（i）＝yo（i）+rmk*ss/3.0
yo（j）＝yo（j）+rmk*ss/3.0
yo（k）＝yo（k）+rmk*ss/3.0
c 计算ip单元各节点的导水矩阵元素d（i，j），（i＝1，2，3）（j＝1，2，3）
do 100 iu＝1，3
do 104 iuu＝1，3
c 当iu＝1时，即ip单元的i节点，分别计算iuu＝1，2，3，即ip单元的i，j，k节点对i节点的ai，ai即公式**的没求和部分
i＝in（ip，iu）
j＝in（ip，iuu）
ai＝txy*（bc（iu，1）*bc（iuu，1）+bc（iu，2）*bc（iuu，2））/ss/4.0
c 在ip单元的三个节点中，排除第3点，只让第1，2点（i，j点）的ai加入i点的导水矩阵元素d（i，j）中，第243行的j可分别轮到i，j，k三点，但第245句的1～9个中，仅有i，j点在i点的idyh中，此句排除了第3点加入i点的导水矩阵元素d（i，j）中随着ip＝ms1+1，ms的循环，到251句时承压水导水矩阵已完全形成
do 106 k＝1，9
if（j.eq.idyh（i，k））d（i，k）＝ai+d（i，k）
106 continue
104 continue
100 continue
80 continue
c 至此，中、深层几何量计算完，以下开始时段循环；时段循环中把潜水也加进来
c 在时段循环中，潜水的导水矩阵，释水矩阵，越流矩阵与含水层厚度有关，所以在时段循环中完成，另外，承压水与时段有关的部分也在时段循环中完成
do 280 ikv＝1，ihv
write（*，*）’ikv＝’，ikv
c 该时段步长赋给dt
dt＝tl（ikv）
c 以下开始计算潜水末刻流场
c 潜水已知水头节点的已知水头赋给本时段初、末刻
do 32 i＝1，nn1
hed（n1+i）＝h0nn1（i，ikv）
h0（n1+i）＝h0nn1（i，ikv）
32 continue
c 先计算潜水几何量，因为潜水的含水层厚度M 随时段而变，所以放在时段循环内
do 33 i＝1，ids1
h0（i）＝hed（i）
c 本时段潜水各节点的含水层厚度ep0（i）用时段初刻水位＝hed（i）减去底板标高zzz（i）求得
ep0（i）＝hed（i）-zzz（i）
c本时段潜水各节点的灌溉回归水量加降水回归水量w（i，ikv）先赋零，以备求和之用
w（i，ikv）＝0.0
33 continue
c 该时段潜水各单元的含水层厚度用三节点之均值
do 299 i＝1，ms1
299 ep1（i）＝（ep0（in（i，1））+ep0（in（i，2））+ep0（in（i，3）））/3.
do 2050 i＝1，ids1
c 同前，对潜水导水矩阵工作单元idyh，释水矩阵e，越流矩阵yo，导水矩阵d，先赋0
do 2050 j＝1，9
idyh（i，j）＝0
e（i）＝0.0
yo（i）＝0.0
2050 d（i，j）＝0.0
do 8015 ip＝1，ms1
i＝in（ip，1）
j＝in（ip，2）
k＝in（ip，3）
i1＝i
j1＝j
k1＝k
do 9015 iv＝1，3
idyh（i1，1）＝i1
do 9415 iu＝j1，k1，k1-j1
do 9615 i2＝2，9
if（idyh（i1，i2）.eq.iu）goto 9415
if（idyh（i1，i2）.ne.0）goto 9615
if（idyh（i1，i2）.eq.0）idyh（i1，i2）＝iu
goto 9415
9615 continue
9415 continue
iuu＝i1
i1＝j1
j1＝k1
k1＝iuu
9015 continue
xi＝zb（i，1）
xj＝zb（j，1）
xk＝zb（k，1）
yi＝zb（i，2）
yj＝zb（j，2）
yk＝zb（k，2）
ss＝abs（（xj-xi）*（yk-yi）-（xk-xi）*（yj-yi））*0.5
jvv＝mqh（ip）
dyk＝fqc（jvv，4）
dyk1＝gh（jvv）
do 776 ig＝1，3
ngh＝in（ip，ig）
do 777 jb＝1，jbn
if（ngh.eq.jbiao（jb））goto 778
777 continue
w（ngh，ikv）＝w（ngh，ikv）-qxia（ikv）*ss/3
778 w（ngh，ikv）＝w（ngh，ikv）+x（ikv）/dt*dyk*ss/3+gq（ikv）*dyk1*ss/3
776 continue
bc（1，1）＝yj-yk
bc（1，2）＝xj-xk
bc（2，1）＝yk-yi
bc（2，2）＝xk-xi
bc（3，1）＝yi-yj
bc（3，2）＝xi-xj
jv＝mqh（ip）
txy＝fqc（jv，1）*ep1（ip）
ts＝fqc（jv，2）
rmk＝fqc（jv，3）
e（i）＝-ts*ss/3.0+e（i）
e（j）＝-ts*ss/3.0+e（j）
e（k）＝-ts*ss/3.0+e（k）
yo（i）＝yo（i）+rmk*ss/3.0
yo（j）＝yo（j）+rmk*ss/3.0
yo（k）＝yo（k）+rmk*ss/3.0
do 1001 iu＝1，3
do 1041 iuu＝1，3
i＝in（ip，iu）
j＝in（ip，iuu）
ai＝txy*（bc（iu，1）*bc（iuu，1）+bc（iu，2）*bc（iuu，2））/ss/4.0
do 1061 k＝1，9
if（j.eq.idyh（i，k））d（i，k）＝ai+d（i，k）
1061 continue
1041 continue
1001 continue
8015 continue
c 至此，该时段潜水几何量计算完
c 时段末水头迭代计数器ikv2，最大误差记录amax
ikv2＝0
9881 amax＝0.0
ikv2＝ikv2+1
write（*，34）ikv，ikv2
34 format（3x，i4，20x，’ikv2＝’，i4）
c 开采井点录用计数器iq
iq＝1
c 对潜水n1个未知水头节点逐点计算该时段末刻水头
do 4002 i＝1，n1
c i节点常数项res：减该时段i点降水灌溉回归量，减i点该时段来自承压水的越流量（越流量计算采用时段初承压水水头，时段末潜水水位），减i点该时段储存量的增加量
res＝-w（i，ikv）-（hed（i+ids1）-hed（i））*yo（i）-e（i）*（hed（i）-h0（i））/dt
if（i.eq.ifdh（iq））then
c 如果i节点恰是开采节点，则i节点常数项res中再加上该时段i点的开采量，然后开采井点录用计数器iq加1，开采井点只能被1个节点录用
res＝res+q（iq，ikv）
iq＝iq+1
endif
c 把与i节点共单元的k节点的导水矩阵元素（k＝1，2，…，最多可9，当k＝1时即i点自己）依次加进i节点常数项res中，k节点的导水矩阵元素采用上一轮迭代出的hed（k）计算
do 3002 j＝1，9
k＝idyh（i，j）
if（k.eq.0）goto 3002
res＝res+d（i，j）*hed（k）
3002 continue
c 把常数项除以：（i节点对i节点自身导水矩阵元素d（i，1）的负数，加，i节点储水阵元素e（i）/dt）
res＝res/（-d（i，1）+e（i）/dt）
c 把i节点的res乘以亚松弛系数zz1，加给上一轮迭代出的hed（i）中，作为这一轮迭代出的hed（i）：
hed（i）＝res*zz1+hed（i）
c 把1-n1个节点中最大的res挑出，并把其点号记到imax1中：
if（abs（res）.le.amax）goto 4002
amax＝abs（res）
imax1＝i
c 循环返回，继续下一个节点：
4002 continue
988 amax＝0.0
232 continue
write（*，*）’amax1＝’，amax，imax1
if（amax.gt.9.999999e-03）goto 9881
c 至此，潜水本时段末刻流场计算完，以下开始中层水本时段末刻流场的计算：
do 232 i＝1，nn1
hed（n2+i）＝h0nn2（i，ikv）
h0（n2+i）＝h0nn2（i，ikv）
iqq＝mi1+1
do 400 i＝ids1+1，n2
c 其中的越流流入量采用（潜水本时段末水位hed（i-ids1）减承压水本时段末水位hed（i））之差hedd计算
hedd＝hed（i-ids1）-hed（i）
heddd＝hed（i）-hed（i+ids1）
res＝-hedd*yo（i）+heddd*yo（i+ids1）-e（i）*（hed（i）-h0（i））/dt
50 if（i.eq.ifdh（iqq））then
res＝res+q（iqq，ikv）
iqq＝iqq+1
endif
do 300 j＝1，9
k＝idyh（i，j）
if（k.eq.0）goto 300
res＝res+d（i，j）*hed（k）
300 continue
res＝res/（-d（i，1）+e（i）/dt）
hed（i）＝res*zz2+hed（i）
if（abs（res）.le.amax）goto 400
amax＝abs（res）
imax2＝i
400 continue
write（*，*）’amax2＝’，amax，imax2
if（amax.gt.9.999999e-03）goto 988
c 至此，中层水全部节点末刻水头迭代完一轮，节点误差最大者如果 ＞0.01m，则988 句进行下一轮迭代，如果≤0.01m，则迭代结束，中层水本时段末刻流场计算完，以下开始深层水本时段末刻流场的计算：
989 amax＝0
iqqq＝mi2+1
do 332 i＝1，nn1
hed（n+i）＝h0nn3（i，ikv）
h0（n+i）＝h0nn3（i，ikv）
332 continue
do 401 i＝ids2+1，n
hedd＝hed（i-ids1）-hed（i）
res＝-hedd*yo（i）-e（i）*（hed（i）-h0（i））/dt
5011 if（i.eq.ifdh（iqqq））then
res＝res+q（iqqq，ikv）
iqqq＝iqqq+1
endif
do 301 j＝1，9
k＝idyh（i，j）
if（k.eq.0）goto 301
res＝res+d（i，j）*hed（k）
301 continue
res＝res/（-d（i，1）+e（i）/dt）
hed（i）＝res*zz2+hed（i）
if（abs（res）.le.amax）goto 401
amax＝abs（res）
imax3＝i
401 continue
write（*，*）’amax3＝’，amax，imax3
if（amax.gt.9.999999e-03）goto 989
c至此，深层水全部节点本时段末刻水头迭代完一轮，节点误差最大者如果＞0.01m，则返回989句进行下一轮迭代，如果≤0.01m，则迭代结束，深层水本时段末刻流场计算完
c以下开始把本时段末刻水位潜，中，深层水全部节点本时段末刻水头hed赋给下时段初刻水位h0，三层拟合点本时段末刻水位hed存入fh
189 do 281 i＝1，ids
h0（i）＝hed（i）
281 continue
do 190 j＝1，igs
mp＝igdh（j）
fh（j，ikv）＝hed（mp）
190 continue
280 continue
c 至此，三层全包括的时段循环完
c 屏幕输出计算区总面积sum2/2，（m2）；sum2是中、深层两层的累加和：
write（*，*）’sum＝’，sum2/2，’m2’
do 1926 i＝1，ids
1926 sss（i）＝hed（i）-h（i）
c’ok1’输出潜水的最终流场：节点号i、坐标zb、最终流场hed、该节点0-ihv时间（最终流场比初始流场）的水位总变幅：
open（1，file＝’ok1’）
write（1，1927）（i，zb（i，1），zb（i，2），hed（i），sss（i），i＝1，ids1）
close（1）
c’ok2’输出中层水的最终流场：节点号i、坐标zb、最终流场hed、该节点0-ihv时间（最终流场比初始流场）的水位总变幅：
open（1，file＝’ok2’）
write（1，1927）（i，zb（i，1），zb（i，2），hed（i），sss（i），i＝ids1+1，ids2）
close（1）
c’ok3’输出深层水的最终流场：节点号i、坐标zb、最终流场hed、该节点0-ihv时间（最终流场比初始流场）的水位总变幅：
open（1，file＝’ok3’）
write（1，1927）（i，zb（i，1），zb（i，2），hed（i），sss（i），i＝ids2+1，ids）
close（1）
1927 format（i4，2f7.0，2f7.2）
c 计算出拟合点各时段误差（计算水位-实测水位）s（i，iuv）
do 8 iuv＝1，ihv
do 2003 i＝1，igs
s（i，iuv）＝fh（i，iuv）-sh（i，iuv）
2003 continue
8 continue
c 输出拟合点各时段误差
open（1，file＝’gan.dat’）
do 3511 i＝1，igs，5
write（1，2993）file1（i），file1（i+1），file1（i+2），file1（i+3），file1（i+4）
2993 format（5（3x，a8，3x））
do 3778 iv＝1，ihv
f0＝fh（i，iv）
f1＝fh（i+1，iv）
f2＝fh（i+2，iv）
f3＝fh（i+3，iv）
f4＝fh（i+4，iv）
s0＝s（i，iv）
s1＝s（i+1，iv）
s2＝s（i+2，iv）
s3＝s（i+3，iv）
s4＝s（i+4，iv）
write（1，3556）iv，f0，s0，f1，s1，f2，s2，f3，s3，f4，s4
3778 continue
3556 format（i3，10f7.2）
3511 continue
close（1）
c 输出拟合节点计算的历时水位fh（i，iuv）、实测的历时水位sh（i，iuv）、拟合误差s（i，iuv），拟合节点名称file1（i）do 2013 i＝1，igs
tt＝0
filee＝file2（i）
open（1，file＝filee）
do 2018 iuv＝1，ihv-1
tt＝tt+tl（iuv）
write（1，2020）tt，fh（i，iuv），sh（i，iuv），s（i，iuv）
2018 continue
write（1，2030）tt，fh（i，iuv），sh（i，iuv），s（i，iuv），file1（i）
close（1）
2013 continue
2020 format（1x，f10.3，3f15.3，2x）
2030 format（1x，f10.3，3f15.3，2x，’＂’，a8，’＂’）
stop
end
对于你说的这个问题，那个3384分指的是你唱这首歌的评分高低，也就是你唱的好不好听。
《sss988》的内容今天先分享到这里了，因为炒股有风险的，大家不要盲目去投资，需要在炒股期间不断学习与积累经验。
评论
游客
回复看帖不回帖都是耍流氓！http://www.dnf70.com/3166.html
国产无套视频在线观看香蕉
回复网站做得不错http://www.zhijian.me/