/[MITgcm]/MITgcm_contrib/ocean_inversion_project/README

Diff of /MITgcm_contrib/ocean_inversion_project/README

Parent Directory | Revision Log | View Revision Graph Revision Graph | View Patch Patch

-revision 1.14 by dimitri,
Tue Oct 28 00:27:48 2003 UTC
+revision 1.21 by dimitri,
Thu Sep 30 23:45:07 2004 UTC
 Line 20 
 salinity initial and boundary conditions
 ===> get MITgcm code from cvs repository
-  CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/gcmpack
+  CVSROOT=:pserver:cvsanon@mitgcm.org:/u/gcmpack
   cvs login ( CVS password: cvsanon )
-  cvs co -P MITgcm
+  cvs co -r checkpoint51n_branch MITgcm_code
+  cvs co -r checkpoint51n_branch MITgcm/verification/testreport
+  cvs co -r checkpoint51n_branch MITgcm/verification/global_with_exf
 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory
-Line 31 
 salinity initial and boundary conditions
+Line 33 
 salinity initial and boundary conditions
 ===> compile and link
+  mkdir bin exe
   cd bin
   cp ../verification/global_with_exf/code/* .
   cp ../ocean_inversion_project/code/* .
-Line 86 
 integration using the global_ocean.90x40
+Line 89 
 integration using the global_ocean.90x40
   cp ../ocean_inversion_project/input/* .
   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
-  mpirun -np 24 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog &
+  cp data.stationary.1 data
+  mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
+    (on orion: qsub job.orion)
  ==========================================================
- Instructions for carrying out time-dependent, 1765-2005,
+ MPI instructions for carrying out time-dependent, 1765-2005,
  anthropogenic carbon perturbation tracer experiments using
  the global_ocean.90x40x15 configuration.
  ==========================================================
-Line 101 
 the global_ocean.90x40x15 configuration.
+Line 105 
 the global_ocean.90x40x15 configuration.
   rm *
   cp ../verification/global_with_exf/code/* .
   cp ../ocean_inversion_project/code/* .
+  cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
+  cp SIZE.h_mpi SIZE.h
   ../tools/genmake2
     (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
   make depend
-Line 115 
 the global_ocean.90x40x15 configuration.
+Line 121 
 the global_ocean.90x40x15 configuration.
   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
-  mitgcmuv > output .txt &
+  mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
+    (on orion: qsub job.orion.td)
  ==========================================================
-Line 134 
 and that it be specified in the genmake2
+Line 141 
 and that it be specified in the genmake2
   cp ../ocean_inversion_project/code/* .
   cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h
   cp ../ocean_inversion_project/write_netCDF/*  .
+  \rm mk_output.F
  ### need customized option file; example below is for nireas
   ../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
   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
+  cp data.1765-2005.1 data
   ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
   ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
   ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
-Line 155 
 and that it be specified in the genmake2
+Line 163 
 and that it be specified in the genmake2
  ======================================================
+ ==========================================================
+ Instructions for generating netcdf output files on orion
+ ==========================================================
+ cd ocean_inversion_project/write_netCDF
+ setenv F_UFMTENDIAN big
+ (check that defaults in  mk_output.F are appropriate)
+ efc -W0 -WB mk_output.F write_nc_phys.F nc_util.F \
+      handle_errors.F write_nc_basisfnctns.F \
+      write_nc_diag_0D.F write_nc_diag_2D.F \
+      -I/u2/dmenem/software/netcdf-3.5.0/include \
+      -L/u2/dmenem/software/netcdf-3.5.0/lib -lnetcdf
+ a.out
+ ======================================================
  % some matlab code for looking at fort.10 debug files
  load fort.10
  tak=zeros(90,40,12);
-Line 195 
 for i=1:30, mydisp(i)
+Line 222 
 for i=1:30, mydisp(i)
   end
  end
  plot(1:30,0*sumtracer,1:30,sumtracer)
+ % check that tracer flux is 1e18 mols / year
+ nb_seconds_per_year=31556880;
+ lon=2:4:360; lat=-78:4:78;
+ mask=readbin('hFacC.data',[90 40 15],1);
+ tracer=zeros(90,40,30); sumtracer=zeros(30,1);
+ for i=1:30, mydisp(i)
+  fn=['PtrFlux' myint2str(i) '.0000000180.data'];
+  tracer(:,:,i)=readbin(fn,[90 40],1);
+  for j=1:length(lat)
+     sumtracer(i) = sumtracer(i) + sum(tracer(:,j,i).*mask(:,j,1)) * ...
+         nb_seconds_per_year* (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
+   end
+ end
+ plot(1:30,0*sumtracer,1:30,sumtracer-mmean(sumtracer))

 Legend:



Removed from v.1.14
 


changed lines


 
Added in v.1.21
 Legend:



Removed from v.1.14
 


changed lines


 
Added in v.1.21
-Removed from v.1.14
+Added in v.1.21

	ViewVC Help
Powered by ViewVC 1.1.22