中国气象辐射资料台站信息表
基于Angstrom-Prescott模型和空间插值的山东省月太阳辐射再估算
收稿日期:２０２０￣０２￣１０
基金项目:国家自然科学基金(４１４０６０１０)
作者简介:王仁政(１９９５—) ꎬ男ꎬ硕士研究生ꎬ研究方向为太阳辐射ꎮ Ｅ￣ｍａｉｌ:２８１８６２２５７７＠ ｑｑ.ｃｏｍ
∗通信作者ꎬ宫响(１９７７—) ꎬ博士ꎬ副教授ꎬ研究方向为太阳能资源估算及神经网络应用ꎮ Ｔｅｌ:１３２１００８６６５５ꎬＥ￣ｍａｉｌ:ｇｏｎｇｘｉａｎｇ＠ ｑｕｓｔ.ｅｄｕ.ｃｎ
在获得太阳能资源分布特征的前提下才能更好地实现对太阳能资源的高效科学利用 [１￣２] ꎮ 但在全球范围内ꎬ
对太阳辐射进行常规观测的台站密度不大ꎬ无法满足分布规律研究和实际应用的需要ꎮ 因此ꎬ如何对某一地
域的太阳辐射进行准确的评估成为了重中之重ꎮ
１９１９ 年ꎬＫｉｍｂａｌｌ [３] 提 出 了 关 于 太 阳 总 辐 射 与 日 照 百 分 率 之 间 存 在 良 好 相 关 性 的 分 析ꎮ １９２４ 年ꎬ
９４
山 东 科 学
２０２０ 年
ｔｈｅ ｔｈｒｅｅ ｎａｔｉｏｎａｌ ｓｔａｔｉｏｎｓꎬ ｔｈｅ Ａｎｇｓｔｒｏｍ￣Ｐｒｅｓｃｏｔｔ ｍｏｄｅｌ ｃｏｅｆｆｉｃｉｅｎｔｓ ( Ａ￣Ｐ ｃｏｅｆｆｉｃｉｅｎｔｓ) ｉｎ ｄｉｆｆｅｒｅｎｔ ｍｏｎｔｈｓ ａｒｅ ｏｂｔａｉｎｅｄ
ｕｓｉｎｇ ｔｈｅ ｌｅａｓｔ ｓｑｕａｒｅ ｆｉｔｔｉｎｇ ｍｅｔｈｏｄ ａｎｄ ｖｅｒｉｆｉｅｄ ｕｓｉｎｇ ｔｈｅ ｏｂｓｅｒｖａｔｉｏｎ ｄａｔａ ｉｎ ２０１６. Ｔｈｅｎꎬ ｔｈｅ ｓｏｌａｒ ｒａｄｉａｔｉｏｎ ｉｎ Ｈｅｚｅꎬ
Ａｎｇｓｔｒｏｍ [４] 首次给出了太阳辐射的气候学计算公式ꎮ １９４０ 年ꎬＰｒｅｓｃｏｔｔ [５] 利用天文总辐射代替 Ａｎｇｓｔｒｏｍ 提出
基于GMI辐射计资料分析南海风速分布特征
195数据库技术Database Technology电子技术与软件工程Electronic Technology & Software Engineering1 引言海面风场是海洋动力学的基本参数，是海洋和大气动量交换、热量交换的重要承载体，直接影响大洋表层物理量的变化。
传统海洋风场主要是靠台站、浮标、船舶等传统方式进行常规观测，由于观测频次和探测范围无法实现对广阔海域的大范围、连续观测，因而利用这些传统观测资料绘画模拟整个海面的风场具有一定的难度[1-4]。
近年来，空间遥感技术的发展，使得高时空分辨率、覆盖全球、连续观测的卫星遥感资料的应用成为可能，而这些遥感资料的大量应用，有利于推动海洋气象事业的发展。
卫星遥感监测方法的逐步成熟，为海面风的监测和研究提供了新的资料来源。
目前可以观测海面风的传感器有微波散射计、微波高度计和微波辐射计[5-7]。
在这些传感器中微波散射计可以测风速、风向，高度计和辐射计可以测风速要素，但在南海研究海风现状中前两者研究较多，在微波辐射计这块研究较少，正因为如此本文将利用辐射计资料研究，弥补南海多种资料分析海面风的缺失。
目前，国内外主要进行了散射计反演风场资料与浮标和船舶观测资料、岛屿站资料、海洋站资料等的检验。
研究结果表明散射计矢量风风速在南海区域具有可信性。
刘春霞[8]等主要分析了QuikSCAT 散射计矢量风资料的统计特征，得出QuikSCAT 矢量风在南海具有可信性，且南海主要盛行两种风，即冬季东北风和夏季西南风，东北风最大中心在巴士海峡、台湾海峡。
利用GMI 辐射计资料和最新的观测数据定量分析误差，而这些误差特征也是GMI 风场数据同化、融合所必须的。
本文重点分析了GMI 辐射计风资料在南海海域的可信性和统计特征，并在此基础上利用辐射计月平均风场资料分析南海月平均风场特征。
综上所述，学者针对南海近海面风场做了较为系统的研究。
但以往的研究时间跨度较短，数据资料空间分辨率较低，且仅风速标量值为统计对象，对于风场经纬向分量长期趋势研究相对较少。
太阳辐射经验系数插值方法的比较
太阳辐射经验系数插值方法的比较买苗;火焰;俞亚勋【摘要】利用江苏省南京、吕泗、淮安3个日射站以及周边省市共10个日射站逐月日照百分率资料和太阳总辐射资料,利用最小二乘法拟合经验系数a、b,并利用除南京、吕泗2个日射站以外的8个日射站的经验系数a、b,采用多种插值方法计算江苏省70个站的经验系数a、b,并据此求出70个站点逐月太阳总辐射值并分析其分布特征。
结果表明：采用反距离权重插值法得到的经验系数a、b,误差最小。
结论可为求解江苏省各地太阳总辐射提供科学参考。
%Based on the monthly sunshine percentage and solar radiation data from 3 solar radiation stations（Nanjing,Lvsi and Huai′an stations） in Jiangsu province and from 7 other solar radiation stations in surrounding provinces,the empirical coefficient a and b were calculated by a least-squares method.The empirical coefficients a and b for about 70 stations in the whole province were obtained by various interpolation schemes.Then the characters of the monthly total solar radiation and its distribution were obtained.The results indicate that the empirical coefficients a and b determined by the inverse distance weighted interpolation method have the least errors among various interpolation schemes.The conclusion could provide scientific references for calculation of solar radiation of Jiangsu province.【期刊名称】《气象与环境学报》【年(卷),期】2011(027)005【总页数】4页(P42-45)【关键词】太阳总辐射;经验系数;最小二乘法;日射站【作者】买苗;火焰;俞亚勋【作者单位】江苏省气候中心,江苏南京210008;江苏省气象科学研究所,江苏南京210008;江苏省气候中心,江苏南京210008【正文语种】中文【中图分类】P422.1引言太阳能是地球上最基本、最重要的能源，影响到地球上所有的物理、生物和化学过程，因此，地表辐射研究在国内外一些重大的研究计划中一直备受重视。
气象数据文件和记录簿表格式
数据文件和记录簿表格式全部试题一、填空：1、自动气象站采集数据文件是指由数据采集器处理后，通过通讯接口直接存储到计算机硬盘中的数据文件。
它是自动气象站与地面气象测报业务软件的接口数据文件2、自动气象站采集数据文件由正点地面气象要素数据文件分钟地面气象要素数据文件、实时地面气象要素数据文件、大风数据文件、正点气象辐射数据文件、实时气象辐射数据文件、正点气象辐射补充数据文件等文件组成3、Z文件（ZIIiiiMM.YYY）中，Z为指示符；IIiii为区站号；MM 为月份，不足两位时，前面补“0”；YYY为年份的后3位。
4、Z文件为随机文件，每月一个，记录采用定长类型，每一条记录218个字节，记录尾用回车换行结束，ASCII字符存盘，每个要素值高位不足补空格。
5、Z文件第一次生成时应进行初始化，初始化的过程是：首先检测Z 文件是否存在，如无当月Z文件，则生成该文件，将全月逐日逐时各要素的位置一律存入相应字长的“-”字符。
6、Z文件按北京时计时，以北京时的00分数据作为正点定时数据。
7、Z文件的第1条记录为本站当月基本参数，每项参数长为5个字节。
8、Z文件中的自动站型号标识：milos系列自动站存入“3”、I型自动站存入“1”、Ⅱ型自动站存入“2”。
9、Z文件中的干湿表通风系数Ai值：扩大107倍后存入。
10、Z文件中每一时次为一条记录，每日24条记录。
记录号的计算方法：N=D 24+T-19 。
上式中，N：记录号；D：北京时日期（月末一天21～23时的日期D 取0）；T：北京时。
11、Z文件中第1条后的每一条记录，存54个要素的正点值，以ASCII 字符写入，除能见度和最小能见度为5个字节外，其他每个要素长度为4字节，最后两位为回车换行符。
12、Z文件中的气压，当气压值≥1000.0hPa时，先减去1000.0，再乘以10后存入。
13、Z文件中若要素缺测，除有特殊规定外，则均应按约定的字长，每个字节位均存入一个“/”字符。
各种遥感数据,地理信息数据网站地址汇总
各种遥感数据，地理信息数据共享网站（至少一百）∙Online Global Satellite Image and Atlas:h ttp://library.gmu.e du/resources/sci/Geog579.htm ∙可以下载Aster,Quic kBird,IKonos,OrbView,LandSat, SRTM,MODIS数据的网站:∙http://glcf.umi /data/∙/data/ NASA数据站...推荐∙National Geophysical Data Cent er GSHHS - A Global Self-consistent, Hierarchical, High-resolution Shoreline Database v1.3 (13 Oct 2005)∙National Geophy sical Data Center Coastline Extractor∙U.S. Geologica l Survey GTOPO30 is a global digital elevation model (DEM) with a horizontal grid spacing of 30 arc seconds (approximately1 kilometer)HYDRO1k Elevation Derivative Database | Australasia∙National Geophy sical Data Center GLOBE Project∙Scripps Institu tion of Oceanography Global Topography | FTP site∙/ *国家MODIS数据共享平台/* 美国国家海洋和大气局/ *GIS公园-很多GIS知识和信息*可以下很多的免费数据http://www.vgt.vito.be/*用法国SPOT卫星生产的植被数据分辨率1Km/*国家基础地理信息系统网站/ *中国卫星遥感数据服务网-可订购多种遥感数据http://www.naturalresources.cs/ *中国自然资源数据库/ *地理信息系统论坛-学习交流好去处/*提供很多有用的GIS扩展模块和工具/ *中国气象科学数据共享服务网以下部分转贴自数据的天空综合资料、降水、SST、地面覆盖资料、风场/OLR/指数资料综合资料库：/cgi-bin/forum/topic.cgi?forum=2&topic=484&start=0Noaa资料库：NCEP资料介绍：/cgi-bin/forum/topic.cgi?forum=13&topic=3/ncep_data/欧洲气象中心资料(grib和NC格式的)：http://www.ecmwf.int/http://data.ecmwf.int/data/d/era40_daily/Levitus资料：/SOURCES/.LEVITUS94/.MONTHLY/Ucar资料/cas/guide/Atmos/Surface/data.htmlNASA资料：ftp:///seasurfaceheight/以前某天全国的天气情况/feature/hi301100.shtml1度×1度资料/datasets/ds083.2/ARGO资料/web/NCEP系统资料：NCEP real-analyses and forecasts /data/NCEP/NCAR REANALYSIS /pub/reanalysis/NCEP Eta /mmb/research/meso.products.htmlNCEP AVN /modelinfo/ftp:///pub/data/nccf/com/gfs/prod/grib格式的ncep资料下载/cgi-bin/forum/topic.cgi?forum=13&topic=561netCDF format NNRP1: 6 hourly, 2.5 degrees, from 1948 to presentftp:///pub/Datasets/ncep.reanalysis/NNRP2: 6 hourly, 2.5 degrees, from 1979 to 2002ftp:///pub/Datasets/ncep.reanalysis2/降水资料CMAP资料:/data_sets/cmap_precip/NCEP资料: /cgi-bin/forum/topic.cgi?forum=2&;topic=484地面资料: /cgi-bin/forum/topic.cgi?forum=2&;topic=445细分辨率的全球径流指标场：/cgi-bin/forum/topic.cgi?forum=2&topic=365全球土壤资料：/CAMPAIGN_DOCS/FTP_SITE/INT_DIS/readmes全国160个站的降水资料/cgi-bin/forum/topic.cgi?forum=13&topic=426&start=0/Website/index.php?ChannelID=43&WCHID=5风场资料：/cdc/data.ncep.reanalysis.derived.html#surface_gauss NCEP-QSCAT混合风场资料/海洋再分析资料：/SOURCES/.UMD/.Carton/.goa/.beta7/海表高度： ftp:///seasurfaceheight/ftp:///pub/sea_surface_height/topex_poseidon/mgdrb/data/MGB_423/MGB423.129.Z海面风场：/海浪要素/cgi-bin/forum/topic.cgi?forum=2&topic=3736中国常规气象观测资料：国气象局气象中心资料室就能拿到，只要一张资料50元的光盘刻录费而已。
AMSU-A全空辐射率资料同化对台风“天鸽”的预报影响研究
AMSU-A全空辐射率资料同化对台风“天鸽”的预报影响研究张涛; 姜立鹏; 师春香; 周自江【期刊名称】《《大气科学学报》》【年(卷),期】2019(042)005【总页数】10页(P705-714)【关键词】资料同化; 全空辐射率; AMSU-A; 台风; 天鸽【作者】张涛; 姜立鹏; 师春香; 周自江【作者单位】国家气象信息中心北京100081【正文语种】中文近年来,随着卫星遥感技术、同化方法、辐射传输模式、数值模式等方面的不断发展,不同平台不同探测方式的卫星资料同化得到了越来越深入的研究(李刚等,2016a;钟亦鸣等,2016;夏宇等,2018)。
在全球大多数的数值天气预报系统中,卫星资料占据了所有同化资料的90%～95%(Bauer et al.,2010;钟亦鸣等,2016)。
然而由于云和降水的影响,超过75%的卫星资料在质量控制时被剔除掉(官莉和陆文婧,2016;Wang et al.,2018a)。
但这并不意味着这些卫星辐射率资料不重要,而是因为它们在数值模式中的使用存在着一定的难度(李刚等,2016b)。
对于卫星观测资料来说,受到云和降水影响的资料信息很难从主要的观测量(如温度和湿度)中剥离出来(Bauer et al.,2011)。
云和降水的发生发展过程往往显示出一些天气过程的重要动力及热力特征,并且模式在预报过程中也对初始场中云水信息较为敏感。
开展对受云和降水影响的卫星资料的同化研究,实现卫星资料的全空同化,对提高资料同化水平,改善数值天气预报效果尤为重要。
2009年3月以来,SSM/I和地球观测系统高级微波扫描辐射计(AMSR-E)的全空微波辐射数据已经被直接同化到ECMWF的业务4DVar资料同化系统(Bauer et al.,2010;Geer et al.,2010)。
随着通用辐射传输模式(CRTM)和预报模式的改进,NCEP在全空微波辐射同化方面取得很大进展,AMSU-A全空辐射率同化自2016年5月12日开始业务化运行(Zhu et al.,2016)。
中国气象局第7号令《气象探测环境和设施保护办法》
中国气象局第7号令《气象探测环境和设施保护办法》第一条为了保护气象探测环境和设施，保证气象探测工作的顺利进行，确保获取的气象探测信息具有代表性、准确性、比较性，提高气候变化的监测能力、气象预报准确率和气象服务水平，为国民经济和人民生活提供可靠保障，根据《中华人民共和国气象法》，制定本办法。
第二条本办法适用于中华人民共和国领域和中华人民共和国管辖的其他海域内气象探测环境和设施的保护。
第三条本办法所称气象探测环境，是指为避开各种干扰保证气象探测设施准确获得气象探测信息所必需的最小距离构成的环境空间。
本办法所称气象探测设施，是指用于各类气象探测的场地、仪器、设备及其附属设施。
第四条国务院气象主管机构负责管理全国气象探测环境和设施的保护工作。
地方各级气象主管机构在上级气象主管机构和同级人民政府的领导下，负责管理本行政区域内气象探测环境和设施的保护工作。
设有气象台站的国务院其他有关部门和省、自治区、直辖市人民政府其他有关部门按照职责，做好本部门气象台站的探测环境和设施的保护工作，并接受同级气象主管机构的指导、监督和行业管理。
其他有关部门按照职责，配合气象主管机构做好气象探测环境和设施的保护工作。
第五条国家依法保护气象探测环境和设施。
任何组织和个人都有保护气象探测环境和设施的义务，有权检举侵占、损毁和擅自移动气象探测设施和破坏气象探测环境的行为。
第六条各级人民政府及有关部门应当加强对气象探测环境和设施保护的宣传教育，树立全民保护气象探测环境和设施的意识。
对在保护气象探测环境和设施工作中作出贡献的单位和个人，给予奖励。
第七条本办法保护以下气象探测环境和设施：(一)国家基准气候站、国家基本气象站、国家一般气象站、自动气象站、太阳辐射观测站、酸雨监测站、生态气象监测站（含农业气象站）的探测环境和设施；(二)高空气象探测站（包括风廓线仪、声雷达、激光雷达等）的探测环境和设施；(三)天气雷达站的探测环境和设施；(四)气象卫星地面接收站（含静止气象卫星地面接收站、极轨气象卫星地面接收站）、卫星测控站、卫星测距站的探测环境和设施；(五)大气本底台站、沙尘暴监测站、污染气象监测站等环境气象监测站的探测环境和设施；(六)遥感卫星辐射校正场的探测环境和设施；(七)闪电探测站的探测环境和设施；(八)GPS气象探测站外场环境；(九)气象专用频道、频率、线路、网络及相应的设施；(十)其他需要保护的气象探测环境和设施。
6.国内气象辐射小时观测数据BUFR编码格式(V1.0)
国内气象辐射小时观测数据BUFR编码格式（V1.0）
（V1.0版）
1 范围
本格式规定了国内气象辐射站小时观测数据及其质量控制码的编码格式、编报规则和代码，适用于国内气象辐射站小时观测数据的编报和传输。
2 格式
编码数据由指示段、标识段、数据描述段、数据段和结束段构成。
2.1 0段——指示段
指示段包括BUFR编码数据的起始标志、BUFR编码数据的长度和BUFR的版本号。
表1 指示段编码说明
2.2 1段——标识段
标识段指示数据编码的主表标识、数据源中心、数据类型、数据子类型、表格版本号、数据的生产时间等信息。
2.3 3段——数据描述段
数据描述段主要指示BUFR资料的数据子集数目、是否压缩以及数据段中所编数据的要素描述符。
2.4 4段——数据段
数据段包括本段段长、保留字段以及数据描述段中的描述符（3 07 196）展开后的所有要素描述符对应数据的编码值。
2.5 5段——结束段
结束段编码说明见表5。
3 自定义描述符和标志表
3.2 C-13-L 本地数据子类型
3.3 自定义代码表0 02 201本地地面传感器标识
3.4 自定义代码表0 20 209作用层情况
3.5 自定义代码表0 20 210作用层状况
4 WMO代码表
4.2 代码表0 02 001 测站类型
4.3 代码表0 08 023 一阶统计
说明：所有一阶统计都取原始数据描述符中定义的单位。
4.4 代码表0 31 021 附加字段意义。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1960.4-1961.2 1965.9-12 1966.8 1967.3-1969.12 1974.1-1975.12
54161
长春
吉林
2
4354
12513
2368
1959.1—2000.12
1967.2 1982.1-12
54292
延吉
吉林
3
4253
12928
1768
45070
1961.1—2000.12
1967.6 1968.7-1972.4 1993.1-4
55591
拉萨
西藏
1
2940
9108
36487
1961.1—2000.12
1967.9 1968.4-1970.12 1985.9
56029
玉树
青海
3
3301
9701
36812
1960.4—2000.12
58847
福州
福建
2
2605
11917
840
1961.1—2000.12
1968.12
59082
韶关
广东
2448
11335
693
1961.1—1990.12
59287
广州
广东
1
2310
11320
410
1973.3-1978.7
56043
果洛
青海
3
3428
10015
37190
1993.1—2000.12
56137
昌都
西藏
2
3109
9710
33060
1961.1—2000.12
1968.12 1969.1-12 1970.2
56146
甘孜
四川
3
3137
10000
33935
1994.1—2000.12
11109
10203
1970.8—1991.12
53487
大同
山西
3
4006
11320
10672
1961.1—2000.12
53543
东胜
内蒙古
3
3950
10959
14604
1992.1—2000.12
53545
伊金霍洛旗
内蒙古
3934
10944
13293
1961.1—1991.12
1970.2-4
安康
陕西
3
3243
10902
2908
1990.1—2000.12
57411
南充
四川
3047
10606
3093
1973.6—1990.12
1989.12
57432
万县
四川
3046
10824
1867
1961.1—1990.12
1967.6-12
57461
宜昌
湖北
3
3042
11118
1331
1961.1—2000.12
中国气象辐射资料台站信息表
区站号
站名
省名
级别
纬度
经度
拔海高度
资料年代
缺测年月
50136
漠河
黑龙江
1
5328
12231
4330
1993.1—2000.12
1993.5-12,1994.6-1996.12
50468
黑河
黑龙江
2
5015
12727
1664
1961.5—2000.12
50527
海拉尔
内蒙古
2
3
4649
13017
812
1961.1—2000.12
1979.7 1980.6-1981.5 1982.10-11
50953
哈尔滨
黑龙江
1
4545
12646
1423
1961.1—2000.12
51076
阿勒泰
新疆
2
4744
8805
7353
1960.3—2000.12
51133
塔城
新疆
2
4644
8300
5349
1993.1—2000.12
51379
奇台*
新疆
4401
8934
7935
1987.1—1990.12
51431
伊宁
新疆
2
4357
8120
6625
1960.2—2000.12
51463
乌鲁木齐
新疆
1
4347
8739
9350
1959.1—2000.12
51567
焉耆
新疆
3
4205
8634
10553
53614
银川
宁夏
2
3829
10613
11114
1961.1—2000.12
1967.8-1972.12
53772
太原
山西
2
3747
11233
7783
1961.1—2000.12
53817
固原
宁夏
3
3600
10616
17530
1985.1—2000.12
53845
延安
陕西
3
3636
10930
9585
1959.3 1967.10
51777
若羌
新疆
3
3902
8810
8883
1957.6—2000.12
51828
和田
新疆
2
3708
7956
13745
1957.6—2000.12
52203
哈密
新疆
2
4249
9331
7372
1961.1—2000.12
1969.8
52267
额济纳旗
内蒙古
1
4157
10104
56173
红原
四川
3
3248
10233
34916
1994.1—2000.12
56196
绵阳
四川
3
3127
10445
4863
1978.1—2000.12
56294
成都
四川
1
3040
10401
5061
1961.1—2000.12
56385
峨眉山
四川
3
2931
10320
(30474)
1961.1—2000.12
1993.1—2000.12
51573
吐鲁番
新疆
3
4256
8912
345
1960.2—2000.12
51628
阿克苏
新疆
3
4110
8014
11038
1993.1—2000.12
51644
库车
新疆
4143
8257
10990
1957.7—1990.12
51709
喀什
新疆
1
3928
7559
12887
1957.7—2000.12
9405
1992.1—2000.12
52418
敦煌
甘肃
2
4009
9441
11390
1957.7—2000.12
1965.1-4
52533
酒泉
甘肃
3
3946
9829
14772
1993.1—2000.12
52681
民勤
甘肃
3
3838
10305
13670
1961.1—2000.12
52754
刚察
青海
3
3720
58265
吕泗
江苏
3
3204
12136
55
1992.1—2000.12
58321
合肥
安徽
2
3152
11714
279
1961.1—2000.12
1962.7-12 1967.7-1968.2 1972.9 1973.8 1974.12 1975.3 1977.3
58362
上海
上海
1
3124
12129
60
10008
(33015)
1993.1—2000.12
52818
格尔木
青海
1
3625
9454
28076
1957.7—2000.12
1969.9
52866
西宁
青海
2
3643
10145
22952
1959.1—2000.12
52889
兰州
甘肃
1
3603
10353
15172
1959.5—2000.12
53068
57494
武汉
湖北
1
3037
11408
231
1961.1—2000.12
1984.3 1984.6
57515
沙坪坝
重庆
2935
10628
2591
1961.1—1987.12
1967.9 1968.5 1968.7-10
57516
重庆
重庆
2