将行列型年月日数据转为单列格式并可提取逐年最大、最小值的自编程序（带实例）

hillside

      现上传本人自编将行列型年月日多年原始数据转为单列格式并可另外提取逐年最大、最小值的程序。对于高手而言似乎不值一提，对于新手而言尚有入门价值。编写目的是为制备其他软件所需的数据格式（主要是美国制）。
      因此，本帖主要针对初学人群。
                                                                      转换、提取程序（fortran）使用效果说明

        编译方式同于常规fortran编译。
       1  数据预处理        转换前原始数据格式如（行为年月，列为日数。为整齐起见，列为31，不足31列根据物理量含义可用0或其他符号代替）：

stnm	stcd	yr	mo	q1	q2	q3	q4	q5	q6	q7	q8	q9	q10	q11	q12	q13	q14	q15	q16	q17	q18	q19	q20	q21	q22	q23	q24	q25	q26	q27	q28	q29	q30	q31
站名	12302250	1965	1	380	380	380	330	393	412	419	426	439	439	419	406	400	386	380	355	342	330	324	317	317	305	300	294	306	311	317	305	330	348	361
站名	12302250	1965	2	361	374	386	393	439	549	651	762	944	1050	1100	1050	953	869	852	877	902	936	987	970	1030	1050	978	844	737	659	651	787	0	0	0
站名	12302250	1965	3	962	894	860	787	659	592	556	549	528	480	459	486	542	549	542	500	473	507	480	446	459	621	996	1540	1950	1900	1790	1760	1980	2480	3100
站名	12302250	1965	4	3720	4540	5340	6030	6300	5700	4720	3750	2900	2230	1730	1360	1120	996	1120	1760	2210	2190	2160	2190	2150	1970	1760	1540	1350	1180	1060	953	852	811	0
站名	12302250	1965	5	852	1110	1520	1990	2060	1860	1810	1990	2730	3020	2630	2460	2580	2390	2000	1680	1450	1250	1180	1550	2460	2950	2740	2640	4160	5560	6190	6420	6300	7190	8050
站名	12302250	1965	6	8450	8260	7560	6840	5820	4840	4160	4540	5520	6350	6600	6390	5990	5840	5450	4960	4590	4700	5620	7080	9230	10800	11000	10300	8790	7000	4790	3740	3200	2850	0

......

           2  预处理后的数据（用作fortran输入数据）

380	380	380	330	393	412	419	426	439
361	374	386	393	439	549	651	762	944
962	894	860	787	659	592	556	549	528
3720	4540	5340	6030	6300	5700	4720	3750	2900
852	1110	1520	1990	2060	1860	1810	1990	2730
8450	8260	7560	6840	5820	4840	4160	4540	5520
2520	2200	2000	1650	1410	1360	1330	1400	1320
2490	1970	1640	1370	1150	1040	952	902	903
3280	2830	2350	2090	1960	1860	1820	1800	1740
647	620	592	565	550	530	575	530	470
398	407	416	393	380	380	367	355	336
406	380	367	355	342	342	330	336	317
380	380	367	355	408	507	685	734	759
340	347	340	340	340	401	436	422	422
1070	1360	1540	1640	2040	2450	2450	2380	2340

......

        3.  转换后结果格式为：
1/ 1/1955,  380.000
1/ 2/1955,  380.000
1/ 3/1955,  380.000
1/ 4/1955,  330.000
1/ 5/1955,  393.000
1/ 6/1955,  412.000
1/ 7/1955,  419.000
1/ 8/1955,  426.000
1/ 9/1955,  439.000
1/10/1955,  439.000
1/11/1955,  419.000
1/12/1955,  406.000
1/13/1955,  400.000
1/14/1955,  386.000
1/15/1955,  380.000
1/16/1955,  355.000
1/17/1955,  342.000
1/18/1955,  330.000
1/19/1955,  324.000
1/20/1955,  317.000
1/21/1955,  317.000
1/22/1955,  305.000
1/23/1955,  300.000
1/24/1955,  294.000
1/25/1955,  306.000
1/26/1955,  311.000
......   
                                                                                               预处理说明                                                                                                                                                                              
              本帖原始数据为水文流量，缺测值用0表示。不同场合， 缺测值的处理有不同的方法。
        因形式上的2月29（非闰年）、30、31日及4月31日不存在，而程序中为便利起见，取整齐的数据排列(行为年，列为31)。以流量而言，在较大流域不会小于0，故将无值设为0。如计算其他缺测不宜用0代替的值，需要对程序略作修改。
                                                                                    将无值替换为0的处理方法
      在EXCEL中，使用替换功能，先选中空白格，复制到待替换位置，再用0值替换。

附件内容说明：
1 将行列型数据（年月日值）转为单列逐日数据(用于IHA输入等).F90
2 从行列型逐日数据中提取年最小、最大值(用于MK输入等).F90
3 实例数据：hydro_station
4 原始数据已作技术处理，不代表原来的含义

fanlaoda999

有列为日行为所 测 值 的 程 序 没