Read optical depth output

[maclion-zhao]

Hi I am running LBLRTM12.17 and get layer optical depth (Dint_***) with double precision, could you provide an IDL code to read the output?

Thanks

[pernak18]

@maclion-zhao try the solution in #59 (comment). it should work with the ODInt files as well.

note that's a Python solution. @jrmascio do you have the IDL script on hand that you could just upload into this issue?

please and thanks

[maclion-zhao]

@maclion-zhao try the solution in #59 (comment). it should work with the ODInt files as well.

note that's a Python solution. @jrmascio do you have the IDL script on hand that you could just upload into this issue?

please and thanks

Thanks, I really need IDL code.

[jrmascio]

@maclion-zhao I will email you the IDL script since GitHub doesn't allow for .pro files to be attached here, but the calling statement to run it is:

read_lbl_file, OD_fileName, od, wv, /double, file_type=0

[zhangxuehai2012]

Hi I am running LBLRTM12.17 and get layer optical depth (ODExact_) files. And I found that you have a program for reading OD files, could you provide a code to read the ODExact_ files?
Thanks very much.

[zhangxuehai2012]

@maclion-zhao I will email you the IDL script since GitHub doesn't allow for .pro files to be attached here, but the calling statement to run it is:

read_lbl_file, OD_fileName, od, wv, /double, file_type=0

Hi I am running LBLRTM12.17 and get layer optical depth (ODExact_) files. And I found that you have a program for reading OD files, could you provide a code to read the ODExact_ files?
Thanks very much.

[Xuly56]

@maclion-zhao I will email you the IDL script since GitHub doesn't allow for .pro files to be attached here, but the calling statement to run it is:

read_lbl_file, OD_fileName, od, wv, /double, file_type=0

Hi I came across the same issue. I'm running LBLRTM12.17 and get layer optical depth (ODint_) files with double precision, but I'm unable to read them. Could you please also provide an IDL script to read these outputs? Since the solution mentioned in 59 seems to work only for TAPE 10,11,12,13. Thanks in advance!

[maclion-zhao]

@maclion-zhao I will email you the IDL script since GitHub doesn't allow for .pro files to be attached here, but the calling statement to run it is:
read_lbl_file, OD_fileName, od, wv, /double, file_type=0

Hi I came across the same issue. I'm running LBLRTM12.17 and get layer optical depth (ODint_) files with double precision, but I'm unable to read them. Could you please also provide an IDL script to read these outputs? Since the solution mentioned in 59 seems to work only for TAPE 10,11,12,13. Thanks in advance!

Hi Xuly,

Here attached is the code Jeana sent to me:

