/[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.20 by dimitri, Thu Dec 18 06:11:14 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_*.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     cp data.stationary.1 data
90     mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
91       (on orion: qsub job.orion)
92    
93    ==========================================================
94    MPI 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   cd MITgcm/exe
100   rm *   rm *
101   cd ../bin   cd ../bin
102   rm *   rm *
103   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
  cp ../ocean_inversion_project/code/.genmakerc .  
104   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
105   ../tools/genmake   cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
106     cp SIZE.h_mpi SIZE.h
107     ../tools/genmake2
108       (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
109   make depend   make depend
110   make   make
111   cd ../exe   cd ../exe
# Line 81  integration using the global_ocean.90x40 Line 114  integration using the global_ocean.90x40
114   cp ../verification/global_with_exf/input/POLY3.COEFFS .   cp ../verification/global_with_exf/input/POLY3.COEFFS .
115   ln -sf ../verification/global_with_exf/input/*.bin .   ln -sf ../verification/global_with_exf/input/*.bin .
116   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
117     cp data.1765-2005 data
118   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
119   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
120     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
121     mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
122       (on orion: qsub job.orion.td)
123    
124    
125    ==========================================================
126    Instructions for carrying out time-dependent, 1765-2005,
127    anthropogenic carbon perturbation tracer experiments using
128    global_ocean.90x40x15 configuration and netcdf output.
129    This assumes that the appropriate libnetcdf.a exists
130    and that it be specified in the genmake2 option file.
131    ==========================================================
132    
133     cd MITgcm/exe
134     rm *
135     cd ../bin
136     rm *
137     cp ../verification/global_with_exf/code/* .
138     cp ../ocean_inversion_project/code/* .
139     cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h
140     cp ../ocean_inversion_project/write_netCDF/*  .
141     \rm mk_output.F
142    
143    ### need customized option file; example below is for nireas
144     ../tools/genmake2 -of linux_ia32_g77_netcdf
145    
146     make depend
147     make
148     cd ../exe
149     cp ../verification/global_with_exf/input/eedata .
150     cp ../verification/global_with_exf/input/data.* .
151     cp ../verification/global_with_exf/input/POLY3.COEFFS .
152     ln -sf ../verification/global_with_exf/input/*.bin .
153     cp ../ocean_inversion_project/input/* .
154     cp data.1765-2005.1 data
155     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
156     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
157     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
158   mitgcmuv > output .txt &   mitgcmuv > output .txt &
159    
160    
161  ======================================================  ======================================================
162    
163    ==========================================================
164    Instructions for generating netcdf output files on orion
165    ==========================================================
166    
167    cd ocean_inversion_project/write_netCDF
168    setenv F_UFMTENDIAN big
169    
170    (check that defaults in  mk_output.F are appropriate)
171    
172    efc -W0 -WB mk_output.F write_nc_phys.F nc_util.F \
173         handle_errors.F write_nc_basisfnctns.F \
174         write_nc_diag_0D.F write_nc_diag_2D.F \
175         -I/u2/dmenem/software/netcdf-3.5.0/include \
176         -L/u2/dmenem/software/netcdf-3.5.0/lib -lnetcdf
177    a.out
178    
179    
180    ======================================================
181    
182  % some matlab code for looking at fort.10 debug files  % some matlab code for looking at fort.10 debug files
183  load fort.10  load fort.10
184  tak=zeros(90,40,12);  tak=zeros(90,40,12);
# Line 115  mypcolor( Line 206  mypcolor(
206  % tracer uptake is approximately 1e18 mols.  % tracer uptake is approximately 1e18 mols.
207  lon=2:4:360; lat=-78:4:78;  lon=2:4:360; lat=-78:4:78;
208  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];
209  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];  
210  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
211  for i=1:30, mydisp(i)  for i=1:30, mydisp(i)
212   fn=['PTRACER' myint2str(i) '.0000003600.001.001.data'];   fn=['PTRACER' myint2str(i) '.0000001800.001.001.data'];
213   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
214   for j=1:length(lat)   for j=1:length(lat)
215    for k=1:length(thk)    for k=1:length(thk)
216     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i)) * ...     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i).*mask(:,j,k)) * ...
217                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
218    end    end
219   end   end
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22