20 |
|
|
21 |
1 ===> get MITgcm code from cvs repository |
1 ===> get MITgcm code from cvs repository |
22 |
|
|
23 |
CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/gcmpack |
cvs co -r checkpoint51n_branch MITgcm_code |
24 |
cvs login ( CVS password: cvsanon ) |
cvs co -r checkpoint51n_branch MITgcm/verification/testreport |
25 |
cvs co -P MITgcm |
cvs co -r checkpoint51n_branch MITgcm/verification/global_with_exf |
26 |
|
|
27 |
2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory |
2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory |
28 |
|
|
31 |
|
|
32 |
3 ===> compile and link |
3 ===> compile and link |
33 |
|
|
34 |
|
mkdir bin exe |
35 |
cd bin |
cd bin |
36 |
cp ../verification/global_with_exf/code/* . |
\cp ../verification/global_with_exf/code/* . |
37 |
cp ../ocean_inversion_project/code/* . |
\cp ../ocean_inversion_project/code/* . |
38 |
rm ptracers_*.F PTRACERS_OPTIONS.h |
\rm ptracers_*.F PTRACERS_OPTIONS.h |
39 |
../tools/genmake2 |
../tools/genmake2 |
40 |
(on orion: ../tools/genmake2 -of linux_ia64_efc+mpi) |
(on orion: ../tools/genmake2 -of linux_ia64_efc+mpi ) |
41 |
|
(on columbia: ../tools/genmake2 -of linux_ia64_efc+mpi_altix ) |
42 |
make depend |
make depend |
43 |
make |
make -j 16 |
44 |
|
|
45 |
4 ===> execute |
4 ===> execute |
46 |
|
|
47 |
cd ../exe |
cd ../exe |
48 |
cp ../verification/global_with_exf/input/eedata . |
\cp ../verification/global_with_exf/input/eedata . |
49 |
cp ../verification/global_with_exf/input/data.* . |
\cp ../verification/global_with_exf/input/data.* . |
50 |
cp ../verification/global_with_exf/input/POLY3.COEFFS . |
\cp ../verification/global_with_exf/input/POLY3.COEFFS . |
51 |
ln -sf ../verification/global_with_exf/input/*.bin . |
\ln -sf ../verification/global_with_exf/input/*.bin . |
52 |
cp ../ocean_inversion_project/input/* . |
\cp ../ocean_inversion_project/input/* . |
53 |
cp data.test data |
\cp data.test data |
54 |
cp data.ptracers.test data.ptracers |
\cp data.ptracers.test data.ptracers |
55 |
mitgcmuv > output.txt |
./mitgcmuv >! output.txt |
56 |
|
|
57 |
5 ===> check that PTRACER output and salinity output are identical. |
5 ===> check that PTRACER output and 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 |
88 |
cp ../ocean_inversion_project/input/* . |
cp ../ocean_inversion_project/input/* . |
89 |
ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin . |
ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin . |
90 |
ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin . |
ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin . |
91 |
mpirun -np 24 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog & |
cp data.stationary.1 data |
92 |
|
mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog & |
93 |
|
(on orion: qsub job.orion) |
94 |
|
|
95 |
|
|
96 |
========================================================== |
========================================================== |
97 |
Instructions for carrying out time-dependent, 1765-2005, |
MPI instructions for carrying out time-dependent, 1765-2005, |
98 |
anthropogenic carbon perturbation tracer experiments using |
anthropogenic carbon perturbation tracer experiments using |
99 |
the global_ocean.90x40x15 configuration. |
the global_ocean.90x40x15 configuration. |
100 |
========================================================== |
========================================================== |
105 |
rm * |
rm * |
106 |
cp ../verification/global_with_exf/code/* . |
cp ../verification/global_with_exf/code/* . |
107 |
cp ../ocean_inversion_project/code/* . |
cp ../ocean_inversion_project/code/* . |
108 |
|
cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h |
109 |
|
cp SIZE.h_mpi SIZE.h |
110 |
../tools/genmake2 |
../tools/genmake2 |
111 |
(on orion: ../tools/genmake2 -of linux_ia64_efc+mpi) |
(on orion: ../tools/genmake2 -of linux_ia64_efc+mpi) |
112 |
make depend |
make depend |
121 |
ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin . |
ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin . |
122 |
ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin . |
ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin . |
123 |
ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat . |
ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat . |
124 |
mitgcmuv > output .txt & |
mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog & |
125 |
|
(on orion: qsub job.orion.td) |
126 |
|
|
127 |
|
|
128 |
========================================================== |
========================================================== |
141 |
cp ../ocean_inversion_project/code/* . |
cp ../ocean_inversion_project/code/* . |
142 |
cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h |
cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h |
143 |
cp ../ocean_inversion_project/write_netCDF/* . |
cp ../ocean_inversion_project/write_netCDF/* . |
144 |
|
\rm mk_output.F |
145 |
|
|
146 |
### need customized option file; example below is for nireas |
### need customized option file; example below is for nireas |
147 |
../tools/genmake2 -of linux_ia32_g77_netcdf |
../tools/genmake2 -of linux_ia32_g77_netcdf |
154 |
cp ../verification/global_with_exf/input/POLY3.COEFFS . |
cp ../verification/global_with_exf/input/POLY3.COEFFS . |
155 |
ln -sf ../verification/global_with_exf/input/*.bin . |
ln -sf ../verification/global_with_exf/input/*.bin . |
156 |
cp ../ocean_inversion_project/input/* . |
cp ../ocean_inversion_project/input/* . |
157 |
cp data.1765-2005 data |
cp data.1765-2005.1 data |
158 |
ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin . |
ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin . |
159 |
ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin . |
ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin . |
160 |
ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat . |
ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat . |
161 |
mitgcmuv > output .txt & |
mitgcmuv > output .txt & |
162 |
|
|
163 |
|
|
164 |
|
========================================================== |
165 |
|
Instructions for carrying out time-dependent, 1765-2005, |
166 |
|
anthropogenic carbon perturbation tracer experiments using |
167 |
|
the ecco1x1 configuration and netcdf output on columbia. |
168 |
|
This assumes that the appropriate libnetcdf.a exists |
169 |
|
and that it be specified in the genmake2 option file. |
170 |
|
It also assumes that surface forcing files are available |
171 |
|
in /nobackup2/menemenl/ocmip/iter69 |
172 |
|
========================================================== |
173 |
|
|
174 |
|
cd MITgcm/exe |
175 |
|
\rm * |
176 |
|
cd ../bin |
177 |
|
\rm * |
178 |
|
\cp ../ocean_inversion_project/code_ecco1x1/* . |
179 |
|
\cp ../ocean_inversion_project/write_netCDF/* . |
180 |
|
\rm mk_output.F |
181 |
|
\cp SIZE.h.96 SIZE.h |
182 |
|
../tools/genmake2 -of linux_ia64_efc+mpi_altix |
183 |
|
make depend |
184 |
|
make -j 16 |
185 |
|
cd ../exe |
186 |
|
\cp ../ocean_inversion_project/input_ecco1x1/* . |
187 |
|
ln -sf ../../iter69/* . |
188 |
|
ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin . |
189 |
|
ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin . |
190 |
|
ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat . |
191 |
|
qsub job01 |
192 |
|
|
193 |
|
|
194 |
|
========================================================== |
195 |
|
Instructions for generating netcdf output files on columbia |
196 |
|
========================================================== |
197 |
|
|
198 |
|
cd ocean_inversion_project/write_netCDF |
199 |
|
setenv F_UFMTENDIAN big |
200 |
|
|
201 |
|
(check that defaults in mk_output.F are appropriate) |
202 |
|
|
203 |
|
efc -W0 -WB mk_output.F write_nc_phys.F nc_util.F \ |
204 |
|
handle_errors.F write_nc_basisfnctns.F \ |
205 |
|
write_nc_diag_0D.F write_nc_diag_2D.F \ |
206 |
|
-I/u/menemenl/software/netcdf-3.5.0/include \ |
207 |
|
-L/u/menemenl/software/netcdf-3.5.0/lib -lnetcdf |
208 |
|
a.out |
209 |
|
|
210 |
|
|
211 |
|
========================================================== |
212 |
|
Instructions for generating netcdf output files on orion |
213 |
|
========================================================== |
214 |
|
|
215 |
|
cd ocean_inversion_project/write_netCDF |
216 |
|
setenv F_UFMTENDIAN big |
217 |
|
|
218 |
|
(check that defaults in mk_output.F are appropriate) |
219 |
|
|
220 |
|
efc -W0 -WB mk_output.F write_nc_phys.F nc_util.F \ |
221 |
|
handle_errors.F write_nc_basisfnctns.F \ |
222 |
|
write_nc_diag_0D.F write_nc_diag_2D.F \ |
223 |
|
-I/u2/dmenem/software/netcdf-3.5.0/include \ |
224 |
|
-L/u2/dmenem/software/netcdf-3.5.0/lib -lnetcdf |
225 |
|
a.out |
226 |
|
|
227 |
|
|
228 |
====================================================== |
====================================================== |
229 |
|
|
230 |
% some matlab code for looking at fort.10 debug files |
% some matlab code for looking at fort.10 debug files |
267 |
end |
end |
268 |
end |
end |
269 |
plot(1:30,0*sumtracer,1:30,sumtracer) |
plot(1:30,0*sumtracer,1:30,sumtracer) |
270 |
|
|
271 |
|
|
272 |
|
% check that tracer flux is 1e18 mols / year |
273 |
|
nb_seconds_per_year=31556880; |
274 |
|
lon=2:4:360; lat=-78:4:78; |
275 |
|
mask=readbin('hFacC.data',[90 40 15],1); |
276 |
|
tracer=zeros(90,40,30); sumtracer=zeros(30,1); |
277 |
|
for i=1:30, mydisp(i) |
278 |
|
fn=['PtrFlux' myint2str(i) '.0000000180.data']; |
279 |
|
tracer(:,:,i)=readbin(fn,[90 40],1); |
280 |
|
for j=1:length(lat) |
281 |
|
sumtracer(i) = sumtracer(i) + sum(tracer(:,j,i).*mask(:,j,1)) * ... |
282 |
|
nb_seconds_per_year* (4*1.113195e+05)^2 * cos(pi*lat(j)/180); |
283 |
|
end |
284 |
|
end |
285 |
|
plot(1:30,0*sumtracer,1:30,sumtracer-mmean(sumtracer)) |