/[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.20 by dimitri, Thu Dec 18 06:11:14 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 -r checkpoint51n_branch 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 219  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.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22