Contents of /MITgcm_contrib/ocean_inversion_project/README

Ocean Inversion Project
=======================

Instructions for using pkg/ptracers to compute tracer
Green's functions for Gruber's ocean inversion project
(http://quercus.igpp.ucla.edu/OceanInversion/).

=====================================

To run a pkg/ptracers test case using salt initial and boundary conditions:

 
 cd MITgcm/bin
 cp ../verification/global_with_exf/code/* .
 cp ../code/.genmakerc .
 cp ../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
 ../tools/genmake
 make depend
 make
 cd ../exe
 cp ../verification/global_with_exf/input/* .
 cp ../input/* .
 cp ../input/data.ptracers.test data.ptracers
 mitgcmuv > output .txt
 
===> check that PTRACER output is identical to salinity output

 diff PTRACER30.0000000020.001.001.data PTRACER01.0000000020.001.001.data
 diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data

===> examine interpolation output for 30-region mask

 matlab
 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

=====================================


Global ocean experiment with ptracers
=====================================

===> set up and integrate 30-tracer computation

 cd MITgcm/bin
 cp ../verification/global_with_exf/code/* .
 cp ../code/.genmakerc .
 cp ../code/* .
 ../tools/genmake
 make depend
 make
 cd ../exe
 cp ../verification/global_with_exf/input/* .
 cp ../input/* .
 mitgcmuv > output .txt

===> look at output

 matlab
 for i=1:30
  tr=readbin(['PTRACER' myint2str(i) ...
              '.0000000020.001.001.data'],[90 40 15],1);
  mypcolor(tr(:,:,2)'), colorbar, pause(1)
 end
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	(http://quercus.igpp.ucla.edu/OceanInversion/).
7
8	=====================================
9
10	To run a pkg/ptracers test case using salt initial and boundary conditions:
11
12
13
14	cd MITgcm/bin
15	cp ../verification/global_with_exf/code/* .
16	cp ../code/.genmakerc .
17	cp ../code/* .
18	cp ../code/ptracers_forcing.F.test ptracers_forcing.F
19	cp ../code/ptracers_init.F.test ptracers_init.F
20	cp ../code/ptracers_read_mask.F.test ptracers_read_mask.F
21	../tools/genmake
22	make depend
23	make
24	cd ../exe
25	cp ../verification/global_with_exf/input/* .
26	cp ../input/* .
27	cp ../input/data.ptracers.test data.ptracers
28	mitgcmuv > output .txt
29
30	===> check that PTRACER output is identical to salinity output
31
32	diff PTRACER30.0000000020.001.001.data PTRACER01.0000000020.001.001.data
33	diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data
34
35	===> examine interpolation output for 30-region mask
36
37	matlab
38	mask=zeros(94,44,30);
39	load fort.10
40	for i=1:length(fort)
41	mask(fort(i,2)+2,fort(i,3)+2,fort(i,1))=fort(i,4);
42	end
43	for i=1:30
44	mypcolor(mask(3:92,3:42,i)');
45	colorbar,title(i),pause
46	end
47	mypcolor(sum(mask,3)'); colorbar
48
49	=====================================
50
51
52	Global ocean experiment with ptracers
53	=====================================
54
55	===> set up and integrate 30-tracer computation
56
57	cd MITgcm/bin
58	cp ../verification/global_with_exf/code/* .
59	cp ../code/.genmakerc .
60	cp ../code/* .
61	../tools/genmake
62	make depend
63	make
64	cd ../exe
65	cp ../verification/global_with_exf/input/* .
66	cp ../input/* .
67	mitgcmuv > output .txt
68
69	===> look at output
70
71	matlab
72	for i=1:30
73	tr=readbin(['PTRACER' myint2str(i) ...
74	'.0000000020.001.001.data'],[90 40 15],1);
75	mypcolor(tr(:,:,2)'), colorbar, pause(1)
76	end