/[MITgcm]/MITgcm_contrib/ocean_inversion_project/README
ViewVC logotype

Diff of /MITgcm_contrib/ocean_inversion_project/README

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

revision 1.3 by dimitri, Tue Sep 23 04:34:23 2003 UTC revision 1.4 by dimitri, Wed Sep 24 04:52:38 2003 UTC
# Line 35  salinity initial and boundary conditions Line 35  salinity initial and boundary conditions
35   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
36   cp ../ocean_inversion_project/code/.genmakerc .   cp ../ocean_inversion_project/code/.genmakerc .
37   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
38   cp ptracers_forcing.F.test ptracers_forcing.F   rm ptracers_init.F ptracers_forcing_surf.F PTRACERS.h
39   rm ptracers_init.F ptracers_forcing.F   rm ptracers_read_mask.F ptracers_read_takahashi.F
40   ../tools/genmake   ../tools/genmake
41   make depend   make depend
42   make   make
# Line 51  salinity initial and boundary conditions Line 51  salinity initial and boundary conditions
51   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
52   cp data.test data   cp data.test data
53   cp data.ptracers.test data.ptracers   cp data.ptracers.test data.ptracers
54   mitgcmuv > output .txt   mitgcmuv > output.txt
55        
56  5 ===> check that PTRACER output for tracer 1 and tracer 30  5 ===> check that PTRACER output and salinity output
57    ===> is identical to salinity output    ===> are identical.
58    
59   diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data   diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data
60   diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data   diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data
# Line 102  clf, contourf(lon,lat,mean(tak(:,:,1),3) Line 102  clf, contourf(lon,lat,mean(tak(:,:,1),3)
102  caxis([-6 6]), colorbar, plotland  caxis([-6 6]), colorbar, plotland
103    
104  % some matlab code for looking at PTRACER output files  % some matlab code for looking at PTRACER output files
105    salt=readbin('S.0000000020.001.001.data',[90 40 15],1);
106  tracer=zeros(90,40,15,30);  tracer=zeros(90,40,15,30);
107  for i=1:30  for i=1:30
108   fn=['PTRACER' myint2str(i) '.0000000020.001.001.data'];   fn=['PTRACER' myint2str(i) '.0000000020.001.001.data'];
# Line 109  for i=1:30 Line 110  for i=1:30
110   clf, mypcolor(tracer(:,:,1,i)'); pause(1)   clf, mypcolor(tracer(:,:,1,i)'); pause(1)
111  end  end
112  tmp=sum(tracer,4); clf, mypcolor(tmp(:,:,1)')  tmp=sum(tracer,4); clf, mypcolor(tmp(:,:,1)')
113    mypcolor(
114    
115    
116  % some matlab code for checking that one year's worth of  % some matlab code for checking that one year's worth of
117  % tracer uptake is approximately 1e18 mols.  % tracer uptake is approximately 1e18 mols.
# Line 134  lon=2:4:360; lat=-78:4:78; Line 137  lon=2:4:360; lat=-78:4:78;
137  thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];  thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];
138  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
139  for i=1:30, mydisp(i)  for i=1:30, mydisp(i)
140   fn=['PTRACER' myint2str(i) '.0000000180.001.001.data'];   fn=['PTRACER' myint2str(i) '.0000003600.001.001.data'];
141   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
142   for j=1:length(lat)   for j=1:length(lat)
143    for k=1:length(thk)    for k=1:length(thk)
# Line 143  for i=1:30, mydisp(i) Line 146  for i=1:30, mydisp(i)
146    end    end
147   end   end
148  end  end
149    plot(1:30,0*sumtracer,1:30,sumtracer)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22