/[MITgcm]/MITgcm_contrib/high_res_cube/matlab/v_integrate.m
ViewVC logotype

Annotation of /MITgcm_contrib/high_res_cube/matlab/v_integrate.m

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


Revision 1.1 - (hide annotations) (download)
Wed Mar 28 18:07:27 2007 UTC (18 years, 3 months ago) by dimitri
Branch: MAIN
CVS Tags: HEAD
Files contributed by JM for computation of barotropic streamfunction on cs510.

1 dimitri 1.1 krd=3; kgr=0;
2    
3     % $Header: /u/gcmpack/MITgcm/utils/matlab/cs_grid/bk_line/use_psiLine.m,v 1.1 2005/09/15 16:46:28 jmc Exp $
4     % $Name: $
5     listI=[02232 04320 06552 08712 10944 13104 15336 17568 19728 21960 24120 26352 ...
6     28584 30600 32832 34992 37224 39384 41616 43848 46008 48240 50400 52632 ...
7     54864 56880 59112 61272 63504 65664 67896];
8     listR=listI(end-11:end);
9     deltaT=listI(end-12:end-1); deltaT=listR-deltaT;
10    
11     Rac='grid_files/';
12     drF=squeeze(rdmds([Rac,'DRF']))'; nr=length(drF);
13    
14     if rem(krd,2) == 1,
15     %- load the grid dx,dy , convert to 10^6 m :
16     dyg=rdmds([Rac,'DYG']); dyg=dyg*1.e-6;
17     fprintf(' load dyg file\n');
18     nc=size(dyg,2); nPg=nc*nc*6;
19     dyg=reshape(dyg,[nPg 1]);
20     end
21    
22     if rem(krd,2) == 1,
23     %-- load velocity :
24     rac='output_files/';
25     %nit=28584;
26     %uu=rdmds([rac,'UVELMASS'],nit);
27     ut=zeros(nPg,1);
28     for n=1:length(listR),
29     nit=listR(n);
30     namf=sprintf([rac,'UVELMASS.%10.10i.data'],nit);
31     fprintf([' process file: ',namf,' ...']);
32     for k=1:nr,
33     fprintf(' %i',k);
34     uk=rdda(namf,[nPg 1],k,'real*4','b');
35     uk=dyg.*uk;
36     ut=ut+drF(k)*deltaT(n)*uk;
37     end
38     fprintf(' done\n');
39     end
40     ut=ut/sum(deltaT);
41     namfil='UintTrsYr.bin';
42     fprintf([' write to file: ',namfil,' ...']);
43     fid=fopen(namfil,'w','b'); fwrite(fid,ut,'real*8'); fclose(fid);
44     fprintf(' done\n');
45     end
46    
47     if krd > 1,
48     clear dyg uk ut
49     dxg=rdmds([Rac,'DXG']); dxg=dxg*1.e-6;
50     fprintf(' load dxg file\n');
51     nc=size(dxg,2); nPg=nc*nc*6;
52     dxg=reshape(dxg,[nPg 1]);
53     end
54    
55     if krd > 1,
56     %-- load velocity :
57     rac='output_files/';
58     %nit=28584;
59     %vv=rdmds([rac,'VVELMASS'],nit);
60     vt=zeros(nPg,1);
61     for n=1:length(listR),
62     nit=listR(n);
63     namf=sprintf([rac,'VVELMASS.%10.10i.data'],nit);
64     fprintf([' process file: ',namf,' ...']);
65     for k=1:nr,
66     fprintf(' %i',k);
67     vk=rdda(namf,[nPg 1],k,'real*4','b');
68     vk=dxg.*vk;
69     vt=vt+drF(k)*deltaT(n)*vk;
70     end
71     fprintf(' done\n');
72     end
73     vt=vt/sum(deltaT);
74     namfil='VintTrsYr.bin';
75     fprintf([' write to file: ',namfil,' ...']);
76     fid=fopen(namfil,'w','b'); fwrite(fid,vt,'real*8'); fclose(fid);
77     fprintf(' done\n');
78     end
79    
80     return

  ViewVC Help
Powered by ViewVC 1.1.22