modtran操作步骤及输入参数说明_第1页
modtran操作步骤及输入参数说明_第2页
modtran操作步骤及输入参数说明_第3页
modtran操作步骤及输入参数说明_第4页
modtran操作步骤及输入参数说明_第5页
免费预览已结束,剩余42页可下载查看

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、本文分两部分,分别是第一部分MODTRAN操作步骤图解和第二部分 MODTRAN参数输入的中文说明 第一部分MODTRAN操作步骤图解 一:计算透过率 1:首先打开MODTRAN软件,然后具体步骤见图解: 点击完后出现如下对话框: KodelKodelAlsosphere(1)Alsosphere(1)RunRun Nuaber1ofNuaber1of1 1 Z CalculationOption HQDTRAN Correlation-Kspeed 5low J ModelAtmosphere 1976USStandard J TypeofAtmosphericPath HDnzont4lP

2、ath ModeofExecution Transmittance 二1 ExecuteWithMultipleScattering NnMultipleScattering TemperatureandPressureAltitudeProfile DefaulttoModel 色 WaterVaporAltitudePinfi加 DefaulttoModel OzoneAltitudeProtile DefaulttoModel 色 MethaneAltitudeProfile DefaulttoModel 二J NitrousOxideAltitudeProtilt DefaulttoM

3、odel CarbonMonoxideAltitudeProfile DefaulttoModel OtherGase5AltitudeProfile Nf5 OutputFileOptions IncludeATMProtiles 1则 UK Help 1 开始对上图中的各个参数进行设置。 ModelAtomosphere ModelAtmosphere 1976USStandard TypeofAtmopheiicPath TropicalModelMidLatitudeSummer A ModeotExecution MidLatitudeWinler SubArcticSummer E

4、xecuteWithMultipleScattering SubArcticWinter 1976USStandard TemperatureandPressureAltitudeProfile NewModelAtmosphereRadiosondedata WaterVaporAltitudeProfile CLIMATdata V 根据自己的情况选择其大气模式。 我们根据其过程是计算大气透过率,其具体设置如下(如果想知道每个参数的详细说明还是需要看说明书,本说明只是对某个操作作具体的步骤分析图解,以中纬度夏季为例,计算其透过率): ModelAtmosphere(1)ModelAtmos

5、phere(1)- -RunNumber1of1I7CRunNumber1of1I7C CalculationOption MODTRAN 三 Correlation-Kspeed Slow ModelAtmosphere MidLatitudeSummer TypeofAtmospheriaPath SlantPath ModeofExecution Transmittance ExecuteWithMultipleScattering NaMultipleScattering TemperatureandPressureAltitudeProfile MidLatitudeSummer W

6、aterVaporAltitudeProfile MidLatitudeSummer jd OzoneAltitudeProfile MidLatitudeSummer jd MethaneAltitudeProfile MidLatitudeSummer NitrousQxideAltitudeProfile MidLatitudeSummer CarbonMonoxideAltitudeProfile MidLatitudeSummer jd OtherGasesAltitudeProfile NotUsed .d OutputFileDptinns IncludeATMProFiles

7、Next|OK Help 设置好上述参数以后,点击NEXT继续; 根据自己的需要选择水汽含量,然后点击NEXT,继续; Aerosols(2)Aerosols(2)RunHimber1ofRunHimber1of AerosolModelU 善 ed Rural-VIS=23km T SeasonalModificatioristoAerosol DeterminedbyModel I StratosphericAerosols(2-30km) BackgroundStratospheric AirMassCharacterforNavyMaritimeAerosol 0 UseCloud/R

