/[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.11 by dimitri, Mon Oct 20 06:12:24 2003 UTC
# Line 22  salinity initial and boundary conditions Line 22  salinity initial and boundary conditions
22    
23   CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/gcmpack   CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/gcmpack
24   cvs login ( CVS password: cvsanon )   cvs login ( CVS password: cvsanon )
25   cvs co MITgcm   cvs co -P MITgcm
26    
27  2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory  2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory
28    
# Line 33  salinity initial and boundary conditions Line 33  salinity initial and boundary conditions
33    
34   cd bin   cd bin
35   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
  cp ../ocean_inversion_project/code/.genmakerc .  
36   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
37   cp ptracers_forcing.F.test ptracers_forcing.F   rm ptracers_*.F
38   rm ptracers_init.F ptracers_forcing.F   ../tools/genmake2
39   ../tools/genmake     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
40   make depend   make depend
41   make   make
42    
# Line 51  salinity initial and boundary conditions Line 50  salinity initial and boundary conditions
50   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
51   cp data.test data   cp data.test data
52   cp data.ptracers.test data.ptracers   cp data.ptracers.test data.ptracers
53   mitgcmuv > output .txt   mitgcmuv > output.txt
54        
55  5 ===> check that PTRACER output for tracer 1 and tracer 30  5 ===> check that PTRACER output and salinity output are identical.
   ===> is identical to salinity output  
56    
57   diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data   diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data
58   diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data   diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data
# Line 74  integration using the global_ocean.90x40 Line 72  integration using the global_ocean.90x40
72   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
73   cp ../ocean_inversion_project/code/.genmakerc .   cp ../ocean_inversion_project/code/.genmakerc .
74   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
75   ../tools/genmake   ../tools/genmake2
76       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
77   make depend   make depend
78   make   make
79   cd ../exe   cd ../exe
# Line 88  integration using the global_ocean.90x40 Line 87  integration using the global_ocean.90x40
87   mitgcmuv > output .txt &   mitgcmuv > output .txt &
88    
89    
90    ==========================================================
91    MPI instructions for carrying out a 3000-year quasi-stationary
92    integration using the global_ocean.90x40x15 configuration.
93    ==========================================================
94    
95    1 ===> compile, link, and execute
96    
97     cd MITgcm/exe
98     rm *
99     cd ../bin
100     rm *
101     cp ../verification/global_with_exf/code/* .
102     cp ../ocean_inversion_project/code/.genmakerc .
103     cp ../ocean_inversion_project/code/* .
104     cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
105     cp SIZE.h_mpi SIZE.h
106     ../tools/genmake2
107       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
108     make depend
109     make
110     cd ../exe
111     cp ../verification/global_with_exf/input/eedata .
112     cp ../verification/global_with_exf/input/data.* .
113     cp ../verification/global_with_exf/input/POLY3.COEFFS .
114     ln -sf ../verification/global_with_exf/input/*.bin .
115     cp ../ocean_inversion_project/input/* .
116     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
117     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
118     mpirun -np 36 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog &
119    
120    
121  ======================================================  ======================================================
122    
123  % some matlab code for looking at fort.10 debug files  % some matlab code for looking at fort.10 debug files
# Line 102  clf, contourf(lon,lat,mean(tak(:,:,1),3) Line 132  clf, contourf(lon,lat,mean(tak(:,:,1),3)
132  caxis([-6 6]), colorbar, plotland  caxis([-6 6]), colorbar, plotland
133    
134  % some matlab code for looking at PTRACER output files  % some matlab code for looking at PTRACER output files
135    salt=readbin('S.0000000020.001.001.data',[90 40 15],1);
136  tracer=zeros(90,40,15,30);  tracer=zeros(90,40,15,30);
137  for i=1:30  for i=1:30
138   fn=['PTRACER' myint2str(i) '.0000000020.001.001.data'];   fn=['PTRACER' myint2str(i) '.0000000020.001.001.data'];
# Line 109  for i=1:30 Line 140  for i=1:30
140   clf, mypcolor(tracer(:,:,1,i)'); pause(1)   clf, mypcolor(tracer(:,:,1,i)'); pause(1)
141  end  end
142  tmp=sum(tracer,4); clf, mypcolor(tmp(:,:,1)')  tmp=sum(tracer,4); clf, mypcolor(tmp(:,:,1)')
143    mypcolor(
144    
 % some matlab code for checking that one year's worth of  
 % tracer uptake is approximately 1e18 mols.  
 lon=2:4:360; lat=-78:4:78;  
 thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];  
 tracer=zeros(90,40,15,30); sumtracer1=zeros(30,1);  
 for i=1:30, mydisp(i)  
  fn=['PTRACER' myint2str(i) '.0000000180.001.001.data'];  
  tracer(:,:,:,i)=readbin(fn,[90 40 15],1);  
  for x=1:length(lon)  
   for y=1:length(lat)  
    for z=1:length(thk)  
     sumtracer1(i) = sumtracer1(i) + tracer(x,y,z,i) * ...  
                     thk(z) * (4*1.113195e+05)^2 * cos(pi*lat(y)/180);  
    end  
   end  
  end  
 end  
145    
146  % some matlab code for checking that one year's worth of  % some matlab code for checking that one year's worth of
147  % tracer uptake is approximately 1e18 mols.  % tracer uptake is approximately 1e18 mols.
148  lon=2:4:360; lat=-78:4:78;  lon=2:4:360; lat=-78:4:78;
149  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];
150    mask=readbin('hFacC.001.001.data',[90 40 15],1);
151  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
152  for i=1:30, mydisp(i)  for i=1:30, mydisp(i)
153   fn=['PTRACER' myint2str(i) '.0000000180.001.001.data'];   fn=['PTRACER' myint2str(i) '.0000001800.001.001.data'];
154   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
155   for j=1:length(lat)   for j=1:length(lat)
156    for k=1:length(thk)    for k=1:length(thk)
157     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i)) * ...     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i).*mask(:,j,k)) * ...
158                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
159    end    end
160   end   end
161  end  end
162    plot(1:30,0*sumtracer,1:30,sumtracer)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22