主要内容

atmoscira

使用COSPAR国际参考大气1986模式

语法

[T高度纬向风]=大气(纬度,ctype,coord,mtype,,行动)

描述

[T高度纬向风]=大气(纬度,ctype,coord,mtype,,行动)实现空间研究委员会(COSPAR)国际参考大气(CIRA)1986年模型的数学表示。CIRA 1986年模型提供了平均气候学。平均气候学包括温度、纬向风和位势高度或压力。它几乎涵盖极间覆盖(80度S至80度N)0公里至120公里。该规定还包括对流层、中层大气和低热层。使用此数学表示作为压力或位势高度的函数。

此函数使用J.Barnett于1990年7月以ASCII格式提供的CIRA数据文件的更正版本。

此函数具有CIRA 1986模型的限制,并限制CIRA 1986模型的值。

CIRA 1986模型将数值限制在地球上80°S到80°N的区域。它还将位势高度限制在0 km到120 km之间。在每个月平均数据集中,该模型忽略了101300帕斯卡或0 m时80°S的数值。它忽略了这些数值,因为这些水平在南极大陆范围内。对于纬向平均压力i在恒定高度坐标下,20公里以下的气压数据不可用。因此,该值是CIRA气候学的最底层。

输入参数

纬度

的数组M大地纬度,以度为单位,其中北纬为正,南纬为负。

ctype

数据类型字符串的坐标类型的表示形式。请指定:

“压力”

压力在帕斯卡

“GPHeight”

重力势高度,单位为米

coord

取决于值ctype,此参数指定以下数组之一:

M

压力帕斯卡

M

重力势高度,单位为米

mtype

数据类型字符串的平均值类型:

“每月”(默认)

每月的值。

“年”

年度值。时有效的ctype价值“压力”.

选择模型取平均值的月份的标量值。这个论点只适用于mtype价值“每月”.

1.(默认)

一月

2.

2月

3.

前进

4.

4月

5.

五月

6.

六月

7.

七月

8.

8月

9

九月

10

十月

11

十一月

12

12月

行动

数据类型字符串的超出范围输入的操作。指定一个:

“错误”
“警告”(默认)
“没有”

输出参数

T

数组的温度:

如果M“每月”

的数组M在开尔文温度,

如果mtype“年”

的数组M7值:

  • 年平均温度(开尔文)

  • 以开尔文为单位的年温度周期振幅

  • 最高月份的年度温度循环阶段

  • 每半年开氏温度周期振幅

  • 最高月份的半年温度循环阶段

  • 以开尔文为单位的三年期温度周期振幅

  • 最高月份的三年生温度循环阶段

高度

如果mtype“每月”,数组M位势高度或M空气压力:

如果ctype“压力”

数组M位势高度

如果ctype“GPHeight”

数组M空气压力

如果mtype“年”,数组M-位势高度的by-7值。该函数仅为北半球定义此阵列(纬度大于0).

  • 年平均位势高度(米)

  • 年位势高度周期振幅,以米为单位

  • 年位势高度周期相位(以最大月份为单位)

  • 每半年位势高度周期振幅以米计

  • 以最大月为单位的半年位势高度周期相位

  • 每三年位势高度周期振幅以米计

  • 年轮位势高度周期在最大值月的阶段

纬向风

纬向风阵列:

如果mtype“每月”

以米/秒为单位的数组。

如果mtype“年”

的数组M7值:

  • 年平均纬向风,单位为米/秒

  • 年纬向风周期振幅(米/秒)

  • 年最大月的纬向风循环阶段

  • 每半年一次的纬向风的周期振幅,单位为米每秒

  • 以最大月份为单位的半年纬向风周期阶段

  • 三年生纬向风的周期振幅,单位为米每秒

  • 以最大月份为单位的太年纬向风周期相位

例子

使用CIRA 1986模式,在纬度45度和1月101,300帕斯卡范围外的行动产生警告,计算月平均值。计算温度值(T),位势高度(alt),以及纬向风(zwind).

[T, alt, zwind] = atmoscira(45, '压力',101300)T = 280.6000 alt = -18 zwind = 3.3000

使用CIRA 1986模式,在纬度45度和20000米的10月,在超出范围的行动产生警告,计算月平均值。计算温度值(T)、压力(压力),以及纬向风(zwind).

[T,压力,zwind]=atmoscira(45,'GPHeight',20000,'Monthly',10)T=215.8500压力=5.5227e+003 zwind=9.5000

10月份在纬度45度和-30度和20000米处使用CIRA 1986模型,超出范围的动作会产生误差。计算温度值(T)、压力(压力),以及纬向风(zwind).

[T,pres,zwind]=atmoscira([45-30],'GPHeight',20000,10,'error')T=215.8500 213.9000 pres=1.0e+003*5.5227 5.6550 zwind=9.5000 4.3000

对于9月,如果超出范围的行动产生警告,则在45度和–30度纬度使用CIRA 1986模型。同时在2000帕斯卡和101300帕斯卡的纬度使用该模型。计算温度的月平均值(T),位势高度(alt),以及纬向风(zwind).

[T, alt, zwind] =大气([45 -30],'压力',[2000 101300],9)T = 223.5395 290.9000 alt = 1.0e+004 * 2.6692 0.0058 zwind = 0.6300 -1.1000

使用CIRA 1986模式在45度纬度和2000帕斯卡范围外的行动产生警告,计算年值。计算温度值(T),位势高度(alt),以及纬向风(zwind).

[T,alt,zwind]=大气(45,'压力',2000,'年')T=221.9596 5.0998 6.5300 1.9499 1.3000 1.0499 1.3000 alt=1.0e+004*2.6465 0.0417 0.0007 0.0087 0.0001 0.0015 0.0002 zwind=4.6099 14.7496 0.6000 1.6499 4.6000 0.5300 1.4000

参考文献

弗莱明,钱德拉,S,舒伯,M. R.,巴奈特,J. J.,0-120 k气温、风、位势高度和气压的月平均全球气候学1988年2月,NASA TM100697

https://ccmc.gsfc.nasa.gov/modelweb/atmos/cospar1.html

在R2007b中引入