8、ainAerosolExtension|NoCloudsorRain jd UseAivny(VSAforAerosolExtension r SurfaceMeteorologicalRange(VIS) 20.00 km WindSpeed=NavyMaritimeAerosol meters/sec 24-hrAverageWindSpeed-NavyMaritime meters/sec RainRate 0.000 min/hr GroundAltitudeAboveSeaLevel 1.000 km UseA+card厂 UseNOVAMoutput r UserSuppliedS

9、pectra Next Piev OK Help 此页面中有的参数需要计算才行,计算其方位角需要用到6s里面的子程序,各个参数设置好以后,点击NEXT,就完成了其参数的设置过程;然后点击运行; FcHadVi.nFcHadVi.nffevCareSlantIaihHIZenith132teoBetryffevCareSlantIaihHIZenith132teoBetry -hfx L L!E&iAEtrE&iAEtr皿InputwsPiltfrInputwsPiltfr&看皿 口。幅叶口曹匕4 4 EQirMlpEQirMlp ObserverHeight1QQ.QQ

10、Q 口。蚓nDandpass DnwnRangeDpstanceKm 然后就可以查看其透过率的图像,其过程如下: Database InteractiveInteractive MultiplefilesMultiplefilesDiforeiLCPlotPlotSpecFluxDiforeiLCPlotPlotSpecFlux EditPlot H2 1.000 Alm:MldLefiiudeSumimer71AMm1i250.0onM E)d:Rur3l-VIS=23fcm 0.020.(10.020.(140.050.040.050.0ED.DED.D - Title 7tans.Tot

11、ail 1 Solid(Red zTrant.H20(BandModel) NotPlatted z Tvans.UniformMixedGd 单 e* NatPlolted 二 Trani.Ozoniei NotPlolted 1rans.Tiracebases HutPlolted Tran .NitrogenContinuum NotPloUjtedl工 Tvan.WaterConHinuym NatIPlDlted Tnant.Mol 之匕 uldfScaUefing HotPlotted Irans.Aerosoh NotPlotted Tran .NitricAdd HoiPlol

12、ted AbtDfpiancBAeroiglr MolPlotted LogToUlTtariinilCaniEe NotPlotted More 0K CWCBI AtJIaScale r XAKISType Micrometers | YAxisType Linear二 MinX R MaxX 14 MmY 0 MaxY 1 XHeiaht 7,5 YHeight 6.5 XDivisions 5 YDivisions 4 XDecimal 0 YDecimal 2 Grki r ExpcwtOutputTdExcel 然后点击ExportOutputToExcel就可以看到其每个波数的透

13、过率与图像 :二 三为肿下曲却楣孙中#3 括忒W WXRHXRH就口&!ttatta:- - 1np-nEnftlEnrtI15M* PcEnidViiiK&vrCajiEPcEnidViiiK&vrCajiE _Liftdj11tdtrri_Liftdj11tdtrri C C V V4 4A A舐皿J.LBktJ.LBkt工iiuhifiiuhif IXTIIXTIHod.41EJ*tQnluHod.41EJ*tQnlu- -KlpKlp 髀(IEIE?6?6&7)(0&7)(04879448794 到此其计算大气透过率的过程就全部完成了。 计算Ra

14、diance步骤 It)delMauphtre(1)-RunBtmber CalculationOpihon MODIRAN二| ConrelaNon-Kspeed Slow| ModelAtmcKSphefe MkdLaliitudeSummei二| TjpeolAtmosphericPath SlantPath| ModeofExeculian ThermalRadiance二| ExecuteWishMultipleScttening TempeiatureandPr&ssuieAltitudeProfile DefaulttoModel| WaterVaporAltitudeP

15、rolile DefaulttoModel二| OzoneAltitudeProfile DefaulttoModel MethaneAltitudeProfile DefaulttoModel,| NdiousOxideAllitudePicfile DefaulttoModel| CaibonMonoxideAltitudeProfile DefaulttoModel-| OtherGasetAltitudeProlile NotUsed二 OutputFileOpbansIncludeATMProfilei Next OK Help 口同区 TransT&sdTransT&

16、;sd L口一5l.n-ELnUELJ_ Aer/笈口1京(2)-RunJhihet1 1pf1 1区 AerosolModelUsed Rural-VI5=2北m SeasonalModiicahons:IoAerosol DeteiminedbyModel 5tialo&pheidcAerosols(2-3 口 km| BackgroundStratospheric AnMassChavactefloiNavyMaritimeAerosol0 UseCloudl/RamAerosolExtensionHoCloudsorRam UseAimj(VSAJforAerasolExlen

17、sfeom SurfaceMetearologucalHangeVIS WindSpeed-NavyMaritimeAerosol 24-hrAvenageWindSpeed 目 NavjMarilimB RainRato GroundAltitudeAboveSc iLevdll UseA+caid Help UserSuppliedSpectra Prev OK PathTypeObscrvarHiighlLZanithAngleandFinalHleighl StrtEndInc bseiveipHeigHht|Km FinalHeightICm ZenithAngiedeg 10000

18、 0nnoo ocnnn LOODO Di0000 D.OOOD 180.00 onnoo o.ooon RadiusofEarth PMhLengthType 0.0000000 km ihoit InitialFrequency 714.0000 【115115 e.Dooooao um 8000OODQ mm FinalFrequency 125DIOUO0 (1/cmJ(1/cmJ I4.005G02 umum 14UO5BD2 niii FrequencyIncfement idooooooidoooooo- - (1/cm|(1/cm| n n01120450112045 uvnu

19、vn iT204431iT204431 mm FWHMofSlitFunction 1.OODO0OO (1/cm (1.113G157 limlim 19.G15G99 nm SlitFunctionType THIANGULAH FWHMTypeABSOLUTE z TjpeofPllaulOutput IRADIANCE PkoutFileUniU|WAVENUMHER . DelinitiBf DgrmdeType|OnlyTbEalIRadZTrans Spectralfluxible|Omit*peclralfluxble三Vi*ualInput Next 0,H2是路径末端的海拔

20、高度) 只有IEMSCT=1或2时,才可以选择多次散射。 7. M1、M2、M3、M4、M5、M6和MDEF是用于改进或提供使用者指定的温度、压力、大气分子 和剖面。当MODEL=0或7时,M1到M6以及MDEF都是零,用户提供剖面数据。 M1=1-6指定大气模式的温度和压力剖面 M2=1-6指定大气模式的H2O剖面 M3=1-6指定大气模式的O3剖面 M4=1-6指定大气模式的CH4剖面 M5=1-6指定大气模式的N2O剖面 M6=1-6指定大气模式的CO剖面 MDDF=0指定大气模式下的其他一些分子剖面 =11976美国标准剖面 =2用户输入剖面 8. IM=0读取MODEL最后输入的类型

21、 =1读取使用者输入的剖面数据 9. NOPRNT=0程序正常运行,生成tape6输出文件 =1将tape6文件内的透射率或辐亮度表格数据,大气剖面的数据最少 =-1生成tape8输出文件 =-2除了tape8文件,还生成光谱冷却比率数据clrates文件 如果NOPRNT设置为-1并进行多次散射计算,沿着视场方向上的光谱漫射和总辐射值将写入tape8文件内。 10. TBOUND0边界温度(K)。在辐射模式IEMSCT=1或2时,对于倾斜路径的 地表或是云表面。 地表反照率(如果TBOUND0,可以是在H2处反照率),SALB等与1减去表面发射率,具值在0到1之间。 负的整数值允许使用者读取

22、DATA/refbkg文件已保存的表面光 谱反照率数据 SALB负值包括:-1=新雪,-2=森林,-3=农田,-4=沙漠,-5=海洋,-6=云盖,-7=干草 CARD1A:DIS,DISAZM,NSTR,LSUN,ISUN,CO2MX,H2OSTR,O3STR,LSUNFL,LBMNAM,LFLTNM,SOLCON FORMAT(2L1,I3,L1,I4,F10.5,2A10,3(1X,A1),4X,F10.3)(MOD4.0) DIS=t,forblank在CARD1中的IMULT= 1时此参数使用,DIS设置为True(t)激活 DISORT多次散射算法,DIS设置为False(域空格)使

23、用IsaaC勺二流 算法计算多次散射 DISORT算法的方位角相关标识。DISAZM设置为True(t)将考虑太阳多次散射视场内的方位角影响,如果仅需要垂直方向辐射,太阳或观测天顶角近似垂直,或太阳的多次散射是较小辐亮度成分,那 么将DISAZM设置为False(恒E空格) DISORT算法的计算流数,推荐此参数设置为8 LSUN设置为f(或空格)将使用缺省的5cm-1光谱分辨率辐照度值; 11.SALB0 地,-8=湿草地, -9=枫叶地,-10=烧过草地。 DISAZM=t,forblank NSTR=2,4,8或16 LSUN=t,forblank LSUN设置为t(即选择FileSun

24、选项)将读取1cm-1间隔的太阳辐照度数据文件且需要ISUN参数的确定 设置三角扫描函数的FWHIW于TOA*阳辐照度的平滑 CO2体积混合比(ppmv),缺省值是330ppmv;当前推荐为360ppmv(1997) 垂直水汽柱特性。如果设置为空格或0,使用缺省的水汽柱含量值。如果第一个非空字符为g,g字符后紧跟水汽含量且单位为gm/cm2(例如g2.0)。如果第一个非空字符为a,a字符后紧跟水汽含量且单位为atm-cm(但J如a3000),或者是一个绝对值来表示水汽柱比例因子 (例如2.0表示两倍于水汽柱含量的缺省 值)。 垂直臭氧柱含量特性。如果设置为空格或0,使用缺省的臭氧柱含量值。如果

25、第一个非空字符为g,g字符后紧跟臭氧含量且单位为gm/cm2(例如g0.0001)。如果第一个非空字符为a,a字符后紧跟臭氧含量且单位为atm-cm(但J如a0.2),或者是一个绝对值来表示臭氧柱比例因子(例如2.0表示两倍于臭氧柱含量的缺省值)。 如果选择t,读取CARD1A1指定的太阳辐照度数据文件。如果LSUN选择t,仅能使用此文件, 否则使用缺省的太阳辐照度数据。LSUNFL可以设置为1,2,3,4;查看CARD1A1。如果LSUNFL设置为f,而LSUN设置为t,那么则使用DATA/newkur.dat文件。如果设置为t,从CARD1A2上读取波段模型参数文件。 否则使用缺省的(1c

26、m-1间隔) 波段模型数据库 (DATA/BMP97_01.BIN)。如果设置为t,从CARD1A3上读取用户指定仪器滤光片函数文件。SOLCON的绝对值(可能接近于+1)是作为TOA太阳辐照度比例 因子来使用。已有的数据文件(在DATA/directory目录内)newkur.dat11368.00W/m2数值相对应,chkur.dat与1359.75W/m2数值相对应,cebchkur.dat与1362.12W/m2数值相对应,thkur.dat与1376.73W/m2数值相对应。 TOA太阳辐照度不进行比例调整 输入太阳常数W/m2ISUN2的整数 CO2MX H2OSTR O3STR

27、LSUNFL=t,forblank. LBMNAM=t,forblank. LFLTNM=t,forblank. SOLCON0) 二 0 在计算中没有气溶胶或云消光 二 1 乡村气溶胶模式,缺省VIS=23km 二 2 乡村气溶胶模式,缺省VIS=5km 二 3 海军海洋型,VIS的设置依赖于风速和相对湿度 二 4 海洋型,VIS=23km(LOWTRAN模式) 二 5 城市型,VIS=5km 二 6 对流层模式,VIS=50km =7使用者定义气溶胶消光系数。读取cards2D,2Dl和2D2,使用者定 义4个海拔区域的消光系数,吸收吸收和不对称因子参数(ARUSS=USS选项对于使用者定

28、义气溶胶提供了更好的适用性) =8水平对流雾,VIS=0.2km =9发散型雾,VIS=0.5km =10沙漠型,根据风速WSS决定能见度 CNOVAM;空格缺省 =N使用海军海洋垂直气溶胶模式输出 ISEASN针对对流层和平流层来选择季节适当的气溶胶剖面。仅在2到10km剖面内的对流层气溶 胶消光系数使用 ISEASN=0季节由MODEL来决定:春一夏季节(MODEL=0,1,2,4,6,7), 秋一冬季节(MODEL=3,5) =1春一夏季节 秋一冬季节 缺省 使用者提供气溶胶光学特性(将附录A) IVULCN参数即可以控制气溶胶剖面的选择也可以确定平流层气溶胶消光类型 IVULCN=0,

29、1背景平流层气溶胶剖面和消光系数 =2中等火山气溶胶剖面和陈年火山气溶胶消光系数 =3强火山气溶胶剖面和新近火山气溶胶消光系数 =4强火山气溶胶剖面和陈年火山气溶胶消光系数 =5中等火山气溶胶剖面和新近火山气溶胶消光系数 =6中等火山气溶胶剖面和背景平流层气溶胶消光系数 =7强火山气溶胶剖面和背景平流层气溶胶消光系数 =8极强火山气溶胶剖面和新近火山气溶胶消光系数 ICSTL是用来描述气团的性质(1到10),其在海军海洋型(IHAZE=3)选择时使用,缺省值为3 ICSTL=1Openocean二2 ARUSS=空格 =USS 10Strongcontinentalinfluence ICST

30、L的数值越大,陆地的影响就越大。 ICLD指定云和雨模式的使用,从地面到相对应的云模式顶部之间的雨剖面是线性递减的。此程序 忽略了云顶的雨 无云,无降雨 积云,云底高度0.66km,云顶高度3.0km 高层云,云底高度2.4km,云顶高度3.0km 层云,云底高度0.33km,云顶高度1.0km 层/积云,云底高度0.66km,云顶高度2.0km 雨层云,云底高度0.16km,云顶高度0.66km 用层云模式(ICLD=3云模式)模拟毛毛雨,地面降雨量2.0mm/hr,在1.0km处降雨量为0.86mm/hr。 用雨层云模式(ICLD=5云模式)模拟小雨,地面降雨量5.0mm/hr,在0.66

31、km处降雨量为2.6mm/hr。 =8用雨层云模式(ICLD=5云模式)模拟中雨,地面降雨量12.5mm/hr, 在0.66km处降雨量为6.0mm/hr。 =9用积云模式(ICLD=1云模式)模拟大雨,地面降雨量25.0mm/hr, 在3.0km处降雨量为0.2mm/hr。 =10用积云模式(ICLD=1云模式)模拟暴雨,地面降雨量75.0mm/hr, 在3.0km处降雨量为1.0mm/hr。 =11读入用户定义的云消光和吸收系数。读取CARDs2D,2Dl和2D2, 用户定义4个海拔区域的消光系数,吸收吸收和不对称因子参数。 =18标准卷云模式(64m半径的冰粒) =19不可见卷云模式(4

32、仙m半径的冰粒) IVSA用于地面边界层气溶胶的陆军垂直结构算法的使用。 IVSA=0不使用 =1垂直结构算法(需要输入CARD2B) VIS指定边界层气象距离,如果设置为0,那么VIS缺省值由IHAZE指定。能见度与地表550nmICLD=0 二1 二2 二3 二4 =5 二6 二7 气溶胶消光系数关系如下, -ln(0.02) EXT550kmJ0.01159km 公式中0.01159km。是在550nm处Rayleigh散射系数。气象距离与能见度的近似关系式为:气象距 离=1.3*能见度 用户指定地表气象距离 由IHAZE来设置气象距离的缺省值 WSS指定现场的风速(在海军海洋型和沙漠型

33、气溶胶模式下) WSS=现场的风速(m/s)。此参数的使用是在IHAZE=3或IHAZE=10 WHH是在海军海洋型气溶胶模式下24小时平均风速。 WHH=24小时平均风速(IHAZE=3)。 对于海军海洋型气溶胶模式(IHAZE=3),如果WSS=WHH=0,那么风速由MODEL参数指定的缺省值来指定,见表3。对于沙漠气溶胶类型(IHAZE=10),如果WSS0 二0 GNDALT 地面海拔高度(km),可以是负值 CARD2A+(如果在CARD2内APLUS=A+)允许用户移动MODTRAN已有的气溶胶,由他们以前的位置到随机的海拔区域(可以交叠)且可以减小和增加它们的范围。重要的一点是能

34、够改变对流层的高度。CARD2A+选项不能够与NOVAM联合使用。 CARD2A+:ZAER11,ZAER12,SCALE1,ZAER21,ZAER22,SCALE2,ZAER31, ZAER32,SCALE3,ZAER41,ZAER42,SCALE4(IfAPLUS=A+) FORMAT(3(1X,F9.0),20X,3(1X,F9.0) 后三个参数是气溶胶区域2,再后三个参数是气溶胶区域3,最后三 气溶胶第i区域底部海拔高度(km) 气溶胶第i区域顶部海拔高度(km) 根据新的底部海拔高度ZAERi1来改变以前的剖面 设置为缺省值,忽略SCALEi0(ZAERi1和ZAERi2都为空格时也

35、为缺省值) 用SCALEi值乘以此区域气溶胶柱密度 用1值乘以此区域气溶胶柱密度(柱密度不变) =1。 除了ICLD=11或0,对于所有云模式CARD2A者隔!要 CARD2A标准结构(卷云模式ICLD=18或19是需要) CARD2A:CTHIK,CALT,CEXT FORMAT(3F8.3) CTHIK是卷云的厚度 CTHIK=0使用统计的厚度1km 0用户定义厚度(km) CALT是卷云的底部海拔高度(km) CALT=0使用已计算的值 0用户定义高度(相对于海平面) 前三个参数是气溶胶区域1 个参数是气溶胶区域4。 ZAERil ZAERi2ZAERil 0.0 =0或空格 在使用缺省

36、值时,SCALEi 一-1 CEXT是0.55微米处消光系数(km) CEXT=0使用0.14*CTHIK 0用户定义消光系数 CARD2A备用结构(云模式ICLD=1到10) CARD2A:CTHIK,CALT,CEXT,NCRALT,NCRSPC,CWAVLN,CCOLWD, CCOLIP,CHUMID,ASYMWD,ASYMIP FORMAT(3F8.3,2I4,6F8.3) CARD2A的这种结构改进不同于卷云的云参数。使用这个CARD会引出CARDs2E1和 CARD2E2。 将CARD2A所有变量设置为-9后,所有变量为缺省值。 CTHIK是云的垂直厚度 CTHIK0用户定义厚度(

37、km) CALT是云的底部海拔高度(km) CALT=0用户定义高度(km),是相对于地面 0使用缺省值 CEXT是云的液态水滴和冰晶的垂直相光系数 CEXT0消光系数 =0缺省值 如果用户定义云/雨剖面,NCRALT是层边界数。 NCRALT=3层边界数(CARD2E1),用户定义云/雨剖面 3使用缺省的云剖面(ICLD) NCRALT最大值为16,NCRALT必须至少为3才能定义云底、云顶和最高边界海拔,在这些地方水滴和冰晶的密度必须等于零。 NCRSPC=2云光谱数据的波长数(CARD2E2) 2使用缺省的光谱数据(ICLD) NCRSPC的最大数据为788。 CWAVLN用于定义云垂直

38、消光的参考波长 如果CWAVLN超出这个范围,则缺省值为0.55微米。CWAVLN只有在CEXT输入时使用。如果 CWAVLN超出了用户定义云光谱数据的光谱范围,那么错误信息会显示且输入中断。 CCOLWD是水滴垂直柱密度 .、)、一,、一3 CCOLWD=0玄液态水滴垂直枉伤度或数重(kmgm/m) =0云冰晶垂直柱密度或数量(kmgm/m3) 0且=105云/雨相对湿度(%) 105假设100%相对湿度 ASYMWD是云液态水滴散射的HenyeyGreenstein相函数不对称因子 |ASYMWD|=1使用用户定义或模式光谱不对称因子 ASYMIP是云冰晶散射的HenyeyGreenste

39、in相函数不对称因子 |ASYMIP|=1使用用户定义或模式光谱不对称因子 CWAVLN =0.2且=200.0参考波长(微米) CARD2B是陆军垂直结构算法(VSA)子卡片输入(DARD2内IVSA=1) CARD2B:ZCVSA,ZTVSA,ZINVSA FORMAT(3F10.3) 此条件由参数VIS、ZCVSA、ZTVSA和ZINVSA来决定 条件1:云/雾在表面;从云/雾底到云/雾顶消光随高度增加而增强。VIS=0条件2:雾/轻雾;随高度到云底其消光增强。0.5VIS=0 条件2:晴/薄雾;随高度到云底其消光增强(但消光小于条件2)。10km=0 条件3:无云幕高度(低层云与地面的

40、距离),但辐射雾层或逆温层或边界层的存在;随雾或边界 层高度增加其消光减弱。ZCVS是A=0 条件4:无云幕高度或逆温层;消光是常数。ZCVSA0且ZINVSA0输入值代表云幕高度 =0高度未知:程序将计算处一个符合条件2的值,对于条件2其缺省 值为1.8km 0输入值即云厚度 =0未知云厚度;缺省值为0.2km ZINVSA是逆温层或边界层高度(km) ZINVSA0输入即高度值 =0使用缺省值 ) incLhane(t14) 9 WM0U7) gent(b) 10 WMOL(第 nitricoxide(NO) 11 WMOL9) sulphurdioxide1SO:) 12 WMOL(IO

41、) nitrogendioxide(N()2) 13 WVKH.tll) ammonia(11 14 WMOL(12) nilricacid(1IN0,那么其他两个参数忽略 如果IHA1=0且ICLD10,使用ICLD1参数 如果IHA1=0且ICLD1=0,使用IVUL1参数 如果AHAZE和EQLWCZ者B是0,那么从IHA1、ICLD1与IVUL1参数读取缺省剖面 ISEA1在ZMDL海拔处的气溶胶的季节选择,查看ISEASN(CARD2) ICHR1用于说明在ZMDL海拔处的两区域之间或用户定义的气溶胶或云区域比较邻近的边界变化 (在IHAZE=7或ICLD=11时需要此参数) ICH

42、R1=0用户定义的气溶胶或云区域无边界变化(区域不邻近) =1有边界变化 注意:ICHR1缺省值为0,如果IHAZEW7或ICLD金11。 CARDS2D,2D1,2D2CARDS2D,2D1,2D2 用于指定所有四个气溶胶海拔区域内的气溶胶和云参数(消光系数、吸收系数和不对称因子)。只有IHAZE=7或ICLD=11(CARD2)时可读取。当ARUSS参数 (CARD2)不设置的时候,其格式如下,但当ARUSS参数设置时,其格式查看附录A。 CARD2D:(IREG(N),N=1,4)(ifIHAZE=7orICLD=11) FORMAT(415) IREG用于指定四个区域内那个区域定义气溶

43、胶和云模式(IHAZE=7orICLD=11)。它用于控制 CARDs2D1和2D2的数目。 注意:区域边界海拔高度的缺省值是02,310,1130,35100km,但此数据可以(在 MODEL=7条件下)被IHA1参数(CARD2c3)覆盖 IREG(N)=0各区域N使用缺省值,N=1,2,3和4 =1输入消光系数、吸收系数和不对称因子 CARD2D1:AWCCON,TITLE FORMAT(E10.3,18A4) AWCCON是由等效的液态水含量(gm/m3)转化为消光系数(km-1)的转化因子。它的数值是 等效的液态水1.0km-1上的消光系数(0.55微米波长上),单位为/1-3 (k

44、mgmm) TITLE用户指定气溶胶或云区域文件 CARD2D2:(VX(I),EXTC(N,I),ABSC(N,I),ASYM(N,I),I=1,47) FORMAT(3(F6.2,2F7.5,F6.4) 注意:如果IREG(N)=1,CARD2D2包括47组(每组4个数);如果IREG(N)1(查看附录时,此CARD输入用户定义气溶胶或云的消光系数、吸收系数 气溶胶或云系数的波长(如果USS选项不使用,见附录A)。如果 IREG(N)=1,VX(I)见表格10。 气溶胶或云消光系数,其归一化在0.55微米处为1km-1 气溶胶或云吸收系数,其归一化在0.55微米处为1km-1 气溶胶或云不

45、对称因子参数A)。在IHAZE=7orICLD=11 VX(I)= EXTC(N,I)= ABSC(N,I)= ASYM(N,I)= I;ibk-IdTheXAmivuth(heRequiredWaxelennihslorthelultinh-ReadCARP212. Index WinvlenjTlh Index elength 1 .2000 75 90000 2 .3000 26 92000 呻j a3J/1 27 10.0000 4 ,5500 2K 10.5910 5 .6943 29 11,0000 6 1.0600 30 115000 7 1.5360 31 12.5000 *

46、2,0000 32 14.8000 9 2.2500 33 15,0000 10 Z5000 34 l(U(K)0 11 2.7000 3? 17.2000 12 3.0000 36 18.5000 13 33923 37 213000 14 3.7500 3X 25.0000 )5 4.5000 39 30.0000 16 5川000 1口 4口.0000 17 55000 41 50.00(0 IK 6.0000 42 6t)P0000 19 6.2000 43 国)川000 20 6.5000 44 100.()000 21 7.2000 45 150.0000 22 7.9000 46

47、 200,000() 23 KJOOO 47 300.0000 24 &7。仆0 CARDS2E1和2E2允许用户来控制非卷云的云参数。如果NCRALT=3(CARD2A),那么读取 CARD2E1;如果NCRSPC=2(CARD2A),那么读取CARD2E2。 CARD2E1:(ZCLD(I,0),CLD(I,0),CLDICE(I,0),RR(I,0),I=1,NCRALT) (ifICLD=1-10,NCRALT-3) FORMAT(4F10.5) ZCLD(I,0)层边界I的高度(相对地面),km,ZCLD必须是单调性的递增,不能高于 MODTRAN大气高度100kmo CLD

48、(I,0)在ZCLD(I,0)高度处的液态水滴密度(g/m3),密度不能是负的。 CLDICE(I,0)在ZCLD(I,0)高度处的冰晶密度(g/m3),密度不能是负的。 内的变量RAINRT输入降雨量时,那么降雨量将替换此参数。因此,如果用户定义降雨量剖面,那么变量RAINRT不能为正值。 CARD2E2:(WAVLEN(I),EXTC(6,I),ABSC(6,I),ASYM(6,I),EXTC(7,I),ABSC(7,I),ASYM(7,I),I=1,NCRSPC)(ifICLD=1-10,NCRSPC_2)FORMAT(7F10.5) CARD2E2内变量用于输入用户定义云光谱数据。 如

49、果CARD2A的CEXT,ASYMWD和ASYMIP参数使用缺省值,那么MODTRAN将使用这些光谱数据作为输入。如果正的垂直云消光系数CEXT有输入,那么消光和吸收系数曲线就会根据CEXT数据进行比例缩放。同样地,如果 CARD2A的不对称因子ASYMWD和ASYMIP小于1,那么他们将取代ASYM(6,I)和ASYM(7,I)数值。 WAVLEN(I)光谱波长(微米)。第一个波长WAVLEN(1)可以低于0微米,且波长 必须必须是递增的。 如果正的垂直云消光系数CEXT有输入(CARD2A),那么参考波长(CWAVLN或0.55微米)必须在WAVLEN(1)与WAVLEN(NCRSPC)之

50、间,否贝UMODTRAN将报错。 EXTC(6,I)WAVLEN(I)波长处液体水滴消光系数(km-1m3/g),如果输入EXTC(6, I)是负值,那么其数值就通过云模式ICLD的缺省值内插获得。 ABSC(6,I)如果大于零:WAVLEN(I)波长处液态水滴吸收系数(km-1m3/g) 如果小于零:WAVLEN(I)波长处液态水滴散射反照率,-1 如果ABSC(6,I)输入值小于-1或其值超过了WAVLEN(I)波长处消光系数,那么ABSC(6,I)由最初确定的云模式ICLD的缺省吸收与消光系数比值来计算, 即用这个比来乘以EXTC(6,I)。ABSC(6,I)输入值大于-1的负值则使用1

51、减去散射反照率后取负,-(1-必)。 WAVLEN(I)波长处液态水滴Henyey-Greensteir射相函数不对称因 如果CARD2A的ASYMWD的大小小于1,那么这些输入就可以忽略。如果ASYM(6,I) 不在-1与1之间,那么ASYM(6,I)数值就通过云模式ICLD的缺省值内插获得。ASYM(6,I) EXTC(7,I)WAVLEN(I)波长处冰晶消光系数(km-1m3/g) 如果输入负值,EXTC(7,I)数值就通过标准的卷云模式(ICLD=18)的消光系数内插获得。 ABSC(7,I)如果大于零:WAVLEN(I)波长处冰晶吸收系数(km-1m3/g) 如果小于零:WAVLEN

52、(I)波长处冰晶散射反照率,-1如果ABSC(7,I)输入值小于-1或其值超过了WAVLEN(I)波长处消光系数,那么ABSC(7,I)由最初确定的标准卷云模式(ICLD=18)的缺省吸收与消光系数比值来计算, 即用这个比来乘以EXTC(7,I)。ABSC(7,I)输入值大于-1的负值则使用1减去散射反照率后取负,-(1-Cl)o ASYM(7,I)WAVLEN(I)波长处冰晶Henyey-Greensteir射相函数不对称因子 如果CARD2A的ASYMIP的大小小于1,那么这些输入就可以忽略。如果ASYM(7,I)不在-1与1之间,那么ASYM(7,I)数值就通过标准卷云模式(ICLD=1

53、8)的波长内插获得。 CARD3:H1,H2,ANGLE,RANGE,BETA,RO,LENN,PHI FORMAT(6F10.3,I5,5X,F10.3) CARD3是用于定义几何路径参数 H1=初始海拔高度(km) H2=最终海拔高度(km)(ITYPE=2) =正切高度(km)(ITYPE=3) 这里强调在辐射消光模式IEMSCT=1或2下,H1(初始海拔高度)通常定义为观测者(或传感器)的位置。H1与H2在透过率模式下(transmittancemode不可以交换。 ANGLE=从H1测量的天顶角(度) RANGE=路径长度(km) BETA=针对H1与H2的地球中心张角(度) 地球半

54、径(km) 如果RO最左边为空格,程序将使用适中的值6371.23(MODEL=7),另外,适于标准模式大 气(MODEL参数)的地球半径见表格11 RO Ixhle1LDLRUIL卜心I山RudiiLMiMDilRmModelALinospheres. Mndd ModelUmosphere EarthRidiiu,RO(km) 0 Iscr-defincd(1orizonialPut!】) X(uusi;d 1 1mpical 637g39 2 MiJlulilLiilcsiiitiiticr 6371.23 个3 Midlatitude.19%11 6371,23 4 SLibarcti

55、csumtner 6356.91 q LIbarelicwinter 6356.91 6 U.S.Standard 637123 7 Usefulefilied 6371.23 LENN=由以下所描述的2a和2e条件下来确定短路径和长路径 如果LENN=1,那么路径将是“长”, 延长正切高度; 如果LENN=0(缺省值), 那么路径将是“短”。 PHI=相对于H1(传感器或初始海拔)的H2(目标或最终海拔)上天顶角 没有必要指定CARD3上的每个变量;根据ITYPE参数来说明这些问题,描述如下。 (1)水平路径(ITYPE=1) (a)指定H1,RANGE (b)如果没有标准大气模式数据可用即

56、MODEL=0(CARD1),那么由CARD2C来详细说明。 (2)两个海拔高度间的垂直或倾斜路径(ITYPE=2) (a)指定H1,H2,ANGLE,和LENN(如果H2H1,使用LENN) (b)指定H1,ANGLE,和RANGE (c)指定H1,H2,和RANGE (d)指定H1,H2,和BETA (e)指定H2,H1,PHI,和LENN(如果H10andRANGE0)THEN CASE2f ELSEIF(PHI0)THEN CASE2e ELSEIF(BETA0)THEN CASE2d ELSEIF(RANGE0ANDANGLE0)THEN CASE2b ELSEIF(RANGE0)T

57、HEN CASE2c ELSE CASE2a ENDIF ITYPE=3,下面的程序用于确定几何输入 IF(PHI0)THEN CASE3c ELSEIF(H2=0)THEN CASE3a ELSE CASE3b ENDIF 表格12给用户列出了CARD3上不同类型大气路径的选择情况 对于太阳直射或月亮辐照度,ITYPE=3的路径(CARD3)结构如下 ALTCARD3:H1,H2,ANGLE,IDAY,RO,ISOURC,ANGLEM FORMAT(3F10.3,15,5X,F10.3,15,F10.3) 观测者海拔高度 顺序天数(ase III 112 Anule BlIA LIAS(Op

58、tional! Pill 2a 革 .: 卓 (*) 2b * 2c * : * 2d : 2e * (*) * * * 3a 东 :; 3b 家 * JC 琳 * * LENNhwedonlywhenIIIH2sndCaw2a,ocH2H1andC ZLSC-COliicrwisv.I.IASisMlmLlicprogmi. *Requiredfnpuls Inhk12AllowedCombmationsofSleuv PathEiraincicrs. H2 太阳或月亮路径正切高度 ANGLE 在H1处的太阳或月亮天顶角 IDAY RO 地球半径(根据MODEL缺省值) ISOURC 辐射源

59、是太阳 辐射源是月亮 ANGLEM 月亮相位角,此角度由太阳、月亮和地球来组成(如果ISOURC=1,此 参数需要;ISOURC=0,由CARD3A1提供太阳参数) 不是H2就是ANGLE应该指定。如果两者都是0,那么就假设为垂直路径(ANGLE=0)。如果 两者都大于零,那么ITYPE=3的程序被调用。如果IDAY没有设置,那么就假设为平均日地距离。 CARDS3A1AND3A2CARDS3A1AND3A2 用于太阳/月亮的散射几何条件 当CARD1内的IEMSCT=2,这些可选的CARDs来控制太阳/月亮的散射几何条件和气溶 胶散射相函数的选择。 CARD3A1:IPARM,IPH,IDA

60、Y,ISOURC(IfIEMSCT=2) FORMAT(415) IPARM= 0,1,2,10,11,12将指定CARD3A2内太阳/月亮的散射几何的方法 选择HenyeyGreenstein气溶胶相函数(查看CARD3A2) 选择用户定义气溶月$相函数(查看CARD3B) 选才?Mie理论生成的MODTRAN模式下气溶胶相函数内部数据库 顺序天数,从1到365用于指定日地距离和指定太阳在天空中的 位置(IPARM=1) 辐射源是太阳 辐射源是月亮 CARD3A2:PARM1,PARM2,PARM3,PARM4,TIME,PSIPO,ANGLEM,G (IfIEMSCT=2) FORMAT(8F10.3) PARM1,PARM2,PARM3,PARM4参数由CARD3A1内IPARM来决定其定义。查看表1

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论