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

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

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

Revision 1.8 - (hide 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 dimitri 1.8 C $Header: /u/gcmpack/MITgcm/pkg/kpp/kpp_readparms.F,v 1.4.6.4 2003/02/05 09:04:35 dimitri Exp $
2 jmc 1.5 C $Name: $
3 heimbach 1.1 #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 dimitri 1.7 INTEGER errIO, iUnit
32 heimbach 1.1
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 dimitri 1.8 kpp_dumpFreq = dumpFreq
64     kpp_taveFreq = taveFreq
65 heimbach 1.1 KPPmixingMaps = .FALSE.
66     KPPwriteState = .FALSE.
67 jmc 1.5 minKPPhbl = drF(1)*0.5
68 heimbach 1.1
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