/[MITgcm]/MITgcm/pkg/sbo/README.sbo
ViewVC logotype

Contents of /MITgcm/pkg/sbo/README.sbo

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

Revision 1.2 - (show annotations) (download)
Tue Feb 18 05:33:55 2003 UTC (21 years, 2 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint51j_post, checkpoint50e_post, checkpoint50c_post, checkpoint48h_post, checkpoint50c_pre, branchpoint-genmake2, checkpoint50d_pre, checkpoint48f_post, checkpoint51e_post, checkpoint51b_post, checkpoint51c_post, checkpoint49, checkpoint48i_post, checkpoint51f_pre, checkpoint51, checkpoint50, checkpoint50d_post, checkpoint51b_pre, checkpoint51h_pre, checkpoint50g_post, checkpoint50b_pre, checkpoint51g_post, checkpoint51f_post, checkpoint50b_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint51d_post, checkpoint50h_post, checkpoint51a_post, checkpoint50e_pre, checkpoint50i_post, checkpoint48g_post, checkpoint51i_pre
Branch point for: branch-genmake2, ecco-branch
Changes since 1.1: +87 -0 lines 
Merging from release1_p12:
o Modifications for using pkg/exf with pkg/seaice
  - improved description of the various forcing configurations
  - added basic radiation bulk formulae to pkg/exf
  - units/sign fix for evap computation in exf_getffields.F
  - updated verification/global_with_exf/results/output.txt
o Added pkg/sbo for computing IERS Special Bureau for the Oceans
  (SBO) core products, including oceanic mass, center-of-mass,
  angular, and bottom pressure (see pkg/sbo/README.sbo).
o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
  to avoid negative values in shallow regions.
  - updated verification/natl_box/results/output.txt
  - updated verification/lab_sea/results/output.txt
o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
o Added useSingleCpuIO option (see PARAMS.h).
o Updated useSingleCpuIO option in mdsio_writefield.F to
  work with multi-field files, e.g., for single-file pickup.
o pkg/seaice:
  - bug fix in growth.F: QNET for no shortwave case
  - added HeffFile for specifying initial sea-ice thickness
  - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
  kpp_transport_t.F, forward_step.F, and the_main_loop.F
o pkg/seaice:
  - adjoint-friendly modifications
  - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
1 Online computation of angular momentume terms
2 =============================================
3
4 Subdirectory pkg/sbo contains a small diagnostic package that computes IERS
5 Special Bureau for the Oceans (SBO) core products, including oceanic mass,
6 center-of-mass, angular, and bottom pressure. The package is enabled by
7 setting "#define ALLOW_SBO" in CPP_OPTIONS.h at compile time and
8 "useSBO=.TRUE." in data.pkg during model execution. There are two separate
9 output products, SBO_DIAG.* and OBPtave*, which are described below.
10
11 ===========
12
13 SBO_DIAG.data is a single file that contains the following fields at every
14 model time step:
15
16 myCurrentTime - model integration time (s)
17 xoamc - x-comp oam due to currents (kg-m**2/s)
18 yoamc - y-comp oam due to currents (kg-m**2/s)
19 zoamc - z-comp oam due to currents (kg-m**2/s)
20 xoamp - x-comp oam due to pressure (kg-m**2/s)
21 yoamp - y-comp oam due to pressure (kg-m**2/s)
22 zoamp - z-comp oam due to pressure (kg-m**2/s)
23 mass - mass of oceans (kg)
24 xcom - x-comp of center-of-mass of oceans (m)
25 ycom - y-comp of center-of-mass of oceans (m)
26 zcom - z-comp of center-of-mass of oceans (m)
27
28 The file is direct access, real*8, binary. Its dimensions are 11 X nt, where
29 nt is number of records. For example, in matlab you can use
30
31 sbo=readbin('SBO_DIAG.data',[11 nt],1,'real*8',0,'native');
32
33 to obtain the 11 X nt matlab array sbo. A copy of readbin.m resides in
34 subdirectory verification/lab_sea/matlab.
35
36 ===========
37
38 OBPtave* files contain ocean-bottom pressure in Pascals (N/m^2). Ocean bottom
39 pressure is cumulated and averaged for sbo_taveFreq seconds. Parameter
40 sbo_taveFreq is specified in data.sbo during model execution. For example:
41
42 # SBO parameters
43 &SBO_PARM01
44 sbo_taveFreq = 86400.
45 &
46
47 will average and dump bottom pressure once per day. The bottom pressure files
48 are 2-dimensional fields with dimension nx X ny equal to those of the model
49 domain. The files have standard MITgcmUV mdsio output format. The files are
50 direct access, real*4, binary. They can be global or tiled (depending on
51 whether "globalFiles=.TRUE." is set in parameter file data). For example, in
52 matlab you can use
53
54 obp=readbin('OBPtave0000000020.data',[nx ny],1);
55
56 to read a global bottom pressure file output at model time step 20. Land is
57 flagged with -999.99 values.
58
59 ===========
60
61 SBO package file description
62
63 Added:
64 pkg/sbo/sbo_calc.F - main SBO routine, lots of comments
65 pkg/sbo/sbo_check.F - check for package consistency
66 pkg/sbo/sbo_diags.F - SBO diagnostics
67 pkg/sbo/sbo_readparms.F - read parameter file data.sbo
68 pkg/sbo/sbo_rho.F - in-situ density computation
69 pkg/sbo/sbo_writevector.F - write direct access binary vector
70 pkg/sbo/SBO.h - main header file
71 pkg/sbo/SBO_OPTIONS.h - SBO compiler time options
72 doc/README.sbo - description of SBO package
73 utils/matlab/read_record.m - matlab routine for reading binary files
74 utils/matlab/readbin.m - matlab routine for reading binary files
75 verification/natl_box/input/data.sbo - example parameter file
76
77 Modified:
78 model/inc/CPP_OPTIONS.h - added ALLOW_SBO
79 model/inc/PARAMS.h - added useSBO
80 model/src/do_the_model_io.F - added call to SBO package
81 model/src/packages_boot.F - added call to SBO package
82 model/src/packages_check.F - added call to SBO package
83 model/src/packages_readparms.F - added call to SBO package
84
85 ===========
86
87 Last modified: 02/04/03 by Dimitris Menemenlis
  ViewVC Help
Powered by ViewVC 1.1.22