/[MITgcm]/MITgcm/pkg/kpp/kpp_readparms.F
ViewVC logotype

Contents of /MITgcm/pkg/kpp/kpp_readparms.F

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

Revision 1.8 - (show annotations) (download)
Tue Feb 18 05:33:55 2003 UTC (21 years, 4 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint50c_post, checkpoint52d_pre, checkpoint50c_pre, checkpoint52j_pre, checkpoint51o_pre, checkpoint54d_post, checkpoint54e_post, checkpoint51l_post, checkpoint48i_post, checkpoint52l_post, checkpoint52k_post, checkpoint55, checkpoint54, checkpoint51, checkpoint50, checkpoint53, checkpoint52, checkpoint50d_post, checkpoint52f_post, checkpoint50b_pre, checkpoint54f_post, checkpoint51f_post, checkpoint51d_post, checkpoint51t_post, checkpoint51n_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint55c_post, checkpoint51j_post, checkpoint52e_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint53d_post, checkpoint48f_post, checkpoint52b_pre, checkpoint54b_post, checkpoint51l_pre, checkpoint52m_post, checkpoint48h_post, checkpoint51q_post, checkpoint51b_pre, checkpoint52b_post, checkpoint52c_post, checkpoint51h_pre, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint52f_pre, checkpoint54a_pre, checkpoint53c_post, branchpoint-genmake2, checkpoint54a_post, checkpoint51r_post, checkpoint51i_post, checkpoint55b_post, checkpoint51b_post, checkpoint51c_post, checkpoint53a_post, checkpoint52d_post, checkpoint53g_post, checkpoint50g_post, checkpoint52a_pre, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint53f_post, checkpoint52j_post, checkpoint50e_post, branch-netcdf, checkpoint50d_pre, checkpoint52n_post, checkpoint53b_pre, checkpoint51e_post, checkpoint55a_post, checkpoint49, checkpoint51o_post, checkpoint51f_pre, checkpoint48g_post, checkpoint53b_post, checkpoint52a_post, checkpoint51g_post, ecco_c52_e35, checkpoint50b_post, checkpoint51m_post, checkpoint53d_pre, checkpoint54c_post, checkpoint51a_post, checkpoint51p_post, checkpoint51u_post
Branch point for: branch-genmake2, branch-nonh, tg2-branch, netcdf-sm0, checkpoint51n_branch
Changes since 1.7: +3 -3 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 C $Header: /u/gcmpack/MITgcm/pkg/kpp/kpp_readparms.F,v 1.4.6.4 2003/02/05 09:04:35 dimitri Exp $
2 C $Name: $
3 #include "KPP_OPTIONS.h"
4
5 SUBROUTINE KPP_READPARMS( myThid )
6 C /==========================================================\
7 C | SUBROUTINE KPP_READPARMS |
8 C | o Routine to read in file data.kpp |
9 C \==========================================================/
10 IMPLICIT NONE
11
12 C === Global variables ===
13 #include "SIZE.h"
14 #include "EEPARAMS.h"
15 #include "PARAMS.h"
16 #include "GRID.h"
17 #include "KPP_PARAMS.h"
18
19 C === Routine arguments ===
20 C myThid - Number of this instance of KPP_INIT
21 INTEGER myThid
22
23 #ifdef ALLOW_KPP
24
25 C === Local variables ===
26 C msgBuf - Informational/error meesage buffer
27 C errIO - IO error flag
28 C iUnit - Work variable for IO unit number
29
30 CHARACTER*(MAX_LEN_MBUF) msgBuf
31 INTEGER errIO, iUnit
32
33 C-- KPP vertical mixing parameters
34 NAMELIST /KPP_PARM01/
35 & kpp_freq, kpp_dumpFreq, kpp_taveFreq,
36 & KPPmixingMaps, KPPwriteState,
37 & minKPPhbl,
38 & epsln, phepsi, epsilon, vonk, dB_dz,
39 & conc1, conam, concm, conc2, zetam,
40 & conas, concs, conc3, zetas,
41 & Ricr, cekman, cmonob, concv, hbf,
42 & zmin, zmax, umin, umax,
43 & num_v_smooth_Ri, num_v_smooth_BV,
44 & num_z_smooth_sh, num_m_smooth_sh,
45 & Riinfty, BVSQcon, difm0, difs0, dift0,
46 & difmcon, difscon, diftcon,
47 & cstar
48
49 _BEGIN_MASTER(myThid)
50
51 WRITE(msgBuf,'(A)') ' KPP_INIT: opening data.kpp'
52 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
53 & SQUEEZE_RIGHT , 1)
54
55
56 CALL OPEN_COPY_DATA_FILE(
57 I 'data.kpp', 'KPP_INIT',
58 O iUnit,
59 I myThid )
60
61 C-- set default KPP vertical mixing parameters
62 kpp_freq = deltaTClock
63 kpp_dumpFreq = dumpFreq
64 kpp_taveFreq = taveFreq
65 KPPmixingMaps = .FALSE.
66 KPPwriteState = .FALSE.
67 minKPPhbl = drF(1)*0.5
68
69 C-----------------------------------------------------------------------
70 C define some non-dimensional constants and
71 C the vertical mixing coefficients in m-k-s units
72 C-----------------------------------------------------------------------
73
74 epsln = 1.e-20
75 phepsi = 1.e-10
76 epsilon = 0.1
77 vonk = 0.40
78 dB_dz = 5.2e-5
79 conc1 = 5.0
80 conam = 1.257
81 concm = 8.380
82 conc2 = 16.0
83 zetam = -0.2
84 conas = -28.86
85 concs = 98.96
86 conc3 = 16.0
87 zetas = -1.0
88
89 c parameters for subroutine "bldepth"
90
91 Ricr = 0.30
92 cekman = 0.7
93 cmonob = 1.0
94 concv = 1.8
95 hbf = 1.0
96
97 c parameters and common arrays for subroutines
98 c "kmixinit" and "wscale"
99
100 zmin = -4.e-7
101 zmax = 0.0
102 umin = 0.0
103 umax = .04
104
105 c parameters for subroutine "Ri_iwmix"
106
107 num_v_smooth_Ri = 0
108 num_v_smooth_BV = 0
109 num_z_smooth_sh = 0
110 num_m_smooth_sh = 0
111 Riinfty = 0.7
112 BVSQcon = -0.2e-4
113
114 difm0 = 0.005
115 difs0 = 0.005
116 dift0 = 0.005
117
118 difmcon = 0.1
119 difscon = 0.1
120 diftcon = 0.1
121
122 C parameters for subroutine "blmix"
123
124 cstar = 10.
125
126 C-----------------------------------------------------------------------
127
128 C-- Read settings from model parameter file "data.kpp".
129 READ(UNIT=iUnit,NML=KPP_PARM01,IOSTAT=errIO)
130 IF ( errIO .LT. 0 ) THEN
131 WRITE(msgBuf,'(A)')
132 & 'S/R INI_PARMS'
133 CALL PRINT_ERROR( msgBuf , 1)
134 WRITE(msgBuf,'(A)')
135 & 'Error reading numerical model '
136 CALL PRINT_ERROR( msgBuf , 1)
137 WRITE(msgBuf,'(A)')
138 & 'parameter file "data.kpp"'
139 CALL PRINT_ERROR( msgBuf , 1)
140 WRITE(msgBuf,'(A)')
141 & 'Problem in namelist KPP_PARM01'
142 CALL PRINT_ERROR( msgBuf , 1)
143 C CALL MODELDATA_EXAMPLE( myThid )
144 STOP 'ABNORMAL END: S/R KPP_INIT'
145 ENDIF
146
147 CLOSE(iUnit)
148
149 WRITE(msgBuf,'(A)') ' KPP_INIT: finished reading data.kpp'
150 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
151 & SQUEEZE_RIGHT , 1)
152
153 _END_MASTER(myThid)
154
155 C-- Everyone else must wait for the parameters to be loaded
156 _BARRIER
157
158 #endif /* ALLOW_KPP */
159
160 return
161 end
  ViewVC Help
Powered by ViewVC 1.1.22