Groundbased测量

地面测量特征的海洋大气相互作用带出使用几种类型的航空公司作为平台,越来越多的测量设备,如商业、钓鱼和研究血管,锚和漂流浮标、固定基础(特别是,钻井平台),飞行员气球,俘虏气球、飞艇和飞机。最常用的载体是商业和渔船。在1984年大约有7690。因此,每日气象测量的总数进行这些船是7690 x 0.3 * 4 = 9228,在第二个图左边是船舶在海上的时间分数,第三个图是每天收到的报告。布拉克内尔地区气象信息中心在实践中,例如,每天收到约3700条消息,也就是说,40%的预期。

人能够接受这个损失的信息如果其另一部分或多或少地均匀分布于大海。不幸的是,情况并非如此,从表1.2是显而易见的,它列出了5度的百分比信息广场有100多,有超过30个样品每月的热带海洋。可以看到,观测数据的可用性在大西洋,印度洋和太平洋海洋亟待改变。此外,事实证明,这些数据不是绝对独立的,不是随机分布在时间和空间(他们聚集在附近的岛屿和船舶航线)。很明显,甚至在密度略有增加和/或改进质量的测量

表1.2。5度的百分比广场> 100 > 30样品每月的热带大西洋,印度洋和太平洋(泰勒,1985)

> 100

> 30

大西洋

印度

太平洋

大西洋

印度

太平洋

采样区域

海洋

海洋

海洋

海洋

海洋

海洋

1967年9月

25

10

10

70年

45

40

1977年9月

30.

15

15

55

40

35

1978年12月

50

40

40

80年

60

55

1980年12月

80年

90年

60

One hundred.

95年

80年

注:值为1980年12月,包括岛和浮标报告。

注:值为1980年12月,包括岛和浮标报告。

除了这些地区会导致明显的平均值的变化特征。

这些初步的评论后,我们给一个简短的概述的一些方法用于地面测量海洋表面大气层的参数及其精度。

海洋表面温度的海洋表面温度测量,利用non-insolated桶,冷凝器或冷却系统摄入,船体管道或通过船体传感器。这些方法的应用志愿观测船舶来自不同国家的船舶(总数的百分雷竞技手机版app比为每个国家舰队)展示在表1.3。这表明,在1980年代中期最常见的海洋表面温度测量的方法是冷却系统摄入的方法。但是桶方法,采用在工业时代的黎明,仍被广泛使用。

这两种方法产生不同的温度价值观:桶值比摄入值降低0.3°C。如果考虑到质量从第一过渡到第二个方法发生在1930年代末,那么这可能是一个原因,其中,全球海洋表面温度上升从测量数据获得在1940年代。但这不是全部。随后的检查显示(泰勒,1985),桶和non-bucket温度之间的差异值在时间和空间变化。特别是在高热流情况(如墨西哥湾流,黑潮或Agulhas地区在冬季)这可能意味着1°C以上的差异。

现在,50年前,还没有达成一致的方法更有效。每种方法都有自己的缺点。的缺点

表1.3。比例的船只使用每种方法六大国家舰队和温度测量的所有志愿观测船舶(世界气象组织,1984年)

方法

表1.3。比例的船只使用每种方法六大国家舰队和温度测量的所有志愿观测船舶(世界气象组织,1984年)

国家

摄入

船体

苏联

_

One hundred.

_

美国

2

98年

- - - - - -

日本

4

96年

- - - - - -

德意志联邦共和国

94年

4

2

英国

91年

5

4

加拿大

96年

4

- - - - - -

所有船只

26

73年

1

桶方法non-unified相关测量技术和决心,严格来说,不是温度的海洋表面温度,但有限的近地表水层深度。non-bucket方法的缺点是相似的。我们注意到在这个连接的强烈依赖温度测量船舶大小、深度的摄入量和摄入温度计的位置。

因此,我们需要国家船舶海洋表面温度测量的误差更加超过+ 0.1°C。

传统船舶空气温度和湿度的测量空气温度使用温度计安装在屏幕上或从便携式干湿表干球测量位于顶部的前桅或桥顶部的迎风面,也就是说,从一个位置的扭曲影响船的船体是最小化。

