Commit e3b9bc66 authored by Joe Taylor's avatar Joe Taylor
Browse files

hopefully final set of updates before internal release

parent 3edcf3a2
......@@ -39,7 +39,7 @@ function B = ICTradModelRU_L1b_aux(sensor,T_ICT,T_CRIS,ICT_Param,ru_term,ru_valu
% .no_perturb_val
% .perturbed_val
%
% 2020-08-19 JKT
% 2020-09-13 JKT
% University of Wisconsin-Madison Space Science and Engineering Center (UW-SSEC)
wn = sensor.v_onaxis;
......
......@@ -2,7 +2,7 @@ function f_h = cris_ru_plot(RU,pname_out,fname_out);
% function f_h = cris_ru_plot(RU,pname_out,fname_out);
%
% This is a function to plot the granule mean (all FOVs, all FORs, all scans) brightness temperature,
% and 3-sigma (k=3 coverage factor) RU in the same format as figure 1 in the
% and 3-sigma (k=3 coverage factor) RU in the same 9-panel format as figure 1 in the
% NASA Cross-track Infrared Sounder (CrIS) Level 1B Radiometric Uncertainty Description Document
%
% Inputs
......@@ -13,8 +13,7 @@ function f_h = cris_ru_plot(RU,pname_out,fname_out);
% Outputs
% f_h handle to figure
%
% JKT 2020-Sep-13
% University of Wisconsin-Madison Space Science and Engineering Center (UW-SSEC)
% JKT, University of Wisconsin-Madison Space Science and Engineering Center (UW-SSEC)
rucol(1,:) = [0.8431,0.1882,0.1529];
rucol(2,:) = [0.9961,0.8784,0.5451];
......@@ -43,13 +42,14 @@ bands{1} = 'lw';
bands{2} = 'mw';
bands{3} = 'sw';
w = 12;
h = 10;
p = 0.01;
f_h = figure;
clf
set(f_h,'Units','inches');
screenpos = get(f_h,'Position');
w = 10;
h = 7.5;
p = 0.01;
set(f_h,...
'Position',[screenpos(1:2) w h],...
'PaperUnits','inches',...
......@@ -79,21 +79,18 @@ for ii = 1:length(bands)
btu.(ru_flds{iflds}) = rad2bt_l1b(v,rad1+radu.(ru_flds{iflds}))-bt1;
end
ru_ax(ii) = subaxis(3,3,ii,'ML',ML,'MR',MR,'MT',MT,'MB',MB,'SV',VS/2,'SH',HS);
ru_ax(ii) = subplot(3,3,ii);
plot(v,nanmean(bt1(:,:),2),'color',[0.1294,0.4431,0.7098],'linewidth',1)
set(gca,'FontSize',fs,'FontWeight',fw,'XLim',[min(v)-20 max(v)+20],'YLim',[180 300])
if ii == 1
ylabel('BT (K)');
else
set(gca,'YTickLabel',[]);
end
if ii == 2
title(sprintf('%s, All FOV mean',strrep(fname_out,'_','\_')));
end
set(gca,'XTickLabel',[]);
grid on
ru_ax(ii+3) = subaxis(3,3,ii+3,'ML',ML,'MR',MR,'MT',MT,'MB',MB,'SV',VS,'SH',HS);
ru_ax(ii+3) = subplot(3,3,ii+3);
hold on
for iflds = 1:length(ru_flds)
dy = btu.(ru_flds{iflds});
......@@ -103,24 +100,21 @@ for ii = 1:length(bands)
xlabel('wavenumber');
if ii == 1
ylabel('3-sigma RU (K)')
else
set(gca,'YTickLabel',[]);
end
grid on
if ii == 2
if ii == 3
pos = get(gca,'position')
h_l = legend(leg_str);
h_pos = get(h_l,'Pos');
set(h_l,'Position',[0.85,h_pos(2:4)]);
h_pos = get(h_l,'Pos')
set(h_l,'Position',[1-h_pos(3),h_pos(2:4)]);
end
ru_ax(ii+6) = subaxis(3,3,ii+6,'ML',ML,'MR',MR,'MT',MT,'MB',MB,'SV',VS,'SH',HS);
ru_ax(ii+6) = subplot(3,3,ii+6);
scatter(nanmean(bt1(:,:),2),nanmean(btu.total(:,:),2),5,v,'filled');
set(gca,'FontSize',fs,'FontWeight',fw,'YLim',[ylims_dbtscatter],'Box','on','XLim',[200 300])
xlabel('BT (K)');
if ii == 1
ylabel('3-sigma RU (K)');
else
set(gca,'YTickLabel',[]);
end
grid on
pos = get(gca,'position');
......@@ -128,7 +122,7 @@ for ii = 1:length(bands)
set(ru_ax(ii+6),'Position',pos);
set(get(hc,'XLabel'),'String','wavenumber','FontSize',8,'FontWeight','n')
cbpos = get(hc,'Pos');
set(hc,'Position',[cbpos(1) cbpos(2)+0.01 cbpos(3)-0 cbpos(4)])
set(hc,'Position',[cbpos(1) cbpos(2)+0.01 cbpos(3) cbpos(4)])
colormap(cmap_btscat)
end
print('-dpng',fullfile(pname_out,sprintf('RU_%s_%02d.png',fname_out,1)));
......
......@@ -11,9 +11,11 @@ function bt = rad2bt_l1b(freq,radiance);
% Outputs:
% bt computed brightnes temp [NCHAN x NFOV x NFOR x NSCAN] units: Kelvin
%
% JKT 2020-Sep-13
% University of Wisconsin-Madison Space Science and Engineering Center (UW-SSEC)
% JKT, University of Wisconsin-Madison Space Science and Engineering Center (UW-SSEC)
% based on rad2bt by DCT 1999-11-11
%
% not required for RU calculation, only for cris_ru_plot.m which is provided as an example of
% plotting the CrIS RU in brightness temperature units [K]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% fundamental constants:
......
......@@ -30,10 +30,12 @@
%
% Refer to the NASA CrIS L1b RU document and ATBD for details on the ICT radiometric model
% Each reflected term in the ICT radiometric model also has a view factor and emissivity, but the uncertainty
% in those quantities produce negligible contributions to the CrIS radiometric uncertainty and are not perturbed for the RU estimate
% J1 has an improved ICT emissivity, and the view factors of some of the reflected terms are set to zero as a result
% in those quantities produce negligible contributions to the CrIS radiometric uncertainty and are not
% perturbed for the RU estimate J1 has an improved ICT emissivity, and the view factors of some of the
% reflected terms are set to zero as a result
%
%
% JKT 2020-Sep-13
% JKT 2022-Mar-11
% University of Wisconsin-Madison Space Science and Engineering Center (UW-SSEC)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
......@@ -41,9 +43,13 @@
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% START: USER INPUT
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Example 1: SNPP
%sat_name = 'SNPP'; % set to 'J1' or 'SNPP'
%l1b_fname = 'SNDR.SNPP.CRIS.20180401T1306.m06.g132.L1B.std.v03_08.G.201231182944.nc'; % set to file name of v3 L1b granule
%l1b_fname = 'SNDR.SNPP.CRIS.20180401T2206.m06.g222.L1B.std.v03_08.G.201231185850.nc'; % set to file name of v3 L1b granule
%% Example 1: NOAA-20/JPSS-1
sat_name = 'J1'; % set to 'J1' or 'SNPP'
l1b_fname = 'SNDR.J1.CRIS.20180401T2254.m06.g230.L1B.std.v03_08.G.201214201329.nc'; % set to file name of v3 L1b granule
l1b_pname = '/Users/joet/git/cris_l1b_user/RU/TEST_INPUT'; % set to path to v3 L1b granule
......
......@@ -11,7 +11,7 @@ function RU = cris_gran_RU_ncparam(sat_name,l1b_fname,saveTotalRU_FLAG,pname_out
% saveTotalRU_FLAG set to 1 if you want to save the RSS RU to a mat file
% pname_out output directory
%
% JKT 2021-Aug-26
% JKT 2022-Mar-02
% University of Wisconsin-Madison Space Science and Engineering Center (UW-SSEC)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
......
......@@ -20,7 +20,7 @@ function nlc_ru = estimate_nlc_ru_L1bonly(sensor,user,l1b_dat,ICT_Param,da2);
% Outputs
% nlc_ru radiometric uncertainty associated with quadratic coefficient (a2) uncertainty
%
% 2020-08-19 JKT
% 2022-03-02 JKT
% University of Wisconsin-Madison Space Science and Engineering Center (UW-SSEC)
PGA_GainALL.lw = l1b_dat.aux.pga_gain_lw;
......
......@@ -55,7 +55,7 @@ function d = read_l1b_for_RU(fname);
% .lat [NFOV x NFOR x NSCAN]
% .lon [NFOV x NFOR x NSCAN]
%
% 2020-08-19 JKT
% 2022-Mar-02 JKT
% University of Wisconsin-Madison Space Science and Engineering Center (UW-SSEC)
%% aux vars
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment