/[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.21 by dimitri, Thu Sep 30 23:45:07 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   CVSROOT=:pserver:cvsanon@mitgcm.org:/u/gcmpack
24   cvs login ( CVS password: cvsanon )   cvs login ( CVS password: cvsanon )
25   cvs co -P MITgcm   cvs co -r checkpoint51n_branch MITgcm_code
26     cvs co -r checkpoint51n_branch MITgcm/verification/testreport
27     cvs co -r checkpoint51n_branch MITgcm/verification/global_with_exf
28    
29  2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory  2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory
30    
# Line 31  salinity initial and boundary conditions Line 33  salinity initial and boundary conditions
33    
34  3 ===> compile and link  3 ===> compile and link
35    
36     mkdir bin exe
37   cd bin   cd bin
38   cp ../verification/global_with_exf/code/* .   cp ../verification/global_with_exf/code/* .
39   cp ../ocean_inversion_project/code/* .   cp ../ocean_inversion_project/code/* .
# Line 56  salinity initial and boundary conditions Line 59  salinity initial and boundary conditions
59    
60   diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data   diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data
61   diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data   diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data
62     diff PTRtave30.0000000020.001.001.data Stave.0000000020.001.001.data
63     diff PTRtave30.0000000020.001.001.data Stave.0000000020.001.001.data
 ==========================================================  
 Instructions for carrying out a 3000-year quasi-stationary  
 integration using the global_ocean.90x40x15 configuration.  
 ==========================================================  
   
  cd MITgcm/exe  
  rm *  
  cd ../bin  
  rm *  
  cp ../verification/global_with_exf/code/* .  
  cp ../ocean_inversion_project/code/* .  
  rm PTRACERS_OPTIONS.h  
  ../tools/genmake2  
    (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)  
  make depend  
  make  
  cd ../exe  
  cp ../verification/global_with_exf/input/eedata .  
  cp ../verification/global_with_exf/input/data.* .  
  cp ../verification/global_with_exf/input/POLY3.COEFFS .  
  ln -sf ../verification/global_with_exf/input/*.bin .  
  cp ../ocean_inversion_project/input/* .  
  ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .  
  ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .  
  mitgcmuv > output .txt &  
64    
65    
66  ==========================================================  ==========================================================
# Line 111  integration using the global_ocean.90x40 Line 89  integration using the global_ocean.90x40
89   cp ../ocean_inversion_project/input/* .   cp ../ocean_inversion_project/input/* .
90   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
91   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
92   mpirun -np 36 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog &   cp data.stationary.1 data
93     mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
94       (on orion: qsub job.orion)
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 126  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 140  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 159  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 171  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 .
# Line 180  and that it be specified in the genmake2 Line 163  and that it be specified in the genmake2
163    
164  ======================================================  ======================================================
165    
166    ==========================================================
167    Instructions for generating netcdf output files on orion
168    ==========================================================
169    
170    cd ocean_inversion_project/write_netCDF
171    setenv F_UFMTENDIAN big
172    
173    (check that defaults in  mk_output.F are appropriate)
174    
175    efc -W0 -WB mk_output.F write_nc_phys.F nc_util.F \
176         handle_errors.F write_nc_basisfnctns.F \
177         write_nc_diag_0D.F write_nc_diag_2D.F \
178         -I/u2/dmenem/software/netcdf-3.5.0/include \
179         -L/u2/dmenem/software/netcdf-3.5.0/lib -lnetcdf
180    a.out
181    
182    
183    ======================================================
184    
185  % some matlab code for looking at fort.10 debug files  % some matlab code for looking at fort.10 debug files
186  load fort.10  load fort.10
187  tak=zeros(90,40,12);  tak=zeros(90,40,12);
# Line 220  for i=1:30, mydisp(i) Line 222  for i=1:30, mydisp(i)
222   end   end
223  end  end
224  plot(1:30,0*sumtracer,1:30,sumtracer)  plot(1:30,0*sumtracer,1:30,sumtracer)
225    
226    
227    % check that tracer flux is 1e18 mols / year
228    nb_seconds_per_year=31556880;
229    lon=2:4:360; lat=-78:4:78;
230    mask=readbin('hFacC.data',[90 40 15],1);
231    tracer=zeros(90,40,30); sumtracer=zeros(30,1);
232    for i=1:30, mydisp(i)
233     fn=['PtrFlux' myint2str(i) '.0000000180.data'];
234     tracer(:,:,i)=readbin(fn,[90 40],1);
235     for j=1:length(lat)
236        sumtracer(i) = sumtracer(i) + sum(tracer(:,j,i).*mask(:,j,1)) * ...
237            nb_seconds_per_year* (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
238      end
239    end
240    plot(1:30,0*sumtracer,1:30,sumtracer-mmean(sumtracer))
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22