屏幕上的空气温度可能有偏见的高值由于百叶箱的辐射加热,或由于船体的热辐射。平均,这样测量的误差±0.5°C和不少于±0.1°C,根据泰勒(1985)。空气温度和通风干湿表受到辐射误差值必然是高于屏幕。此外,湿度测量的准确性主要取决于暴露的位置和方法。在这方面,选址正确屏幕更可取,因为错误可以更容易评估。空气温度测量结果得到的苏联船只门证实这一结论。事实证明(泰勒,1985),当测量产生的2.5繁荣扩展从桥上的船和剖面浮标测量的区别是0.04 + 0.12°C。

近地表湿度数据通常是由干湿球温度的湿度测定法。很明显,扭曲由船舶辐射加热和屏幕应该考虑在这种情况下。但由于增加空气温度和相对湿度下降部分相互补偿,比湿是或多或少地准确地决定的。显然,最严重的错误来自污染湿球温度的灯芯(由盐沉积)说,从雨水的蒸发或从湿喷屏幕。今天这是一个常识问题,船舶测量往往高估空气含湿量和平均偏差是0.5 + 0.2克/公斤(泰勒,1985)。

表面风速1984年根据世界气象组织的数据有关,测定风速风向的测量风速进行只有三分之一的所有船只;在其他船只估计风速的数据可视化观测的海洋国家。这两种方法的风速决心并非没有问题。第二种方法,尤其是基于假设之间存在一一对应风力和海况也假设风等效连接波弗特估计和风速是已知的在所有情况下,至少有两个缺点:主体性和应用局限于白天。错误的风速取决于这个方法不同±2 m / s当风力等于一个点在蒲福风级±5 m / s时风力在蒲福风级= 5分,根据泰勒(1983)。

从第一个方法出现的问题造成的扭曲影响船的船体。通常认为,最具代表性的网站风速计安装前桅的顶部或弓繁荣。深入分析风的速度场附近的船证实这个结论(见图1.3)。但即使在这里仍有相当大的扭曲。此外,事实证明,错误的测定风速风向的测量风速显著依赖之间的比例风速计的高度曝光和上层建筑的高度上桥,以及相对风向从船的航线。数据在表1.4见证。可以看到,对于弓风测定风速风向的测量的准确性的前桅不同风速的1 - 8%。对于其他风,风速度可以测量不同于参考价值50%甚至更多。

表1.4。相对误差(%)测定风速风向的测量风速

相对风向

从45到115°

的位置

风速

从-45年到45°

从255年到

从115年到255°

作者

板凳上测量

(米/秒)

风(鞠躬)

风315°(梁)

斯特恩(风)

一直等啊(1974)

气象浮标

8

0

_

_

10

8

从6到-21

- - - - - -

清(1976)

弓的繁荣

8

从3到5

从45 - 11

- - - - - -

Kidwall和·赛甘·(1978)

弓的繁荣

10

从2到6

从- 1到8

从-30年到-65年

大池塘(1982)

Oceanographie浮标

10

0

4

- - - - - -

罗曼诺夫et al。(1983)

之外的影响

10

从1到4

从4到6

从6到2

域的船

模型

注意:积极的值对应于高估;负值对应于风速的低估。

注意:积极的值对应于高估;负值对应于风速的低估。

风速计船桥
高度的风速计和上层建筑的高度上桥= 1.9:1的情况下风向。
图1.4比较获得的每月平均风速值与标准船舶测量数据和测量研究气象平台,根据泰勒(1983)。

上面的估计,提出了适用于几乎瞬时(平均超过一个小的时间间隔)风的速度值。有趣的是确定多少船从浮标测量不同于适当的测量和固定平台每月的平均时间。比较的结果如图1.4所示,很明显,船的测量系统过高估计每月平均风速。证实了这个事实,如果进一步的研究,可能意味着重要和深远的气候影响。

来确定风速(温度和湿度)从不同高度大气无线电探空仪、飞船、飞机和类似的测量。

在海洋表面净辐射通量

辐射通量(辐射平衡)底层表面代表之间的差额由底层表面太阳辐射通量同化和底层表面的有效排放,定义为向上或向下长波热辐射通量之间的区别。

