/[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.6 by dimitri, Thu Sep 25 03:01:59 2003 UTC revision 1.14 by dimitri, Tue Oct 28 00:27:48 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_*.F   rm ptracers_*.F PTRACERS_OPTIONS.h
38   ../tools/genmake   ../tools/genmake2
39       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
40   make depend   make depend
41   make   make
42    
# Line 56  salinity initial and boundary conditions Line 56  salinity initial and boundary conditions
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     diff PTRtave30.0000000020.001.001.data Stave.0000000020.001.001.data
60     diff PTRtave30.0000000020.001.001.data Stave.0000000020.001.001.data
61    
62    
63  ==========================================================  ==========================================================
64  Instructions for carrying out a 3000-year quasi-stationary  MPI instructions for carrying out a 3000-year quasi-stationary
65  integration using the global_ocean.90x40x15 configuration.  integration using the global_ocean.90x40x15 configuration.
66  ==========================================================  ==========================================================
67    
68  1 ===> compile, link, and execute   cd MITgcm/exe
69     rm *
70     cd ../bin
71     rm *
72     cp ../verification/global_with_exf/code/* .
73     cp ../ocean_inversion_project/code/* .
74     rm PTRACERS_OPTIONS.h
75     cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
76     cp SIZE.h_mpi SIZE.h
77     ../tools/genmake2
78       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
79     make depend
80     make
81     cd ../exe
82     cp ../verification/global_with_exf/input/eedata .
83     cp ../verification/global_with_exf/input/data.* .
84     cp ../verification/global_with_exf/input/POLY3.COEFFS .
85     ln -sf ../verification/global_with_exf/input/*.bin .
86     cp ../ocean_inversion_project/input/* .
87     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
88     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
89     mpirun -np 24 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog &
90    
91    
92    ==========================================================
93    Instructions for carrying out time-dependent, 1765-2005,
94    anthropogenic carbon perturbation tracer experiments using
95    the global_ocean.90x40x15 configuration.
96    ==========================================================
97    
98     cd MITgcm/exe
99     rm *
100     cd ../bin
101     rm *
102     cp ../verification/global_with_exf/code/* .
103     cp ../ocean_inversion_project/code/* .
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     cp data.1765-2005 data
115     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
116     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
117     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
118     mitgcmuv > output .txt &
119    
120    
121    ==========================================================
122    Instructions for carrying out time-dependent, 1765-2005,
123    anthropogenic carbon perturbation tracer experiments using
124    global_ocean.90x40x15 configuration and netcdf output.
125    This assumes that the appropriate libnetcdf.a exists
126    and that it be specified in the genmake2 option file.
127    ==========================================================
128    
129   cd MITgcm/exe   cd MITgcm/exe
130   rm *   rm *
131   cd ../bin   cd ../bin
132   rm *   rm *
133   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
  cp ../ocean_inversion_project/code/.genmakerc .  
134   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
135   ../tools/genmake   cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h
136     cp ../ocean_inversion_project/write_netCDF/*  .
137    
138    ### need customized option file; example below is for nireas
139     ../tools/genmake2 -of linux_ia32_g77_netcdf
140    
141   make depend   make depend
142   make   make
143   cd ../exe   cd ../exe
# Line 81  integration using the global_ocean.90x40 Line 146  integration using the global_ocean.90x40
146   cp ../verification/global_with_exf/input/POLY3.COEFFS .   cp ../verification/global_with_exf/input/POLY3.COEFFS .
147   ln -sf ../verification/global_with_exf/input/*.bin .   ln -sf ../verification/global_with_exf/input/*.bin .
148   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
149     cp data.1765-2005 data
150   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
151   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
152     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
153   mitgcmuv > output .txt &   mitgcmuv > output .txt &
154    
155    
# Line 115  mypcolor( Line 182  mypcolor(
182  % tracer uptake is approximately 1e18 mols.  % tracer uptake is approximately 1e18 mols.
183  lon=2:4:360; lat=-78:4:78;  lon=2:4:360; lat=-78:4:78;
184  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];
185  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];  
186  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
187  for i=1:30, mydisp(i)  for i=1:30, mydisp(i)
188   fn=['PTRACER' myint2str(i) '.0000003600.001.001.data'];   fn=['PTRACER' myint2str(i) '.0000001800.001.001.data'];
189   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
190   for j=1:length(lat)   for j=1:length(lat)
191    for k=1:length(thk)    for k=1:length(thk)
192     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i)) * ...     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i).*mask(:,j,k)) * ...
193                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
194    end    end
195   end   end
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22