/[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.13 by dimitri, Tue Oct 21 06:21:51 2003 UTC revision 1.23 by dimitri, Fri Oct 8 06:21:29 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 -P 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
60   diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data   diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data
61     diff PTRtave30.0000000020.001.001.data Stave.0000000020.001.001.data
62     diff PTRtave30.0000000020.001.001.data Stave.0000000020.001.001.data
63    
64    
65  ==========================================================  ==========================================================
66  Instructions for carrying out a 3000-year quasi-stationary  MPI instructions for carrying out a 3000-year quasi-stationary
67  integration using the global_ocean.90x40x15 configuration.  integration using the global_ocean.90x40x15 configuration.
68  ==========================================================  ==========================================================
69    
# Line 70  integration using the global_ocean.90x40 Line 74  integration using the global_ocean.90x40
74   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
75   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
76   rm PTRACERS_OPTIONS.h   rm PTRACERS_OPTIONS.h
77     cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
78     cp SIZE.h_mpi SIZE.h
79   ../tools/genmake2   ../tools/genmake2
80     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
81   make depend   make depend
# Line 82  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   mitgcmuv > output .txt &   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  MPI instructions for carrying out a 3000-year quasi-stationary  MPI instructions for carrying out time-dependent, 1765-2005,
98  integration using the global_ocean.90x40x15 configuration.  anthropogenic carbon perturbation tracer experiments using
99    the global_ocean.90x40x15 configuration.
100  ==========================================================  ==========================================================
101    
102   cd MITgcm/exe   cd MITgcm/exe
# Line 96  integration using the global_ocean.90x40 Line 105  integration using the global_ocean.90x40
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/* .
  rm PTRACERS_OPTIONS.h  
108   cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h   cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
109   cp SIZE.h_mpi SIZE.h   cp SIZE.h_mpi SIZE.h
110   ../tools/genmake2   ../tools/genmake2
# Line 109  integration using the global_ocean.90x40 Line 117  integration using the global_ocean.90x40
117   cp ../verification/global_with_exf/input/POLY3.COEFFS .   cp ../verification/global_with_exf/input/POLY3.COEFFS .
118   ln -sf ../verification/global_with_exf/input/*.bin .   ln -sf ../verification/global_with_exf/input/*.bin .
119   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
120     cp data.1765-2005 data
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   mpirun -np 36 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog &   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
124     mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
125       (on orion: qsub job.orion.td)
126    
127    
128  ==========================================================  ==========================================================
129  Instructions for carrying out time-dependent, 1765-2005,  Instructions for carrying out time-dependent, 1765-2005,
130  anthropogenic carbon perturbation tracer experiments using  anthropogenic carbon perturbation tracer experiments using
131  the global_ocean.90x40x15 configuration.  global_ocean.90x40x15 configuration and netcdf output.
132    This assumes that the appropriate libnetcdf.a exists
133    and that it be specified in the genmake2 option file.
134  ==========================================================  ==========================================================
135    
136   cd MITgcm/exe   cd MITgcm/exe
# Line 126  the global_ocean.90x40x15 configuration. Line 139  the global_ocean.90x40x15 configuration.
139   rm *   rm *
140   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
141   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
142   ../tools/genmake2   cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h
143     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)   cp ../ocean_inversion_project/write_netCDF/*  .
144     \rm mk_output.F
145    
146    ### need customized option file; example below is for nireas
147     ../tools/genmake2 -of linux_ia32_g77_netcdf
148    
149   make depend   make depend
150   make   make
151   cd ../exe   cd ../exe
# Line 136  the global_ocean.90x40x15 configuration. Line 154  the global_ocean.90x40x15 configuration.
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 .
# Line 146  the global_ocean.90x40x15 configuration. Line 164  the global_ocean.90x40x15 configuration.
164  ==========================================================  ==========================================================
165  Instructions for carrying out time-dependent, 1765-2005,  Instructions for carrying out time-dependent, 1765-2005,
166  anthropogenic carbon perturbation tracer experiments using  anthropogenic carbon perturbation tracer experiments using
167  global_ocean.90x40x15 configuration and netcdf output.  the ecco1x1 configuration and netcdf output on columbia.
168  This assumes that the appropriate libnetcdf.a exists  This assumes that the appropriate libnetcdf.a exists
169  and that it be specified in the genmake2 option file.  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   cd MITgcm/exe
175   rm *   \rm *
176   cd ../bin   cd ../bin
177   rm *   \rm *
178   cp ../verification/global_with_exf/code/* .   \cp ../ocean_inversion_project/code_ecco1x1/* .
179   cp ../ocean_inversion_project/code/* .   \cp ../ocean_inversion_project/write_netCDF/*  .
180   cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h   \rm mk_output.F
181   cp ../ocean_inversion_project/write_netCDF/*  .   \cp SIZE.h.144 SIZE.h
182     ../tools/genmake2 -of linux_ia64_efc+mpi_altix
 ### need customized option file; example below is for nireas  
  ../tools/genmake2 -of linux_ia32_g77_netcdf  
   
183   make depend   make depend
184   make   make -j 16
185   cd ../exe   cd ../exe
186   cp ../verification/global_with_exf/input/eedata .   \cp ../ocean_inversion_project/input_ecco1x1/* .
187   cp ../verification/global_with_exf/input/data.* .   ln -sf ../../iter69/* .
  cp ../verification/global_with_exf/input/POLY3.COEFFS .  
  ln -sf ../verification/global_with_exf/input/*.bin .  
  cp ../ocean_inversion_project/input/* .  
  cp data.1765-2005 data  
188   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
189   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
190   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
191   mitgcmuv > output .txt &   qsub job01
192    
193    
194    ==========================================================
195    Instructions for generating netcdf output files on orion
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/u2/dmenem/software/netcdf-3.5.0/include \
207         -L/u2/dmenem/software/netcdf-3.5.0/lib -lnetcdf
208    a.out
209    
210    
211  ======================================================  ======================================================
# Line 220  for i=1:30, mydisp(i) Line 250  for i=1:30, mydisp(i)
250   end   end
251  end  end
252  plot(1:30,0*sumtracer,1:30,sumtracer)  plot(1:30,0*sumtracer,1:30,sumtracer)
253    
254    
255    % check that tracer flux is 1e18 mols / year
256    nb_seconds_per_year=31556880;
257    lon=2:4:360; lat=-78:4:78;
258    mask=readbin('hFacC.data',[90 40 15],1);
259    tracer=zeros(90,40,30); sumtracer=zeros(30,1);
260    for i=1:30, mydisp(i)
261     fn=['PtrFlux' myint2str(i) '.0000000180.data'];
262     tracer(:,:,i)=readbin(fn,[90 40],1);
263     for j=1:length(lat)
264        sumtracer(i) = sumtracer(i) + sum(tracer(:,j,i).*mask(:,j,1)) * ...
265            nb_seconds_per_year* (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
266      end
267    end
268    plot(1:30,0*sumtracer,1:30,sumtracer-mmean(sumtracer))
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22