Jacchia-77大气模型

Jacchia-77大气模型，可以较为方便的计算不同高度的大气参数，可作为系统的嵌入工具

% P2Z_CIRA86    Converts pressure to altitude and temperature
%
%    Converts pressure to altitude and temperature assuming cira86
%    climataology. Linear interpolation between log（pressure）
%    latitude and day of year is performed.
%
% FORMAT   [ZT] = p2z_cira86（ P  LAT  DOY ）
%        
% OUT   Z     Altitude in meter for all combinations of P LAT DOY 
%             Size is [npnlatndoy]
%       T     Temperature in kelvin for all combinations of P LAT DOY 
%             Size is [nznlatndoy]
% IN    P     Pressure in pascal. Can be a scalar or column vector
%       LAT   Latitude. Can be a scalar or column vector [-90 to 90].
%       DOY   Day of year. Can be a scalar or column vector

% 2007-11-01   Created by Bengt Rydberg

function [zt]=p2z_cira86（platdoy）
                                                                   %&%
rqre_datatype（ p @istensor1 ）;                                    %&%
rqre_datatype（ lat @istensor1 ）;                                  %&%
rqre_datatype（ doy @istensor1 ）;                                  %&%
if any（doy）>367 | any（doy<1）                                       %&%
  error（‘doy range is [1367]‘）                                    %&%
end                                                                %&%
if any（lat）>90 | any（lat<-90）                                      %&%
   error（‘latitude range is -90 to 90‘）                            %&%
end                                                                %&%

% Create a a mjd matching the given doy （ie. the inverse of mjd2doy）
mjd = date2mjd（ 2000 1 1 0 0 0 ） + doy - 1;


cira86 = fullfile（ atmlab（‘ARTS_xmlDATA_PATH‘） ‘climatology‘ ‘cira86‘ ...
                                                               ‘cira86.z.xml‘ ）;

G = gf_artsxml（ cira86 ‘CIRA86‘ ‘z_field‘ ）;

G = atmdata_regrid（ G { p lat 0 mjd } ）;

if length（mjd） > 1
  z = reshape（ G.DATA [size（G.DATA1） size（G.DATA2） size（G.DATA4）] ）;
else
  z = G.DATA;
end

if nargout > 1
  cira86 = fullfile（ atmlab（‘ARTS_xmlDATA_PATH‘） ‘climatology‘ ‘cira86‘ ...
                                                               ‘cira86.t.xml‘ ）;
  G = gf_artsxml（ cira86 ‘CIRA86‘ ‘t_field‘ ）;

  G = atmdata_regrid（ G { p lat 0 mjd } ）;

  if length（mjd） > 1
    t = reshape（ G.DATA [size（G.DATA1） size（G.DATA2） size（G.DATA4）] ）;
  else
    t = G.DATA;
  end
end  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     154112  2017-07-03 20:41  Jacchia-77大气模型\Jacchia-77 Atmospheric Model （en）.exe

     文件      85170  2017-07-03 21:05  Jacchia-77大气模型\Jacchia77.out

     文件       2458  2017-07-03 21:10  Jacchia-77大气模型\新建文本文档.m

     目录          0  2018-06-25 00:15  Jacchia-77大气模型

----------- ---------  ---------- -----  ----

               241740                    4