/[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.17 by dimitri, Tue Oct 28 01:31:02 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 -r checkpoint51n_branch 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 57  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     cp data.stationary.1 data
90     mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
91    
92    
93    ==========================================================
94    Instructions for carrying out time-dependent, 1765-2005,
95    anthropogenic carbon perturbation tracer experiments using
96    the global_ocean.90x40x15 configuration.
97    ==========================================================
98    
99     cd MITgcm/exe
100     rm *
101     cd ../bin
102     rm *
103     cp ../verification/global_with_exf/code/* .
104     cp ../ocean_inversion_project/code/* .
105     ../tools/genmake2
106       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
107     make depend
108     make
109     cd ../exe
110     cp ../verification/global_with_exf/input/eedata .
111     cp ../verification/global_with_exf/input/data.* .
112     cp ../verification/global_with_exf/input/POLY3.COEFFS .
113     ln -sf ../verification/global_with_exf/input/*.bin .
114     cp ../ocean_inversion_project/input/* .
115     cp data.1765-2005 data
116     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
117     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
118     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
119     mitgcmuv > output .txt &
120    
121    
122    ==========================================================
123    Instructions for carrying out time-dependent, 1765-2005,
124    anthropogenic carbon perturbation tracer experiments using
125    global_ocean.90x40x15 configuration and netcdf output.
126    This assumes that the appropriate libnetcdf.a exists
127    and that it be specified in the genmake2 option file.
128    ==========================================================
129    
130   cd MITgcm/exe   cd MITgcm/exe
131   rm *   rm *
132   cd ../bin   cd ../bin
133   rm *   rm *
134   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
  cp ../ocean_inversion_project/code/.genmakerc .  
135   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
136   ../tools/genmake   cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h
137     cp ../ocean_inversion_project/write_netCDF/*  .
138    
139    ### need customized option file; example below is for nireas
140     ../tools/genmake2 -of linux_ia32_g77_netcdf
141    
142   make depend   make depend
143   make   make
144   cd ../exe   cd ../exe
# Line 82  integration using the global_ocean.90x40 Line 147  integration using the global_ocean.90x40
147   cp ../verification/global_with_exf/input/POLY3.COEFFS .   cp ../verification/global_with_exf/input/POLY3.COEFFS .
148   ln -sf ../verification/global_with_exf/input/*.bin .   ln -sf ../verification/global_with_exf/input/*.bin .
149   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
150     cp data.1765-2005 data
151   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
152   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
153     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
154   mitgcmuv > output .txt &   mitgcmuv > output .txt &
155    
156    
# Line 116  mypcolor( Line 183  mypcolor(
183  % tracer uptake is approximately 1e18 mols.  % tracer uptake is approximately 1e18 mols.
184  lon=2:4:360; lat=-78:4:78;  lon=2:4:360; lat=-78:4:78;
185  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];
186  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];  
187  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
188  for i=1:30, mydisp(i)  for i=1:30, mydisp(i)
189   fn=['PTRACER' myint2str(i) '.0000003600.001.001.data'];   fn=['PTRACER' myint2str(i) '.0000001800.001.001.data'];
190   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
191   for j=1:length(lat)   for j=1:length(lat)
192    for k=1:length(thk)    for k=1:length(thk)
193     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i)) * ...     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i).*mask(:,j,k)) * ...
194                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
195    end    end
196   end   end
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22