Eric,
There is no ready ARTS method for setting up backend_channel_response.
But I think the description
arts -d backend_channel_response
is quite clear and you should be able to use Matlab or Python to set it
up. What you need to do is basically to repeat for each channel what
backend_channel_responseFlat produces.
There is one demo in Atmlab doing this. Attached.
Bye,
Patrick
On 2022-08-15 17:12, Eric Sauvageat wrote:
Dear ARTS community,
For my ground-based ozone retrievals, I'm trying to implement a
rectangular channel response for my spectrometer. Actually, I am dealing
with a filter bank spectrometer with refined resolution at the line
center so I am trying to implement individual channel widths for each
spectrometer channels.
It seems that the workspace method "backend_channel_responseFlat" would
be the correct one to use for this purpose but it seems to only use a
single "resolution" value for all channels. Is that correct and if yes,
is there any ways to input an "array of channel width" to get different
response for each channel ?
Thanks in advance for your help,
% QARTS_BACKEND_DEMO Example on how to introduce a passband filter backend
%
% A simple simulation of a ground-based temperature radiometer, having
% passbbands of different widths.
%
% FORMAT [f,y] = qarts_backend_demo( zplat, za )
%
% OUT f Centre frequency of passbands.
% y Simulated brightness temperature (RJTB)
% IN zplat Platform altitude.
% za Zenith angle.
% 2013-04-17 Created, Patrick Eriksson.
function [f,y] = qarts_backend_demo(zplat,za)
%= Atmlab settings
%
arts_xmldata_path = atmlab( 'ARTS_XMLDATA_PATH' );
%
if isnan( arts_xmldata_path )
error('You need to ARTS_XMLDATA_PATH to run this example.');
end
%
fascod = fullfile( arts_xmldata_path, 'planets', 'Earth', 'Fascod' );
%= Init Q structures
%
Q = qarts;
%
Q.INCLUDES = { fullfile( 'ARTS_INCLUDES', 'general.arts' ), ...
fullfile( 'ARTS_INCLUDES', 'agendas.arts' ), ...
fullfile( 'ARTS_INCLUDES', 'continua.arts' ), ...
fullfile( 'ARTS_INCLUDES', 'planet_earth.arts' ) };
%= Define agendas
%
Q.PPATH_AGENDA = { 'ppath_agenda__FollowSensorLosPath' };
Q.PPATH_STEP_AGENDA = { 'ppath_step_agenda__GeometricPath' };
Q.IY_SPACE_AGENDA = { 'iy_space_agenda__CosmicBackground' };
Q.IY_SURFACE_AGENDA = { 'iy_surface_agenda__UseSurfaceRtprop' };
Q.IY_MAIN_AGENDA = { 'iy_main_agenda__Emission' };
%= Species
%
Q.ABS_SPECIES(1).TAG{1} = 'N2-SelfContStandardType';;
Q.ABS_SPECIES(2).TAG{1} = 'H2O-PWR98';
Q.ABS_SPECIES(3).TAG{1} = 'O2-PWR93';
%= Set-up atmosphere and surface
%
Q.ATMOSPHERE_DIM = 1;
%
Q.P_GRID = z2p_simple( [-100:250:50e3]' );
%
atm = 'midlatitude-summer';
Q.RAW_ATMOSPHERE = fullfile( arts_xmldata_path, 'planets', 'Earth', ...
'Fascod', atm, atm );
Q.RAW_ATM_EXPAND_1D = false;
%
Q.REFELLIPSOID = ellipsoidmodels( 'SphericalEarth' );
Q.Z_SURFACE = 0;
%= Absorption
%
Q.ABS_LINES_FORMAT = 'none';
%
Q.ABSORPTION = 'OnTheFly';
%= Various
%
Q.CLOUDBOX_DO = false;
Q.J_DO = false;
%
Q.STOKES_DIM = 1;
Q.PPATH_LMAX = -1;
%
Q.YCALC_WSMS = { 'yCalc' };
%
Q.IY_UNIT = 'RJBT';
%= Define spectrometer passbands
%
f = 1e9 * [ 51.26 52.28 53.86 54.94 56.66 57.3 58.0 ]';
df = 1e6 * [ 230 230 230 230 600 1000 2000 ]';
%= Init sensor definition
%
Q.SENSOR_POS = zplat;
Q.SENSOR_LOS = za;
%
H = qartsSensor;
H.SENSOR_NORM = true;
H.ANTENNA_DO = false;
%= Include spectrometer
%
H.BACKEND_DO = true;
H.F_BACKEND = f;
%
B.name = 'Spectrometer channel response function';
B.gridnames = {'Frequency'};
B.dataname = 'Response';
B.data = [1 1];
%
% Here we need to define one response for each channel, but only the
% width differs
for i = 1 : length(df)
B.grids = { df(i)/2*[-1 1] };
H.BACKEND_CHANNEL_RESPONSE{i} = B;
end
%= Frequency grid
%
Q.F_GRID = linspace( f(1)-df(1)*0.51, f(end)+df(end)*0.51, 100 )';
%= One option to manage with a shorter F_GRID is to use "filling", but make
% very careful tests before starting to use it:
%
%H.FILL_FGRID = [ 4 5 ];
%= Finalise sensor definition
%
Q.SENSOR_RESPONSE = H;
Q.SENSOR_DO = true;
%= Calculate spectrum/spectra
%
y = arts_y( Q );