请注意,太阳辐射能量黑体(身体的状态热力学平衡温度6000 K)。第九其辐射强度是由普朗克公式= (2 hc2 / X5) (exp (hc / kXT) - 1) ~ \ h = 6.62 x 10”20 js是普朗克常数;k = 1.38 x 10 ~ 9 j / k是玻尔兹曼常数;c = 3 x 108 m / s的速度是光在真空中;X是波长;T是绝对温度。从这个公式,第九= 0 1 = 0,X = oo。但由于路不能承担负在波长范围内从X = 0到X = oo,至少有一个最大的功能是必须存在。区分/ A / X和将产生的表达为零得到极值的必要条件,也就是说,[(Xc / 5 kXT) + 1 - exp (hc / kXT)) = 0。

上面的方程有一个根在hc / kXt = 4.9651,从XmT =常数K = 0.2897厘米,其中Xm是波长符合第九的最大。发现二阶导数和确定其签署第九我们向自己保证函数必须有一个最大的在上面的点表示。平等Xm T = const表达了维恩定律表明,温度的上升伴随着变化最大的方向更短的波长。众所周知,下垫面和大气的温度远低于太阳的温度。因此,太阳辐射是集中在短波(从0.3到5 | im,内可见)光谱范围,和底层表面和大气排放集中在长波(从4到100海里)光谱范围。

Actinometric测量只在进行研究和海洋船只。很明显,这些信息不足以重建独立组件的空间结构的辐射平衡,因此没有其他选择,只能利用间接的方法。它们的实际实现寻找经验关系关联组件参数确定它们的辐射平衡(太阳辐射在大气边界上,太阳高度,总数量和低级云,水和空气温度),和顺向重建信息缺少标准的测量数据。与独立的测量显示(见Strokina, 1989),至少每月的平均时间,计算和观测值之间的差异的辐射平衡组件热带和温带纬度相对海洋达2 - 4%的单位。

垂直涡流动量通量,热量和水分在海洋大气界面

垂直涡流动量通量(t)、热(H)和水分(E)是由平等z = pu 'w”、H = pcpw不”和E - pw 'q, w和u '是垂直和水平(面向的方向表面风)速度波动;7“和q是温度和含湿量的波动;p和cp空气密度和定压比热容;overbar指定时间平均。估计涡通量根据这些定义它是必要的,首先,测量风速的波动,温度和湿度;其次,获得的瞬时值共变u 'w”、w不”w 'q ';第三,他们平均的限制范围内适当选择的时间间隔。

速度测量风力波动用热线风速计,声学风速计、风速计与球形米应对压力变化,杯风速计,叶片风速计和旗帜风速计注册风轴上的负载。在这个群测量设备中,声学风速计基于使用的声传播时间之间的线性相关的固定部分路径引导在一个方向,风速是最受欢迎的方向相同。声风速计不同积极与其他类型的风速仪输出信号的线性依赖于风速波动,通过简单的检测的必要的风速分量和小时间差。这些优势弥补他们的缺点——平均以上基础的长度是10厘米的。

测量的温度波动影响热电偶,热敏电阻和电阻温度计。所有这些传感器的特点是精度高和小时间差。第一个试图衡量湿度波动在1961年由戴尔(看到戴尔,1974)。干球和溼球温度计用于此目的。随着时间的推移,取而代之的是更敏感的光学湿度计,一般的工作原理是基于的依赖光吸收在一定

(通常是红外线或紫外线)波长乐队在空气中水蒸气浓度。这种湿度计的例子是一个湿度计的紫外线光源(所谓lightman-alpha湿度计)、微波refractometer-hygrometer、陀螺测试仪与石英振荡器传感器和一个红外窄带湿度计。

主要的需求设计了传感器来测量风速的波动,温度和湿度,协方差如下。首先,他们应该高灵敏度低惯性传感器。否则,艾迪通量测定将参加了大量的错误。例如,叶片风速计的使用时限为0.3年代测量垂直速度波动四米的高度从底层表面导致涡流动量通量和热量的低估了20 - 30%。同样的趋势,但对艾迪的水分通量,发生在使用干球和湿球温度的温度计(看到Kuharets等人工智能,1980)。

第二,获得统计结果稳定所需的平均时间应该足够小,一方面,排除低频振荡的影响天气,和足够大,另一方面,尽可能覆盖整个范围的动荡起伏。根据现代概念的湍流波动风速和湿度在海洋表面大气层空间尺度上从分数几十毫米,甚至数百米,在低频范围(范围的大空间尺度上)光谱本质上不同于彼此:温度和湿度光谱,特别是水平速度谱,扩展到较低的频率比垂直速度谱。这姿势问题时选择光谱最低分离湍流波动范围从天气尺度的扰动。