;+
; NAME: read_lbl_file.pro
;
;
;
; PURPOSE:;-
;
; This file is an IDL (Interactive Data Language, by Research Systems, Inc.)
; procedure to read a LBLRTM format spectral file.
;
; It can handle all types of files (monochromatic or scanned
; optical depth, transmittance, or radiance), monochromatic
; radiance (which contains both transmittance and radiance).
;
;
; CATEGORY:
;
; --------------------------------------------------------------------------
; | |
; | Copyright 2002 - 2005, Atmospheric & Environmental Research, Inc. (AER).|
; | This software may be used, copied, or redistributed as long as it is |
; | not sold and this copyright notice is reproduced on each copy made. |
; | This model is provided as is without any express or implied warranties. |
; | (http://www.rtweb.aer.com/) |
; | |
; --------------------------------------------------------------------------
;
; CALLING SEQUENCE:
;
;
;
; INPUTS:
;
; filename: filename
; file_type: 1: radiance from monochromatic radiance calculation,
; 2: transmittance " " " "
; 3: Aerosol absorbtance from spectral aerosol transmittance
; file (TAPE20)
; 4: Aerosol scattering from spectral aerosol transmittance
; file (TAPE20)
; 5: Aerosol asymmetry from spectral aerosol transmittance
; file (TAPE20)
; 6: transmittance from monochromatic radiance calculation,
; 0: otherwise (scanned transmittance or radiance, optical depth)
;
; OPTIONAL INPUTS:
;
;
;
; KEYWORD PARAMETERS:
;
; double: keyword set if reading in file in double precision
; vstart: return data starting with wavenumber vstart
; n_file: number of files read in
;
;
;
; OUTPUTS:
;
; spectrum: LBLRTM spectral data
; frequency: frequency grid
; npts: Number of points
;
; OPTIONAL OUTPUTS:
;
;
;
; COMMON BLOCKS:
;
;
;
; SIDE EFFECTS:
;
;
;
; RESTRICTIONS:
;
; Where multiple LBLRTM files (say sub-files) exist
; on a single file (e.g. from multiple scans or IMRG <> 0)
; it will read one sub-file.
;
; PROCEDURE:
;
;
;
; EXAMPLE:
;
;
;
; MODIFICATION HISTORY:
;
; ----------------------------------------------------------------------------
; Written by:
;
; William Gallery
;
; Updated by:
; Patrick Brown
;
; Updated by:
; Mark Shephard, June, 2004; mshep@aer.com
; 131 Hartwell Ave
; Lexingtion, MA
; USA
; 02421
;
; ----------------------------------------------------------------------------

pro read_lbl_file,filename,spectrum,frequency,npts,file_type=file_type,double=double,vstart=vstart,n_file=n_file,ierr

ierr=0

unit=1

if keyword_set(file_type) then file_type = file_type else file_type = 0
if keyword_set(vstart) then vstart = vstart else vstart = 0
if keyword_set(n_file) then nf = n_file else nf = 1

openr,unit,filename, /get_lun , /f77_unf ;,/SWAP_ENDIAN

for i = 1,nf do begin ;loop over n_files, finish with the last one

;read header
header = strarr(80,1)
;header = strarr(100,1)
readu,unit, header
    ;print, 'First Panel Header'
print,header
header = transpose(header)

pv1 = 0.0d & pv2 = 0.0d

    IF keyword_set(double) THEN BEGIN         
        pdv = 0.0d 
        np_arr = LONARR(2) 
    ENDIF ELSE BEGIN 
        pdv = 0.0
        np_arr = LONARR(1)
    ENDELSE
    panel = 0

    ;;the integer np is being stored in a real in LBLRTM panel header
    ;;therefore, it can be located in different spots
    ;;on the unix and linux (due to bit swapping)

if (not(eof(unit))) THEN readu,unit,pv1,pv2,pdv,np_arr

print,pv1,pv2,pdv,np_arr

    ;;find where the value is located
    w_p = where(np_arr GT 0)
    np = np_arr(w_p(0))
    ;print, 'np', np

if (np le 0) then begin
	print, 'File may be corrupt, choose new file.'
	ierr=1
	return
endif

v1 = pv1 & v2 = pv2 & dv = pdv    
npts = 0l

; n_max = 6000000l
; n_max = 7000000l
n_max = 90000001

    ;;Maximum of points in spectrum
    IF keyword_set(double) THEN BEGIN 
        spectrum = dblarr(n_max)          ;dbl
    ENDIF ELSE spectrum = fltarr(n_max)   ;sgl

on_ioerror, io_error
error = ''

;loop over panels
while (np gt 0) do BEGIN
        IF keyword_set(double) THEN BEGIN 
            temp1 = dblarr(np)            ;dbl
        ENDIF ELSE temp1 = fltarr(np)     ;sgl
    readu,unit,temp1

        IF keyword_set(double) THEN BEGIN 
            temp2 = dblarr(np)            ;dbl
        ENDIF ELSE temp2 = fltarr(np)     ;sgl
    if file_type ge 1 then begin
        readu,unit,temp2
    ENDIF

        IF keyword_set(double) THEN BEGIN 
            temp3 = dblarr(np)            ;dbl
        ENDIF ELSE temp3 = fltarr(np)     ;sgl
    if file_type gt 2 AND file_type lt 6 then begin
        readu,unit,temp3
    endif

    if npts+np gt n_max then begin
        error = 'Attempting to read more than n_max points of spectrum'
	print, 'npts = ',npts,'  n_max = ',n_max
        print,error
        stop, 'Enter .cont to continue or retall to exit'
	goto,finish
    endif

    if pv2 ge vstart then begin
        panel = panel+1
	if panel eq 1 then  begin 
           v1 = pv1
           ;print,'start accumulating spectrum, v1 = ',v1
        endif
	case 1 of
		(file_type lt 2): spectrum(npts) = temp1
		(file_type eq 3): spectrum(npts) = temp1
		(file_type eq 5): spectrum(npts) = temp3
		(file_type eq 6): spectrum(npts) = temp2
			    else: spectrum(npts) = temp2
	endcase
        npts = npts+np
    endif

    ; read next panel header
    if (not(eof(unit))) then BEGIN 
            readu, unit, pv1, pv2, pdv, np_arr 
            np = np_arr(w_p(0))
            ;print, pv1, pv2, pdv, np
        ENDIF ELSE np = -99
    ;print,pv1,pv2,pdv,np
    if np gt 0 then v2 = pv2     ;Reset v2  if valid panel
endwhile

finish:

spectrum = spectrum(0:npts-1)

;create frequency grid
frequency = v1+findgen(npts)*dv

endfor

free_lun,unit

return

io_error:

print,'IO Error****: ',!err_string

;error = strmessage(!err)
;print,'IO Error****: ',error

if (!err ne -118) then close,unit

print,npts
if npts gt 0 then begin
spectrum = spectrum(0:npts-1)
frequency = v1+findgen(npts)*dv
return
endif else begin
retall
endelse

end

[zhangxuehai2012]

Dear Zhao Thanks very much! We have successfully read the OD file using the RC_utils file on github. Thank you again for your help！ Xuehai Zhang Henan University of Technology E-mail: ***@***.*** At 2025-11-20 23:04:39, "maclion-zhao" ***@***.***> wrote: maclion-zhao left a comment (AER-RC/LBLRTM#80) @maclion-zhao I will email you the IDL script since GitHub doesn't allow for .pro files to be attached here, but the calling statement to run it is: read_lbl_file, OD_fileName, od, wv, /double, file_type=0 Hi I came across the same issue. I'm running LBLRTM12.17 and get layer optical depth (ODint_) files with double precision, but I'm unable to read them. Could you please also provide an IDL script to read these outputs? Since the solution mentioned in 59 seems to work only for TAPE 10,11,12,13. Thanks in advance! Hi Xuly, Here attached is the code Jeana sent to me: ;+ ; NAME: read_lbl_file.pro ; ; ; ; PURPOSE:;- ; ; This file is an IDL (Interactive Data Language, by Research Systems, Inc.) ; procedure to read a LBLRTM format spectral file. ; ; It can handle all types of files (monochromatic or scanned ; optical depth, transmittance, or radiance), monochromatic ; radiance (which contains both transmittance and radiance). ; ; ; CATEGORY: ; ; -------------------------------------------------------------------------- ; | | ; | Copyright 2002 - 2005, Atmospheric & Environmental Research, Inc. (AER).| ; | This software may be used, copied, or redistributed as long as it is | ; | not sold and this copyright notice is reproduced on each copy made. | ; | This model is provided as is without any express or implied warranties. | ; | (http://www.rtweb.aer.com/) | ; | | ; -------------------------------------------------------------------------- ; ; CALLING SEQUENCE: ; ; ; ; INPUTS: ; ; filename: filename ; file_type: 1: radiance from monochromatic radiance calculation, ; 2: transmittance " " " " ; 3: Aerosol absorbtance from spectral aerosol transmittance ; file (TAPE20) ; 4: Aerosol scattering from spectral aerosol transmittance ; file (TAPE20) ; 5: Aerosol asymmetry from spectral aerosol transmittance ; file (TAPE20) ; 6: transmittance from monochromatic radiance calculation, ; 0: otherwise (scanned transmittance or radiance, optical depth) ; ; OPTIONAL INPUTS: ; ; ; ; KEYWORD PARAMETERS: ; ; double: keyword set if reading in file in double precision ; vstart: return data starting with wavenumber vstart ; n_file: number of files read in ; ; ; ; OUTPUTS: ; ; spectrum: LBLRTM spectral data ; frequency: frequency grid ; npts: Number of points ; ; OPTIONAL OUTPUTS: ; ; ; ; COMMON BLOCKS: ; ; ; ; SIDE EFFECTS: ; ; ; ; RESTRICTIONS: ; ; Where multiple LBLRTM files (say sub-files) exist ; on a single file (e.g. from multiple scans or IMRG <> 0) ; it will read one sub-file. ; ; PROCEDURE: ; ; ; ; EXAMPLE: ; ; ; ; MODIFICATION HISTORY: ; ; ---------------------------------------------------------------------------- ; Written by: ; ; William Gallery ; ; Updated by: ; Patrick Brown ; ; Updated by: ; Mark Shephard, June, 2004; ***@***.*** ; 131 Hartwell Ave ; Lexingtion, MA ; USA ; 02421 ; ; ---------------------------------------------------------------------------- pro read_lbl_file,filename,spectrum,frequency,npts,file_type=file_type,double=double,vstart=vstart,n_file=n_file,ierr ierr=0 unit=1 if keyword_set(file_type) then file_type = file_type else file_type = 0 if keyword_set(vstart) then vstart = vstart else vstart = 0 if keyword_set(n_file) then nf = n_file else nf = 1 openr,unit,filename, /get_lun , /f77_unf ;,/SWAP_ENDIAN for i = 1,nf do begin ;loop over n_files, finish with the last one ;read header header = strarr(80,1) ;header = strarr(100,1) readu,unit, header ;print, 'First Panel Header' print,header header = transpose(header) pv1 = 0.0d & pv2 = 0.0d IF keyword_set(double) THEN BEGIN pdv = 0.0d np_arr = LONARR(2) ENDIF ELSE BEGIN pdv = 0.0 np_arr = LONARR(1) ENDELSE panel = 0 ;;the integer np is being stored in a real in LBLRTM panel header ;;therefore, it can be located in different spots ;;on the unix and linux (due to bit swapping) if (not(eof(unit))) THEN readu,unit,pv1,pv2,pdv,np_arr print,pv1,pv2,pdv,np_arr ;;find where the value is located w_p = where(np_arr GT 0) np = np_arr(w_p(0)) ;print, 'np', np if (np le 0) then begin print, 'File may be corrupt, choose new file.' ierr=1 return endif v1 = pv1 & v2 = pv2 & dv = pdv npts = 0l ; n_max = 6000000l ; n_max = 7000000l n_max = 90000001 ;;Maximum of points in spectrum IF keyword_set(double) THEN BEGIN spectrum = dblarr(n_max) ;dbl ENDIF ELSE spectrum = fltarr(n_max) ;sgl on_ioerror, io_error error = '' ;loop over panels while (np gt 0) do BEGIN IF keyword_set(double) THEN BEGIN temp1 = dblarr(np) ;dbl ENDIF ELSE temp1 = fltarr(np) ;sgl readu,unit,temp1 IF keyword_set(double) THEN BEGIN temp2 = dblarr(np) ;dbl ENDIF ELSE temp2 = fltarr(np) ;sgl if file_type ge 1 then begin readu,unit,temp2 ENDIF IF keyword_set(double) THEN BEGIN temp3 = dblarr(np) ;dbl ENDIF ELSE temp3 = fltarr(np) ;sgl if file_type gt 2 AND file_type lt 6 then begin readu,unit,temp3 endif if npts+np gt n_max then begin error = 'Attempting to read more than n_max points of spectrum' print, 'npts = ',npts,' n_max = ',n_max print,error stop, 'Enter .cont to continue or retall to exit' goto,finish endif if pv2 ge vstart then begin panel = panel+1 if panel eq 1 then begin v1 = pv1 ;print,'start accumulating spectrum, v1 = ',v1 endif case 1 of (file_type lt 2): spectrum(npts) = temp1 (file_type eq 3): spectrum(npts) = temp1 (file_type eq 5): spectrum(npts) = temp3 (file_type eq 6): spectrum(npts) = temp2 else: spectrum(npts) = temp2 endcase npts = npts+np endif ; read next panel header if (not(eof(unit))) then BEGIN readu, unit, pv1, pv2, pdv, np_arr np = np_arr(w_p(0)) ;print, pv1, pv2, pdv, np ENDIF ELSE np = -99 ;print,pv1,pv2,pdv,np if np gt 0 then v2 = pv2 ;Reset v2 if valid panel endwhile finish: spectrum = spectrum(0:npts-1) ;create frequency grid frequency = v1+findgen(npts)*dv endfor free_lun,unit return io_error: print,'IO Error****: ',!err_string ;error = strmessage(!err) ;print,'IO Error****: ',error if (!err ne -118) then close,unit print,npts if npts gt 0 then begin spectrum = spectrum(0:npts-1) frequency = v1+findgen(npts)*dv return endif else begin retall endelse end — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you commented.Message ID: ***@***.***>