Documentation of make_annual_heatflux

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

Help text

  MIX data

Cross-Reference Information

Listing of script make_annual_heatflux

clean
lims = [-0.1 360 -90 90];
cd /home/disk/hayes2/dvimont/csiro/data
nc = netcdf('shfl_qm1_0036-0425.nc', 'nowrite');
  [lat, lon, yk, xk] = get_nclatlon2(lims, nc);
  hfl = nc{'hfl', 1}(1:330,:,:,:);
nc = close(nc);
[ntim, nmo, nlat, nlon] = size(hfl);
hfl2 = repmat(NaN, [ntim nlat nlon]);
for i = 1:ntim;
  hfl2(i,:,:) = mean(squeeze(hfl(i,:,:,:)));
end
hfl = hfl2; clear hfl2;

nc = netcdf('sevp_qm1_0036-0425.nc', 'nowrite');
  evp = nc{'evp', 1}(1:330,:,:,:);
nc = close(nc);
evp2 = repmat(NaN, [ntim nlat nlon]);
for i = 1:ntim;
  evp2(i,:,:) = mean(squeeze(evp(i,:,:,:)));
end
evp = evp2; clear evp2;
evp = evp*2.5e6/(24*3600);

nc = netcdf('srgn_qm1_0036-0425.nc', 'nowrite');
  rgn = nc{'rgn', 1}(1:330,:,:,:);
nc = close(nc);
rgn2 = repmat(NaN, [ntim nlat nlon]);
for i = 1:ntim;
  rgn2(i,:,:) = mean(squeeze(rgn(i,:,:,:)));
end
rgn = rgn2; clear rgn2;

nc = netcdf('ssgn_qm1_0036-0425.nc', 'nowrite');
  sgn = nc{'sgn', 1}(1:330,:,:,:);
nc = close(nc);
sgn2 = repmat(NaN, [ntim nlat nlon]);
for i = 1:ntim;
  sgn2(i,:,:) = mean(squeeze(sgn(i,:,:,:)));
end
sgn = sgn2; clear sgn2;

hflx = evp+hfl+rgn-sgn;

cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Flux
save MIX_HFLX_tot_ANN.mat hfl evp rgn sgn hflx lat lon lims




%  Keep monthly resolution

clean

lims = [89 306 -80 80];

%  MIX data

cd /home/disk/hayes2/dvimont/csiro/data
nc = netcdf('shfl_qm1_0036-0425.nc', 'nowrite');
  [lat, lon, yk, xk] = get_nclatlon2(lims, nc);
  hfl = nc{'hfl', 1}(1:330,:,yk,xk);
nc = close(nc);
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Flux
save MIX_HFL_pac_mon.mat hfl lat lon lims
hflx = hfl;
clear hfl;

cd /home/disk/hayes2/dvimont/csiro/data
nc = netcdf('sevp_qm1_0036-0425.nc', 'nowrite');
   evp = nc{'evp', 1}(1:330,:,yk,xk);
nc = close(nc);
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Flux
evp = evp*2.5e6/(24*3600);
save MIX_EVP_pac_mon.mat evp lat lon lims
hflx = hflx+evp;
clear evp;

cd /home/disk/hayes2/dvimont/csiro/data
nc = netcdf('srgn_qm1_0036-0425.nc', 'nowrite');
   rgn = nc{'rgn', 1}(1:330,:,yk,xk);
nc = close(nc);
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Flux
save MIX_RGN_pac_mon.mat rgn lat lon lims
hflx = hflx+rgn;
clear rgn;

cd /home/disk/hayes2/dvimont/csiro/data
nc = netcdf('ssgn_qm1_0036-0425.nc', 'nowrite');
  sgn = nc{'sgn', 1}(1:330,:,yk,xk);
nc = close(nc);
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Flux
save MIX_SGN_pac_mon.mat sgn lat lon lims
hflx = hflx-sgn;
clear sgn;

cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Flux
save MIX_HFLX_pac_mon.mat hflx lat lon lims