/[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.4 by dimitri, Wed Sep 24 04:52:38 2003 UTC revision 1.19 by dimitri, Thu Dec 18 03:17:56 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 53  salinity initial and boundary conditions Line 52  salinity initial and boundary conditions
52   cp data.ptracers.test data.ptracers   cp data.ptracers.test data.ptracers
53   mitgcmuv > output.txt   mitgcmuv > output.txt
54        
55  5 ===> check that PTRACER output and salinity output  5 ===> check that PTRACER output and salinity output are identical.
   ===> are identical.  
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 83  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    
142    ### need customized option file; example below is for nireas
143     ../tools/genmake2 -of linux_ia32_g77_netcdf
144    
145     make depend
146     make
147     cd ../exe
148     cp ../verification/global_with_exf/input/eedata .
149     cp ../verification/global_with_exf/input/data.* .
150     cp ../verification/global_with_exf/input/POLY3.COEFFS .
151     ln -sf ../verification/global_with_exf/input/*.bin .
152     cp ../ocean_inversion_project/input/* .
153     cp data.1765-2005.1 data
154     ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
155     ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
156     ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
157   mitgcmuv > output .txt &   mitgcmuv > output .txt &
158    
159    
160  ======================================================  ======================================================
161    
162    ==========================================================
163    Instructions for generating netcdf output files on orion
164    ==========================================================
165    
166    cd ocean_inversion_project/write_netCDF
167    setenv F_UFMTENDIAN big
168    
169    (check that defaults in  mk_output.F are appropriate)
170    
171    efc -W0 -WB mk_output.F write_nc_phys.F nc_util.F \
172         handle_errors.F write_nc_basisfnctns.F \
173         write_nc_diag_0D.F write_nc_diag_2D.F \
174         -I/u2/dmenem/software/netcdf-3.5.0/include \
175         -L/u2/dmenem/software/netcdf-3.5.0/lib -lnetcdf
176    a.out
177    
178    
179    ======================================================
180    
181  % some matlab code for looking at fort.10 debug files  % some matlab code for looking at fort.10 debug files
182  load fort.10  load fort.10
183  tak=zeros(90,40,12);  tak=zeros(90,40,12);
# Line 117  mypcolor( Line 205  mypcolor(
205  % tracer uptake is approximately 1e18 mols.  % tracer uptake is approximately 1e18 mols.
206  lon=2:4:360; lat=-78:4:78;  lon=2:4:360; lat=-78:4:78;
207  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];
208  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];  
209  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);  tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
210  for i=1:30, mydisp(i)  for i=1:30, mydisp(i)
211   fn=['PTRACER' myint2str(i) '.0000003600.001.001.data'];   fn=['PTRACER' myint2str(i) '.0000001800.001.001.data'];
212   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);   tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
213   for j=1:length(lat)   for j=1:length(lat)
214    for k=1:length(thk)    for k=1:length(thk)
215     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i)) * ...     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i).*mask(:,j,k)) * ...
216                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);                    thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
217    end    end
218   end   end
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22