matlab_class/gcmfaces_misc/check_loop.m


dirIn='release1/nctiles_budg/budgMo/'

for tt=1:238;
  tt
  budgIn.tend=read_nctiles([dirIn 'tend/tend'],'tend',tt);
  budgIn.trU=read_nctiles([dirIn 'trU/trU'],'trU',tt);
  budgIn.trV=read_nctiles([dirIn 'trV/trV'],'trV',tt);
  if dirIn(end-1)=='o';
    budgIn.trW=read_nctiles([dirIn 'trW/trW'],'trW',tt);
    %
    budgIn.trWtop=budgIn.trW;
    budgIn.trWbot=budgIn.trW(:,:,2:50);
    budgIn.trWbot(:,:,50)=0;
    %
    %budgIn.tend(isnan(tend))=0;
    %budgIn.trU(isnan(trU))=0; budgIn.trV(isnan(trV))=0;
    %budgIn.trWtop(isnan(trWtop))=0; budgIn.trWbot(isnan(trWbot))=0;
  else;
    budgIn.trWtop=read_nctiles([dirIn 'trWtop/trWtop'],'trWtop',tt);
    budgIn.trWbot=read_nctiles([dirIn 'trWbot/trWbot'],'trWbot',tt);
  end;
  %
  nr=size(budgIn.tend{1},3);
  for kk=1:nr; prec(kk,:)=check_budg(budgIn,kk); end;
  if 0;
    figureL;
    plot(log10(prec(:,1)),'b');
    hold on;
    plot(log10(prec(:,2)),'r');
    plot(log10(prec(:,3)),'k');
    title(num2str(tt));
  end;%if 0;
  if tt==1;
    store_prec=prec;
  else;
    store_prec(:,:,tt)=prec;
  end;
end;%for tt=[1 119 238];

max(store_prec(:))
eval(['save ' dirIn 'residuals.mat store_prec;']);

1	gforget	1.1
2	gforget	1.2	dirIn='release1/nctiles_budg/budgMo/'
3	gforget	1.1
4			for tt=1:238;
5			tt
6			budgIn.tend=read_nctiles([dirIn 'tend/tend'],'tend',tt);
7			budgIn.trU=read_nctiles([dirIn 'trU/trU'],'trU',tt);
8			budgIn.trV=read_nctiles([dirIn 'trV/trV'],'trV',tt);
9			if dirIn(end-1)=='o';
10			budgIn.trW=read_nctiles([dirIn 'trW/trW'],'trW',tt);
11			%
12			budgIn.trWtop=budgIn.trW;
13			budgIn.trWbot=budgIn.trW(:,:,2:50);
14			budgIn.trWbot(:,:,50)=0;
15			%
16			%budgIn.tend(isnan(tend))=0;
17			%budgIn.trU(isnan(trU))=0; budgIn.trV(isnan(trV))=0;
18			%budgIn.trWtop(isnan(trWtop))=0; budgIn.trWbot(isnan(trWbot))=0;
19			else;
20			budgIn.trWtop=read_nctiles([dirIn 'trWtop/trWtop'],'trWtop',tt);
21			budgIn.trWbot=read_nctiles([dirIn 'trWbot/trWbot'],'trWbot',tt);
22			end;
23			%
24			nr=size(budgIn.tend{1},3);
25			for kk=1:nr; prec(kk,:)=check_budg(budgIn,kk); end;
26			if 0;
27			figureL;
28			plot(log10(prec(:,1)),'b');
29			hold on;
30			plot(log10(prec(:,2)),'r');
31			plot(log10(prec(:,3)),'k');
32			title(num2str(tt));
33			end;%if 0;
34			if tt==1;
35			store_prec=prec;
36			else;
37			store_prec(:,:,tt)=prec;
38			end;
39			end;%for tt=[1 119 238];
40
41			max(store_prec(:))
42			eval(['save ' dirIn 'residuals.mat store_prec;']);
43