/[MITgcm]/MITgcm/verification/lab_sea/README

Diff of /MITgcm/verification/lab_sea/README

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

-revision 1.4 by dimitri,
Sat Dec 28 10:11:11 2002 UTC
+revision 1.6 by dimitri,
Wed Apr 30 07:04:08 2003 UTC
 Line 24 
 NCEP reanalysis (see pkg/seaice/SEAICE_F
    vwindFile      = 'v10m.labsea79'    # 10-m meridional wind
    atempFile      = 'tair.labsea1979'  # 2-m air temperature
    aqhFile        = 'qa.labsea1979'    # 2-m specific humidity
-   lwfluxFile     = 'flo.labsea1979'   # longwave radiation
+   lwdownFile     = 'flo.labsea1979'   # downward longwave radiation
-   swfluxFile     = 'fsh.labsea1979'   # shortwave radiation
+   swdownFile     = 'fsh.labsea1979'   # downward shortwave radiation
    precipFile     = 'prate.labsea1979' # precipitation
    evapFile       = 'evap.labsea1979'  # evaporation
 Line 84 
 Generating 1-CPU executable with SEAICE_
    make
    mv mitgcmuv mitgcmuv_exf_forcing
+-CPU executable with SEAICE_EXTERNAL_FORCING, without READ_EVAP
+   cd ../../../verification/lab_sea/input
+   ln -sf ../code/SIZE.h .
+   ln -sf ../code/CPP_EEOPTIONS.h .
+   ln -sf ../code/CPP_OPTIONS.h .
+   ln -sf ../code/ECCO_CPPOPTIONS.BULK.h ECCO_CPPOPTIONS.h
+   ln -sf ../code/SEAICE_OPTIONS.h .
+   ../../../tools/genmake -makefile
+     ==> on alhena use:
+     ../../../tools/genmake -platform=o2k_noopt -makefile
+   make clean
+   make depend
+   make
+   mv mitgcmuv mitgcmuv_compute_evap
  Generating 1-CPU executable with SEAICE_EXTERNAL_FLUXES:
    cd ../../../verification/lab_sea/input
    ln -sf ../code/SIZE.h .
-Line 180 
 To run Experiment 3:
+Line 195 
 To run Experiment 3:
    mkdir exp3b
    mv *tave.0000000001.data exp3b
- Use matlab script lookat_exp3.m to compare
+ Use matlab script lookat_exp3.m to look at
- the output of exp2 to that of exp1:
+ results from exp3a and exp3b:
    cd ../../../verification/lab_sea/matlab
    matlab
    lookat_exp3
  If the dynamic ice solvers are correct, they should converge to a
- constant solution.  The tests above show that both the ADI and the LSR
+ solution that is constant at each latitude.
- solver are unable to handle periodic domains and therefore that they
- cannot be "correctly" parallelized.
  Instructions for running Experiment 4
-Line 217 
 NPSEUDO in data.seaice at the expense of
+Line 230 
 NPSEUDO in data.seaice at the expense of
  Also as sea-ice ages (thickens) the difference at the tiles becomes
  increasingly small.  For solver pkg/seaice/adi.F a value of NPSEUDO=10
  appears adequate for forward integrations but cannot be used for sea-ice
- adjoint model.  Work is underway to fix this problem.
+ adjoint model.  With NPSEUDO=300 or greater, the tile disappears, but the
+ computational cost is very large (11 times that of teh forward model).
  Instructions for running Experiment 5
-Line 337 
 the output of exp9 to that of exp8:
+Line 351 
 the output of exp9 to that of exp8:
    cd ../../../verification/lab_sea/matlab
    matlab
    lookat_exp9
+ Instructions for running Experiment 10
+ ======================================
+ This is a 1-cpu, 10-hour integration similar to exp8, but
+ with bulk formula computation of evaporation fields.
+ To run Experiment 10:
+   cd ../../../verification/lab_sea/input
+   ln -sf DATA.10hours_exf data
+   ln -sf DATA.SEAICE.adi_exf data.seaice
+   mitgcmuv_compute_evap >&! output.txt
+   mkdir exp10
+   mv *tave.0000000010.data exp10
+ Use matlab script lookat_exp10.m to compare
+ the output of exp10 to that of exp8:
+   cd ../../../verification/lab_sea/matlab
+   matlab
+   lookat_exp10

 Legend:



Removed from v.1.4
 


changed lines


 
Added in v.1.6
 Legend:



Removed from v.1.4
 


changed lines


 
Added in v.1.6
-Removed from v.1.4
+Added in v.1.6

	ViewVC Help
Powered by ViewVC 1.1.22