/[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.16 by dimitri, Tue Oct 28 00:37:28 2003 UTC revision 1.25 by dimitri, Mon Nov 1 20:03:35 2004 UTC
# Line 20  salinity initial and boundary conditions Line 20  salinity initial and boundary conditions
20    
21  1 ===> get MITgcm code from cvs repository  1 ===> get MITgcm code from cvs repository
22    
23   CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/gcmpack   cvs co -r checkpoint51n_branch MITgcm_code
24   cvs login ( CVS password: cvsanon )   cvs co -r checkpoint51n_branch MITgcm/verification/testreport
25   cvs co -r checkpoint51n_branch MITgcm   cvs co -r checkpoint51n_branch MITgcm/verification/global_with_exf
26    
27  2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory  2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory
28    
# Line 31  salinity initial and boundary conditions Line 31  salinity initial and boundary conditions
31    
32  3 ===> compile and link  3 ===> compile and link
33    
34     mkdir bin exe
35   cd bin   cd bin
36   cp ../verification/global_with_exf/code/* .   \cp ../verification/global_with_exf/code/* .
37   cp ../ocean_inversion_project/code/* .   \cp ../ocean_inversion_project/code/* .
38   rm ptracers_*.F PTRACERS_OPTIONS.h   \rm ptracers_*.F PTRACERS_OPTIONS.h
39   ../tools/genmake2   ../tools/genmake2
40     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)     (on orion:    ../tools/genmake2 -of linux_ia64_efc+mpi       )
41       (on columbia: ../tools/genmake2 -of linux_ia64_efc+mpi_altix )
42   make depend   make depend
43   make   make -j 16
44    
45  4 ===> execute  4 ===> execute
46    
47   cd ../exe   cd ../exe
48   cp ../verification/global_with_exf/input/eedata .   \cp ../verification/global_with_exf/input/eedata .
49   cp ../verification/global_with_exf/input/data.* .   \cp ../verification/global_with_exf/input/data.* .
50   cp ../verification/global_with_exf/input/POLY3.COEFFS .   \cp ../verification/global_with_exf/input/POLY3.COEFFS .
51   ln -sf ../verification/global_with_exf/input/*.bin .   \ln -sf ../verification/global_with_exf/input/*.bin .
52   cp ../ocean_inversion_project/input/* .   \cp ../ocean_inversion_project/input/* .
53   cp data.test data   \cp data.test data
54   cp data.ptracers.test data.ptracers   \cp data.ptracers.test data.ptracers
55   mitgcmuv > output.txt   ./mitgcmuv >! output.txt
56      
57  5 ===> check that PTRACER output and salinity output are identical.  5 ===> check that PTRACER output and salinity output are identical.
58    
59   diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data   diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data
# Line 86  integration using the global_ocean.90x40 Line 88  integration using the global_ocean.90x40
88   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
89   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
90   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
91   mpirun -np 24 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog &   cp data.stationary.1 data
92     mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
93       (on orion: qsub job.orion)
94    
95    
96  ==========================================================  ==========================================================
97  Instructions for carrying out time-dependent, 1765-2005,  MPI instructions for carrying out time-dependent, 1765-2005,
98  anthropogenic carbon perturbation tracer experiments using  anthropogenic carbon perturbation tracer experiments using
99  the global_ocean.90x40x15 configuration.  the global_ocean.90x40x15 configuration.
100  ==========================================================  ==========================================================
# Line 101  the global_ocean.90x40x15 configuration. Line 105  the global_ocean.90x40x15 configuration.
105   rm *   rm *
106   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
107   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
108     cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
109     cp SIZE.h_mpi SIZE.h
110   ../tools/genmake2   ../tools/genmake2
111     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
112   make depend   make depend
# Line 115  the global_ocean.90x40x15 configuration. Line 121  the global_ocean.90x40x15 configuration.
121   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
122   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
123   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
124   mitgcmuv > output .txt &   mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
125       (on orion: qsub job.orion.td)
126    
127    
128  ==========================================================  ==========================================================
# Line 134  and that it be specified in the genmake2 Line 141  and that it be specified in the genmake2
141   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
142   cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h   cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h
143   cp ../ocean_inversion_project/write_netCDF/*  .   cp ../ocean_inversion_project/write_netCDF/*  .
144     \rm mk_output.F
145    
146  ### need customized option file; example below is for nireas  ### need customized option file; example below is for nireas
147   ../tools/genmake2 -of linux_ia32_g77_netcdf   ../tools/genmake2 -of linux_ia32_g77_netcdf
# Line 146  and that it be specified in the genmake2 Line 154  and that it be specified in the genmake2
154   cp ../verification/global_with_exf/input/POLY3.COEFFS .   cp ../verification/global_with_exf/input/POLY3.COEFFS .
155   ln -sf ../verification/global_with_exf/input/*.bin .   ln -sf ../verification/global_with_exf/input/*.bin .
156   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
157   cp data.1765-2005 data   cp data.1765-2005.1 data
158   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
159   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
160   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
161   mitgcmuv > output .txt &   mitgcmuv > output .txt &
162    
163    
164    ==========================================================
165    Instructions for carrying out time-dependent, 1765-2005,
166    anthropogenic carbon perturbation tracer experiments using
167    the ecco1x1 configuration and netcdf output on columbia.
168    This assumes that the appropriate libnetcdf.a exists
169    and that it be specified in the genmake2 option file.
170    It also assumes that surface forcing files are available
171    in /nobackup2/menemenl/ocmip/iter69
172    ==========================================================
173    
174     cd MITgcm/exe
175     \rm *
176     cd ../bin
177     \rm *
178     \cp ../ocean_inversion_project/code_ecco1x1/* .
179     \cp ../ocean_inversion_project/write_netCDF/*  .
180     \rm mk_output.F
181     \cp SIZE.h.96 SIZE.h
182     ../tools/genmake2 -of linux_ia64_efc+mpi_altix
183     make depend
184     make -j 16
185     cd ../exe
186     \cp ../ocean_inversion_project/input_ecco1x1/* .
187     ln -sf ../../iter69/* .
188     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
189     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
190     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
191     qsub job01
192    
193    
194    ==========================================================
195    Instructions for generating netcdf output files on columbia
196    ==========================================================
197    
198    cd ocean_inversion_project/write_netCDF
199    setenv F_UFMTENDIAN big
200    
201    (check that defaults in  mk_output.F are appropriate)
202    
203    efc -W0 -WB mk_output.F write_nc_phys.F nc_util.F \
204         handle_errors.F write_nc_basisfnctns.F \
205         write_nc_diag_0D.F write_nc_diag_2D.F \
206         -I/u/menemenl/software/netcdf-3.5.0/include \
207         -L/u/menemenl/software/netcdf-3.5.0/lib -lnetcdf
208    a.out
209    
210    
211    ==========================================================
212    Instructions for generating netcdf output files on orion
213    ==========================================================
214    
215    cd ocean_inversion_project/write_netCDF
216    setenv F_UFMTENDIAN big
217    
218    (check that defaults in  mk_output.F are appropriate)
219    
220    efc -W0 -WB mk_output.F write_nc_phys.F nc_util.F \
221         handle_errors.F write_nc_basisfnctns.F \
222         write_nc_diag_0D.F write_nc_diag_2D.F \
223         -I/u2/dmenem/software/netcdf-3.5.0/include \
224         -L/u2/dmenem/software/netcdf-3.5.0/lib -lnetcdf
225    a.out
226    
227    
228  ======================================================  ======================================================
229    
230  % some matlab code for looking at fort.10 debug files  % some matlab code for looking at fort.10 debug files
# Line 195  for i=1:30, mydisp(i) Line 267  for i=1:30, mydisp(i)
267   end   end
268  end  end
269  plot(1:30,0*sumtracer,1:30,sumtracer)  plot(1:30,0*sumtracer,1:30,sumtracer)
270    
271    
272    % check that tracer flux is 1e18 mols / year
273    nb_seconds_per_year=31556880;
274    lon=2:4:360; lat=-78:4:78;
275    mask=readbin('hFacC.data',[90 40 15],1);
276    tracer=zeros(90,40,30); sumtracer=zeros(30,1);
277    for i=1:30, mydisp(i)
278     fn=['PtrFlux' myint2str(i) '.0000000180.data'];
279     tracer(:,:,i)=readbin(fn,[90 40],1);
280     for j=1:length(lat)
281        sumtracer(i) = sumtracer(i) + sum(tracer(:,j,i).*mask(:,j,1)) * ...
282            nb_seconds_per_year* (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
283      end
284    end
285    plot(1:30,0*sumtracer,1:30,sumtracer-mmean(sumtracer))
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22