三十层全球大洋环流模式使用说明
编写人：周天军、金向泽、张学洪
中国科学院大气物理研究所LASG
2000年5月20日
目 录
一、模式的使用……………………………………………………. （一）模式运行………………………………………………….. （二）模式输入参数…………………………………………….. （三）模式结果输出…………………………………………….. 二、模式的结构………………………………………………….. （一）模式设计原理…………………………………………….. （二）模式子程序结构………………………………………….. （三）硬盘和存储文件管理…………………………………….. （四）连续运行步骤…………………………………………….. 三、模式变量和参数表…………………………………………..
一、 模式的使用
（一）模式运行
1.
源程序维护
三十层全球大洋环流模式的专用目录下，包含三个子目录：
⑴FILE ：用于放置源程序，具体包括24个后缀名为“.f ”的FORTRAN 程序、6个后缀名为“.h ”的参数和数组定义文件，以及makefile 文件。
⑵FORT22：用于放置每年12月第31天的相关物理量场，以便于模式重新启动、继续积分，文件名为“fort.22.????Dec ”，其中的“???? ”表示当前的模式年。
⑶OMMEAN ：用于放置模式的积分结果。
注意模式版本一旦确定后，不要轻易对源程序进行改动。 2.
模式运行描述
在UNIX 操作系统下，在海洋模式专用目录下的FILE 子目录中，运行make 命令，利用makefile 文件对源程序进行优化编译，生成的可执行文件名为L30T63_OGCM，将其在FILE 子目录的上级目录－海洋模式专用目录下运行。
模式开始运行后，屏幕输出信息用于对模式运行情况进行监控，包括各种参数的设定、能量守恒情况等。可将其储存到任一给定的文件中，以便于事后需要时，对模式的运行过程进行分析检查。
此外，专用目录下有一文本文件modeltime ，在模式积分过程中，它随时被改写，用于监测模式积分的当前日期。
（二）模式输入参数
包括三方面内容：模式将要积分的时间长度（以月为单位）、积分结果输出的时间间隔（以年为单位）和模式起始状态（1为从静止海洋状态开始积分，0为接着上次的积分结果继续积分）的设定。上述参数在ccinput 文件中给定。
另外，还要读取海表强迫场，海表强迫场资料存放在无格式数据文件MODEL.FRC 中，它由rdriver.f 子程序从90号通道自动读入，无须模式运控人员干预。
（三）模式结果输出
模式输出量包括逐月的海表高度、海冰面积、海冰厚度、对流频率（即发生对流的次数）、海水位温、海水盐度、海流纬向速度和海流经向速度。结果的输出由SSAVE.f 子程序完成。
需要说明的是，海洋模式的结果输出，资料量是非常巨大的。因此，在其定型版中，为了节省存储空间，对于三维变量，例如海水位温、海水盐度、海流纬向速度和海流经向速度，
每年的前11个月，只输出最上面5层的结果；最后一个月（12月），输出全部30层的结果，资料存储方式为无格式。如果需要输出全部三十层的结果，只需在SSAVE.f 子程序中，将变量KLV 定义为30即可。
二、 模式的结构
（一）模式设计原理
参见技术文本《三十层大洋环流模式的原理》。 （二）模式子程序结构
文件名
功能
子程序名
Included files
param.h grid parameters: IMT=194, JMT=92, KM=30
tracer numbers: NTRA=2
IMM=193, JMM=91, KMP1=31, KMM1=29
comblk.h 16 common blocks for all the model’s variables
dncoef.h Normalized (1972)
isopyc.h common block for isopycnal diffusion variables pconst.h common blocks for land/sea mask
1-d grid-, time-, and physical parameters
pmix.h
& Philander (1981).
Main program
program
Subroutines
accumy.f ACCUMM
use
addps.f barotr.f
compensating the loss of gross mass
solve the barotropic mode of momentum equations
ADDPS BAROTR BCLINIC COMQD
bclinic.f solve the baroclinic mode of momentum equations comqd.f calculate Haney’s Q and D const.f
CONST
finished)
convadj.f GFDL's full convective adjustment
B-C formula
CONV ADJ
density.f DENSITY energy.f ENERGY
mass .
filter.f grids.f ice.f inirun.f intfor.f
perform 1-D zonal fouriour filtering produce all the grid parameters Sea ice model Initialize the model
Interpolates monthly mean fields
FILTER GRIDS ICESNOW INIRUN INTFOR ISOPYI ISOPYC ISOADV ISOPLUX
readyc.f READYC
and baroclinic modes of momentum equations
readyt.f READYT
related to ηmodes of momentum equations
rdriver.f Read in climatological forcing fields smth.f
variables; Perform 1-D zonal smoothing
ssave.f
RDRIVER SMOOT3 SMT
SSA VE year.
tracer.f
solve tracer equations (isopycnal mixing not finished)
TRACER, INVTRI
upwell.f calculate vertical velocity vinteg.f
Perform vertical integration
UPWELL VINTEG
SMOOTH
isopycmix.f isopycnal mixing scheme of Gent_McWilliams (1990)
（三）硬盘和存储文件管理
模式结果输出存放在子目录OMMEAN 下，均为月平均输出量，各月的文件名依次被
命名为MMEAN????Jan 、MMEAN????Feb 、MMEAN????Mar 、MMEAN????Apr 、MMEAN????May 、MMEAN????Jun 、MMEAN????Jul 、MMEAN????Aug 、MMEAN????Sep 、MMEAN????Oct 、MMEAN????Nov 、MMEAN????Dec ，其中的“???? ”代表当前的模式年。
（四）连续运行步骤
模式启动时初始场的设定，可分为两种情况：
⑴模式从静止海洋状态开始积分，此时需读取初始位温和盐度场（观测的年平均气候状况），该资料存储在无格式文件TSinitial 中，它由inirun.f 子程序从81号通道自动读入，无须模式运控人员干预；
⑵模式接着上一次积分的某个中间状态继续积分，此时，须从FORT22子目录中，将模式启动年的海洋状态数据（包括海表高度、位温、盐度、速度、海冰面积、厚度等），拷贝到当前目录下的fort.22文件中，即：cp FORT22/fort.22.????Dec fort.22，其中“???? ”代表确定的某一模式年份，例如2220等。
注意模式启动积分时的初始场，是从静止海洋还是从中间状态继续积分，要在ccinput 文件中设定。需要说明的是，海洋模式的收敛过程，需要长时间的积分（上千年），因此，各种试验的开展，应从完成spinup 过程之后开始（模式的定型版已完成1165年的积分，其中前面1000年采用了加速收敛的非同步积分技术，后面165年为同步积分），亦即在开展试验研究时，模式启动的初始场，应选取第二种情况，如从第1165年的12月底开始。
三、 模式变量和参数表
表3.1 三维模式预报量（在COMMON /NVAR1/中定义）
变量名
说明
单位
UB (IMT,JMT) VB (IMT,JMT) H0 (IMT,JMT) UBP(IMT,JMT) VBP(IMT,JMT) H0P(IMT,JMT) U (IMT,JMT,KM) V (IMT,JMT,KM) WS(IMT,JMT,KMP1) UP(IMT,JMT,KM) VP(IMT,JMT,KM)
Barotropic velocity in zonal direction Barotropic velocity in meridional direction Sea sueface elevation
Same as UB(i,j) but for previous step Same as VB(i,j) but for previous step H0(i,j) at previous step
Total zonal velocity, positive for eastward
ms ms m ms ms m ms sec ms ms
Total meridional velocity, positive for southward ms Vertical velocity in η-coordinate Same as U(i,j,k) but for previous step Same as V(i,j,k) but for previous step
AT(IMT,JMT,KM,NTRA)
ATB(IMT,JMT,KM,NTRA)
表3.2 二维模式预报量（在COMMON /NVAR2/中定义） 变量名 说明
H0L(IMT,JMT) Sea surface elevations deduced from H0(i,j) H0F(IMT,JMT) Sea surface elevations deduced from H0(i,j) H0BL(IMT,JMT) Sea surface elevations deduced from H0(i,j) H0BF(IMT,JMT) Sea surface elevations deduced from H0(i,j) UTL(IMT,JMT,KM) Zonal velocity deduced from U(i,j,k) UTF(IMT,JMT,KM) Zonal velocity deduced from U(i,j,k) VTL(IMT,JMT,KM) Meridional velocity deduced from V(i,j,k) VTF(IMT,JMT,KM) Meridional velocity deduced from V(i,j,k)
单位
m m m m ms ms ms ms
表3.3 海冰预报量（在COMMON /NSICE/中定义） 变量名 说明
ITICE(IMT,JMT) Index field of sea ice(=1 for ice and 0 for water) ALEAD(IMT,JMT) Percentage of area covered by ice TLEAD(IMT,JMT) The increase of SST in the leads
HI(IMT,JMT) Thickness of sea ice
K Meter
单位
变量名
表3.4 海表强迫场（在COMMON /NCYC/中定义）
说明
surface for the m-th month
单位
SU3(I2,J1,12) 10Nm SV3(I2,J1,12) 10Nm
sea surface for the m-th month
PSA3(I2,J1,12) the m-th month
TSA3(I2,J1,12) HPa
for the m-th month
℃ psu Wm
SSS3(I2,J1,12) m-th month
SWV3(I2,J1,12)
the surface for the m-th month
UV A3(I2,J1,12)
the m-th month
the m-th month
ms -
QAR3(I2,J1,12)
CLD3(I2,J1,12)
percentage month
表3.5 模式输出量（在COMMON /NMMN/中定义）
变量名 说明
Monthly mean sea surface elevation
WSMON(I2,J1,K0Monthly mean vertical velocity in z-coordinate
Monthly mean temperature Monthly mean salinity Monthly mean zonal velocity Monthly mean meridional velocity
Monthly mean sea-ice thickness
HDMON(IMT,JMT)Monthly mean average thickness of sea ice in one grid ICMON(IMT,JMT,2
单位
m ms ℃ psu ms ms m m unitless
Frequency of convection within a month
表3.6物理参数之一（在COMMON /NMIX/中定义）
参数名
赋值
说明
AM
0.5×105m 2⋅s −1 for 50°N-60°S, Lateral eddy viscosity
coefficient
2.0×105m 2⋅s −1 for the other region
1.0×10−4m 2⋅s −1 2.0×103m 2⋅s −1 0.3×10−4m 2⋅s −1
=A HV (0)
Vertical viscosity coefficient Lateral diffusion coefficient Vertical diffusion coefficient Diffusivity between ice & water
AMV AH AHV(KM) AHICE
表3.7物理参数之二（在COMMON /NCN1/ 中定义）
参数名 赋值 说明
D0
1029.0 kg ⋅m −3
Reference density of sea water, ρo
CP
3901.0J ⋅kg −1⋅K −1
specific heat of sea water at constant pressure, c p
G C0 TBICE OD0 SAG CAG OD0CP RSD ASEA
9.806 m ⋅s −2 2.6×10−3 -1.8°C
Acceleration of gravity, g bottom friction coefficient frozen point of sea water
ρ−1
sin g cos g
g =10°(bias angle of bottom friction) g =10°(bias angle of bottom friction)
total area of the model’s surface layer
(ρC )
−1
p
(A MV ρ0)−1
Asselin Filter coefficient, B AFB2
1−2×B
Asselin Filter coefficient, C AFC2
1−2×C
Asselin Filter coefficient, αT AFT2
1−2×T
表3.8 积分和时间控制参数（在COMMON /NCMN/中定义）
参数名 赋值 说明
DTB DTC DTS DTB2 DTC2 NBB NCC NSS
∆t B =FLOAT(IDTB) ∆t C =FLOAT(IDTC)*60.0 ∆t T =FLOAT(IDTS)*3600.0
2∆t B 2∆t C
IDTB=120 seconds IDTC=240 minutes IDTS= 8 hrs
Barotropic steps within one baroclinic step
Baroclinic steps within one thermohaline step
steps for thermohaline process per day
N B =∆t C /∆t B N C =∆t T /∆t C
24/IDTS
ONBB ONBC ONCC ISB ISC IST
(N B +1)−1 (N C ⋅N B +1)−1 (N C +1)−1
0 for Euler-forward scheme 0 for Euler-forward scheme 0 for Euler-forward scheme
switch on Euler-forward or leap-frog scheme for barotropic mode switch on Euler-forward or leap-frog scheme for baroclinic mode
switch on Euler-forward or leap-frog scheme for thermohaline process
/31, 28, 31, 30, 31, 30, 31, 31, 30, day number of each calendar month (12) (12) (12)
31, 30, 31/ 16, 15, 16/
'Jun', 'Jul', 'Aug', 'Sep', 'Oct', 'Nov', 'Dec' /
MSTART /16, 14, 16, 15, 16, 15, 16, 16, 15, the day of mid-month
ABMON / 'Jan', 'Feb', 'Mar',' Apr', 'May', Character string for naming the
monthly mean outputs
三十层全球大洋环流模式使用说明
编写人：周天军、金向泽、张学洪
中国科学院大气物理研究所LASG
2000年5月20日
目 录
一、模式的使用……………………………………………………. （一）模式运行………………………………………………….. （二）模式输入参数…………………………………………….. （三）模式结果输出…………………………………………….. 二、模式的结构………………………………………………….. （一）模式设计原理…………………………………………….. （二）模式子程序结构………………………………………….. （三）硬盘和存储文件管理…………………………………….. （四）连续运行步骤…………………………………………….. 三、模式变量和参数表…………………………………………..
一、 模式的使用
（一）模式运行
1.
源程序维护
三十层全球大洋环流模式的专用目录下，包含三个子目录：
⑴FILE ：用于放置源程序，具体包括24个后缀名为“.f ”的FORTRAN 程序、6个后缀名为“.h ”的参数和数组定义文件，以及makefile 文件。
⑵FORT22：用于放置每年12月第31天的相关物理量场，以便于模式重新启动、继续积分，文件名为“fort.22.????Dec ”，其中的“???? ”表示当前的模式年。
⑶OMMEAN ：用于放置模式的积分结果。
注意模式版本一旦确定后，不要轻易对源程序进行改动。 2.
模式运行描述
在UNIX 操作系统下，在海洋模式专用目录下的FILE 子目录中，运行make 命令，利用makefile 文件对源程序进行优化编译，生成的可执行文件名为L30T63_OGCM，将其在FILE 子目录的上级目录－海洋模式专用目录下运行。
模式开始运行后，屏幕输出信息用于对模式运行情况进行监控，包括各种参数的设定、能量守恒情况等。可将其储存到任一给定的文件中，以便于事后需要时，对模式的运行过程进行分析检查。
此外，专用目录下有一文本文件modeltime ，在模式积分过程中，它随时被改写，用于监测模式积分的当前日期。
（二）模式输入参数
包括三方面内容：模式将要积分的时间长度（以月为单位）、积分结果输出的时间间隔（以年为单位）和模式起始状态（1为从静止海洋状态开始积分，0为接着上次的积分结果继续积分）的设定。上述参数在ccinput 文件中给定。
另外，还要读取海表强迫场，海表强迫场资料存放在无格式数据文件MODEL.FRC 中，它由rdriver.f 子程序从90号通道自动读入，无须模式运控人员干预。
（三）模式结果输出
模式输出量包括逐月的海表高度、海冰面积、海冰厚度、对流频率（即发生对流的次数）、海水位温、海水盐度、海流纬向速度和海流经向速度。结果的输出由SSAVE.f 子程序完成。
需要说明的是，海洋模式的结果输出，资料量是非常巨大的。因此，在其定型版中，为了节省存储空间，对于三维变量，例如海水位温、海水盐度、海流纬向速度和海流经向速度，
每年的前11个月，只输出最上面5层的结果；最后一个月（12月），输出全部30层的结果，资料存储方式为无格式。如果需要输出全部三十层的结果，只需在SSAVE.f 子程序中，将变量KLV 定义为30即可。
二、 模式的结构
（一）模式设计原理
参见技术文本《三十层大洋环流模式的原理》。 （二）模式子程序结构
文件名
功能
子程序名
Included files
param.h grid parameters: IMT=194, JMT=92, KM=30
tracer numbers: NTRA=2
IMM=193, JMM=91, KMP1=31, KMM1=29
comblk.h 16 common blocks for all the model’s variables
dncoef.h Normalized (1972)
isopyc.h common block for isopycnal diffusion variables pconst.h common blocks for land/sea mask
1-d grid-, time-, and physical parameters
pmix.h
& Philander (1981).
Main program
program
Subroutines
accumy.f ACCUMM
use
addps.f barotr.f
compensating the loss of gross mass
solve the barotropic mode of momentum equations
ADDPS BAROTR BCLINIC COMQD
bclinic.f solve the baroclinic mode of momentum equations comqd.f calculate Haney’s Q and D const.f
CONST
finished)
convadj.f GFDL's full convective adjustment
B-C formula
CONV ADJ
density.f DENSITY energy.f ENERGY
mass .
filter.f grids.f ice.f inirun.f intfor.f
perform 1-D zonal fouriour filtering produce all the grid parameters Sea ice model Initialize the model
Interpolates monthly mean fields
FILTER GRIDS ICESNOW INIRUN INTFOR ISOPYI ISOPYC ISOADV ISOPLUX
readyc.f READYC
and baroclinic modes of momentum equations
readyt.f READYT
related to ηmodes of momentum equations
rdriver.f Read in climatological forcing fields smth.f
variables; Perform 1-D zonal smoothing
ssave.f
RDRIVER SMOOT3 SMT
SSA VE year.
tracer.f
solve tracer equations (isopycnal mixing not finished)
TRACER, INVTRI
upwell.f calculate vertical velocity vinteg.f
Perform vertical integration
UPWELL VINTEG
SMOOTH
isopycmix.f isopycnal mixing scheme of Gent_McWilliams (1990)
（三）硬盘和存储文件管理
模式结果输出存放在子目录OMMEAN 下，均为月平均输出量，各月的文件名依次被
命名为MMEAN????Jan 、MMEAN????Feb 、MMEAN????Mar 、MMEAN????Apr 、MMEAN????May 、MMEAN????Jun 、MMEAN????Jul 、MMEAN????Aug 、MMEAN????Sep 、MMEAN????Oct 、MMEAN????Nov 、MMEAN????Dec ，其中的“???? ”代表当前的模式年。
（四）连续运行步骤
模式启动时初始场的设定，可分为两种情况：
⑴模式从静止海洋状态开始积分，此时需读取初始位温和盐度场（观测的年平均气候状况），该资料存储在无格式文件TSinitial 中，它由inirun.f 子程序从81号通道自动读入，无须模式运控人员干预；
⑵模式接着上一次积分的某个中间状态继续积分，此时，须从FORT22子目录中，将模式启动年的海洋状态数据（包括海表高度、位温、盐度、速度、海冰面积、厚度等），拷贝到当前目录下的fort.22文件中，即：cp FORT22/fort.22.????Dec fort.22，其中“???? ”代表确定的某一模式年份，例如2220等。
注意模式启动积分时的初始场，是从静止海洋还是从中间状态继续积分，要在ccinput 文件中设定。需要说明的是，海洋模式的收敛过程，需要长时间的积分（上千年），因此，各种试验的开展，应从完成spinup 过程之后开始（模式的定型版已完成1165年的积分，其中前面1000年采用了加速收敛的非同步积分技术，后面165年为同步积分），亦即在开展试验研究时，模式启动的初始场，应选取第二种情况，如从第1165年的12月底开始。
三、 模式变量和参数表
表3.1 三维模式预报量（在COMMON /NVAR1/中定义）
变量名
说明
单位
UB (IMT,JMT) VB (IMT,JMT) H0 (IMT,JMT) UBP(IMT,JMT) VBP(IMT,JMT) H0P(IMT,JMT) U (IMT,JMT,KM) V (IMT,JMT,KM) WS(IMT,JMT,KMP1) UP(IMT,JMT,KM) VP(IMT,JMT,KM)
Barotropic velocity in zonal direction Barotropic velocity in meridional direction Sea sueface elevation
Same as UB(i,j) but for previous step Same as VB(i,j) but for previous step H0(i,j) at previous step
Total zonal velocity, positive for eastward
ms ms m ms ms m ms sec ms ms
Total meridional velocity, positive for southward ms Vertical velocity in η-coordinate Same as U(i,j,k) but for previous step Same as V(i,j,k) but for previous step
AT(IMT,JMT,KM,NTRA)
ATB(IMT,JMT,KM,NTRA)
表3.2 二维模式预报量（在COMMON /NVAR2/中定义） 变量名 说明
H0L(IMT,JMT) Sea surface elevations deduced from H0(i,j) H0F(IMT,JMT) Sea surface elevations deduced from H0(i,j) H0BL(IMT,JMT) Sea surface elevations deduced from H0(i,j) H0BF(IMT,JMT) Sea surface elevations deduced from H0(i,j) UTL(IMT,JMT,KM) Zonal velocity deduced from U(i,j,k) UTF(IMT,JMT,KM) Zonal velocity deduced from U(i,j,k) VTL(IMT,JMT,KM) Meridional velocity deduced from V(i,j,k) VTF(IMT,JMT,KM) Meridional velocity deduced from V(i,j,k)
单位
m m m m ms ms ms ms
表3.3 海冰预报量（在COMMON /NSICE/中定义） 变量名 说明
ITICE(IMT,JMT) Index field of sea ice(=1 for ice and 0 for water) ALEAD(IMT,JMT) Percentage of area covered by ice TLEAD(IMT,JMT) The increase of SST in the leads
HI(IMT,JMT) Thickness of sea ice
K Meter
单位
变量名
表3.4 海表强迫场（在COMMON /NCYC/中定义）
说明
surface for the m-th month
单位
SU3(I2,J1,12) 10Nm SV3(I2,J1,12) 10Nm
sea surface for the m-th month
PSA3(I2,J1,12) the m-th month
TSA3(I2,J1,12) HPa
for the m-th month
℃ psu Wm
SSS3(I2,J1,12) m-th month
SWV3(I2,J1,12)
the surface for the m-th month
UV A3(I2,J1,12)
the m-th month
the m-th month
ms -
QAR3(I2,J1,12)
CLD3(I2,J1,12)
percentage month
表3.5 模式输出量（在COMMON /NMMN/中定义）
变量名 说明
Monthly mean sea surface elevation
WSMON(I2,J1,K0Monthly mean vertical velocity in z-coordinate
Monthly mean temperature Monthly mean salinity Monthly mean zonal velocity Monthly mean meridional velocity
Monthly mean sea-ice thickness
HDMON(IMT,JMT)Monthly mean average thickness of sea ice in one grid ICMON(IMT,JMT,2
单位
m ms ℃ psu ms ms m m unitless
Frequency of convection within a month
表3.6物理参数之一（在COMMON /NMIX/中定义）
参数名
赋值
说明
AM
0.5×105m 2⋅s −1 for 50°N-60°S, Lateral eddy viscosity
coefficient
2.0×105m 2⋅s −1 for the other region
1.0×10−4m 2⋅s −1 2.0×103m 2⋅s −1 0.3×10−4m 2⋅s −1
=A HV (0)
Vertical viscosity coefficient Lateral diffusion coefficient Vertical diffusion coefficient Diffusivity between ice & water
AMV AH AHV(KM) AHICE
表3.7物理参数之二（在COMMON /NCN1/ 中定义）
参数名 赋值 说明
D0
1029.0 kg ⋅m −3
Reference density of sea water, ρo
CP
3901.0J ⋅kg −1⋅K −1
specific heat of sea water at constant pressure, c p
G C0 TBICE OD0 SAG CAG OD0CP RSD ASEA
9.806 m ⋅s −2 2.6×10−3 -1.8°C
Acceleration of gravity, g bottom friction coefficient frozen point of sea water
ρ−1
sin g cos g
g =10°(bias angle of bottom friction) g =10°(bias angle of bottom friction)
total area of the model’s surface layer
(ρC )
−1
p
(A MV ρ0)−1
Asselin Filter coefficient, B AFB2
1−2×B
Asselin Filter coefficient, C AFC2
1−2×C
Asselin Filter coefficient, αT AFT2
1−2×T
表3.8 积分和时间控制参数（在COMMON /NCMN/中定义）
参数名 赋值 说明
DTB DTC DTS DTB2 DTC2 NBB NCC NSS
∆t B =FLOAT(IDTB) ∆t C =FLOAT(IDTC)*60.0 ∆t T =FLOAT(IDTS)*3600.0
2∆t B 2∆t C
IDTB=120 seconds IDTC=240 minutes IDTS= 8 hrs
Barotropic steps within one baroclinic step
Baroclinic steps within one thermohaline step
steps for thermohaline process per day
N B =∆t C /∆t B N C =∆t T /∆t C
24/IDTS
ONBB ONBC ONCC ISB ISC IST
(N B +1)−1 (N C ⋅N B +1)−1 (N C +1)−1
0 for Euler-forward scheme 0 for Euler-forward scheme 0 for Euler-forward scheme
switch on Euler-forward or leap-frog scheme for barotropic mode switch on Euler-forward or leap-frog scheme for baroclinic mode
switch on Euler-forward or leap-frog scheme for thermohaline process
/31, 28, 31, 30, 31, 30, 31, 31, 30, day number of each calendar month (12) (12) (12)
31, 30, 31/ 16, 15, 16/
'Jun', 'Jul', 'Aug', 'Sep', 'Oct', 'Nov', 'Dec' /
MSTART /16, 14, 16, 15, 16, 15, 16, 16, 15, the day of mid-month
ABMON / 'Jan', 'Feb', 'Mar',' Apr', 'May', Character string for naming the
monthly mean outputs