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

Contents of /MITgcm_contrib/ocean_inversion_project/README

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

Revision 1.10 - (show annotations) (download)
Fri Oct 10 16:33:09 2003 UTC (22 years ago) by dimitri
Branch: MAIN
Changes since 1.9: +1 -1 lines 
Modified ocean_inversion_project/README
1 Ocean Inversion Project
2 =======================
3
4 Instructions for using pkg/ptracers to compute tracer
5 Green's functions for Gruber's ocean inversion project
6 (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 ===============================================
17 First check that pkg/ptracers works OK by using
18 salinity initial and boundary conditions
19 ===============================================
20
21 1 ===> get MITgcm code from cvs repository
22
23 CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/gcmpack
24 cvs login ( CVS password: cvsanon )
25 cvs co -P MITgcm
26
27 2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory
28
29 cd MITgcm
30 cvs co -d ocean_inversion_project MITgcm_contrib/ocean_inversion_project
31
32 3 ===> compile and link
33
34 cd bin
35 cp ../verification/global_with_exf/code/* .
36 cp ../ocean_inversion_project/code/.genmakerc .
37 cp ../ocean_inversion_project/code/* .
38 rm ptracers_*.F
39 ../tools/genmake
40 make depend
41 make
42
43 4 ===> execute
44
45 cd ../exe
46 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/* .
51 cp data.test data
52 cp data.ptracers.test data.ptracers
53 mitgcmuv > output.txt
54
55 5 ===> check that PTRACER output and salinity output are identical.
56
57 diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data
58 diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data
59
60
61 ==========================================================
62 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 rm *
70 cd ../bin
71 rm *
72 cp ../verification/global_with_exf/code/* .
73 cp ../ocean_inversion_project/code/.genmakerc .
74 cp ../ocean_inversion_project/code/* .
75 ../tools/genmake
76 make depend
77 make
78 cd ../exe
79 cp ../verification/global_with_exf/input/eedata .
80 cp ../verification/global_with_exf/input/data.* .
81 cp ../verification/global_with_exf/input/POLY3.COEFFS .
82 ln -sf ../verification/global_with_exf/input/*.bin .
83 cp ../ocean_inversion_project/input/* .
84 ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
85 ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
86 mitgcmuv > output .txt &
87
88
89 ==========================================================
90 MPI instructions for carrying out a 3000-year quasi-stationary
91 integration using the global_ocean.90x40x15 configuration.
92 ==========================================================
93
94 1 ===> compile, link, and execute
95
96 cd MITgcm/exe
97 rm *
98 cd ../bin
99 rm *
100 cp ../verification/global_with_exf/code/* .
101 cp ../ocean_inversion_project/code/.genmakerc .
102 cp ../ocean_inversion_project/code/* .
103 cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
104 cp SIZE.h_mpi SIZE.h
105 ../tools/genmake
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
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)
  ViewVC Help
Powered by ViewVC 1.1.22