matlab_class/gcmfaces_misc/check_snap.m


listBudgs={'budgMo','budgHo','budgSo','budgMi','budgHi','budgSi'};
dt=(0000174564-0000000732)*3600;
for ii=1:length(listBudgs);
  budgName=listBudgs{ii};
  fil=['r4it11.c65i/nctiles_budg/' budgName '/tend/tend'];
  ncload([fil '.0001.nc'],'t0'); ncload([fil '.0001.nc'],'t1');
  tend=0*mygrid.hFacC;
  if ii>3; tend=tend(:,:,1); end;

  for tt=1:length(t0);
    tmp1=read_nctiles(['r4it11.c65i/nctiles_budg/' budgName '/tend/tend'],'tend',tt);  
    tend=tend+(t1(tt)-t0(tt))*tmp1;
  end;
  ini=read_nctiles(['r4it11.c65i/nctiles_budg/' budgName '/initial/snapshot'],'snapshot');
  fin=read_nctiles(['r4it11.c65i/nctiles_budg/' budgName '/final/snapshot'],'snapshot');

  test0=nanstd((fin-ini)-tend)/nanstd(tend);
  fprintf('%s : %3.2g \n',budgName,test0);

end;


1	gforget	1.1
2			listBudgs={'budgMo','budgHo','budgSo','budgMi','budgHi','budgSi'};
3			dt=(0000174564-0000000732)*3600;
4			for ii=1:length(listBudgs);
5			budgName=listBudgs{ii};
6			fil=['r4it11.c65i/nctiles_budg/' budgName '/tend/tend'];
7			ncload([fil '.0001.nc'],'t0'); ncload([fil '.0001.nc'],'t1');
8			tend=0*mygrid.hFacC;
9			if ii>3; tend=tend(:,:,1); end;
10
11			for tt=1:length(t0);
12			tmp1=read_nctiles(['r4it11.c65i/nctiles_budg/' budgName '/tend/tend'],'tend',tt);
13			tend=tend+(t1(tt)-t0(tt))*tmp1;
14			end;
15			ini=read_nctiles(['r4it11.c65i/nctiles_budg/' budgName '/initial/snapshot'],'snapshot');
16			fin=read_nctiles(['r4it11.c65i/nctiles_budg/' budgName '/final/snapshot'],'snapshot');
17
18			test0=nanstd((fin-ini)-tend)/nanstd(tend);
19			fprintf('%s : %3.2g \n',budgName,test0);
20
21			end;
22
23