/[MITgcm]/MITgcm_contrib/gael/matlab_class/gcmfaces_misc/check_snap.m
ViewVC logotype

Annotation of /MITgcm_contrib/gael/matlab_class/gcmfaces_misc/check_snap.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (hide annotations) (download)
Wed Feb 4 14:12:43 2015 UTC (10 years, 5 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint65x, checkpoint65r, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, HEAD
- check_budg.m : switch back to upward convention, and clean-up
- check_snap.m : double check initial and final snapshots

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    

  ViewVC Help
Powered by ViewVC 1.1.22