[go: up one dir, main page]
Menu ▾ ▴

[57f2a1]: / anal / GUISPERT.m  Maximize  Restore  History

Download this file

96 lines (90 with data), 3.1 kB 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
% GUISPERT.m: special hacks

% GUISDAP v8.1  03-04-30 Copyright EISCAT

%

% script for GUISPERTs and other power users. Put here any commands by

% which you wish to change the data before the analysis is started

%

% See also: GUIZARD GUIDITOR



% Correct the doppler sign

if (strcmp(name_ant(1:3),'vhf') & d_time(1,1)>1999) | ...

   (strcmp(name_ant(2:3),'2m') & (d_date>datenum(2003,3,21,8,50,0) & d_date<datenum(2003,10,25) | d_date>datenum(2009,10,1))) | ...

   (strcmp(name_ant(1:3),'kir') & d_date>datenum(2009,11,23) & d_date<datenum(2009,11,23,18,42,0)) | ...

   (strcmp(name_ant(1:3),'kir') & d_date>datenum(2013,1,1)) | ...

   (strcmp(name_ant(1:3),'sod') & d_date>datenum(2013,1,1))

  d_data=conj(d_data);

end

if ~exist('Magic_const')

 if strcmp(name_ant(1:3),'vhf')

  if d_date<datenum(2007,8,14)

   a_Magic_const=1.367;

  end

 elseif strcmp(name_ant(1:3),'uhf')

  if d_time(1,1)>2004

   a_Magic_const=0.80;

   if d_date>datenum(2009,10,26) % polariser back

    a_Magic_const=1.14;  %MTR

   end

   if d_date>datenum(2010,02,02) % polariser calibrated against Kiruna

    a_Magic_const=0.70;  %MTR

   end

  end

 elseif strcmp(name_ant(1:3),'42m')

  if d_date<datenum(2007,3,1)

   a_Magic_const=1.5;

  elseif d_date<datenum(2007,3,18)

   a_Magic_const=1.87;

  elseif d_date<datenum(2007,5,21)

   a_Magic_const=2.45; %AW

  elseif d_date<datenum(2007,10,15)

   a_Magic_const=1.59; %after snow clearence IH

  elseif d_date<datenum(2009,05,01)

   a_Magic_const=2.0;

  elseif d_date<datenum(2009,06,01)

   a_Magic_const=4.0;

  else

   a_Magic_const=1.5;

  end

 end

end

if strcmp(name_ant(1:3),'42m') & ~exist('ad_coeff_no_Pt','var')

 d=find(ch_gain~=ch_gain(1));

 % Only transmission affected for tx on 42m and rx on 32m

 ch_gain(d)=sqrt(a_Magic_const)*ch_gain(d);

end



% Range not always recorded

if name_site=='K' | name_site=='S'

 if isempty(ch_range)

  global v_lightspeed

  if exist('ch_height','var')

    ch_range=height_to_range(ch_height,ch_el(1))*(v_lightspeed*p_dtau*1e-6/2/1e3);

  else

    ch_range=height_to_range(300,ch_el(1))*(v_lightspeed*p_dtau*1e-6/2/1e3);

    warning('GUISDAP:default','range is dummy!')

  end

 end

 if d_date>datenum(2013,1,1)

  ch_fradar(:)=radar_freqs(3);

 end

end



% Correct transmitter power when wrongly recorded

%d_secs=tosecs(d_time(1,:));

%if (name_site=='T' & d_time(1,1)==2001 & d_secs>27334800 & d_secs<25660800) 

% After 09:38 UT 9 Nov 2001 

%  Pt=Pt*2

%end



if name_site=='V' & d_date>datenum(2003,1,1)

 [ch_el ch_az ch_gain]=vhf_elaz(ch_el(1),0,10^4.31/2);

end



if strcmp(name_ant(1:3),'uhf')

 if (d_date>=datenum(2000,10,01) & d_date<=datenum(2001,02,08)) | ...

    (d_date>=datenum(2001,06,28) & d_date<=datenum(2001,10,15,13,31,00)) 

     fprintf('GUISPERT: doubling uhf tx power estimate\n')

     ch_Pt=ch_Pt*2;

 elseif d_date>=datenum(2001,4,2) & d_date<=datenum(2001,4,10) 

% 2-10 April 2001, the UHF transmitter ran with only one klystron, but this was

% not reflected in the power values recorded in the parameter block.

% Should only affect the UP run on April 3-4.

  fprintf('GUISPERT: halving uhf tx power estimate\n')

  ch_Pt=ch_Pt/2;

 end

end