/[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.2 by dimitri, Thu Sep 18 02:33:38 2003 UTC revision 1.8 by dimitri, Sat Sep 27 07:45:51 2003 UTC
# Line 3  Ocean Inversion Project Line 3  Ocean Inversion Project
3    
4  Instructions for using pkg/ptracers to compute tracer  Instructions for using pkg/ptracers to compute tracer
5  Green's functions for Gruber's ocean inversion project  Green's functions for Gruber's ocean inversion project
6  (http://quercus.igpp.ucla.edu/OceanInversion/).  (see quercus.igpp.ucla.edu/OceanInversion/ for details).
7    
8    Preprocessed OceanInversion input files are available
9    under directories region_mask, takahashi, and
10    atm_co2 in this package.  All netcdf input files
11    have been converted to binary format in order to avoid
12    having to link netcdf library with MITgcm code.  See
13    respective README files in each directory for details.
14    
15    
16  ===============================================  ===============================================
# Line 22  salinity initial and boundary conditions Line 29  salinity initial and boundary conditions
29   cd MITgcm   cd MITgcm
30   cvs co -d ocean_inversion_project MITgcm_contrib/ocean_inversion_project   cvs co -d ocean_inversion_project MITgcm_contrib/ocean_inversion_project
31    
32  3 ===> compile  3 ===> compile and link
33    
34   cd bin   cd bin
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 ../ocean_inversion_project/code/ptracers_forcing.F.test ptracers_forcing.F   rm ptracers_*.F
  cp ../ocean_inversion_project/code/ptracers_init.F.test ptracers_init.F  
  cp ../ocean_inversion_project/code/ptracers_read_mask.F.test ptracers_read_mask.F  
39   ../tools/genmake   ../tools/genmake
40   make depend   make depend
41   make   make
# Line 38  salinity initial and boundary conditions Line 43  salinity initial and boundary conditions
43  4 ===> execute  4 ===> execute
44    
45   cd ../exe   cd ../exe
46   cp ../verification/global_with_exf/input/* .   cp ../verification/global_with_exf/input/eedata .
47     cp ../verification/global_with_exf/input/data.* .
48     cp ../verification/global_with_exf/input/POLY3.COEFFS .
49     ln -sf ../verification/global_with_exf/input/*.bin .
50   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
51   cp ../ocean_inversion_project/input/data.ptracers.test data.ptracers   cp data.test data
52   mitgcmuv > output .txt   cp data.ptracers.test data.ptracers
53     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
59    
60    
61  ======================================================  ==========================================================
62  Second generate binary input files.  netcdf files are not OK  Instructions for carrying out a 3000-year quasi-stationary
63    integration using the global_ocean.90x40x15 configuration.
64  ======================================================  ==========================================================
65  ======================================================  
66  ======================================================  1 ===> compile, link, and execute
67  ======================================================  
68     cd MITgcm/exe
69   cd MITgcm/bin   rm *
70     cd ../bin
71     rm *
72   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
73   cp ../code/.genmakerc .   cp ../ocean_inversion_project/code/.genmakerc .
74   cp ../code/* .   cp ../ocean_inversion_project/code/* .
  cp ../code/ptracers_forcing.F.test ptracers_forcing.F  
  cp ../code/ptracers_init.F.test ptracers_init.F  
  cp ../code/ptracers_read_mask.F.test ptracers_read_mask.F  
75   ../tools/genmake   ../tools/genmake
76   make depend   make depend
77   make   make
78   cd ../exe   cd ../exe
79   cp ../verification/global_with_exf/input/* .   cp ../verification/global_with_exf/input/eedata .
80   cp ../input/* .   cp ../verification/global_with_exf/input/data.* .
81   cp ../input/data.ptracers.test data.ptracers   cp ../verification/global_with_exf/input/POLY3.COEFFS .
82   mitgcmuv > output .txt   ln -sf ../verification/global_with_exf/input/*.bin .
83     cp ../ocean_inversion_project/input/* .
84  ===> examine interpolation output for 30-region mask   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
85     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
86   matlab   mitgcmuv > output .txt &
87   mask=zeros(94,44,30);  
88   load fort.10  
89   for i=1:length(fort)  ==========================================================
90    mask(fort(i,2)+2,fort(i,3)+2,fort(i,1))=fort(i,4);  MPI instructions for carrying out a 3000-year quasi-stationary
91   end  integration using the global_ocean.90x40x15 configuration.
92   for i=1:30  ==========================================================
93    mypcolor(mask(3:92,3:42,i)');  
94    colorbar,title(i),pause  1 ===> compile, link, and execute
95   end  
96   mypcolor(sum(mask,3)'); colorbar   cd MITgcm/exe
97     rm *
98  =====================================   cd ../bin
99     rm *
   
 Global ocean experiment with ptracers  
 =====================================  
   
 ===> set up and integrate 30-tracer computation  
   
  cd MITgcm/bin  
100   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
101   cp ../code/.genmakerc .   cp ../ocean_inversion_project/code/.genmakerc .
102   cp ../code/* .   cp ../ocean_inversion_project/code/* .
103     cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
104     cp SIZE.h_mpi SIZE.h
105   ../tools/genmake   ../tools/genmake
106   make depend   make depend
107   make   make
108   cd ../exe   cd ../exe
109   cp ../verification/global_with_exf/input/* .   cp ../verification/global_with_exf/input/eedata .
110   cp ../input/* .   cp ../verification/global_with_exf/input/data.* .
111   mitgcmuv > output .txt   cp ../verification/global_with_exf/input/POLY3.COEFFS .
112     ln -sf ../verification/global_with_exf/input/*.bin .
113  ===> look at output   cp ../ocean_inversion_project/input/* .
114     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
115   matlab   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
116   for i=1:30   mpirun -np 36 dplace -s1 -c25-60 mitgcmuv < /dev/null > & ! errlog &
117    tr=readbin(['PTRACER' myint2str(i) ...  
               '.0000000020.001.001.data'],[90 40 15],1);  
   mypcolor(tr(:,:,2)'), colorbar, pause(1)  
  end  
118    
119    ======================================================
120    
121    % some matlab code for looking at fort.10 debug files
122    load fort.10
123    tak=zeros(90,40,12);
124    for n=1:length(fort)
125     m=fort(n,1); i=fort(n,2); j=fort(n,3);
126     if i>0&i<91&j>0&j<41, tak(i,j,m)=fort(n,4); end
127    end
128    lon=2:4:360; lat=-78:4:78;
129    clf, contourf(lon,lat,mean(tak(:,:,1),3)',-10:10)
130    caxis([-6 6]), colorbar, plotland
131    
132    % some matlab code for looking at PTRACER output files
133    salt=readbin('S.0000000020.001.001.data',[90 40 15],1);
134    tracer=zeros(90,40,15,30);
135    for i=1:30
136     fn=['PTRACER' myint2str(i) '.0000000020.001.001.data'];
137     tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
138     clf, mypcolor(tracer(:,:,1,i)'); pause(1)
139    end
140    tmp=sum(tracer,4); clf, mypcolor(tmp(:,:,1)')
141    mypcolor(
142    
143    
144    % some matlab code for checking that one year's worth of
145    % tracer uptake is approximately 1e18 mols.
146    lon=2:4:360; lat=-78:4:78;
147    thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];
148    mask=readbin('hFacC.001.001.data',[90 40 15],1);
149    tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
150    for i=1:30, mydisp(i)
151     fn=['PTRACER' myint2str(i) '.0000001800.001.001.data'];
152     tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
153     for j=1:length(lat)
154      for k=1:length(thk)
155       sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i).*mask(:,j,k)) * ...
156                      thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
157      end
158     end
159    end
160    plot(1:30,0*sumtracer,1:30,sumtracer)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22