From 5c307c0e67dce431d2109daa75342f79ab0f18b7 Mon Sep 17 00:00:00 2001 From: Nichole Barry Date: Wed, 3 Aug 2022 23:08:03 +1000 Subject: [PATCH] param transfer defaults --- .../beam_modeling/fhd_struct_init_antenna.pro | 26 +++++++++++++++++-- 1 file changed, 24 insertions(+), 2 deletions(-) diff --git a/fhd_core/beam_modeling/fhd_struct_init_antenna.pro b/fhd_core/beam_modeling/fhd_struct_init_antenna.pro index e9f242c1..f5e73a31 100644 --- a/fhd_core/beam_modeling/fhd_struct_init_antenna.pro +++ b/fhd_core/beam_modeling/fhd_struct_init_antenna.pro @@ -4,8 +4,8 @@ FUNCTION fhd_struct_init_antenna,obs,beam_model_version=beam_model_version,$ psf_dim=psf_dim,psf_max_dim=psf_max_dim,beam_offset_time=beam_offset_time,$ inst_tile_ptr=inst_tile_ptr,ra_arr=ra_arr,dec_arr=dec_arr,fractional_size=fractional_size,$ kernel_window=kernel_window,beam_per_baseline=beam_per_baseline,$ - beam_gaussian_decomp=beam_gaussian_decomp,conserve_memory=conserve_memory,$ - import_pyuvdata_beam_filepath=import_pyuvdata_beam_filepath,$ + beam_gaussian_decomp=beam_gaussian_decomp,beam_gauss_param_transfer=beam_gauss_param_transfer,$ + conserve_memory=conserve_memory,import_pyuvdata_beam_filepath=import_pyuvdata_beam_filepath,$ use_psf_resolution=use_psf_resolution,_Extra=extra t0=Systime(1) @@ -23,6 +23,28 @@ if keyword_set(beam_gaussian_decomp) and keyword_set(kernel_window) then begin kernel_window=0 endif +;Default the parameter transfer if not set +if keyword_set(beam_gauss_param_transfer) then begin + ;Default to instrumental beam if not set + if (beam_gauss_param_transfer EQ 1) then beam_gauss_param_transfer = 'decomp' + + ;Set transfer to the instrumental beam or gaussian beam. Currently only available for the MWA + if (beam_gauss_param_transfer EQ 'decomp') or (beam_gauss_param_transfer EQ 'gauss') then begin + if instrument EQ 'mwa' then begin + pointing_num = mwa_get_pointing_number(obs,string=1) + beam_gauss_param_transfer = filepath(instrument + '_decomp_params_pointing' + pointing_num + '.sav',$ + root=rootdir('FHD'),subdir='instrument_config') + endif else begin + message, 'Gaussian decomposition parameter defaults not currently set for non-MWA instruments' + endelse + endif + + ;Match the psf size to the parameters from the file + psf_transfer = getvar_savefile(beam_gauss_param_transfer,'psf') + psf_resolution = psf_transfer.resolution + psf_dim = psf_transfer.dim +endif + n_tiles=obs.n_tile n_freq=obs.n_freq n_pol=obs.n_pol