/[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.5 by dimitri, Wed Sep 24 06:50:26 2003 UTC revision 1.13 by dimitri, Tue Oct 21 06:21:51 2003 UTC
# Line 22  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    
# Line 33  salinity initial and boundary conditions Line 33  salinity initial and boundary conditions
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   rm ptracers_init.F ptracers_forcing_surf.F PTRACERS.h   rm ptracers_*.F PTRACERS_OPTIONS.h
38   rm ptracers_read_mask.F ptracers_read_takahashi.F   ../tools/genmake2
39   ../tools/genmake     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
40   make depend   make depend
41   make   make
42    
# Line 64  Instructions for carrying out a 3000-yea Line 63  Instructions for carrying out a 3000-yea
63  integration using the global_ocean.90x40x15 configuration.  integration using the global_ocean.90x40x15 configuration.
64  ==========================================================  ==========================================================
65    
66  1 ===> compile, link, and execute   cd MITgcm/exe
67     rm *
68     cd ../bin
69     rm *
70     cp ../verification/global_with_exf/code/* .
71     cp ../ocean_inversion_project/code/* .
72     rm PTRACERS_OPTIONS.h
73     ../tools/genmake2
74       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
75     make depend
76     make
77     cd ../exe
78     cp ../verification/global_with_exf/input/eedata .
79     cp ../verification/global_with_exf/input/data.* .
80     cp ../verification/global_with_exf/input/POLY3.COEFFS .
81     ln -sf ../verification/global_with_exf/input/*.bin .
82     cp ../ocean_inversion_project/input/* .
83     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
84     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
85     mitgcmuv > output .txt &
86    
87    
88    ==========================================================
89    MPI instructions for carrying out a 3000-year quasi-stationary
90    integration using the global_ocean.90x40x15 configuration.
91    ==========================================================
92    
93   cd MITgcm/exe   cd MITgcm/exe
94   rm *   rm *
95   cd ../bin   cd ../bin
96   rm *   rm *
97   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
  cp ../ocean_inversion_project/code/.genmakerc .  
98   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
99   ../tools/genmake   rm PTRACERS_OPTIONS.h
100     cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
101     cp SIZE.h_mpi SIZE.h
102     ../tools/genmake2
103       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
104   make depend   make depend
105   make   make
106   cd ../exe   cd ../exe
# Line 84  integration using the global_ocean.90x40 Line 111  integration using the global_ocean.90x40
111   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
112   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
113   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
114     mpirun -np 36 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog &
115    
116    
117    ==========================================================
118    Instructions for carrying out time-dependent, 1765-2005,
119    anthropogenic carbon perturbation tracer experiments using
120    the global_ocean.90x40x15 configuration.
121    ==========================================================
122    
123     cd MITgcm/exe
124     rm *
125     cd ../bin
126     rm *
127     cp ../verification/global_with_exf/code/* .
128     cp ../ocean_inversion_project/code/* .
129     ../tools/genmake2
130       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
131     make depend
132     make
133     cd ../exe
134     cp ../verification/global_with_exf/input/eedata .
135     cp ../verification/global_with_exf/input/data.* .
136     cp ../verification/global_with_exf/input/POLY3.COEFFS .
137     ln -sf ../verification/global_with_exf/input/*.bin .
138     cp ../ocean_inversion_project/input/* .
139     cp data.1765-2005 data
140     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
141     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
142     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
143     mitgcmuv > output .txt &
144    
145    
146    ==========================================================
147    Instructions for carrying out time-dependent, 1765-2005,
148    anthropogenic carbon perturbation tracer experiments using
149    global_ocean.90x40x15 configuration and netcdf output.
150    This assumes that the appropriate libnetcdf.a exists
151    and that it be specified in the genmake2 option file.
152    ==========================================================
153    
154     cd MITgcm/exe
155     rm *
156     cd ../bin
157     rm *
158     cp ../verification/global_with_exf/code/* .
159     cp ../ocean_inversion_project/code/* .
160     cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h
161     cp ../ocean_inversion_project/write_netCDF/*  .
162    
163    ### need customized option file; example below is for nireas
164     ../tools/genmake2 -of linux_ia32_g77_netcdf
165    
166     make depend
167     make
168     cd ../exe
169     cp ../verification/global_with_exf/input/eedata .
170     cp ../verification/global_with_exf/input/data.* .
171     cp ../verification/global_with_exf/input/POLY3.COEFFS .
172     ln -sf ../verification/global_with_exf/input/*.bin .
173     cp ../ocean_inversion_project/input/* .
174     cp data.1765-2005 data
175     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
176     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
177     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
178   mitgcmuv > output .txt &   mitgcmuv > output .txt &
179    
180    
# Line 116  mypcolor( Line 207  mypcolor(
207  % tracer uptake is approximately 1e18 mols.  % tracer uptake is approximately 1e18 mols.
208  lon=2:4:360; lat=-78:4:78;  lon=2:4:360; lat=-78:4:78;
209  thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];  thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];
210  tracer=zeros(90,40,15,30); sumtracer1=zeros(30,1);  mask=readbin('hFacC.001.001.data',[90 40 15],1);
 for i=1:30, mydisp(i)  
  fn=['PTRACER' myint2str(i) '.0000000180.001.001.data'];  
  tracer(:,:,:,i)=readbin(fn,[90 40 15],1);  
  for x=1:length(lon)  
   for y=1:length(lat)  
    for z=1:length(thk)  
     sumtracer1(i) = sumtracer1(i) + tracer(x,y,z,i) * ...  
                     thk(z) * (4*1.113195e+05)^2 * cos(pi*lat(y)/180);  
    end  
   end  
  end  
 end  
   
 % some matlab code for checking that one year's worth of  
 % tracer uptake is approximately 1e18 mols.  
 lon=2:4:360; lat=-78:4:78;  
 thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];  
211  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
212  for i=1:30, mydisp(i)  for i=1:30, mydisp(i)
213   fn=['PTRACER' myint2str(i) '.0000003600.001.001.data'];   fn=['PTRACER' myint2str(i) '.0000001800.001.001.data'];
214   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
215   for j=1:length(lat)   for j=1:length(lat)
216    for k=1:length(thk)    for k=1:length(thk)
217     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i)) * ...     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i).*mask(:,j,k)) * ...
218                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
219    end    end
220   end   end
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22