垂直涡通量的动量,热量和水分的情况下更简单,因为函数> Sww(有限公司)(在Sww(>)的谱密度垂直速度波动;有限公司是频率)下降到几乎为零的空间尺度上的几百米表面边界层和几公里的行星边界层。因此,cospectra coSww (co) > SwT (co), coSwq(有限公司),在大约相同的空间尺度上也应该下降为零,因此协方差紫外线= J Suw (co) da >, w不”= J SwT (>) dco w 'q”= f Swq (co) dco明确的物理意义,这些都是积分谱从尺度约10厘米(基线长度的乐器)低频函数的最小coSww(有限公司)。因此,为了提供可靠的估计垂直涡通量的频率特性测量设备必须切断(减少到零)范围内的最低coSww(有限公司)。这个条件是满足如果测量覆盖的范围10”3 < ojz / u < 5或者1 ct3 <因为/ u < 10, z是高度的测量;>是赫兹的频率。它遵循从这个,在z = 5 m和u = 5 m / s的最低频率必须等于10 ~ 3赫兹,和平均的周期必须不少于15分钟。

最后,信息的可靠性是由流动的程度的基础设备安装。如果基础不稳定,它的运动将导致扭曲的垂直风速的波动和涡动通量波频率。这些扭曲可以消除通过注册设备动作和修改数据。

我们认为简单的过程占Yolkov提出的振动效应和Koprov (1974)。

让£代表设备的垂直位移对其平均位置,让w ' 0 = w ' + u罪£6 + d / dt, u ' 0 = u ' 4 -£du / dz, T ' T + 0 =£8 T / 8 z, e”0 =£de / dz是风速的波动,温度和绝对湿度投手所扭曲的效果。这里u, w T和e是风速的平均值,温度和湿度;u ' w ' T ' e '湍流波动风速、温度和绝对湿度没有设备基础运动(£= 0);第二项右边的第一个平等描述底层表面的斜率的影响;6是倾角(纵摇角)。

让我们编译产品u ' 0 w ' 0 w w ' 0 t ' 0, ' 0 e ' 0然后平均。然后考虑到£6不与u ' w ' T和e”,而且£d£/ 8 T = 0和6罪x0 u 'w”= u ' 0 w ' 0 - uOl;du / dz, wT = - 8 t / dz,

估计这些表达式的第二条款在右手边我们假设£= 1 m, 6 = 1°, u = 10 m / s, 8 t / 8 z = 0.01°C / m, 8 u / dz = 0.1 8 e / 8 z = 0.005克/ m4。替换这些值收益u6t 8 u / 8 k 150 cm2 z / s2, u9£8 t / 8 z % 0.20°C cm / s, ud ^ de / dz x 0.75 x 10 ~ 6克/ cm2。典型的u值'w”、w不”,方能“等于400 cm2 / s2, 1°C cm / s和3 x 10”分别为6克/厘米2 s,艾迪通量误差引起的俯仰可高达30%。因此,一个必要条件的可靠性信息表面边界层的湍流特性的精确位移余量或稳定的高度测量。

这些要求很严格。也许这(加上复杂的性质和测量设备的高成本)解释了这样一个事实:如今,缺乏数据艾迪动量通量,感受到热量和湿气多达30年前当海气相互作用实验研究才刚刚开始。

1.4.2卫星测量

在讨论表面大气层的卫星测量的特点,我们的预订设备安装在卫星寄存器而不是大气层特征本身,而是不同的特征(如热辐射强度或强度电磁波散射)间接相关的那些我们感兴趣的。此外,注册是几个参数的函数特征的地球表面和大气层。自然,这使卫星数据的解释,导致需要在多个光谱测量时间间隔以单独的不同参数的影响,,主要刻度上启用备用卫星数据的地面测量数据。很明显,卫星测量的精度是有限的地面数据,因此,利用卫星信息不是质量,而是相反,其数量、提供研究全球尺度海气相互作用的可能性。

我们接下来考虑卫星测量的主要特征的表层特征和他们的能力。

