Instructions are best viewed with a fixed width font, e.g. Lucida Sans Typewriter Date: January 2011 USER INSTRUCTIONS FOR LBLRTM --------------------------------------------------------------------- FILE ASSIGNMENTS FOR LBLRTM TAPE3 UNFORMATTED LINE FILE WITH LBLRTM BLOCKING EXTERNAL FILE NOT REQUIRED FOR IHIRAC=0 IHIRAC=9 IHIRAC NE 0,9 ; ITEST=1 TAPE4 NLTE VIBRATIONAL TEMPERATURES (POPULATIONS) BY LAYER ONLY REQUIRED FOR IHIRAC=4 TAPE5 LBLRTM INPUT FILE TAPE6 LBLRTM OUTPUT FILE TAPE7 FILE OF MOLECULAR COLUMN AMOUNTS FROM LBLATM ONLY FOR IATM=1; IPUNCH=1 (CARD 3.1) TAPE9 FILE OF EFFECTIVE LINES FOR LBLF4 CREATED BY LINF4 TAPE10 OPTICAL DEPTH RESULTS FROM LINE BY LINE CALCULATION LAST LAYER FOR IMRG EQ 0 LAYER BY LAYER FOR IMRG EQ 1 TAPE11 SPECTRAL RESULTS FROM SCANFN AND INTRPL JEMIT=-1: ABSORPTION JEMIT= 0: TRANSMITTANCE JEMIT= 1: RADIANCE TAPE12 MONOCHROMATIC RESULTS IEMIT=0: OPTICAL DEPTH IEMIT=1: RADIANCE/TRANSMITTANCE INCLUDES AEROSOL CONTRIBUTION FOR IAERSL=1; IEMIT=1 TAPE13 MONOCHROMATIC RESULTS FOR WEIGHTING FUNCTIONS IEMIT=0: OPTICAL DEPTH IEMIT=1: RADIANCE/TRANSMITTANCE ONLY CREATED FOR IMRG= 3 TO 18 TAPE14 MONCHROMATIC RESULTS INLUDING AEROSOL CONTRIBUTION IEMIT=0: OPTICAL DEPTH IEMIT=1: RADIANCE/TRANSMITTANCE ONLY CREATED FOR IAERSL=1; IEMIT=0 TAPE20 SPECTRAL AEROSOL TRANSMITTANCES TOTAL AEROSOL CONTRIBUTION FOR IEMIT=0 LAYER BY LAYER CONTRIBUTION FOR IEMIT=1 ONLY CREATED FOR IAERSL=1 TAPE29 FILE CONTAINING VALUES OF Y FOR PLOTTING ONLY FOR IPLOT EQ 1 TAPE39 AFGL PLOT FILE OUTPUT JACOBIAN/DERIVATIVE FILES (in the AJ directory): RDderivDNW_xx_lll Layer downwelling Jacobian/derivative files RDderivUPW_xx_lll Layer upwelling Jacobian/derivative files (includes downwelling and surface components) LEV_RDderivDNW_xx_lll Level downwelling Jacobian/derivative files LEV_RDderivUPW_xx_lll Level upwelling Jacobian/derivative files (includes downwelling and surface components) Files related to surface derivatives LEV_RDderivE-R_-1_000 TOA radiance Jacobian wrt surface emittance/reflectance (with E + R = 1) LEV_RDderivEMI_-1_000 TOA radiance Jacobian wrt surface emittance LEV_RDderivRFL_-1_000 TOA radiance Jacobian wrt surface reflectance LEV_RDderivTSF_-1_000 TOA radiance Jacobian wrt surface temperature Note: -> xx signifies the parameter for which the radiance analytic Jacobian is calculated -> lll is the layer number starting from the lowest altitude -> Layer output Jacobian/derivative files: - are in the same TAPE12 two-panel binary format - first panel containing the Jacobian, and the second the total transmission - scanned layer files (IMRG=42,43) will only have a single panel (like TAPE10) -> Level output Jacobian/derivative files: - are in single-panel binary format (like TAPE10) Users desiring to use the new Analytical Jacobian capability of LBLRTM are strongly urged to run the scripts provided for that purpose. The scripts provide a good mechanism to understand the utilization of LBLRTM for Analytic Jacobians and to compare the results with symmetric finite difference calculations. ------------------------------------------------------------------------------------------ RECORD RECORD 1.1 CXID: 80 characters of user identification (80A1) CXID(1) is the flag which determines program initialization and termination. The actual input data stream for LBLRTM commences with the record containing a '$' in CXID(1). Any record that are read prior to a record containing a '$' in CXID(1) are ignored. Since LBLRTM automatically recycles and reads Record 1.1 at the end of each run (this allows the stacking of multiple runs), each subsequent run of LBLRTM must also have a '$' in CXID(1) on Record 1.1. In order to effect a normal termination of LBLRTM, CXID(1) must be set to '%', on the final record of the input file. Thus a standard LBLRTM input file would begin with a '$' and end with a '%'. Records on the input file following the record with CXID(1) = '%' are ignored by LBLRTM. Note: LBLRTM will internally set CXID(1) = ' ' for all output. RECORD RECORD 1.2 IHIRAC, ILBLF4, ICNTNM, IAERSL, IEMIT, ISCAN, IFILTR, IPLOT, ITEST, IATM, IMRG, ILAS, IOD, IXSECT, MPTS, NPTS 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 54-55, 60, 65, 70, 72-75, 77-80 4X,I1, 4X,I1, 4X,I1, 4X,I1, 4X,I1, 4X,I1, 4X,I1, 4X,I1, 4X,I1, 4X,I1, 3X,A2, 4X,I1, 4X,I1, 4X,I1, 1X,I4, 1X,I4 IHIRAC (0,1,2,3,4,9) selects desired version of HIRAC = 0 HIRAC HIRAC not activated; line-by-line calculation is bypassed (skips to selected function) = 1 HIRAC1 (Voigt profile) = 2 HIRACL (Lorentz profile, not available in LBLRTM) = 3 HIRACD (Doppler profile, not available in LBLRTM) = 4 NLTE Option (Non Local Thermodynamic Equilibrium) -state populations as a function of altitude required on TAPE4 = 9 central line contribution omitted (functions 1-3) ILBLF4 (0,1,2) flag for LBLF4 (Line-By-Line Function 4) (LBLF4 extends bound of line by line calculation beyond 64 halfwidths from line center) = 0 LBLF4 not activated (line by line bound is 64 halfwidths) (Note: if IHARAC>0 and line coupling is used in generating the input TAPE3 then ILBLF4 = 1) = 1 line by line bound is 25 cm-1 for all layer pressures (altitudes) = 2 line by line bound is 25 cm-1 for layers with pressures .GT. 0.5 mb and 5 cm-1 for layers with pressures .LE. 0.5 mb (ILBLF4 = 2 saves computational time above 50 km with minimal loss of accuracy for most cases) ICNTNM (0,1,2,3,4,5,6) flag for continuum (CONTNM) = 0 no continuum calculated = 1 all continua calculated, including Rayleigh extinction where applicable = 2 H2O self not calculated, all other continua/Rayleigh extinction calculated = 3 H2O foreign not calculated, all other continua/Rayleigh extinction calculated = 4 H2O self and foreign not calculated, all other continua/Rayleigh extinction calculated = 5 Rayleigh extinction not calculated, all other continua calculated = 6 Individual continuum scale factors input (Requires Record 1.2a) RECORD RECORD 1.2 (Continued) IAERSL (0,1,5,7,9) flag for aerosols (LOWTRN) = 0 no aerosols used = 1 internal LOWTRAN aerosol models = 5 spectral optical depths by layer from file 'in_lblrtm_cld' = 7 user defined aerosol models = 9 use precalculated aerosols (TAPE20 from a previous aerosol run) IEMIT (0,1,2,3) = 0 optical depth only = 1 radiance and transmittance (Radiance Units: W / cm^2 sr^-1 cm^-1 = 2 solar radiance (requires previously calculated optical depths or transmittances and binary solar radiation file SOLAR.RAD) = 3 radiance analytic Jacobian/derivative - requires IMRG=40,41,42 or 43 - requires a subdirectory named "AJ" to contain the analytic Jacobian/derivative files Note: - requires previously calculated optical depths files (created by setting IEMIT=0, IMRG=1, IOD=3) - units (dR/dx): temperature: { W / (cm^2 sr cm^-1) ] / [K] } molecules: { W / (cm^2 sr cm^-1) ] / [log(volume mixing ratio)] } sfc temp: { W / (cm^2 sr cm^-1) ] / [K] } sfc emiss/refl: { W / (cm^2 sr cm^-1) ] / [K] } ISCAN (0,1,2,3) flag for SCANFN 1 = scanning function or INTRPL 2 = interpolating procedure or FFTSCN 3 = Fast Fourier Transform scan IFILTR (0,1) flag for FILTR 1 = yes IPLOT (0,1) flag for PLTLBL 1 = yes ITEST (0,1) flag for TEST 1 = yes IATM (0,1) flag for LBLATM 1 = yes -------------------------------------------------------- | | | INPUT FOR RECORD 1.2 IS CONTINUED ON NEXT PAGE | | | -------------------------------------------------------- RECORD 1.2 (Continued) IMRG (0 to 9; 12-18; 22-28; 32; 35-36; 40-43; 45-46) selects merge option IMRG = 0 KFILE rewound after each layer > 1 KFILE not rewound after each layer = 0 normal merge; result on MFILE; last layer optical depth on KFILE = 1 optical depths only; results for each layer on different file = 2 optical depths from precalculated KFILE merged sequentially onto MFILE ** uses merge path as defined by precalculated KFILE ** = 9 radiance merge with aerosols; precalculated KFILE ** uses merge path as defined by precalculated KFILE ** --------------------------------------------------------- | | | NOTE: KFILE usually maps to TAPE10 | | MFILE " " " TAPE12 | | LFILE alternates with MFILE, mapping to TAPE11 | | NFILE usually maps to TAPE13 | | KODFIL usually maps to | | ODdeflt_ for IMRG=1 and IOD=0 | | ODint_ for IMRG=1 and IOD=1 | | ODint_ for IMRG=1 and IOD=3 | | ODint_ for IMRG=1 and IOD=4 | | ODexact_ for IMRG=1 and IOD=2 | | | | layer optical depth file = "ODint_" | | upwelling total optical depth file = "ODtoupw_" | | downwelling total optical depth file = "ODtodnw_"| --------------------------------------------------------- IMRG options are continued on next page. ---------------------------------------- C --------------------------------------------------------------- C C - If IOD = 1 or 4 then calculate optical depths for each C layer with DV = DVOUT (using DVSET if IOD=4) and maintain C separately. Use PTHODI as the name of the optical depth files. C This requires the format HFMODL, which is produced by C calling the SUBROUTINE QNTIFY. C C - If IOD = 2 and IMERGE = 1 then calculate optical depths C for each layer using the exact DV of each layer C Use PTHODE as the name of the optical depth files. C This requires the format HFMODL, which is produced by C calling the SUBROUTINE QNTIFY. C C - If IOD=3 and IMRG=1 then calculate layer optical depths and C and interpolate all layers to the dv of the final layer C (used for analytic derivative calculation) C C - If calculating optical depths using the default procedure, C sending output to a different file for each layer (IEMIT=0, C IOD=0, and IMRG=1), then use PTHODI for the optical depth C pathnames. C C - Otherwise, use TAPE10. For IOD=1, calculate optical depths C for each layer with DV = DVOUT (from DVSET in TAPE5, carried C in by COMMON BLOCK /IODFLG/ (interpolation in PNLINT). C RECORD RECORD 1.2 (Continued) IMRG (Continued) IMRG OPTIONS FOR MERGED OUTPUT ============================= Monochromatic Input, on precalculated KFILE. Output is determined by IEMIT. ____________________________________________________________________________ IMRG = ----------------------------------------------- MERGE PATH: space to ground ground to space tangent Final RADIANCE CALCULATION: downwelling upwelling Results OBSERVER PATH: uplooking downlooking limb on file --------------- --------------- ------- ------- monochromatic (total merged) A C B MFILE (12) (32) (22) Monochromatic Input, on precalculated KODFILs. Output is determined by IEMIT. ____________________________________________________________________________ IMRG = ----------------------------------------------- MERGE PATH: space to ground ground to space tangent Final RADIANCE CALCULATION: downwelling upwelling Results OBSERVER PATH: uplooking downlooking limb on file --------------- --------------- ------- ------- monochromatic (total merged) 40 41 MFILE scanned results 42 43 For radiance analytic Jacobians/derivatives (IMRG.eq.40,41,42 or 43), a previous run with IMRG.eq.1 and IOD.eq.3 is required to provide the required optical depth files. IMRG OPTIONS FOR MERGED OUTPUT BY LAYER ======================================= Monochromatic Input, on precalculated KODFILs. Output is determined by IEMIT. ____________________________________________________________________________ IMRG = Final ----------------------------------------------- Results MERGE PATH: space to ground ground to space tangent by layer RADIANCE CALCULATION: downwelling upwelling on files --------------- --------------- ------- -------- monochromatic (by layer, 1 file/layer) 45 46 MFILE IMRG options are continued on next page. ---------------------------------------- RECORD RECORD 1.2 (Continued) IMRG (Continued) IMRG OPTIONS FOR SEQUENTIAL OUTPUT ================================== Monochromatic Input, calculated in current run. _______________________________________________ IMRG = ----------------------------------------------- MERGE PATH: space to ground ground to space tangent Results RADIANCE CALCULATION: downwelling upwelling on file --------------- --------------- ------- ------- monochromatic 3 4 7 NFILE scanned results 13 14 17 NFILE filtered results 23 24 27 NFILE Monochromatic Input, on precalculated KFILE. ____________________________________________ IMRG = ----------------------------------------------- MERGE PATH: space to ground ground to space tangent Results RADIANCE CALCULATION: downwelling upwelling on file --------------- --------------- ------- ------- monochromatic 5 6 8 NFILE scanned results 15 16 18 NFILE filtered results 25 26 28 NFILE Monochromatic Input, on precalculated KODFILs. ____________________________________________ IMRG = ----------------------------------------------- Results MERGE PATH: space to ground ground to space tangent Results RADIANCE CALCULATION: downwelling upwelling on file --------------- --------------- ------- ------- scanned results (flux calculations) 35 36 NFILE -------------------------------------------------------- | | | INPUT FOR RECORD 1.2 IS CONTINUED ON NEXT PAGE | | | -------------------------------------------------------- RECORD RECORD 1.2 (Continued) The actual results produced on MFILE or NFILE are determined by both IEMIT on Record 1.2 and JEMIT on either Record 6 (scanned results) or Record 7.1 (filtered results). The following table describes the options: RESULTS PRODUCED ON MFILE OR NFILE ------------------------------------------ TYPE OF RESULT IEMIT=0 IEMIT=1 ------------------ ----------------- ----------------- monochromatic optical depths transmittance & radiance scanned (JEMIT=0) transmittance transmittance scanned (JEMIT=1) n/a radiance filtered (JEMIT=0) transmittance transmittance filtered (JEMIT=1) n/a radiance In the context of LBLRTM calculations, weighting functions are readily obtained from the NFILE sequential transmittances by "subtracting" adjacent files (dT or dR), and employing the stored layer boundary definitions for P or z, i.e. (dlnP or dz). ( Use Option IOPT = 2 in PLTLBL for subtraction of LBLRTM files ) For the filter option (IMRG.ge.23 .and. IMRG.le.28), the results, including the differences of the total transmittances or radiances, are also printed on unit IPR (TAPE6) For the flux computations (IMRG.eq.35 .or. IMRG.eq.36), the results are merged from the top of the atmosphere down for IMRG.eq.35 to provide the downward radiance at each level, and from the ground up for IMRG.eq.36 to provide upward radiance at each level. "SEQUENTIAL" is defined at the total optical depth, transmittance, or radiance from H1 to each layer boundary (Ln), always beginning at H1: H1 L1 L2 L3 * *----> *----|----> *----|----|----> etc... -------------------------------------------------------- | | | INPUT FOR RECORD 1.2 IS CONTINUED ON NEXT PAGE | | | -------------------------------------------------------- RECORD RECORD 1.2 (Continued) ******* ILAS (0,1,2) flag for laser options ************ not currently available in LBLRTM = 0 no = 1 single laser frequency selected (V1 on record 1.3) = 2 multiple laser frequency selected (first laser frequency is V1 on record 1.3; subsequent laser frequencies are on records 5.1 ... 5.N) IOD (0,1,2,3,4) flag for layering control in optical depth calculations = 0 normal calculation (default) = 1 uses exact calculated dv for each layer and interpolates to the spacing set by DVOUT = 2 uses exact calculated dv for each layer (requires IMRG = 1) = 3 sets DVSET to DV for last layer and interpolates all layers to the DVSET grid (requires IMRG = 1 or IMRG=40-43) = 4 uses DVSET for last layer and interpolates to the spacing set by DVOUT IXSECT (0,1) flag for cross-sections = 0 no cross-sections included in calculation = 1 cross-sections included in calculation (input file FSCDXS containing the cross-section file information must be provided. See Table II following Record 3.7.1 for a description of file FSCDXS) MPTS number of optical depth values printed for the beginning and ending of each panel as a result of convolution for current layer (for laser option ILAS > 1, MPTS should be set to 0) (for MPTS < O, output printing is suppressed) NPTS number of values printed for the beginning and ending of each panel as result of merge of current layer with previous layers (optical depth for IEMIT=0; radiance and transmission for IEMIT=1) (for laser option ILAS > 1, NPTS should be set to 0) RECORD RECORD 1.2a (required if ICNTNM = 6) XSELF, XFRGN, XCO2C, XO3CN, XO2CN, XN2CN, XRAYL free format XSELF H2O self broadened continuum absorption multiplicative factor XFRGN H2O foreign broadened continuum absorption multiplicative factor XCO2C CO2 continuum absorption multiplicative factor XO3CN O3 continuum absorption multiplicative factor XO2CN O2 continuum absorption multiplicative factor XN2CN N2 continuum absorption multiplicative factor XRAYL Rayleigh extinction multiplicative factor RECORD RECORD 1.2.1 (required if IEMIT = 2; otherwise omit) NOTE: IEMIT = 2 requires binary solar radiance file SOLAR.RAD, containing extraterrestrial source spectra in units W/(m2 cm-1). INFLAG, IOTFLG, JULDAT 1-5, 6-10, 13-15 I5, I5, 2X, I3 INFLAG (0,1,2,3) input flag for solar radiance calculation = 0 input previously calculated radiance and transmittance from TAPE12 (default) = 1 input previously calculated optical depth from TAPE12 = 2 input previously calculated upwelling radiance and transmittance from TAPE12, downwelling radiance and transmittance from 'SOL.PATH.T2', and solar reflection function from 'SOL.REFLECTANCE' (all binary files). = 3 input previously calculated transmittance from TAPE12 (CHARTS-type output) IOTFLG (0,1,2) output flag for solar radiance calculation (output to TAPE13) = 0 output attenuated solar radiance (default) = 1 output total radiance (requires INFLAG = 0) = 2 output total radiance, including effects of downwelling thermal and reflection from solar path (requires INFLAG = 2) JULDAT Julian day associated with calculation (1-365/366 starting January 1). Used to calculate Earth distance from sun. A value of 0 (default) indicates no scaling of solar source function using earth-sun distance. Note 1: For total downwelling radiance use INGLAG=0 and IOTFLG=1. Requires previous LBLRTM run with IEMIT=1 to obtain TAPE12. Note 2: For total upwelling radiance use INFLAG=2 and IOTFLG=2. Requires two previous LBLRTM runs with IEMIT=1. Total radiance given by R1+S2*T2*r1*T1 _ | Observer |-O-| | - \|/ Sun /|\ --O-- | /|\ | S2,/ / | T2/ / R1,T1| / |/_ R3 | / / | / / r1||/_|/_r3 Ground --------------------------- /////////////////////////// Run 1: Downwelling radiance and transmittance (R3,T2): set ANGLE to solar zenith angle. After run is complete copy TAPE12 to SOL.PATH.T2. Run 2: Upwelling radiance (R1,T1): set angle to observer angle. Make sure H1, H2, TBOUND, SREMISS and SRREFL are correct for this geometry. RECORD RECORD 1.3 (required if IHIRAC > 0; IAERSL > 0; IEMIT = 1; IATM = 1; or ILAS > 0; otherwise omit) V1, V2, SAMPLE, DVSET, ALFAL0, AVMASS, DPTMIN, DPTFAC, ILNFLG, DVOUT, NMOL_SCAL 1-10, 11-20, 21-30, 31-40, 41-50, 51-60, 61-70, 71-80, 85, 90-100, 105 E10.3, E10.3, E10.3, E10.3, E10.3, E10.3, E10.3, E10.3, 4X,I1, 5X,E10.3, 3x,I2 V1 beginning wavenumber value for the calculation (VLAS = V1 for ILAS = 1,2) V2 ending wavenumber value for the calculation ( (V2-V1) must be less than 2020 cm - 1 ) SAMPLE number of sample points per mean halfwidth (between 1. and 4.) (default = 4.) DVSET selected DV for the final monochromatic calculation if positive, must be within 20% of DV for final monochromatic calculation determined by LBLRTM If negative, final DV = abs(DVSET). However, this option causes line width readjustment which may effect the validity of the results. Use with extreme caution. ALFAL0 average collision broadened halfwidth (cm - 1/atm) (default = 0.04) AVMASS average molecular mass (amu) for Doppler halfwidth (default = 36) DPTMIN minimum molecular optical depth below which lines will be rejected (negative value defaults to DPTMIN = 0.0002) DPTFAC factor multiplying molecular continuum optical depth to determine optical depth below which lines will be rejected (negative value defaults to DPTFAC = 0.001) ILNFLG flag for binary record of line rejection information = 0 line rejection information not recorded (default) = 1 write line rejection information to REJ1, REJ4 = 2 read line rejection information from REJ1, REJ4 DVOUT selected DV grid for the optical depth "monochromatic" output spacing (must be less than or equal to default spacing or DVSET if nonzero). NMOL_SCAL Enables the scaling of the atmospheric profile for selected species NMOL_SCAL is the highest molecule number for which scaling will be applied See Record(s) 1.3.a/1.3.b.n RECORD RECORD 1.3.a (required if NMOL_SCAL > 0; otherwise omit) HMOL_SCAL(M=1,39) 1-39 39A1 HMOL_SCAL(M) for molecule M (column M in Record 1.3.a: HMOL_SCAL(M) XMOL_SCAL(M) from Record 1.3.b interpreted and used as follows: ' ' no scaling applied (overrides any value on Record 1.3.b) '0' scaling factor is zero (layer column amounts are set to 0. for all levels) '1' scaling factor used directly to scale profile 'c' or 'C' column amount to which the profile is to be scaled (molec/cm^2) 'd' or 'D' column amount in Dobson units to which the profile is to be scaled 'm' or 'M' volume mixing ratio (ppv) wrt dry air for the total column to which the profile will be scaled 'p' or 'P' value of Precipitable Water Vapor (cm) to which the profile will be scaled (water vapor only) The scaling factor is the same at all levels M Molecule number as indicated in the following Table: 1: H2O 2: CO2 3: O3 4: N2O 5: CO 6: CH4 7: O2 8: NO 9: SO2 10: NO2 11: NH3 12: HNO3 13: OH 14: HF 15: HCL 16: HBR 17: HI 18: CLO 19: OCS 20: H2CO 21: HOCL 22: N2 23: HCN 24: CH3CL 25: H2O2 26: C2H2 27: C2H6 28: PH3 29: COF2 30: SF6 31: H2S 32: HCOOH 33: HO2 34: O 35:CLONO2 36: NO+ 37: HOBR 38: C2H4 39: CH3OH RECORD 1.3.b.(1..n) (only if NMOL_SCAL > 0; otherwise omit) # records depends on NMOL_SCAL (XMOL_SCAL, M=1,mol_max) (8E15.7) XMOL_SCAL(M) value interpreted from HMOL_SCAL(M) to be used for scaling the profile of molecule M RECORD RECORD 1.4 (required if IEMIT = 1, or both IEMIT=2 and IOTFLG=2; otherwise omit) TBOUND, SREMIS(1), SREMIS(2), SREMIS(3), SRREFL(1), SRREFL(2), SRREFL(3), surf_refl 1-10, 11-20, 21-30, 31-40, 41-50, 51-60, 61-70, 75 E10.3, E10.3, E10.3, E10.3, E10.3, E10.3, E10.3 4X,1A temperature and emissivity parameters for boundary at H2 (end of path) (applies to downlooking, uplooking and tangent paths) for down looking case, the effect of reflected atmospheric radiance above lower boundary included for SRREFL(I) > 0. (reflected atmosphere from H1 to H2 included) TBOUND temperature of boundary (K) SREMIS(I) frequency dependent boundary emissivity coefficients (I = 1,2,3) EMISSIVITY = SREMIS(1) + SREMIS(2)*V + SREMIS(3)*(V**2) *** NOTE: Entering a value for SREMIS(1) < 0 allows for direct input of boundary emissivities from file 'EMISSIVITY' SRREFL(I) frequency dependent boundary reflectivity coefficients (I = 1,2,3) REFLECTIVITY = SRREFL(1) + SRREFL(2)*V + SRREFL(3)*(V**2) *** NOTE: Entering a value for SRREFL(1) < 0 allows for direct input of boundary reflectivities from file 'REFLECTIVITY' surf_refl specifies the surface type used in computing the reflected downward radiance 's' or ' ' - is for a specular surface 'l' - is for a lambertian surface Note: For the surf_refl = 'l' option If IATM=0 the appropriate angle is specified on Record 2.1; Otherwise the angle is determined from the geometry in lblatm (IATM>0) RECORD RECORD 1.4 (Continued) --------------------------------------------------------------------------------------------------- | FORMAT FOR 'EMISSION' and 'REFLECTION' FILES | | -------------------------------------------- | | | | Record 1: | | V1, V2, DV, NLIM | | | | 1-10, 11-20, 21-30, 35-40 | | | | E10.3, E10.3, E10.3, 5X,I5 | | | | | | V1 beginning wavenumber of input (should be less than V1 on RECORD 1.3 from TAPE5) | | | | V2 ending wavenumber of input (should be greater than V2 on RECORD 1.3 from TAPE5) | | | | DV spectral spacing of input (cannot be less than 0.5 cm-1) | | | | NLIM number of points in input (must be less than or equal to 4040) | | | | | | Record 2(i), i=1,NLIM: | | | | ZDATA | | | | 1-10 | | | | E10.3 | | | | | | ZDATA Spectrally dependent boundary emission or reflection value | | (one spectral value per card) | | | --------------------------------------------------------------------------------------------------- RECORD RECORD 1.5 (required for Analytic Jacobian calculation: IEMIT=3 and IMRG=40,41,42 or 43) NSPCRT 1-5 I5 NSPCRT (-1,0,1,...N) species number for radiance analytic Jacobian/derivative -1 = surface emissivity and surface skin temperature 0 = atmospheric temperature profile 1 -> Nmol = reference number in TABLE 1 for molecular species RECORD RECORD 1.6a (required if IMRG = 35-36, 40-41, 45-46; otherwise omit) NOTE: IMRG = 35-36, 40-41, 45-46 require separate optical depth files for each layer PTHODL, LAYTOT 1-55, 57-60 A55, 1X, I4 PTHODL pathname for precalculated optical depth files Example: For PTHODL="ODdeflt_", code will look for files ODdeflt_01, ... ------- LAYTOT total number of layers used in radiative transfer *** NOTE *** LAYTOT ignored for IEMIT = 3 ***************************************************************************** ****** these records applicable only if LBLATM not selected (IATM=0) ****** LAYER INPUT (MOLECULES ONLY) RECORD RECORD 2.1 IFORM, NLAYRS, NMOL, SECNTO, ZH1, ZH2, ZANGLE 2 3-5, 6-10, 11-20, 41-48, 53-60, 66-73 1X,I1 I3, I5, F10.2, 20X, F8.2, 4X, F8.2, 5X, F8.3 IFORM (0,1) column amount format flag = 0 read PAVE(L), WKL(M,L), WBROADL(L) in F10.4, E10.3, E10.3 formats (default) = 1 read PAVE(L), WKL(M,L), WBROADL(L) in E15.7 format NLAYRS number of layers (maximum of 200) NMOL value of highest molecule number used (default = 7; maximum of 35) See Table I for molecule numbers. SECNTO user entered scale factor for the column amount for the layers defined by NLAYRS if positive, looking up if negative, looking down normal value = 1.0 ZH1 observer altitude ZH2 end point altitude ZANGLE mean zenith angle for path calculation (degrees) NOTE: With the surf_refl = 'l' option on Record 1.4 ZANGLE must be set. (ZANGLE .gt.90 and .le.180 ) RECORD RECORD 2.1.1 PAVE(L), TAVE(L), SECNTK(L), ITYL(L), IPATH, ALTZ(L-1), PZ(L-1), TZ(L-1), ATLZ(L), PZ(L), TZ(L) 1-10, 11-20, 21-30, 31-33, 34-35, 37-43, 44-51, 52-58, 59-65, 66-73, 74-80 F10.4, F10.4, F10.4, A3, I2, 1X,F7.2, F8.3, F7.2, F7.2, F8.3, F7.2 IF IFORM=1 from record 2.1, then the following format is used: 1-15, 16-25, 26-35, 36-38, 39-40, 42-48, 49-56, 57-63, 64-70, 71-78, 79-85 E15.7, F10.4, F10.4, A3, I2, 1X,F7.2, F8.3, F7.2, F7.2, F8.3, F7.2 PAVE(L) average pressure of layer (millibars) (**If IFORM=1, then PAVE in E15.7 format**) TAVE(L) average temperature of layer (K) SECNTK(L) if SECNTK .ne. 0, it is the value of the column amount scale factor that will be used in the calculation for this layer a non-zero value for SECNTK will override the scale factor value, SECNTO . RECORD RECORD 2.1.1 (Continued) ITYL(L) overrides the LBLRTM internal calulation of ITYPE for the layer L. ITYPE controls the DV ratio of the previous layer to the current layer. A blank defaults to the internal calulation. = 0 specifies 1/1 1 specifies 2/1 2 specifies 3/2 4 specifies 5/4 5 specifies 6/5 99 specifies first layer (i.e. no ratio) IPATH = 0 determined by sign of SECNTO 1 looking down merge (upwelling radiance calculated) 2 tangent layer merge 3 looking up merge (downwelling radiance calculated) ALTZ(L-1) altitude for bottom of current layer (information only) PZ(L-1) pressure at ALTZ(L-1) (information only) TZ(L-1) temperature at ALTZ(L-1) - used by LBLRTM for Planck Function Calculation ** NOTE ** ALTZ(L-1), PZ(L-1) and TZ(L-1) are only required for the first layer. LBLRTM assumes that these quantites are equal to the top of the previous layer for L > 1. ALTZ(L) altitude for top of current layer (information only) PZ(L) pressure at ALTZ(L) (information only) TZ(L) temperature at ALTZ(L) - used by LBLRTM for Planck Function Calculation RECORD RECORD 2.1.2 (WKL(M,L), M=1, 7), WBROADL(L) (8E10.3) WKL(M,L) column densities (molecules/cm**2) or mixing ratios (ppv) for 7 molecular species WBROADL(L) column density for broadening gases (molecules/cm**2) **NOTE** If IFORM=1, then WKL(M,L) and WBROADL(L) are in 8E15.7 format RECORD 2.1.3 only if (NMOL .GT . 7) # records depends on NMOL (WKL(M,L), M=8, NMOL) (8E10.3) NMOL is set from LINFIL (TAPE3) (NMOL limited to 35 in LBLRTM) **NOTE: If IFORM=1 then WKL(M,L) in 8E15.7 format** REPEAT RECORDS 2.1.1 through 2.1.3 for the remaining layers (up to NLAYRS) ***************************************************************************** ****** these records applicable only if LBLATM not selected (IATM=0) ****** **************** and cross-sections ARE selected (IXSECT=1) *************** LAYER INPUT (CROSS-SECTIONS ONLY) RECORD RECORD 2.2 IXMOLS, IXSBIN 5, 15 I5, 5X, I5 IXMOLS number of cross-section molecules to be inputed (maximum of 35) IXSBIN (0,1) flag to deselect pressure convolution of cross-sections = 0 cross-sections convolved with pressure = 1 cross-sections not convolved with pressure RECORD RECORD 2.2.1 XSNAME(I), I=1, IXMOLS (7A10,(/,8A10)) XSNAME name of the cross-section molecules to be used (Note: the cross-sections specified here must also appear in the cross-section input file FSCDXS. See Table II after RECORD 3.7.1) RECORD RECORD 2.2.2 IFRMX, NLAYXS, IXMOL, SECNTX, HEDXS 2 3-5, 6-10, 11-20, 21-80 1X,I1 I5, I5, F10.2, 15A4 IFRMX (0,1) column amount format flag = 0 read PAVX, XAMNT(M,L), WBROADX(L) in F10.4, E10.3, E10.3 formats (default) = 1 read PAVX, XAMNT(M,L), WBROADX(L) in E15.7 format NLAYXS number of layers (maximum of 200), should be the same as NLAYRS on Record 2.1 IXMOL number of cross-sections (maximum of 35), should be the same as IXMOLS on Record 2.2 SECNTX user entered scale factor for the column amount for the layers defined by NLAYXS if positive, looking up if negative, looking down normal value = 1.0 *** Note: this quantity is for information only, SECNTX should be the same as SECNT0 on Record 2.1, which is used instead of this quantity. HEDXS 60 character user identification **************************************************************************************************************** **** Note: all the quantites on this Record, are for information only, LBLRTM will use the corresponding **** **** values from Record 2.1.1. This record is retained to maintain compatibility with 2.1.1 and **** **** is consistent with the output from LBLATM with IPU = 1; this provides the flexibility of **** **** independent utilization of this data with external programs. **** **************************************************************************************************************** RECORD RECORD 2.2.3 PAVX, TAVX, SECKXS, ITYX(L), IPATX, ALTX(L-1), PX(L-1), TX(L-1), ATLX(L), PX(L), TX(L) 1-10, 11-20, 21-30, 31-33, 34-35, 37-43, 44-51, 52-58, 59-65, 66-73, 74-80 F10.4, F10.4, F10.4, A3, I2, 1X,F7.2, F8.3, F7.2, F7.2, F8.3, F7.2 PAVX ave. pressure of layer (millibars) (**If IFRMX=1, then PAVX in E15.7 format**) TAVX ave. temperature of layer (K) SECKXS if SECKXS .ne. 0, it is the value of the column amount scale factor that will be used in the calculation for this layer a non-zero value for SECKXS will override the scale factor value, SECNTX . ITYX(L) overrides the LBLRTM internal calulation of ITYPE for the layer L. ITYPE controls the DV ratio of the previous layer to the current layer. A blank defaults to the internal calulation. = 0 specifies 1/1 1 specifies 2/1 2 specifies 3/2 4 specifies 5/4 5 specifies 6/5 99 specifies first layer (i.e. no ratio) IPATX = 0 determined by sign of SECNTX 1 looking down merge (upwelling radiance calculated) 2 tangent layer merge 3 looking up merge (downwelling radiance calculated) ALTX(L-1) altitude for bottom of current layer PX(L-1) pressure at ALTX(L-1) TX(L-1) temperature at ALTX(L-1) ** NOTE ** ALTX(L-1), PX(L-1) and TX(L-1) are only required for the first layer. LBLRTM assumes that these quantites are equal to the top of the previous layer for L > 1. ALTX(L) altitude for top of current layer PX(L) pressure at ALTZ(L) TX(L) temperature at ALTZ(L) RECORD RECORD 2.2.4 (XAMNT(I,L), I=1, 7), WBROADX(L) (8E10.3) XAMNT(I,L) column densities (molecules/cm**2) or mixing ratios (ppv) for 7 cross-section molecules WBROADX(L) column density for broadening gases (molecules/cm**2), ** information only ** should be the same as WBROADL, on Record 2.1.2 **NOTE** If IFRMX=1, then XAMNT(I,L) and WBROADX(L) are in 8E15.7 format RECORD 2.2.5 only if (IXMOL .GT . 7) # records depends on IXMOL (XAMNT(I,L), I=8, IXMOL) (8E10.3) **NOTE: If IFRMX=1 then XAMNT(I,L) in 8E15.7 format** REPEAT RECORDS 2.2.3 through 2.2.5 for the remaining layers (up to NLAYXS) ***************************************************************************** ******** these records applicable if LBLATM selected (IATM=1) ******** LBLATM RECORD RECORD 3.1 MODEL, ITYPE, IBMAX, ZERO, NOPRNT, NMOL, IPUNCH, IFXTYP, MUNITS, RE, HSPACE, VBAR, REF_LAT 5, 10, 15, 20, 25, 30, 35, 36-37, 39-40, 41-50, 51-60, 61-70, 81-90 I5, I5, I5, I5, I5, I5, I5, I2, 1X, I2, F10.3, F10.3, F10.3, 10x, F10.3 MODEL selects atmospheric profile = 0 user supplied atmospheric profile = 1 tropical model = 2 midlatitude summer model = 3 midlatitude winter model = 4 subarctic summer model = 5 subarctic winter model = 6 U.S. standard 1976 ITYPE selects type of path = 1 horizontal path (constant pressure, temperature), use RECORD 3.2H = 2 slant path from H1 to H2, use RECORD 3.2 = 3 slant path from H1 to space (see HSPACE), use RECORD 3.2 IBMAX selects layering for LBLRTM = 0 LBLRTM layers are generated internally (default) > 0 IBMAX is the number of layer boundaries read in on Record 3.3B which are used to define the layers used in LBLRTM calculation If IBMAX is set to a negative value, the layer boundaries are specified in PRESSURE (mbars). ZERO = 2 zeroes absorber amounts which are less than 0.1 percent of total NOPRNT = 0 full printout = 1 selects short printout NMOL number of molecular species (default = 7; maximum value is 35) IPUNCH = 0 layer data not written (default) = 1 layer data written to unit ITAPE7)PU (TAPE7) = 2 atmospheric data written to AJ_atmosphere and AJ_xs_amnts for Analytic Jacobian IFXTYP = 0 leave ITYL(L) and ITYX(L) fields blank on TAPE7 (if IPUNCH = 1, default) = 1 calculate and write ITYL(L) and ITYX(L) to TAPE7 (if IPUNCH = 1) = 2 calculate and write ITYL(L) and ITYX(L) to file "IFIXTYPE" (if IPUNCH = 1) = -2 use values of ITYL(L) and ITYX(L) from file "IFIXTYPE" NOTE: See RECORD 2.1.1 and RECORD 2.2.3 for ITYL(L) and ITYX(L) MUNITS = 0 write molecular column amounts to TAPE7 (if IPUNCH = 1, default) = 1 write molecular mixing ratios to TAPE7 (if IPUNCH = 1) RECORD RECORD 3.1 (Continued) RE radius of earth (km) defaults for RE=0: a) MODEL 0,2,3,6 RE = 6371.23 km b) 1 RE = 6378.39 km c) 4,5 RE = 6356.91 km HSPACE altitude definition for space (default = 100 km) internal models defined to 120 km VBAR frequency for refractive geometry calculation (default: VBAR = (V1+V2) / 2 ) (V1,V2 from Record 1.3) REF_LAT latitude of location of calculation (degrees) defaults for REF_LAT = 0: a) MODEL 0,2,3,6 REF_LAT = 45.0 degrees b) MODEL 1 REF_LAT = 15.0 c) MODEL 4,5 REF_LAT = 60.0 ----------------------------------------------------------------------------- ------------------ SLANT PATH (ITYPE = 2,3) (MODEL = 0-6) ------------------ RECORD RECORD 3.2 H1, H2, ANGLE, RANGE, BETA, LEN, HOBS 1-10, 11-20, 21-30, 31-40, 41-50, 51-55, 61-70 F10.3, F10.3, F10.3, F10.3, F10.3, I5, 5X,F10.3 H1 observer altitude (km) If IBMAX is negative, H1 is provided in pressure units (mbars) H2 for ITYPE = 2, H2 is the end point altitude (km) ITYPE = 3, H2 is the tangent height (km) for H2 .GT. 0. if H2 = 0. ANGLE determines tangent height If IBMAX is negative, H2 is provided in pressure units (mbars) ANGLE zenith angle at H1 (degrees) RANGE length of a straight path from H1 to H2 (km) BETA earth centered angle from H1 to H2 (degrees) LEN = 0 short path (default) = 1 long path through a tangent height LEN is only used for H1 > H2 (ANGLE > 90`) for ITYPE = 2, only 3 of the first 5 parameters are required to specify the path, e.g., H1, H2, ANGLE or H1, H2 and RANGE for ITYPE = 3, H1 = observer altitude must be specified. Either H2 = tangent height or ANGLE must be specified. Other parameters are ignored. HOBS Height of observer, used only for informational purposes in satellite-type simulations when computing output geometry above 120 km. RECORD 3.3 options RECORD RECORD 3.3A For IBMAX = 0 (from RECORD 3.1) AVTRAT, TDIFF1, TDIFF2, ALTD1, ALTD2 1-10, 11-20, 21-30, 31-40, 41-50 F10.3, F10.3, F10.3, F10.3, F10.3 AVTRAT maximum Voigt width ratio across a layer (if zero, default = 1.5) TDIFF1 maximum layer temperature difference at ALTD1 (if zero, default = 5 K) TDIFF2 maximum layer temperature difference at ALTD2 (if zero, default = 8 K) ALTD1 altitude of TDIFF1 (if zero, default = 0 Km) ALTD2 altitude of TDIFF2 (if zero, default = 100 Km) RECORD RECORD 3.3B For IBMAX > 0 (from RECORD 3.1) ZBND(I), I=1, IBMAX altitudes of LBLRTM layer boundaries (8F10.3) If IBMAX < 0 PBND(I), I=1, ABS(IBMAX) pressures of LBLRTM layer boundaries (8F10.3) -------------------------------------------------------------------------------- ----------------------------------------------------------------------------- User Defined Atmospheric Profile -------------------------------- (MODEL = 0) -------------------------------- RECORD RECORD 3.4 IMMAX, HMOD 5, 6-29 I5, 3A8 IMMAX number of atmospheric profile boundaries If IMMAX is set to a negative value, the level boundaries are specified in PRESSURE (mbars). HMOD 24 character description of profile RECORD RECORD 3.5 ZM, PM, TM, JCHARP, JCHART, JLONG, (JCHAR(M),M =1,39) 1-10, 11-20, 21-30, 36, 37, 39, 41 through 80 E10.3, E10.3, E10.3, 5x, A1, A1, 1x, A1, 1x, 39A1 ZM boundary altitude (km). If IMMAX < 0, altitude levels are computed from pressure levels PM. If any altitude levels are provided, they are ignored if IMMAX < 0 (exception: The first input level must have an accompanying ZM for input into the hydrostatic equation) PM pressure (units and input options set by JCHARP) TM temperature (units and input options set by JCHART) JCHARP flag for units and input options for pressure (see Table I) JCHART flag for units and input options for temperature (see Table I) JLONG flag for reading long record for molecular information = L read VMOL(M) in 8E15.8 format JCHAR(K) flag for units and input options for the K'th molecule (see Table I) RECORD RECORD 3.6.1 ... 3.6.N VMOL(M), M=1, NMOL 8E10.3 VMOL(M) density of the M'th molecule in units set by JCHAR(K) **NOTE** If JLONG=L, then VMOL(M) is in 8E15.8 format REPEAT records 3.5 and 3.6.1 to 3.6.N for each of the remaining IMMAX boundaries ---------------------------------------------------------------------------- **************************************************************************** TABLE I. Units and input options for the K'th molecule TABLE I USER OPTIONS FOR PRESSURE, TEMPERATURE, AND MOLECULAR DENSITY JCHARP PRESSURE 1-6 default to value for specified model atmosphere (JCHARP) " ",A pressure in (mb) B " " (atm) C " " (torr) JCHART TEMPERATURE 1-6 default to value for specified model atmosphere (JCHART) " ",A ambient temperature in deg (K) B " " " " " " (C) JCHAR(M) (M): AVAILABLE MOLECULAR SPECIES: 1: H2O 2: CO2 3: O3 4: N2O 5: CO 6: CH4 7: O2 8: NO 9: SO2 10: NO2 11: NH3 12: HNO3 13: OH 14: HF 15: HCL 16: HBR 17: HI 18: CLO 19: OCS 20: H2CO 21: HOCL 22: N2 23: HCN 24: CH3CL 25: H2O2 26: C2H2 27: C2H6 28: PH3 29: COF2 30: SF6 31: H2S 32: HCOOH 33: HO2 34: O 35:CLONO2 36: NO+ 37: HOBR 38: C2H4 39: CH3OH potential choice of units for above species: JCHAR = 1-6 - default to value for specified model atmosphere = " ",A - volume mixing ratio (ppmv): = B - number density (cm-3) = C - mass mixing ratio (gm/kg) = D - mass density (gm m-3) = E - partial pressure (mb) = F - dew point temp (K) *H2O only* = G - dew point temp (C) *H2O only* = H - relative humidity (percent) *H2O only* = I - available for user definition JCHAR must be less than "J" ---------------------------------------------------------------------------- ---------------- HORIZONTAL PATH (ITYPE = 1) (MODEL = 0-6) --------------- RECORD RECORD 3.2H H1, , , RANGEF 1-10, 31-40 F10.3, 10X, 10X, F10.3 H1 = observer altitude (km) RANGEF = path length (km) RECORD RECORD 3.3 not used for this case RECORD RECORD 3.4 records included for MODEL = 0 only RECORD 3.5 (see MODEL = 0 input records) RECORD 3.6 ----------------------------------------------------------------------------- ----------------------------------------------------------------------------- ***************************************************************************** ********** these records applicable if LBLATM selected (IATM=1) *********** **************** and cross-sections ARE selected (IXSECT=1) *************** XAMNTS RECORD RECORD 3.7 IXMOLS, IPRFL, IXSBIN 5, 10, 15 I5, I5, I5 IXMOLS number of cross-section molecules to be inputed (maximum of 35) IPRFL (0,1) = 0 user input profile = 1 standard profile as determined by LBLATM IXSBIN (0,1) flag to deselect pressure convolution of cross-sections = 0 cross-sections convolved with pressure (default) = 1 cross-sections not convolved with pressure RECORD RECORD 3.7.1 XSNAME(I), I=1, IXMOLS (7A10,(/,8A10)) XSNAME name of the cross-section molecules to be used (Note: the cross-sections specified here must also appear in the cross-section input file FSCDXS. See Table II) TABLE II. Structure of file FSCDXS. TABLE II Below is a sample file for FSCDXS. The format for the fields of file FSCDXS is given below: XNAME, V1X, V2X, DVX, NTEMP, IFRM, CFRM, (XFILS(I),I=1,NTEMP) FORMAT (A10,2F10.4,F10.7,I5,5X,I5,A1,4X,6A10) Sample file FSCDXS: ------------------------------------------------------------------------------------------------------------------------ MOLECULE V1 V2 DV NTEMP FORMAT FILET1 FILET2 FILET3 FILET4 FILET5 FILET6 CCL4 786.0010 805.9980 .01095126 1 86N CCL4 F11 830.0090 859.9990 .01095325 4 90N F11AT1 F11AT2 F11AT3 F11AT4 F11 1060.0040 1099.9980 .01095726 1 86N F11BT1 F12 860.0080 939.9960 .01095576 6 90 F12AT1 F12AT2 F12AT3 F12AT4 F12AT5 F12AT6 F12 1070.0050 1179.9940 .01095726 4 90N F12BT1 F12BT2 F12BT3 F12BT4 %%% ------------------------------------------------------------------------------------------------------------------------ XNAME 10 character cross-section name, left-justified, all CAPS, which matches an alias from the following table: Cross-Section Molecules and Aliases ------------------------------------------------- Alias(1) Alias(2) Alias(3) Alias(4) ---------- ---------- ---------- ---------- CLONO2 CLNO3 HNO4 CHCL2F CFC21 F21 CCL4 CCL3F CFCL3 CFC11 F11 CCL2F2 CF2CL2 CFC12 F12 C2CL2F4 C2F4CL2 CFC114 F114 C2CL3F3 C2F3CL3 CFC113 F113 N2O5 HNO3 CF4 CFC14 F14 CHCLF2 CHF2CL CFC22 F22 CCLF3 CFC13 F13 C2CLF5 CFC115 F115 NO2 PAN ACET CH3COCH3 ACETONE CH3C(O)CH3 CH3CN Note: If a molecule has cross-section data in more than one spectral region, each region must have a separate entry in FSCDXS. A '%' in column 1 will terminate the read of FSCDXS TABLE II. Structure of file FSCDXS. (Continued) V1X Beginning wavenumber for the cross-section data (cm-1) V2X Ending wavenumber for the cross-section data (cm-1) DVX Wavenumber spacing for the cross-section data (cm-1) ** for information only ** NTEMP number of temperature dependent cross-section files to be read in. '1' indicates no temperature dependence. Maximum of '6' files (i.e. 6 temperatures) is allowed. ***** no default ***** IFRM is the format specifier = 86 - HITRAN 86 format = 90 - HITRAN 90 format ** default ** Note: the headers are different between HITRAN 86 and HITRAN 90 CFRM is the blocking specifier = 'N' - unblocked data (100 characters/record) = ' ' - blocked data (51*100 characters/record) ** default ** XFILS(I) 10 character cross-section file name (may include directory name if length is adequate) ----------------------------------------------------------------------------- User Defined Atmospheric Profile -------------------------------- (IPRFL = 0) -------------------------------- RECORD RECORD 3.8 LAYX, IZORP, XTITLE 5, 10, 11-60 I5, I5 A50 LAYX number of atmospheric profile boundaries IZORP (0,1) flag which determines value of ZORP on Record 3.8.1 = 0 ZORP is an altitude in KM = 1 ZORP is a pressure in millibars XTITLE 50 character description of profile RECORD RECORD 3.8.1 ZORP, (JCHAR(K),K =1,28) 1-10, 16 through 50 F10.3, 5X, 35A1 ZORP boundary altitude (km) or pressure (millibars) as determined by IZORP on Record 3.8 JCHAR(K) flag for units and input options for the K'th cross-section JCHAR = 1 - default to value for global model atmosphere = " ",A - volume mixing ratio (ppmv) RECORD RECORD 3.8.2 ... 3.8.N DENX(K), K=1, IXMOLS 8E10.3 DENX(K) density of the K'th cross-section in units set by JCHAR(K) REPEAT records 3.8.1 to 3.8.N for each of the remaining LAYX boundaries ---------------------------------------------------------------------------- ****************************************************************************************** Spectral Absorber (IAERSL = 5) ------------------------------------------------------------------------------------------ ****************************************************************************************** Spectral Absorption by Layer: This option exists to enable the user to provide slowly varying optical depths by layer This absorption may be due to clouds or aerosols Data must be on file 'in_lblrtm_cld' All Abs formats are free form. space delimited. RECORD RECORD Abs.1 n_freq number of frequency points RECORD RECORD Abs.2 v_cloud_freq(l) n_freq frequencies at which layer spectral optical depths are provided RECORD RECORD Abs.3 n_layer number of layers for which the spectral optical depths are provided must correlate with the layers being used in the radiative transfer RECORD RECORD Abs.4.1 ilay layer index press layer pressure (mbar) (for information only) RECORD RECORD Abs.4.2 cloudodlayer(n_freq) n_layer records of optical depths with n_freq values in each record Repeat Records Abs.4.1 and Abs.4.2 for each of the remaining n_layer layers. ****************************************************************************************** AEROSOLS (IAERSL = 1,7) ------------------------------------------------------------------------------------------ Subroutine LOWTRN Stripped down version of LOWTRAN 7 run as a subroutine to supply LBLRTM with the attenuation due to aerosols, clouds, fogs, and rain. Further details of most of the aerosol parameters are contained in the LOWTRAN7 User Guide (AFGL-TR-88-0177), the LOWTRAN6 Report (AFGL-TR-83-0187) the LOWTRAN5 Report (AFGL-TR-80-0167) and the Millimeter Aerosol Report (AFGL-TR-79-0253) **************************** WARNING TO USERS! **************************** The LOWTRAN options have not been updated or tested in recent versions of LBLRTM. Support is no longer available from AER for these options. Userswho wish to model an absorber with slowly varying spectral dependence have the option of seting IAERSL=5 and inputting spectrally dependent optical depths at each level. ****************************************************************************************** Activated by IAERSL = 1 or 7 (on record 1.2) ** Note: 1) CTHIK and CALT should be place on specified ** ** model levels (ZBND or PBND from Record 3.3B) ** ** 2) MODEL is read in LBLATM ** ** 3) M1, M2, and M3 are set depending on MODEL ** RECORD sequence as follows 4.1 IHAZE,ISEASN,IVULCN,ICSTL,ICLD,IVSA,VIS,WSS,WHH,RAINRT,GNDALT Format (6I5,5F10.3) ***Optional RECORDS*** (if ICLD = 18,19, or 20) RECORD 4.2 CTHIK,CALT,CEXT,ISEED Format (3F10.3,I10) (if IVSA = 1) RECORD 4.3 ZCVSA,ZTVSA,ZINVSA Format (3F10.3) (if IAERSL = 7) RECORD 4.4 ML,TITLE Format (I5,18A4) ** RECORD 4.5 is repeated ML times ** RECORD 4.5 ZMDL,AHAZE,EQLWCZ,RRATZ,IHA1,ICLD1,IVUL1,ISEA1,ICHR1 Format (4F10.3,5I5) (if IHAZE = 7 RECORD 4.6.1 (IREG(I),I=1,4) or ICLD = 11) Format (4I5) RECORD 4.6.2 AWCCON(N),TITLE(N) Format (E10.3,18A4) ** Note: N = IREG(1)+IREG(2)+IREG(3)+IREG(4) RECORD 4.6.3 (VX(I),EXTC(N,I),ABSC(N,I),ASYM(N,I),I=1,47) Format (3(F6.2,2F7.5,F6.4)) ** 16 records ** ****************************************************************************************** RECORD 4.1 IHAZE,ISEASN,IVULCN,ICSTL,ICLD,IVSA,VIS,WSS,WHH,RAINRT,GNDALT FORMAT (6I5,5F10.3) IHAZE selects the type of extinction and a default meteorological range for the boundary-layer aerosol model (0 to 2 km altitude) IHAZE = 0 no aerosol attenuation included in calculation = 1 Rural extinction, (23 km VIS. default profile) = 2 Rural extinction, (5 km VIS. default profile) = 3 Navy Maritime extinction, sets own VIS. = 4 Maritime extinction, 23 km VIS. (LOWTRAN 5 Model) = 5 Urban extinction,(5 km VIS. default profile) = 6 Tropospheric extinction, (50 km VIS. default profile) = 7 User defined aerosol extinction coefficients see Records 4.6.1 - 4.6.3 = 8 FOG1 (advection fog) extinction, 0.2 km VIS. = 9 FOG2 (radiation fog) extinction, 0.5 km VIS. = 10 DESERT extinction sets own visibility from wind speed (WSS). If VIS is also specified on Record 4.1 it will override the default IHAZE value of VIS ISEASN selects the seasonal dependence of the vertical profiles for both the tropospheric (2 to 10 km) and stratospheric (10 to 30 km) aerosols. ISEASN = 0 Defaults to season of MODEL (see RECORD 3.1) (MODEL 0,1,2,4,6,7) summer (MODEL 3,5) winter = 1 spring-summer = 2 fall-winter IVULCN selects both the profile and extinction type for the stratospheric aerosols and determines transition profiles above the stratosphere to 100 km. IVULCN = 0 default to stratospheric background = 1 stratospheric background = 2 aged volcanic extinction; moderate volcanic profile = 3 fresh volcanic extinction; high volcanic profile = 4 aged volcanic extinction; high volcanic profile = 5 fresh volcanic extinction; moderate volcanic profile = 6 background stratospheric extinction; moderate volcanic profile = 7 background stratospheric extinction; high volcanic profile = 8 fresh volcanic extinction; extreme volcanic profile RECORD RECORD 4.1 (Continued) ICSTL is the air mass character (1 to 10); only used with Navy maritime model (IHAZE=3) (default value = 3) ICSTL = 1 open ocean . . . 10 strong continental influence ICLD determines the inclusion of cirrus cloud attenuation or gives a choice of five cloud models and 5 rain models ICLD for cloud and or rain ICLD = 0 no clouds or rain = 1 Cumulus cloud; base .66 km; top 3.0 km = 2 Altostratus cloud; base 2.4 km; top 3.0 km = 3 Stratus cloud; base .33 km; top 1.0 km = 4 Stratus; Strato Cumulus; base .66 km; top 2.0 km = 5 Nimbostratus cloud; base .16 km; top .66 km = 6 2.0 mm/hr drizzle (modeled with cloud 3) rain 2.0 mm/hr at 0.0 km to 0.22 mm/hr at 1.5 km = 7 5.0 mm/hr light rain (modeled with cloud 5) rain 5.0 mm/hr at 0.0 km to 0.2 mm/hr at 2.0 km = 8 12.5 mm/hr moderate rain (modeled with cloud 5) rain 12.5 mm/hr at 0.0 km to 0.2 mm/hr at 2.0 km = 9 25.0 mm/hr heavy rain (modeled with cloud 1) rain 25.0 mm/hr at 0.0 km to 0.2 mm/hr at 3.0 km =10 75.0 mm/hr extreme rain (modeled with cloud 1) rain 75.0 mm/hr at 0.0 km to 0.2 mm/hr at 3.5 km =11 read in user defined cloud extinction and absorption can only be used with IAERSL = 7 =18 Standard Cirrus Model =19 Sub-visual Cirrus Model =20 NOAA Cirrus model (LOWTRAN 6 Model) RECORD RECORD 4.1 (Continued) IVSA determines the use of the Army Vertical Structure Algorithm for aerosols in the boundary layer. IVSA = 0 not used = 1 vertical structure algorithm VIS = meterological range (km) (when specified, supersedes default value set by IHAZE) WSS = current wind speed (m/s) only with (IHAZE=3, or IHAZE=10) WHH = 24 hour average wind speed (m/s) only with (IHAZE=3) *** Note: For the Navy Maritime Model (IHAZE=3), if WSS=WHH=0.0, default wind speeds are set as follows: WSS and WHH MODEL Model Atmosphere Default Wind Speed (m/s) 0 User-Defined (Horiz Path) 6.9 1 Tropical 4.1 2 Midlatitude Summer 4.1 3 Midlatitude Winter 10.29 4 Subarctic Summer 6.69 5 Subarctic Winter 12.35 6 U.S. Standard 7.2 7 User-Defined 6.9 For the DESERT aerosol model (IHAZE=10), if WSS<0, the default wind speed is 10.0 m/s. RAINRT = rain rate (mm/hr) Default value is zero. Used to top of cloud when cloud is present, when no clouds, rain rate used to 6 km. GNDALT = altitude of surface relative to sea level (km) used to modify aerosol profiles below 6 km altitude ************************************************************************************ Optional input records after Record 4.1 Selected by parameters ICLD, IVSA, and IHAZE on Record 4.1 ********************************************************************************** RECORD RECORD 4.2 CTHIK,CALT,CEXT,ISEED (if ICLD = 18, 19, or 20) FORMAT (3F10.3,I10) Input record for cirrus altitude profile subroutine when ICLD = 18, 19, or 20 CTHIK = cirrus thickness (km) 0 use thickness statistics > 0 user defined thickness CALT = cirrus base altitude (km) 0 use calculated value > 0 user defined base altitude CEXT = extinction coefficient (km-1) at 0.55 microns 0 use 0.14*CTHIK > 0 user defined extinction coefficient ISEED = random number initialization flag 0 use default mean values for cirrus > 0 initial value of SEED for RANDM function (option not implemented in current version) ********************************************************************************** RECORD RECORD 4.3 ZCVSA, ZTVSA, ZINVSA (if IVSA = 1) FORMAT (3F10.3) Input record for Army vertical structure algorithm subroutine when IVSA = 1 ZCVSA = cloud ceiling height (km) ZCVSA > 0 known cloud ceiling ZCVSA = 0 unknown cloud ceiling height program calculates cloud height (default = 1.8 km) ZCVSA < 0 no cloud ceiling ZTVSA = thickness of cloud or fog (km) = 0 defaults to 0.2 km ZINVSA = height of the inversion layer (km) = 0 defaults to 2 km (0.2 km for fog) < 0 no inversion layer ********************************************************************************** ********************************************************************************** RECORD RECORD 4.4 ML,TITLE (if IAERSL = 7) FORMAT (I5,18A4) Additional aerosol profile ML = number of aerosol profiles levels to be inserted (maximum of 34) TITLE = identification of new model aerosol profile ** Note: RECORD 4.5 is repeated ML times ** RECORD 4.5 (if IAERSL = 7) ZMDL,AHAZE,EQLWCZ,RRATZ,IHA1,ICLD1,IVUL1,ISEA1,ICHR1 FORMAT (4F10.3,5I5) ZMDL = Altitude of layer boundary (km) AHAZE = aerosol visible extinction coef (km-1) at altitude ZMDL EQLWCZ = equivalent liquid water content (gm m-3) at alt ZMDL for the aerosol, cloud or fog models **** either AHAZE or EQLWCZ is allowed **** RRATZ = rain rate (mm/hr) at altitude ZMDL IHA1 = aerosol model used for spectral dependence of extinction IVUL1 = stratospheric aerosol model used for spectral dependence of extinction at altitude ZMDL ICLD1 = cloud model to be used for spectral dependence of extinction at altitude Z **** only one of IHA1, ICLD1 or IVUL1 is allowed **** IHA1.NE.0 others ignored IHA1.EQ.0 and ICLD1.NE.0 use ICLD1 **** if AHAZE and EQLWCZ are both zero, default profile loaded **** according to IHAZ1,ICLD1,IVUL1 ISEA1 = aerosol season control for the altitude ZMDL (only used with IHA1 or IVUL1.NE.0) ICHR1 = indicates a boundary change between two or more adjacent user defined aerosol or cloud regions at altitude ZMDL (required for IHAZE=7 or ICLD=11) ** Note: defaults to 0 for IHAZE.ne.7 or ICLD.ne.11 ** = 0 no boundary change = 1 signifies boundary change *********************************************************************************** RECORDS 4.6.1 - 4.6.3 read in the user defined cloud extinction and absorption (IHAZE = 7 or ICLD = 11) RECORD RECORD 4.6.1 (IREG(I),I=1,4) FORMAT (4I5) IREG = specifies which of the four altitude regions a user defined aerosol or cloud model will use ** Note: regions default to 1) 0 - 2 km 2) 3 - 10 km 3) 11 - 30 km 4) 35 - 100 km but can be overridden with 'IHA1' settings on Record 4.5 (IAERSL=7) ***** The following records (4.6.2 - 4.6.3) are repeated N times where N has the following value: ***** ***** ***** ***** N = IREG(1) + IREG(2) + IREG(3) + IREG(4) ***** ***** ***** RECORD RECORD 4.6.2 AWCCON(N),TITLE(N) FORMAT (E10.3,18A4) AWCCON(N) = conversion factor from equivalent liquid water content (gm/m3) to extinction coefficient (km-1). It is numerically equal to the equivalent liquid water content corresponding to an extinction coefficient of 1.0 km-1, at a wavelength of 0.55 microns. AWCCON has units of (km-gm/m3) TITLE(N) = for an aerosol or cloud region (up to 72 characters) RECORD RECORD 4.6.3 (VX(I),EXTC(N,I),ABSC(N,I),ASYM(N,I),I=1,47) FORMAT (3(F6.2,2F7.5,F6.4)) *** Note: there are 16 records per region VX(I) = wavelength of aerosol coefficient (not used by program but must be the wavelengths defined in array VX2 in subroutine EXTDTA) EXTC(N,I) = aerosol extinction coefficient, normalized so that EXTC for a wavelength of 0.55 microns (I=4) is 1.0 km-1. ABSC(N,I) = aerosol absorption coefficient, normalized so that EXTC for a wavelength of 0.55 microns (I=4) is 1.0 km-1. ASYM(N,I) = aerosol asymmetry factor ************************************************************************************ LASER (ILAS > 2) RECORD RECORD 5.1 ... 5.N VLAS 1-20 F20.8 VLAS laser wavenumber value (cm-1) (negative value terminates multiple LASER option) REPEAT RECORD 5.1 SCNMRG ( for scanned sequential results; IMRG between 13 and 18; 35-36 ) RECORD RECORD 6 HWHM, V1, V2, JEMIT, JFN, JVAR, SAMPL, NNFILE, NPTS 1-10, 11-20, 21-30, 34-35, 39-40, 44-45, 46-55, 71-75, 76-80 10.3, F10.3, F10.3, 3X,I2, 3X,I2, 3X,I2, F10.4, 15X, I5, I5 HWHM (Half Width Half Maximum) Note: HWHM is first zero crossing of periodic functions for JFN < 0. HWHM is redefined as HWHM=(FIRST ZERO)/(PI/SCALE) V1 beginning wavenumber value for performing SCAN V2 ending wavenumber value for performing SCAN JEMIT = 0 SCAN convolved with transmission = 1 SCAN convolved with radiance JFN selects choice of scanning function JFN Function Half-width Bound Default Sample --- -------- ---------------- -------------- = 0 rectangular 1.0 0.5 = 1 triangular 2.0 2.0 = 2 gaussian 4.0 4.0 = 3 sinc squared 54.1826 4.0 = 4 sinc 119.332818 4.0 Note: For each function BOUND = HWHM*(Half-width Bound). BOUND represents the overlap required before V1 and after V2 in the monochromatic calculation in order to obtain a correct result from the scanning function. ( see Record 1.3 ) ------------------------------------------------------------------------------------------- Note: If JFN < 0, HWHM contains the value for the first zero crossing of periodic function. JVAR flag for variable HWHM = 0 no variation = 1 HWHM(vi) = HWHM(v1) * (vi / v1) SCNMRG (Continued) RECORD RECORD 6 (Continued) SAMPL number of sample points per half width = 0 gives default value for each function < 0 this variable specifies the output spectral spacing (DELVO cm-1) The value of SAMPL is calculated internally as SAMPL = HWHM/DELVO NNFILE unit number for scanned sequential output defaults to NFILE (= 13) or previous value of NNFILE if doing multiple LBLRTM runs NPTS number of values to be printed for the beginning and ending of each panel for current scanned file ***** Note: if the DV of the data to be scanned is larger than DVINT = HWHM/12., ***** ***** the data is first interpolated to DVINT before it is scanned. ***** ***** Four-point interpolation is used for this procedure. ***** RECORD RECORD 6.1 (required for IMRG = 35, 36, 45, or 46; otherwise omit) DIRCOS 1-10 F10.8 DIRCOS direction cosine of radiance computation for external calculation of fluxes from Gaussian quadrature summation over one or more angles. FLTMRG ( for sequential results with filter; IMRG between 23 and 28 ) RECORD RECORD 7.1 V1F, DVF, NPTS, JEMIT, NNFILE, HEDDR 1-10, 11-20, 21-25, 26-30, 31-35, 46-80 F10.3, F10.4, I5, I5, I5, 10X, 8A4,A3 V1F wavenumber of initial filter value DVF wavenumber increment between filter values NPTS number of filter values (if NPTS.LT.0, previous values of XF utilized) JEMIT = 0 transmittance integrated over filter function = 1 radiance integrated over filter function NNFILE unit number for filtered sequential intermediate output defaults to NFILE (= 13) or previous value of NNFILE if doing multiple LBLRTM runs HEDDR 35 characters of user identification RECORD RECORD 7.2 IVAR 1-80 20A4 IVAR is format specification for reading filter values, XF RECORD RECORD 7.3 . . . . . . . . . 7.N XF format specification given by record 7.2 XF NPTS values of FILTRFN SCANFN (ISCAN = 1) RECORD RECORD 8.1 HWHM, V1, V2, JEMIT, JFN, JVAR, SAMPL, IUNIT, IFILST, NIFILS, JUNIT, NPTS, param 1-10, 11-20, 21-30, 34-35, 39-40, 44-45, 46-55, 59-60, 64-65, 69-70, 74-75, 76-80, 81-90 F10.3, F10.3, F10.3, 3X,I2, 3X,I2, 3X,I2, F10.4, 3X,I2, 3X,I2, 3X,I2, 3X,I2, I5, f10.5 HWHM (Half Width Half Maximum) negative value terminates SCANFN option Notes: 1. HWHM is first zero crossing of periodic functions for JFN < 0. HWHM is redefined as HWHM=(FIRST ZERO)/(PI/SCALE) 2. HWHM is instrument field of view half angle (in degrees) for JFN=5,6 (e.g., for FOV of 10.0 degrees, HWHM=5.0) V1 beginning wavenumber value for performing SCAN V2 ending wavenumber value for performing SCAN JEMIT = -1 SCANFN convolved with absorption (1.0 - transmission) = 0 SCANFN convolved with transmission = 1 SCANFN convolved with radiance JFN selects choice of scanning function JFN Function HWHM Default HWHM Bound Default Sample --- -------- ------ ------------------ -------------- = 0 rectangular 1.0 0.5 = 1 triangular 2.0 2.0 = 2 gaussian param 4.0 4.0 = 3 sinc squared 54.1826 4.0 = 4 sinc 119.332818 4.0 = 5 FOV correction 1.0 0.5 (uses center of box) = 6 FOV correction 1.0 0.5 (uses left edge of box) SCANFN (Continued) RECORD RECORD 8.1 (Continued) Notes: 1. For each function BOUND = HWHM*(Half-width Bound). BOUND represents the overlap required before V1 and after V2 in the monochromatic calculation. (Record 1.3) 2. If JFN < 0, HWHM contains the value for the first zero crossing of periodic function. 3. JFN cannot hold the values of -5 or -6. 4. JFN = 5,6 cannot be used with IMRG = 13-18,35,36 on Record 1.2. ------------------------------------------------------------------------------------------- JVAR flag for variable HWHM = 0 no variation = 1 HWHM(vi) = HWHM(v1) * (vi / v1) SAMPL number of sample points output per half width = 0 gives default value for each function < 0 this variable specifies the output spectral spacing (DELVO cm-1) The value of SAMPL is calculated internally as SAMPL = HWHM/DELVO Note: SAMPL must be < 0 and carry the value of the output spectral spacing for JFN=5,6. IUNIT unit designation of file to be scanned (default is MFILE) IFILST initial file from IUNIT to be scanned NIFILS number of files to be scanned starting at IFILST JUNIT file containing scanned results (default is JUNIT, file 11) NPTS number of values to be printed for the beginning and ending of each panel for current scanned file param additional parameter for selected functions, e.g. HWHM for JFN=2 ***** Note: if the DV of the data to be scanned is larger than DVINT = HWHM/12., ***** ***** the data is first interpolated to DVINT before it is scanned. ***** ***** Interpolated data is written to named file 'SCNINTF'. ***** ***** Four-point interpolation is used for this procedure. ***** REPEAT RECORD 8.1 A '-1.' within columns 1-10 will terminate scanning. INTRPL (ISCAN = 2) RECORD RECORD 9.1 DVO, V1, V2, JEMIT, I4PT, IUNIT, IFILST, NIFILS, JUNIT, NPTS 1-10, 11-20, 21-30, 31-35, 36-40, 56-60, 61-65, 66-70, 71-75, 76-80 F10.3, F10.3, F10.3, I5 I5, 15X,I5, I5, I5, I5, I5 DVO wavenumber spacing for interpolated result negative value terminates INTRPL option V1 beginning wavenumber value for performing INTRPL V2 ending wavenumber value for performing INTRPL JEMIT = -1 interpolation of absorption (1 - transmission) = 0 interpolation of transmission = 1 interpolation of radiance I4PT = 0 linear interpolation of data points = 1 four-point interpolation of data points IUNIT unit designation of file to be interpolated (default is MFILE) IFILST initial file from IUNIT to be interpolated NIFILS number of files to be interpolated starting at IFILST JUNIT file containing interpolated results (default is JUNIT, file 11) NPTS number of values to be printed for the beginning and ending of each panel for current interpolated file REPEAT RECORD 9.1 A '-1.' within columns 1-10 will terminate interpolation. FFTSCN (ISCAN = 3) RECORD RECORD 10.1 HWHM, V1, V2, JEMIT, JFNin, MRATin, DVOUT, IUNIT, IFILST, NIFILS, JUNIT, IVX, NOFIX 1-10, 11-20, 21-30, 31-35, 36-40, 41-45, 46-55, 56-60, 61-65, 66-70, 71-75, 76-78, 79-80 F10.3, F10.3, F10.3, I5, I5, I5, F10.3, I5, I5, I5, I5, I3, I2 HWHM (Half Width Half Maximum) negative value terminates FFTSCN option Note: HWHM is the maximum optical path difference of an equivalent interferometer for JFNin < 0. V1 beginning wavenumber value for performing FFTSCN V2 ending wavenumber value for performing FFTSCN JEMIT = 0 convolve with transmittance = 1 convolve with radiance FFTSCN (Continued) RECORD RECORD 10.1 (Continued) JFNin selects choice of scanning function JFNin Scanning Function Apodization Function a/HWHM a/FZ CR ----- ----------------- -------------------- ------ ---- -- = 0 boxcar = 1 1-v/a, |v|a (sinc squared) (triangle) = 2 exp(-0.5*(v/a)**2) exp(-2*pi*(a*x)**2) 0.849322 (NA) 10 (gauss) (gauss) = 3 (sin(pi*x*a)/pi*x*a))**2 1-x*a, |x|=<1/a 2.257609 1.0 40 (sinc squared) 0, |x|>1/a (triangle) = 4 sin(u)/u 1, |x|=<1/a 3.314800 2.0 160 (sinc) 0, |x|>1/a (rectangle) = 5 J(5/2,u)/(u**(5/2)) (1-(x*a)**2)**2 2.100669 0.91728 20 (Beer) = 6 sinc(u)+c1*(sinc(u+pi)+ (1+2*c1*cos(pi*x*a))/ 2.195676 1.0 20 sinc(u-pi)) (1+2*c1) (Hamming) = 7 sinc(u)+0.5*(sinc(u+pi)+ (1+cos(pi*x*a))/2 2.0 1.0 20 sinc(u-pi)) (Hanning) NORTON-BEER FUNCTIONS: _____________________ sum{Ci*(1-(x*a)**2)**i} from i=0 to 4, for 0 =< x =< 1/a and 0 for x > 1/a = 8 weak Ci = 0.384093,-0.087577, 2.57027 40 0.703484, 0., 0. = 9 moderate Ci = 0.152442,-0.136176, 2.36771 40 0.983734, 0., 0. = 10 strong Ci = 0.045335, 0., 2.07176 20 0.554883,0.,0.399782 FFTSCN (Continued) RECORD RECORD 10.1 (Continued) JFNin Scanning Function Apodization Function a/HWHM a/FZ CR ----- ----------------- -------------------- ------ ---- -- OTHER FUNCTIONS: _______________ = 11 NA 1 for x < p variable 100 (Brault) and (1+cos**2(u'))/2, for p/a < x < 1/a where (0 < p < 1) = 12 NA Io(pi*p*sqrt(1-(x*a)**2)/ variable 10 (Kaiser-Bessel) Io(pi*p), for 0 =< x =< 1/a and 0 for x > 1/a where (2 < p < 4) = 13 c1*sinc(u)+ NA variable 160 c2*sinc(u-2*Pi* v_offset/a) (Kiruna, assymetric) _______________________________________________________________________________________________________________________ LEGEND: - v = frequency (cm-1) - x = optical path difference (cm) - u = (2*pi*v)/a - u' = pi*(x*a-p)/(1-p) - a = 1/L, where L is the maximum optical path difference of an equivalent interferometer. "a" determines the resolution, or the width of the scanning function. - FZ = distance from the center of the scanning function to the first zero. - CR = critical value for the ratio of the extent of the spectrum (v2-v1) and the HWHM. When the ratio is less than CR, the apodization function is calculated as the FFT of the scanning function. When the ratio is greater than CR, the apodization function is calculated analytically. - J(5/2,u)/(u**(5/2)) = ((3-u**2)*sin(u)-3*u*cos(u))/u**5, where J(n,u) is the Bessel function of order n. - Io = modified Bessel function of the first kind, order 0. - c1 = 0.428752 _______________________________________________________________________________________________________________________ Notes: 1. If JFNin < 0, HWHM contains the value for the maximum optical path difference of an equivalent interferometer, apodized to give the scanning function given by |JFNin|. 2. The Brault and Kaiser-Bessel functions depend upon the parameter p. The values of a/HWHM and a/FZ for these functions also depend upon p. FFTSCN (Continued) RECORD RECORD 10.1 (Continued) MRATin ratio of HWHM of the scanning function to the halfwidth of the boxcar; for prescanning with a boxcar (default = 12.). If MRATin < 0., no boxcaring is performed. DVOUT output grid spacing. The scanned spectral file is interpolated onto the grid defined by V1, V2, and DVOUT. IUNIT unit designation of file to be scanned (default is 12) IFILST initial file from IUNIT to be scanned NIFILS number of files to be scanned starting at IFILST JUNIT file containing scanned results (default is 11) IVX (-1,0,1) flag for calculation of scanning function = -1 scanning function is calculated as the FFT of the apodization function = 0 program determines how to calculate the scanning function, using CR = 1 scanning function is calculated analytically NOFIX flag for deconvolution (used when prescanning with a boxcar) = 0 deconvolve the scanned spectrum with the boxcar = nonzero, do not deconvolve with the boxcar RECORD RECORD 10.2 (required if ABS(JFNin) > 10; otherwise omit) PARM1, PARM2, PARM3 1-10, 11-20, 21-30 F10.4, F10.4, F10.4 PARM1 For JFNin = 11, value of p, where (0 < p < 1) For JFNin = 12, value of p, where (2 < p < 4) For JFNin = 13, value of v_offset PARM2 For JFNin = 13, value of c1 PARM3 For JFNin = 13, value of c2 REPEAT RECORDS 10.1, 10.2 A '-1.' within columns 1-10 will terminate FFT scanning. FLTRFN (IFILTR = 1) RECORD RECORD 11.1 V1F, DVF, NPTS, JEMIT, IUNIT, IFILST, NIFILS, JUNIT, HEDDR 1-10, 11-20, 21-25, 26-30, 31-35, 36-40, 41-45, 46-50, 51-85 F10.3, F10.4, I5, I5, I5, I5, I5, I5, 8A4,A3 V1F wavenumber of initial filter value negative value terminates the FLTRFN option DVF wavenumber increment between filter values negative value indicates V1F value is center wavenumber of filter function (NOTE: requires NPTS to be an odd number) NPTS number of filter values (if NPTS.LT.0, previous values of XF utilized) JEMIT = -1 absorption integrated over filter function = 0 transmittance integrated over filter function = 1 radiance integrated over filter function IUNIT file designation for integration IFILST initial file to be integrated NIFILS number of files to be integrated starting at IFILST JUNIT file designation for output file (named "FLT_OUT") HEDDR 35 characters of user identification FLTRFN (Continued) RECORD RECORD 11.2 IVAR 1-80 20A4 IVAR is format specification for reading filter values, XF FLTRFN (Continued) RECORD RECORD 11.3 . . . . . . . . . 11.N XF format specification given by record 11.2 XF NPTS values of FILTRFN REPEAT RECORD 11.1 A '-1.' within columns 1-10 will terminate FLTRFN. PLTLBL RECORD RECORD 12.1 CPRGID, CEX 1-60, 79-80 A60, 18X, A2 CPRGID identifying header for plot files CEX parameter to indicate external file name to be provided by user value of 'EX' indicates external file names to be supplied **************************************************************************************************** **************************************************************************************************** **************** Use Records 12.2A and 12.3A for PLOT (IOPT = 0,1) ****************** **************** Use Record 12.2B for File Difference or Ratio (IOPT = 2,3) ****************** **************************************************************************************************** **************************************************************************************************** RECORD RECORD 12.2A (IOPT = 0,1) V1, V2, XSIZE, DELV,NUMSBX,NOENDX, LFILE,LSKIPF, SCALE, IOPT, I4P, IXDEC 1-10, 11-20, 21-30, 31-40, 45, 50, 55, 60, 61-70, 72, 75, 80 F10.4, F10.4, F10.4, F10.4, I5, I5, I5, I5, F10.3, I2, I3, I5 V1 initial wavenumber of the plot, negative value terminates the PLTLBL option V2 the final wavenumber of the plot XSIZE number of inches of the X-axis. DELV number of wavenumbers (cm - 1) per major division. NUMSBX * IOPT ne 1 NUMSBX controls the number of subdivisions per major division of x-axis. NUMSBX = 1 gives no tic-marks per major division = 2 gives one tic-mark per major division * IOPT eq 1 NUMSBX controls the use of symbols for overlayed plots NUMSBX = 0 produces a line plot without symbols > 0 produces a line plot with symbols at every NUMSBX point < 0 produces a point plot with a symbol at every NUMSBX point NOENDX * IOPT =0,2,3 NOENDX controls the numbers at either end of the x-axis NENDX = 0 includes the numbers at both ends of the axis 1 suppresses the numbers at either end of the axis 2 suppresses the beginning number 3 suppresses the ending number * IOPT eq 1 NOENDX selects the symbol to be used on the overlayed plot. For a list of symbols see the applicable NCAR plotter information. PLTLBL (Continued) RECORD RECORD 12.2A (Continued) LFILE tape number of file to be read from LBLRTM. LSKIPF number of files to be skipped in tape LFILE number of file to be plotted will be LSKIPF + 1 for plotting MFILE (TAPE12) LSKIPF = 0 SCALE enables one to enlarge or reduce a plot. IOPT = 0 for plot, USE instructions for records 12.2A and 12.3A = 1 for overlay plot, USE instructions for records 12.2A and 12.3A * will suppress the plotting of the header and axes, * and will overlay the line on the previous plot axes. * See also NUMSBX, NOENDX and YMIN. = 2 for file difference, USE instructions for record 12.2B instead of 12.2A = 3 for file ratio, USE instructions for record 12.2B instead of 12.2A I4P = 0 for linear connection of points, = 1 for four point interpolation IXDEC number of figures after decimal point on x-axis --------------------------------------------------------------------------------------------------------- RECORD RECORD 12.2.1A (IOPT = 0,1; CEX = 'EX') CFILEN(1) 1-25 A25 CFILEN(1) is the external name of the data file to be plotted or used for the overlay data (i.e. LFILE) --------------------------------------------------------------------------------------------------------- PLTLBL (Continued) RECORD RECORD 12.3A (IOPT = 0,1) YMIN, YMAX, YSIZE, DELY,NUMSBY,NOENDY, IDEC, JEMIT, JPLOT,LOGPLT, JHDR, JOUT,JPLTFL 1-10, 11-20, 21-30, 31-40, 45, 50, 55, 60, 65, 70, 72, 77, 80 G10.4, G10.4, F10.3, F10.3, I5, I5, I5, I5, I5, I5, I2,3X, I2, I3 YMIN is Y value at bottom of Y-axis. * IOPT = 1 will cause YMIN to determine the vertical offset (in Y) * for the overlayed plot. YMIN is in the units of the previous plot * axes and any offset is applied with respect to the previously * defined axes. (i.e. YMIN = 10. will offset the new plot 10. units up) YMAX is Y value at top of Y-axis. (if log plot selected, YMIN and YMAX are the exponent values at the bottom and top of the plot). YSIZE is the number of inches for the y-axis. DELY is the number of y units per major division (= 1. when log plot is selected). NUMSBY is the number of subdivisions per major division of y-axis. NOENDY controls the plotting of the values at either end of y-axis (= 1 supresses plotting). IDEC is number of figures after decimal point on linear y-axis. JEMIT = 0 for transmission, = 1 for radiance. JPLOT = 0 plots transmission for JEMIT = 0 = 1 plots optical depth for JEMIT = 0 = 2 plots attenuation in decibels for JEMIT = 0 = 0 plots radiance in watts/(cm^2 sr cm^-1) for JEMIT = 1 (if YMIN and YMAX.GT.1, these values are interpreted as brightness temperatures from which minimum and maximum values of radiance are assigned). = 1 plots radiance in equivalent brightness temperature for JEMIT = 1 LOGPLT = 0 for linear y-axis, = 1 for log y-axis. JHDR = 0 for plot of header data = 1 for suppression of header data plot PLTLBL (Continued) RECORD RECORD 12.3A (Continued) JOUT = 0 for plot to system plotfile/screen = 1 for plot values written to JPLTFL (unformatted) = 2 for plot to system plotfile/screen and plot values written to JPLTFL (unformatted) = 3 for plot values written to JPLTFL (formatted) = 4 for plot to system plotfile/screen and plot values written to JPLTFL (formatted) Note: For JOUT=3 & 4 the program writes out a description of the header followed by both the frequency and data arrays. Format is below. FORMAT (3X,1P,E15.8,4X,E15.8) JPLTFL file to which plot values are written for JOUT .gt. 0 (default = 29) the structure of JPLTFL (unformatted) is standard LBLRTM file structure (see Table III), plot values are written out in SUBROUTINE FSCLIN REPEAT RECORD 12.2A or RECORD 12.2B A '-1.' within columns 1-10 will terminate plotting. PLTLBL (Continued) *************************************************************************************************** ****************************** USE RECORD 12.2B FOR IOPT > 1 ************************************** *************************************************************************************************** RECORD RECORD 12.2B (IOPT = 2,3) V1, V2, JFILE,JSKIPF, LFILE,LSKIPF, IOPT, MFILE 1-10, 11-20, 45, 50, 55, 60, 72, 80 F10.4, F10.4, 20X,I5, I5, I5, I5, 10X,I2, 3X,I5 V1 initial wavenumber of the difference or ratio, negative value terminates the PLTLBL option V2 the final wavenumber of the difference or ratio JFILE tape number of file to be read from LBLRTM (no default) JSKIPF number of files to be skipped in tape JFILE number of file to be used will be JSKIPF + 1 LFILE tape number of file to be read from LBLRTM (no default) LSKIPF number of files to be skipped in tape LFILE number of file to be used will be LSKIPF + 1 IOPT = 0 for plot, USE Records 12.2A and 12.3A = 1 for overlay plot, USE Records 12.2A and 12.3A = 2 for file difference (JFILE-LFILE) = 3 for file ratio (JFILE/LFILE) MFILE tape number of file for difference/ratio output (no default) Subroutine FILOPT will difference or ratio two LBLRTM output files. Files must be 1) UNFORMATTED LBLRTM files 2) Single quantity files i.e. contain only ONE of the following: a) Optical Depths b) Transmittance c) Radiance e) Temperature 3) The difference in DV and V1 between the files must be less than 1.E-8 *** NOTE: Standard LBLRTM output files usually contain BOTH Transmittance and Radiance. The user must create either scanned or plot files for use as input to this routine. --------------------------------------------------------------------------------------------------------- (IOPT = 2,3; CEX = 'EX') RECORD RECORD 12.2.1B CFILEN(1) 1-25 A25 CFILEN(1) is the external name for the data file associated with JFILE. RECORD RECORD 12.2.2B CFILEN(2) 1-25 A25 CFILEN(2) is the external name for the data file associated with LFILE. RECORD RECORD 12.2.3B CFILEN(3) 1-25 A25 CFILEN(3) is the external name for the data file associated with the difference/ratio file, MFILE. --------------------------------------------------------------------------------------------------------- REPEAT RECORD 12.2A or RECORD 12.2B A '-1.' within columns 1-10 will terminate plotting. TABLE III. Structure of file JPLTFL (standard LBLRTM structure). TABLE III record 1 /PLTHDR/ XID(10),SEC,P0,T0,HMOL(64),W(64),WBROAD,DVT,V1V, V2V,TBOUND,EMISIV,FSCDID(17),NMOL,NLAYER, YID1,YID(10) -------------------------------------------------------------------------- record 2 /JPLTFL/ V1P,V2P,DVP,NLIM record 3 NLIM values of Y from V1P to V2P with increment DVP ---------- repeats records 2 and 3 until all Y values written --------- record 4 LBLRTM endfile --------------- repeats for each plot request ------------------------------- NOTES: 1. for machines with 32 bit words the following variables may have been set to double precision: XID(10),SEC,HMOL(64),V1V,V2V,YID(10) V1P,V2P 2. fftscn.f has parameters that are hardware dependent. Before compiling, check to make sure the appropriate statements are uncommented. 3. the placement of the decimal point in the input supercedes that of a given format. 4. if WKL(M,L) or XAMNT(I,L) is greater than one, then column density is assumed for that value of M or I and L; if either is less than one, then mixing ratio is assumed for that value of M or I and L. This allows for combining density and mixing ratio information in the input. 5. the character '%' signals the end of TAPE5 and terminates LBLRTM. _______________________________________________________________________________________________________________________ If the NLTE (IHIRAC = 4) option is selected a separate file (TAPE4) must be created. The following instructions describe the contents of a TAPE4 file. Note : TAPE4 permits selection of vibrational states; only the states listed in Table IV are allowed at present, with up to 9 isotopologues for each state. Table IV -------- Molecule State Energy Degeneracies -------- ----- -----(cm-1)----- ------------ H2O 000 0. 1 --- 010 1594.750 1 020 3151.630 1 100 3657.053 1 001 3755.930 1 030 4666.793 1 110 5234.977 1 011 5333.269 1 CO2 --- 00001 0. 1 01101 667.380 2 10002 1285.409 1 02201 1335.132 2 10001 1388.185 1 11102 1932.470 2 03301 2003.246 2 11101 2076.856 2 00011 2349.143 1 20003 2548.366 1 12202 2585.022 2 20002 2671.143 1 04401 2671.717 2 12201 2760.725 2 20001 2797.135 1 01111 3004.012 2 10012 3612.842 1 02211 3659.273 2 10011 3714.783 1 11112 4247.706 2 03311 4314.914 2 11111 4390.629 2 20013 4853.623 1 04411 4970.931 1 20012 4977.834 1 20011 5099.660 1 O3 -- 000 0. 1 010 700.931 1 001 1042.084 1 100 1103.140 1 020 1399.275 1 011 1726.528 1 110 1796.261 1 002 2057.892 1 101 2110.785 1 200 2201.157 1 111 2785.245 1 003 3041.200 1 004 3988. 1 005 4910. 1 006 5803. 1 007 6665. 1 008 7497. 1 009 8299. 1 CO -- 0 0. 1 1 2143.272 1 2 4260.063 1 NO -- 0 0. 1 1 1878.077 1 2 3724.067 1 TAPE4(RECORD 1.1) TIT 1-80 A80 TIT: 80 characters of user indentification REPEAT TAPE4-RECORD 1.1 up to 20 times. TAPE4(RECORD 1.2) IVIB, MOLNEQ 5, 10 I5 I5 IVIB: (0,1) selects input for vibrational data: = 1 (vibrational temperatures (K)) = 0 (population) MOLNEQ: place holder for future input -------------------------------------------------------------------------------- ****** these records applicable only if IVIB eq 1 *** -------------------------------------------------------------------------------- TAPE4(RECORD 2.1) TEXTLINE 1-80 A80 TEXTLINE: character string beginning with two or more hyphens ¡®-¡®and including the molecule name directly after the hyphens and optional spaces. The word ¡®VIBRATIONAL¡¯ should be included with the molecule name for vibrational data. Example record: ----- H2O VIBRATIONAL STATE DATA TAPE4(RECORD 2.2) NUMSTATE, ISOSTATE, IDSTATE, EESTATE, NDGSTATE Free Format NUMSTATE: numerical index corresponding to the AFGL mapping for the vibrational state ISOSTATE: isotopologue number of the vibrational state IDSTATE: state identifier enclosed in single quotes EESTATE: energy level (cm^-1) of the state NDGSTATE: number of degeneracies of the state Example records for H2O: 1 1 ¡®000¡¯ 0. 1 2 1 ¡®010¡¯ 1594.750 1 3 1 ¡®020¡¯ 3151.630 1 4 1 ¡®100¡¯ 3657.053 1 5 1 ¡®001¡¯ 3755.930 1 6 1 ¡®030¡¯ 4666.793 1 7 1 ¡®110¡¯ 5234.977 1 8 1 ¡®011¡¯ 5333.269 1 TAPE4(RECORD 2.3) TEXTLINE 1-80 A80 TEXTLINE: character string beginning with one or more hypens ¡®-¡®and including the molecule name directly after the hyphens and optional spaces. Example record: ----- H2O TEMPERATURE PROFILES TAPE4(RECORD 2.4) ALT, TMP(1), TMP(2:N) Free Format ALT: altitude (km) TMP(1): kinetic temperature (K) for this altitude TMP(2:N): two state vibrational temperature (K) with respect to the ground state for the N vibrational modes listed in Record 2.2. These values should be set to zero if the mode is not desired. Example records for H2O: 7. 243.200 243.203 243.200 243.205 243.214 or 7. 243.200 0, 0, 0, 0, 0, 0. 0. Repeat records 2.1, 2.2, 2.3, and 2.4 for each non-LTE molecules (H2O, CO2, O3, CO, and NO) desired. All molecules need not be input. TAPE4(RECORD 3.1) TEXTLINE 1-80 A80 TEXTLINE: character string beginning with one or more hypens ¡®-¡®and including the string ¡°END OF DATA¡¯ directly after the hyphens and optional spaces. Example record: ---- END OF DATA ________________________________________________________________________________________________________________________