# v8 --> similar settings as in Bernard et al. (2016), i.e., microwindows, species fitted as profile (or scaled); # The apriori has been estimated from HIPPO + ACE-FTS (50-90N) # The Sa diagonal has been estimated from HIPPO + ACE-FTS (50-90N) and scaled by the layer thickness # Using HITRAN 2012 (According with Bernard et al some weak OCS lines are missing in HITRAN 2012) file.in.stalayers = /data/Campaign/TAB/local/station.layers file.in.refprofile = reference.prf file.in.spectrum = t15asc.4 file.in.modulation_fcn = /data/Campaign/TAB/ilsFiles/ils_poly2.mf file.in.phase_fcn = /data/Campaign/TAB/ilsFiles/ils_poly2.phs file.in.sa_matrix = sa.input file.in.isotope = /data1/ebaumer/tab/ocs/x.ocs/isotope.input file.in.solarlines = /data/ebaumer/Code/sfit-linelist/linelist/solar/120621/solar.dat file.in.linelist = /data1/ebaumer/tab/ocs/x.ocs/02043.721622-02058.728378.hbin #HIT2012 # Output file names file.out.pbpfile = pbpfile file.out.statevec = statevec file.out.k_matrix = k.output file.out.kb_matrix = kb.output file.out.sa_matrix = sa.complete file.out.retprofiles = rprfs.table file.out.aprprofiles = aprfs.table file.out.ab_matrix = ab.output file.out.summary = summary file.out.parm_vectors = parm.output file.out.seinv_vector = seinv.output file.out.sainv_matrix = sainv.complete file.out.smeas_matrix = smeas.target file.out.shat_matrix = shat.complete file.out.ssmooth_matrix = ssmooth.target file.out.ak_matrix = ak.target file.out.solarspectrum = solar.output file.out.g_matrix = d.complete # Definition for retrieval gases gas.layers = 47 gas.profile.list = OCS O3 CO gas.profile.OCS.correlation = T gas.profile.OCS.correlation.type = 2 gas.profile.OCS.correlation.width = 4.0 gas.profile.OCS.correlation.minalt = 0.000 gas.profile.OCS.correlation.maxalt = 100.000 gas.profile.OCS.logstate = F gas.profile.OCS.scale = 1.0 gas.profile.OCS.sigma = 0.04599566 0.06929772 0.08158473 0.09226264 0.10457941 0.11894357 0.13557441 0.15520003 0.17862274 0.20600774 0.23766667 0.27436473 0.3135012 0.35426525 0.39184109 0.42232827 0.44795542 0.46631773 0.47708333 0.47606848 0.45776483 0.42424578 0.38336206 0.34022356 0.29700998 0.25465884 0.2111356 0.17093443 0.13577354 0.10715383 0.08776228 0.0761655 0.07031301 0.06810526 0.06824017 0.06978619 0.07222712 0.07560497 0.07954985 0.08364822 0.08983117 0.09532826 0.10157186 0.10714341 0.11174826 0.1170233 0.12291308 gas.profile.O3.correlation = T gas.profile.O3.correlation.type = 2 #1 gas.profile.O3.correlation.minalt = 0.0 #10 (ORIG) gas.profile.O3.correlation.maxalt = 100.0 #150.0 gas.profile.O3.correlation.width = 4.0 #3.0 gas.profile.O3.logstate = F gas.profile.O3.scale = 1.0d0 gas.profile.O3.sigma = 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 gas.profile.H2O.correlation = T gas.profile.H2O.correlation.type = 2 gas.profile.H2O.correlation.width = 4.0 gas.profile.H2O.correlation.minalt = 0.000 gas.profile.H2O.correlation.maxalt = 120.000 gas.profile.H2O.logstate = T gas.profile.H2O.scale = 1.0 gas.profile.H2O.sigma = 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 gas.profile.CO.correlation = T gas.profile.CO.correlation.type = 2 gas.profile.CO.correlation.width = 8.0 gas.profile.CO.correlation.minalt = 0.000 gas.profile.CO.correlation.maxalt = 120.000 gas.profile.CO.logstate = F gas.profile.CO.scale = 1.0 gas.profile.CO.sigma = 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 gas.column.list = CO2 H218O H2O #H218O #CO23 CO22 #CO2 #CO21 H2O gas.column.H2O.scale = 1.0 gas.column.H2O.sigma = 1.0 gas.column.H218O.scale = 1.0 gas.column.H218O.sigma = 1.0 gas.column.CO2.scale = 1.0 gas.column.CO2.sigma = 1.0 gas.column.CO21.scale = 1.0 gas.column.CO21.sigma = 1.0 gas.column.CO22.scale = 1.0 gas.column.CO22.sigma = 0.05 gas.column.CO23.scale = 1.0 gas.column.CO23.sigma = 0.05 #gas.column.CO2.scale = 1.0 #gas.column.CO2.sigma = 0.02 gas.column.CO.scale = 1.0 gas.column.CO.sigma = 1.0 # Forward model parameters fw.raytonly = F fw.isotope_separation = T fw.delnu = 0.1000 fw.lshapemodel = 0 fw.linemixing = T #T (ORIG) fw.linemixing.gas = CO2 fw.solar_spectrum = T fw.pressure_shift = T fw.apod_fcn = T #T (ORIG) fw.apod_fcn.type = 2 fw.apod_fcn.order = 2 fw.phase_fcn = F fw.phase_fcn.type = 2 fw.phase_fcn.order = 2 fw.emission = F # Retrieval parameter rt = T rt.lm = F rt.lm.gamma_start = 1.0e5 rt.lm.gamma_inc = 10.0 rt.lm.gamma_dec = 10.0 rt.convergence = 0.1 #rt.tolerance = 0.05 rt.max_iteration = 25 rt.ifcalc_se = T rt.dwshift = F rt.wshift = T rt.wshift.type = 2 rt.wshift.apriori = 0.000 rt.wshift.sigma = 0.100 rt.slope = T rt.slope.apriori = 0.000 rt.slope.sigma = 0.010 rt.curvature = F rt.curvature.apriori = 0.000 rt.curvature.sigma = 0.100 rt.phase = T rt.phase.apriori = 0.000 rt.phase.sigma = 0.100 rt.apod_fcn = T rt.apod_fcn.apriori = 0.000 rt.apod_fcn.sigma = 0.050 rt.phase_fcn = F rt.phase_fcn.apriori = 0.000 rt.phase_fcn.sigma = 0.050 rt.solshift = T rt.solshift.apriori = 0.000 rt.solshift.sigma = 0.10 rt.solstrnth = T # Originally set as F rt.solstrnth.apriori = 0.000 rt.solstrnth.sigma = 0.1 rt.temperature = F # Microwindows and their parameters band = 1 2 3 4 band.1.nu_start = 2047.85 band.1.nu_stop = 2048.24 band.1.zshift = F band.1.calc_point_space = 4.0e-4 band.1.wave_factor = 1.000 band.1.max_opd = 180.0000 band.1.omega = 2.3923 band.1.apodization_code = 0 band.1.gasb = OCS O3 CO2 band.2.nu_start = 2049.77 band.2.nu_stop = 2050.18 band.2.zshift = F band.2.calc_point_space = 4.0e-4 band.2.wave_factor = 1.000 band.2.max_opd = 180.0000 band.2.omega = 2.3923 band.2.apodization_code = 0 band.2.gasb = OCS O3 CO2 CO H218O band.3.nu_start = 2051.18 band.3.nu_stop = 2051.46 band.3.zshift = F band.3.calc_point_space = 4.0e-4 band.3.wave_factor = 1.000 band.3.max_opd = 180.0000 band.3.omega = 2.3923 band.3.apodization_code = 0 band.3.gasb = OCS O3 H2O band.4.nu_start = 2054.33 band.4.nu_stop = 2054.67 band.4.zshift = F band.4.calc_point_space = 4.0e-4 band.4.wave_factor = 1.000 band.4.max_opd = 180.0000 band.4.omega = 2.3923 band.4.apodization_code = 0 band.4.gasb = OCS H2O O3 H218O CO2 #-------- band.5.nu_start = 2030.68 band.5.nu_stop = 2031.06 band.5.zshift = F band.5.calc_point_space = 4.0-4 band.5.wave_factor = 1.000 band.5.max_opd = 180.0000 band.5.omega = 2.3923 band.5.apodization_code = 0 band.5.gasb = O3 CO21 CO22 band.5.nu_start = 2035.35 band.5.nu_stop = 2035.45 band.5.zshift = F band.5.calc_point_space = 4.0-4 band.5.wave_factor = 1.000 band.5.max_opd = 180.0000 band.5.omega = 2.3923 band.5.apodization_code = 0 band.5.gasb = O3 CO2 #CO21 band.5.nu_start = 2045.30 band.5.nu_stop = 2045.80 band.5.zshift = F band.5.calc_point_space = 4.0e-4 band.5.wave_factor = 1.000 band.5.max_opd = 180.0000 band.5.omega = 3.1100 band.5.apodization_code = 0 band.5.gasb = OCS CO2 O3 H2O CO CO22 #CO21 band.8.nu_start = 2055.65 band.8.nu_stop = 2055.95 band.8.zshift = F band.8.calc_point_space = 4.0e-4 band.8.wave_factor = 1.000 band.8.max_opd = 180.0000 band.8.omega = 3.1100 band.8.apodization_code = 0 band.8.gasb = OCS H2O O3 H218O #CO21 CO # Original file sp.snr = 1 sp.snr.1.nu_start = 2035.35 sp.snr.1.nu_stop = 2055.95 sp.snr.1.snr = 400.0 # Kb derivative calculations kb = T kb.temperature = T kb.slope = T #ORIG WAS F kb.curvature = F #VERY VERY SMALL kb.solshft = T kb.solstrnth = T #ORIG WAS F kb.phase = T #ORIG WAS F kb.dwshift = F #ORIG WAS F kb.wshift = T #ORIG WAS F kb.apod_fcn = T #ORIG WAS F kb.phase_fcn = T #VERY VERY SMALL kb.zshift = T #ORIG WAS F kb.sza = T kb.omega = T #ORIG WAS F kb.max_opd = F #VERY VERY SMALL kb.line = T kb.line.type = 1 kb.line.gas = retrieval # Output Files Section out.level = 1 out.gas_spectra = T out.gas_spectra.type = 1 out.sa_matrix = T out.statevec = T out.g_matrix = T out.k_matrix = T out.shat_matrix = F out.retprofiles = T out.aprprofiles = T out.ab_matrix = F out.ak_matrix = T out.summary = T out.pbpfile = T out.channel = F out.parm_vectors = F out.seinv_vector = T out.sainv_matrix = F out.smeas_matrix = F out.ssmooth_matrix = T out.raytrace = T out.raytrace.type = 1 out.solarspectrum = F out.levmardet = F out.xscdetail = F