海洋表面温度众所周知,地面发射的光谱代表乐队的不规则交替传输(所谓的透明光谱窗口)和热辐射的吸收。在前者,海洋表面的发射,在通过大气层,减到最小的程度,后者在最大程度上。因此,来表示海洋的温度表面看来不足以测量发射的强度在一个或另一个光谱窗口,然后确定适当的温度排放强度作为注册辐射计。

让我们检查,例如,黑体的发射。根据普朗克公式(见26页),这样一个机构的排放强度是第九= (2 hc2 / A5) (exp (hc / kAT) - 1) _1,从它小波长时,摘要(- hc / kXT)»1普朗克公式减少第九Wienn公式= (2 hc2 / X5) exp (- hc / kXT)描述突然下降的Ix紫光谱范围;另一方面,对于较大的波长,当摘要(hc / kXT)«1 + (hc / kXT),第九普朗克公式转换成Rayleigh-Jeans公式= (2 c / a4) k T描述下降慢得多的Ix在红外和无线电波的光谱范围。因此,在其他条件相同的情况下,使用的透明窗似乎红外范围(波长范围从3.4到4.2点)会比那些来自middle-infra-red范围(波长范围从8.1到14.0点),特别是从微波range1(波长范围从0.1至21.0厘米),如果没有效果的阳光,导致增加的面积镜面反射的太阳辐射,因此,排放通量和可比性反射太阳辐射。因为这个近红外线的透明窗的选择范围意味着能够接收无声的信息仅限于晚上时间,而对于一个透明窗等medium-IR范围没有限制。

在红外辐射测量范围的主要缺点(波长范围从1到15点)是海洋排放由水蒸气的影响明显减弱,当通过大气气溶胶和云。他们的影响力可以被排除在外,而不是热发射测量红外光谱范围内,一个无线电热发射的微波测量范围。这里,吸收微波发射只发生在云含有大的雨滴。但事情并非如此简单的在这种情况下,要么因为灵敏度相对较低的排放强度变化在海洋表面温度和发射率的强烈依赖当地的属性的海浪和泡沫的存在,冰和表面活性的物质。

最后,更多的辐射测量的不确定性来源之一是所谓的集肤效应(逆温层一层近地表海洋层厚度为几毫米)。注意,数据在红外和微波辐射测量的范围描述辐射和广播亮度温度分别定义为温度的黑体发射等于实际的发射。因为后者是指层厚度一毫米的顺序,反演的温度在这一层,创建主要通过蒸发,会导致低估的海洋表面温度与它的价值由传统方法相比在一层厚度从几几十厘米到几米。如果我们包括近地表温度差异的事实逆温层不是常数(见表1.5),那么建立一个通用的可能性之间的依赖辐射或射电亮度温度和海洋表面温度测量的传统方法是不可能的因为这引入额外的错误。

如今,卫星测量的均方根误差的海洋表面温度不低于0.6°C在红外光谱范围内,和1.0°C在微波范围。

1这个范围内根据无线电波的分类定义。波属于这个范围有更长的波长比长波辐射的范围。

表1.5。在近地表逆温层温差

温度

区别

源数据

作者

(°C)

评论

实验室

保尔森和帕克

1.14 - -1.81

小取回,没有波浪

测量

(1972)

山(1972)

1.0

动态风速更少

比35 cm / s

0.1

动态的风速

超过35 cm / s,波

出现和发展

Katsaros等等。(1977)

0.4 - -2.4

没有风

自然

伍德考克和

0.5 - -1.0

夜间,平静

测量

Stommel (1947)

尤因和去

0.6

夜间,风速

(1960)

小于0.5 m / s

Malevskii-Malevich

0.3 - -0.5

风速小于

(1970)

10 m / s,最有可能的

估计

克劳斯et al。(1970)

0.35

夜间

0.5

白天

Bortkovskii et al。

0.22 - -0.50

开放的海洋;风速

(1974)

小于4 m / s

0.20 - -0.30

风速随

6到8 m / s

Grassel (1976)

0.17 - -0.21

风速变化从1

10 m / s

Katsaros等等。(1977)

3.4

风速小于

1米/秒,没有波浪

斯库利(1967)

0.2

风速对

2.5 m / s,万里无云的天空

0.0

相同,但与云

Hunjua和

0.1 - -0.5

风速高达

Andreyev (1974)

6 - 8米/秒,海浪

3点

保尔森和

