clear,
cd /home/disk/rocinante/DATA/WRF/GFDL_CM3_RCP85/sfc_vars/6hrly
plotdir='~/WRF/PLOTs/diagnostics/GFDL_CM3_RCP85';

filelist=dir('wrfout*');
varlist={'T2','PSFC','RAINC','RAINNC','RAINSH','U10','V10'};

ncload('wrfout_d02_1970-03-14_00:00:00','XLAT','XLONG'),


for vi=1:length(varlist)
  tmpdat=NaN([size(XLAT) length(filelist)]);

  for fi=1:length(filelist)

    ncload(filelist(fi).name,varlist{vi}),

    tmpdat(:,:,fi) = eval(varlist{vi});
    eval(['clear ' varlist{vi}]),
  end
  clear fi,

  eval([varlist{vi} ' = tmpdat;']),
  clear tmpdat,
end
clear vi ans,

%
load coast
coast_lon=long;
coast_lat=lat;
clear long lat ans,

ti=1;
for vi=1:length(varlist)
  eval(['tmpdat=squeeze(' varlist{vi} '(:,:,ti));']),

  if length(findstr('RAIN',varlist{vi}))
    eval(['tmpdat2=squeeze(' varlist{vi} '(:,:,ti+1));']),
    tmpdat=tmpdat2-tmpdat;
    clear tmpdat2
  figure(vi), clf,
hpc=pcolor(XLONG,XLAT,squeeze(PSFC(:,:,1))/100);
set(hpc,'edgealpha',0),
hold on,
hpp=plot(coast_lon,coast_lat,'k-');
set(hpp,'color',[1 1 1]*.5);
hold off,
ht=title(varlist{vi});
hc=colorbar('h');

print('-dpng',[plotdir '/map_PSFC_1970-03-14_00:00:00.png']),

