shown in Table 1. Several transitions are concatenated and separated by the UNIX end-of-line character “\n”.e commands to select the parameters of the O2 line selected for SMILES-2 (isotopologue tag=71) and those of nearby lines are shown in Table 2. e HITRAN database [16] is used and the data are stored per isotopologue into Python dictionaries (hitranZip). We also create the string variable (outCabs_lines) to be used with AMATERASU-V1.3.2Radiative transfer with no polarizationBoth versions of AMATERASU are used to compute outgoing radiation near the SMILES-2 O2 line without considering the Zeeman eects. Figure 6 shows spectra computed for a LOS with a tangent height of 52 km and FiF6SMILES-2 O2 line without and with Zeeman effect. A line from O3 is seen near 115 MHz. The Jacobian is computed without Zeeman effect and 1.5 km vertical resolution for the profiles in the vector x (Eq. 2). TabT3Definition of the atmospheric state with AMATERASU as well as the frequencies (10000 points near the O2 line) and the LOS tangent heights (50 points between 30 and 120 km). We also create a magnetic field profile B with spherical coordinates shown in Fig. 5 (right-most panel).import amaterasu.atm.atm as AAdicAtm = AA.Atm().getProf( zkm= np.linspace(20,130,151))print (dicAtm.keys()) #’p’, ‘t’, ‘z’, 1, 2, 3, …nalt = dicAtm[“z”].size #150, Number of altitudesB = np.array( [[60.e-6,0.,0.]]*nalt ) #array [nalt,3]: |B| (T), theta and phi (deg)FGHZ=np.linspace(773.834-3, 773.834+3, 10000) #Frequencies (GHz), nF = 10000TANH= np.linspace(30,120,50) #Tangent-heights (km), nTanh=50TabT4 Unpolarized radiative transfer with AMATERASU-V1import amaterasu.common.main as AMama = AM.Amaterasu()ama.inCatm[“zabs”] = dicAtm[“z”]ama.inCabs[“f”] = FGHZama.inCrt[“flag”]=1 #Radiative transfer computationama.inCrt[“unit”] = “rj” #Radiance unit in Rayleigh-Jeans brightness temperatureama.inCrt[“ObsHeight”] = [600.] #Observation altitude (needed for LOS elevation angles)ama.inCrt[“RadiusEarth”] = 6370. #Earth radius (km)ama.inCrt[“LimbStep_km”] = 3. #Integration segment size along LOS ama.inCrt[“ObsTanH”] = TANHama.setCout() #Initiatialize the model with the given settingama.outCabs[“lines”] = outCabs_amalines.replace(“\n”, “”) ama.outCatm = dicAtm #replace the default atmosphere with that created previouslyama.setAabs() #Absorption coefficient -> ama.outAabs[“abscoef”], array [nF,nalt,1]ama.setArt() #Radiative transfer radv1 =ama.outArd[“rad”][:,:,0] #Numpy array [nTanh,nF,1]112 情報通信研究機構研究報告 Vol. 65 No. 1 （2019）4 衛星センサによる宇宙からの地球環境観測
元のページ ../index.html#116