Documentation of flux_gr_ice

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

[swl1, lwl1, lhl1, shl1] = getnc('clim_5093.nc', v1, v2, v3, v4);

Cross-Reference Information

Listing of script flux_gr_ice

clear
v1 = 'FLNS'; v2 = 'FSNS'; v3 = 'LHFLX'; v4 = 'SHFLX'; v5 = 'TS';
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/cdtemp
[sww1, lww1, lhw1, shw1] = getnc('time_wgr_5093.nc', v1, v2, v3, v4);
[swc1, lwc1, lhc1, shc1] = getnc('time_cgr_5093.nc', v1, v2, v3, v4);
[lat, lon] = getll('clim_5093.nc');
get_global; default_global;

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[sww2, lww2, lhw2, shw2] = getnc('time_wgr_hard.nc', v1, v2, v3, v4);
[swc2, lwc2, lhc2, shc2] = getnc('time_cgr_hard.nc', v1, v2, v3, v4);
%[swl2, lwl2, lhl2, shl2] = getnc('time_clim_hard.nc', v1, v2, v3, v4);

tim = sort([1:12:120 2:12:120 3:12:120 4:12:120 5:12:120 6:12:120]);
FRAME = [30 120 -45 30];

dw = lww1 + lhw1 + shw1 - sww1;
dc = lwc1 + lhc1 + shc1 - swc1;
sw = lww2 + lhw2 + shw2 - sww2;
sc = lwc2 + lhc2 + shc2 - swc2;

lab = 'Total Surface Flux';
cint = 2.5; clev = [-60:cint:60];
i = 2;

if i == 1;
  var1 = squeeze(dw(tim,:,:));
  var2 = squeeze(dc(tim,:,:));
elseif i == 2;
  var1 = squeeze(sw(tim,:,:));
  var2 = squeeze(dw(tim,:,:));
end

[tstat, dof] = ttest(var1, var2);
score = tscore(dof, 2.5);
top = squeeze(mean(var1 - var2))/(3-i);
sp(1);
  clma;
  mcont(top, clev, 'giso', [0 180]);
  greyshd(abs(tstat), score);
  if i == 1;
    title(['NDJFMA ' lab ':  (WGR\_DOM - CGR\_DOM) / 2']);
  elseif i == 2;
    title(['NDJFMA ' lab ':  (WGR\_SOM - WGR\_DOM)']);
  end
  xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);

if i == 1;
  var1 = squeeze(sw(tim,:,:));
  var2 = squeeze(sc(tim,:,:));
elseif i == 2;
  var1 = squeeze(sc(tim,:,:));
  var2 = squeeze(dc(tim,:,:));
end

[tstat, dof] = ttest(var1, var2);
score = tscore(dof, 2.5);
top = squeeze(mean(var1 - var2))/(3-i);
sp(2);
  clma;
  mcont(top, clev, 'giso', [0 180]);
  greyshd(abs(tstat), score);
  if i == 1;
    title(['NDJFMA ' lab ':  (WGR\_SOM - CGR\_SOM) / 2']);
  elseif i == 2;
    title(['NDJFMA ' lab ':  (CGR\_SOM - CGR\_DOM)']);
  end
  xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/DOM_SOM
% print -dps2