/[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.12 by dimitri, Tue Oct 21 03:25:30 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 15  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    
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/* .
  cp ../ocean_inversion_project/code/.genmakerc .  
36   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
37   cp ../ocean_inversion_project/code/ptracers_forcing.F.test ptracers_forcing.F   rm ptracers_*.F PTRACERS_OPTIONS.h
38   cp ../ocean_inversion_project/code/ptracers_init.F.test ptracers_init.F   ../tools/genmake2
39   cp ../ocean_inversion_project/code/ptracers_read_mask.F.test ptracers_read_mask.F     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
  ../tools/genmake  
40   make depend   make depend
41   make   make
42    
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  ======================================================   cd MITgcm/exe
67  ======================================================   rm *
68     cd ../bin
69   cd MITgcm/bin   rm *
70   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
71   cp ../code/.genmakerc .   cp ../ocean_inversion_project/code/.genmakerc .
72   cp ../code/* .   cp ../ocean_inversion_project/code/* .
73   cp ../code/ptracers_forcing.F.test ptracers_forcing.F   rm PTRACERS_OPTIONS.h
74   cp ../code/ptracers_init.F.test ptracers_init.F   ../tools/genmake2
75   cp ../code/ptracers_read_mask.F.test ptracers_read_mask.F     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
  ../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 &
  mask=zeros(94,44,30);  
  load fort.10  
  for i=1:length(fort)  
   mask(fort(i,2)+2,fort(i,3)+2,fort(i,1))=fort(i,4);  
  end  
  for i=1:30  
   mypcolor(mask(3:92,3:42,i)');  
   colorbar,title(i),pause  
  end  
  mypcolor(sum(mask,3)'); colorbar  
   
 =====================================  
87    
88    
89  Global ocean experiment with ptracers  ==========================================================
90  =====================================  MPI instructions for carrying out a 3000-year quasi-stationary
91    integration using the global_ocean.90x40x15 configuration.
92    ==========================================================
93    
94     cd MITgcm/exe
95     rm *
96     cd ../bin
97     rm *
98     cp ../verification/global_with_exf/code/* .
99     cp ../ocean_inversion_project/code/.genmakerc .
100     cp ../ocean_inversion_project/code/* .
101     rm PTRACERS_OPTIONS.h
102     cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
103     cp SIZE.h_mpi SIZE.h
104     ../tools/genmake2
105       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
106     make depend
107     make
108     cd ../exe
109     cp ../verification/global_with_exf/input/eedata .
110     cp ../verification/global_with_exf/input/data.* .
111     cp ../verification/global_with_exf/input/POLY3.COEFFS .
112     ln -sf ../verification/global_with_exf/input/*.bin .
113     cp ../ocean_inversion_project/input/* .
114     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
115     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
116     mpirun -np 36 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog &
117    
 ===> set up and integrate 30-tracer computation  
118    
119   cd MITgcm/bin  ==========================================================
120    Instructions for carrying out time-dependent, 1765-2005,
121    anthropogenic carbon perturbation tracer experiments using
122    the global_ocean.90x40x15 configuration.
123    ==========================================================
124    
125     cd MITgcm/exe
126     rm *
127     cd ../bin
128     rm *
129   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
130   cp ../code/.genmakerc .   cp ../ocean_inversion_project/code/.genmakerc .
131   cp ../code/* .   cp ../ocean_inversion_project/code/* .
132   ../tools/genmake   ../tools/genmake2
133       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
134   make depend   make depend
135   make   make
136   cd ../exe   cd ../exe
137   cp ../verification/global_with_exf/input/* .   cp ../verification/global_with_exf/input/eedata .
138   cp ../input/* .   cp ../verification/global_with_exf/input/data.* .
139   mitgcmuv > output .txt   cp ../verification/global_with_exf/input/POLY3.COEFFS .
140     ln -sf ../verification/global_with_exf/input/*.bin .
141  ===> look at output   cp ../ocean_inversion_project/input/* .
142     cp data.1765-2005 data
143   matlab   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
144   for i=1:30   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
145    tr=readbin(['PTRACER' myint2str(i) ...   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
146                '.0000000020.001.001.data'],[90 40 15],1);   mitgcmuv > output .txt &
147    mypcolor(tr(:,:,2)'), colorbar, pause(1)  
  end  
148    
149    ======================================================
150    
151    % some matlab code for looking at fort.10 debug files
152    load fort.10
153    tak=zeros(90,40,12);
154    for n=1:length(fort)
155     m=fort(n,1); i=fort(n,2); j=fort(n,3);
156     if i>0&i<91&j>0&j<41, tak(i,j,m)=fort(n,4); end
157    end
158    lon=2:4:360; lat=-78:4:78;
159    clf, contourf(lon,lat,mean(tak(:,:,1),3)',-10:10)
160    caxis([-6 6]), colorbar, plotland
161    
162    % some matlab code for looking at PTRACER output files
163    salt=readbin('S.0000000020.001.001.data',[90 40 15],1);
164    tracer=zeros(90,40,15,30);
165    for i=1:30
166     fn=['PTRACER' myint2str(i) '.0000000020.001.001.data'];
167     tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
168     clf, mypcolor(tracer(:,:,1,i)'); pause(1)
169    end
170    tmp=sum(tracer,4); clf, mypcolor(tmp(:,:,1)')
171    mypcolor(
172    
173    
174    % some matlab code for checking that one year's worth of
175    % tracer uptake is approximately 1e18 mols.
176    lon=2:4:360; lat=-78:4:78;
177    thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];
178    mask=readbin('hFacC.001.001.data',[90 40 15],1);
179    tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
180    for i=1:30, mydisp(i)
181     fn=['PTRACER' myint2str(i) '.0000001800.001.001.data'];
182     tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
183     for j=1:length(lat)
184      for k=1:length(thk)
185       sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i).*mask(:,j,k)) * ...
186                      thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
187      end
188     end
189    end
190    plot(1:30,0*sumtracer,1:30,sumtracer)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22