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

[73c35b]: / exps / cp1l / guispert.m  Maximize  Restore  History

Download this file

46 lines (42 with data), 1.4 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
% guispert.m: special experiment specific hacks

% GUISDAP v1.70   01-11-30 Copyright EISCAT

%

% See also: GUISPERT GUIZARD

%

global p_rep

ch_Pt=ch_Pt(1);

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

  p_rep=312480;

  ch_height=292.9;

end

if d_date>datenum(2000,09,06,12,0,0) & d_date<datenum(2000,09,06,23,0,0)

 polhv=[1.7980 -178.1795 4508.4];

 %read hv from txlog

 if ~exist('txlog','var')

  txlog=get_tx('U0900.txt',[52]);

  d=find(txlog(:,2)<25); txlog(d,2)=25; %ch_Pt<0

 end

 hv=interp1(txlog(:,1),txlog(:,2),mean(datenum(d_time)),'nearest');

 ch_Pt=polyval(polhv,hv)*1000/2;

 fprintf('GUISPERT: HV=%g kV -> TX=%.1f MW\n',hv,ch_Pt/1e6)

 d=find(analysis_altit<270); analysis_altit(d)=[]; analysis_maxwidth(d)=[];

end



if name_site=='K' && d_date>datenum(2001,04,02) & d_date<datenum(2001,06,25)

 ch_Pt=ch_Pt(1)/2.0;  % one klystron but 2 assumed.

end



if (d_date>=datenum(2012,10,31) & d_date<datenum(2012,10,31,23,59,0))

  f_ch_Pt=1020e3;

  fprintf('GUISPERT: fixing uhf tx power estimate 1')

  if (d_date>=datenum(2012,10,31,12,50,0)) 

    f_ch_Pt=430e3;

    fprintf('GUISPERT: fixing uhf tx power estimate 2')

  end

  if (d_date>=datenum(2012,10,31,13,44,0))

    f_ch_Pt=1030e3;

    fprintf('GUISPERT: fixing uhf tx power estimate 3')

  end

  if (d_date>=datenum(2012,10,31,15,03,0))

    f_ch_Pt=1150e3;

    fprintf('GUISPERT: fixing uhf tx power estimate 4')

  end  

end