0.15 - -0.30

风速为5.5

辛普森(1981)

9.2米/秒

Panin (1985)

~ 2.0

风速小于

3米/秒

空气温度和湿度的测定大气中温度在任何高度(包括表面边界层)减少登记的地球辐射在二氧化碳的不同部分吸收带(波长4.3和15 | im),也就是说,在发射的光谱区间的强度主要由控制温度的垂直分布和二氧化碳的浓度,但不是由湿度和气溶胶。主要先决条件由著名的物理事实,表示各个部分的排放二氧化碳吸收带是形成于不同大气层的规定和不变的垂直分布,因此二氧化碳浓度排放之间的对应关系,这些层的温度必须单值。

在万里无云的天空的依赖第九排放强度符合波长k温度T的气氛应该形式(见Malkevich, 1973)

在9天顶角;5 x是底层表面的发射率;BX (T) = (2 hc / k3) [exp (hc / kT) - l] 1是普朗克函数;Px(£9)的函数辐射传输(比辐射通过层与厚度£入射辐射的边界层)在垂直方向形成一个角9;Px(£)是传输函数平均值9上楼梯的方向;£= p / ps等压坐标系是一个垂直变量;p是大气压力;c = 3 x 108 m / s是光速;h = 6.62 x 10”20 J s是普朗克常数;k = 1.38 x 10”9 J / k是玻尔兹曼常数;下标s表示属于底层表面。第一项的右边(1.4.1)描述底层表面的发射,发射的第二项描述大气层之间包含底层表面(£= 1)和水平&和第三项描述了从底层表面热发射反映。

如果我们假设大海第一次近似是一个完美的黑体(< 5 = 1),和第九的测量卫星的最低点(9 = 0,< 5 = 0),然后(1.4.1)的基础上

这个表达式代表弗雷德霍姆一阶方程对T(0允许温度的垂直剖面的确定如果第九辐射强度的值和传输函数Px。,相同的故事也有可能在原则上如果九(9)的角分布。但在

9)= dxBx (Ts) Px (l, 9) - Bx \ _Tm f1 bp xpx(我,9)B.inm - ^ dz,

第九= BX (TS) PX (1,0) - BxlT (m d£。

这将需要很大的错误是因为现实水平不均匀的温度场和云的存在。

因此,我们意识到解决辐射反问题理论的必要性,以确定环境排放数据的属性。这个问题是不适定的,小信号测量中的错误可以导致大的错误估计环境参数。

让我们通过一个简单的例子来说明。让我们假设一个方程的解x f (x) -一个已经被发现。很明显,x = f ~ 1 (a),在那里/ _1的逆函数f使用定理对衍生品的逆函数dx / x = (p (a) da / a。这个表达式集之间的关系的错误输入(da / a)和输出(dx / x)的系统。< p (a) = a / f (f ~ l (a)的放大系数(a)]是一个错误。基于前面的表达式dx / x,我们得出这样的结论:输出错误将在大型大(p{)或,否则,在小f”(a),在这些值接近的一个衍生的根源f (x)。

实现可接受的精度必须要么减少输入da /错误,或者,如果这是不可能的,调整问题。后者意味着函数/ (x)是被近距离函数英尺(x),这样,在感兴趣的范围,新的误差放大系数(p ^)不超过一定的限制规定的值。

因此,要获得一个不适定问题的解决方案是必要的先验信息的解决方案,允许排除大量实验数据中的错误。这些问题是翻倍的,应用于大气中,要求和测量参数是随机函数,所以在寻找适当的解决方案类型的积分方程(1.4.2),一个应用统计正规化与地球的辐射场的统计特性和大气参数作为先验信息的字段。它的使用经验表明,垂直温度分布之间的标准差24-kilometre越低层大气的重建数据的卫星和无线电探空仪的测量通常是大约1.5°C。

与垂直剖面的温度,湿度是由测量的垂直剖面的发射/水蒸气吸收带的不同部分(波长6.3和20 - 25 | im)。在这种情况下,所有其他条件相同的情况下,方程(1.4.2)写成,佤邦(£)= (pjg) jo-1“df的含水率是常压塔之间的上层大气边界(^ = 0)和四级佤邦是大气的整体含水率;q是空气含湿量。

继续阅读:热量平衡

这篇文章有用吗?

0 0