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
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