atmoscira - Use COSPAR International Reference Atmosphere 1986 model

Syntax

[T altitude zonalWind] = atmoscira(latitude, ctype, coord, mtype, month, action)

Description

[T altitude zonalWind] = atmoscira(latitude, ctype, coord, mtype, month, action) implements the mathematical representation of the Committee on Space Research (COSPAR) International Reference Atmosphere (CIRA) from 1986 model. The CIRA 1986 model provides a mean climatology. The mean climatology consists of temperature, zonal wind, and geopotential height or pressure. It encompasses nearly pole-to-pole coverage (80 degrees S to 80 degrees N) for 0 km to 120 km. This provision also encompasses the troposphere, middle atmosphere, and lower thermosphere. Use this mathematical representation as a function of pressure or geopotential height.

This function uses a corrected version of the CIRA data files provided by J. Barnett in July 1990 in ASCII format.

This function has the limitations of the CIRA 1986 model and limits the values for the CIRA 1986 model.

The CIRA 1986 model limits values to the regions of 80 degrees S to 80 degrees N on Earth. It also limits geopotential heights from 0 km to 120 km. In each monthly mean data set, the model omits values at 80 degrees S for 101,300 pascal or 0 m. It omits these values because these levels are within the Antarctic land mass. For zonal mean pressure in constant altitude coordinates, pressure data is not available below 20 km. Therefore, this value is the bottom level of the CIRA climatology.

Inputs

latitude

Array of m geopotential heights, in meters.

ctype

String that defines representation of coordinate type. Specify:

`'Pressure'`	Pressure in pascal
`'GPHeight'`	Geopotential height in meters

coord

Depending on the value of ctype, this argument specifies one of the following arrays:

`m`	Pressures in pascal
`m`	Geopotential height in meters

mtype

String that selects one of the following mean value types:

`'Monthly'` (default)	Monthly values.
`'Annual'`	Annual values. Valid when `ctype` has a value of `'Pressure'`.

month

Scalar value that selects the month in which the model takes the mean values. This argument applies only when mtype has a value of 'Monthly'.

`1` (default)	January
`2`	February
`3`	March
`4`	April
`5`	May
`6`	June
`7`	July
`8`	August
`9`	September
`10`	October
`11`	November
`12`	December

action

String that defines action for out-of-range input. Specify one:

'Error'

'Warning' (default)

'None'

Outputs

T

Array of temperatures:

If m is 'Monthly'

Array of m temperatures, in kelvin

If mtype is 'Annual'

Array of m-by-7 values:

Annual mean temperature in kelvin
Annual temperature cycle amplitude in kelvin
Annual temperature cycle phase in month of maximum
Semiannual temperature cycle amplitude in kelvin
Semiannual temperature cycle phase in month of maximum
Terannual temperature cycle amplitude in kelvin
Terannual temperature cycle phase in month of maximum

altitude

If mtype is 'Monthly', an array of m geopotential heights or m air pressures:

If `ctype` is `'Pressure'`	Array `m` geopotential heights
If `ctype` is `'GPHeight'`	Array `m` air pressure

If mtype is 'Annual', an array of m-by-7 values for geopotential heights. The function defines this array only for the northern hemisphere (latitude is greater than 0).

Annual mean geopotential heights in meters
Annual geopotential heights cycle amplitude in meters
Annual geopotential heights cycle phase in month of maximum
Semiannual geopotential heights cycle amplitude in meters
Semiannual geopotential heights cycle phase in month of maximum
Terannual geopotential heights cycle amplitude in meters
Terannual geopotential heights cycle phase in month of maximum

zonalWind

Array of zonal winds:

If mtype is 'Monthly'

Array in meters per second.

If mtype is 'Annual'

Array of m-by-7 values:

Annual mean zonal winds in meters per second
Annual zonal winds cycle amplitude in meters per second
Annual zonal winds cycle phase in month of maximum
Semiannual zonal winds cycle amplitude in meters per second
Semiannual zonal winds cycle phase in month of maximum
Terannual zonal winds cycle amplitude in meters per second
Terannual zonal winds cycle phase in month of maximum

Examples

Using the CIRA 1986 model at 45 degrees latitude and 101,300 pascal for January with out-of-range actions generating warnings, calculate the mean monthly values. Calculate values for temperature (T), geopotential height (alt), and zonal wind (zwind).

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

Using the CIRA 1986 model at 45 degrees latitude and 20,000 m for October with out-of-range actions generating warnings, calculate the mean monthly values. Calculate values for temperature (T), pressure (pres), and zonal wind (zwind).

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

Use the CIRA 1986 model at 45 and –30 degrees latitude and 20,000 m for October with out-of-range actions generating errors. Calculate values for temperature (T), pressure (pres), and zonal wind (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

For September, with out-of-range actions generating warnings, use the CIRA 1986 model at 45 degrees latitude and –30 degrees latitude. Also use the model at 2000 pascal and 101,300 pascal. Calculate mean monthly values for temperature (T), geopotential height (alt), and zonal wind (zwind).

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

Using the CIRA 1986 model at 45 degrees latitude and 2000 pascal with out-of-range actions generating warnings, calculate annual values. Calculate values for temperature (T), geopotential height (alt), and zonal wind (zwind).

[T alt zwind] = atmoscira( 45, 'Pressure', 2000, 'Annual' )
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

References

Fleming, E. L., Chandra, S., Shoeberl, M. R., Barnett, J. J., Monthly Mean Global Climatology of Temperature, Wind, Geopotential Height and Pressure for 0-120 km, NASA TM100697, February 1988

http://modelweb.gsfc.nasa.gov/atmos/cospar1.html