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 |
adding sbo to packages.conf and "useSBO=.TRUE." in data.pkg during model |
8 |
execution. There are two separate output products, SBO_DIAG.* and |
9 |
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 |