/[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.21 - (hide annotations) (download)
Fri Sep 18 11:40:22 2009 UTC (14 years, 9 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62c, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62w, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint63g, checkpoint64, checkpoint62, checkpoint63, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint62b, checkpoint64x, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint61v, checkpoint61w, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.20: +15 -8 lines 
add double diffusive contributions as a hack, turned off by default
for now, and the code can be excluded with a CPP-flag
EXCLUDE_KPP_DOUBLEDIFF just as EXCLUDE_KPP_SHEAR_MIX
1 mlosch 1.21 C $Header: /u/gcmpack/MITgcm/pkg/kpp/kpp_readparms.F,v 1.20 2009/07/08 21:57:51 jmc Exp $
2 jmc 1.5 C $Name: $
3 heimbach 1.1 #include "KPP_OPTIONS.h"
4    
5 edhill 1.9 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
6     CBOP 0
7     C !ROUTINE: KPP_READPARMS
8    
9     C !INTERFACE:
10 heimbach 1.1 SUBROUTINE KPP_READPARMS( myThid )
11 edhill 1.9
12     C !DESCRIPTION:
13     C Routine to read in file data.kpp
14    
15     C !USES:
16 heimbach 1.1 IMPLICIT NONE
17     #include "SIZE.h"
18     #include "EEPARAMS.h"
19     #include "PARAMS.h"
20     #include "KPP_PARAMS.h"
21    
22 edhill 1.9 C !INPUT PARAMETERS:
23 heimbach 1.1 INTEGER myThid
24 edhill 1.9 CEOP
25 heimbach 1.1
26 edhill 1.9 C !LOCAL VARIABLES:
27 heimbach 1.1 #ifdef ALLOW_KPP
28     C msgBuf - Informational/error meesage buffer
29     C errIO - IO error flag
30     C iUnit - Work variable for IO unit number
31    
32     CHARACTER*(MAX_LEN_MBUF) msgBuf
33 dimitri 1.7 INTEGER errIO, iUnit
34 heimbach 1.1
35     C-- KPP vertical mixing parameters
36     NAMELIST /KPP_PARM01/
37     & kpp_freq, kpp_dumpFreq, kpp_taveFreq,
38 jmc 1.12 & KPPmixingMaps, KPPwriteState, KPP_ghatUseTotalDiffus,
39 mlosch 1.21 & KPPuseDoubleDiff, LimitHblStable,
40 heimbach 1.1 & minKPPhbl,
41     & epsln, phepsi, epsilon, vonk, dB_dz,
42     & conc1, conam, concm, conc2, zetam,
43     & conas, concs, conc3, zetas,
44     & Ricr, cekman, cmonob, concv, hbf,
45 jmc 1.16 & zmin, zmax, umin, umax,
46 heimbach 1.1 & num_v_smooth_Ri, num_v_smooth_BV,
47     & num_z_smooth_sh, num_m_smooth_sh,
48     & Riinfty, BVSQcon, difm0, difs0, dift0,
49     & difmcon, difscon, diftcon,
50 mlosch 1.21 & Rrho0, dsfmax,
51 heimbach 1.1 & cstar
52    
53     _BEGIN_MASTER(myThid)
54    
55     WRITE(msgBuf,'(A)') ' KPP_INIT: opening data.kpp'
56     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
57     & SQUEEZE_RIGHT , 1)
58    
59    
60     CALL OPEN_COPY_DATA_FILE(
61     I 'data.kpp', 'KPP_INIT',
62     O iUnit,
63     I myThid )
64    
65     C-- set default KPP vertical mixing parameters
66 mlosch 1.21 kpp_freq = deltaTClock
67     kpp_dumpFreq = dumpFreq
68     kpp_taveFreq = taveFreq
69     KPPmixingMaps = .FALSE.
70     KPPwriteState = .FALSE.
71     KPPuseDoubleDiff = .FALSE.
72     LimitHblStable = .TRUE.
73 jmc 1.12 KPP_ghatUseTotalDiffus = .FALSE.
74 dfer 1.18 minKPPhbl = UNSET_RL
75 heimbach 1.1
76     C-----------------------------------------------------------------------
77     C define some non-dimensional constants and
78     C the vertical mixing coefficients in m-k-s units
79     C-----------------------------------------------------------------------
80    
81 jmc 1.16 epsln = 1. _d -20
82     phepsi = 1. _d -10
83     epsilon = 1. _d -1
84     vonk = 0.4 _d 0
85     dB_dz = 5.2 _d -5
86     conc1 = 5. _d 0
87     conam = 1.257 _d 0
88     concm = 8.380 _d 0
89     conc2 = 16. _d 0
90     zetam = -0.2 _d 0
91     conas = -28.86 _d 0
92     concs = 98.96 _d 0
93     conc3 = 16. _d 0
94     zetas = -1. _d 0
95    
96     C parameters for subroutine "bldepth"
97    
98     Ricr = 0.3 _d 0
99     cekman = 0.7 _d 0
100     cmonob = 1. _d 0
101     concv = 1.8 _d 0
102     hbf = 1. _d 0
103    
104     C parameters and common arrays for subroutines
105     C "kmixinit" and "wscale"
106    
107     zmin = -4. _d -7
108     zmax = 0. _d 0
109     umin = 0. _d 0
110     umax = 4. _d -2
111    
112     C parameters for subroutine "Ri_iwmix"
113 heimbach 1.1
114     num_v_smooth_Ri = 0
115     num_v_smooth_BV = 0
116     num_z_smooth_sh = 0
117     num_m_smooth_sh = 0
118 jmc 1.16 Riinfty = 0.7 _d 0
119     BVSQcon = -0.2 _d -4
120 heimbach 1.1
121 jmc 1.16 difm0 = 5. _d -3
122     difs0 = 5. _d -3
123     dift0 = 5. _d -3
124    
125     difmcon = 0.1 _d 0
126     difscon = 0.1 _d 0
127     diftcon = 0.1 _d 0
128 heimbach 1.1
129 mlosch 1.21 C parameters for double diffusion routine "KPP_DOUBLEDIFF"
130    
131     Rrho0 = 1.9 _d 0
132     dsfmax = 10. _d -3
133    
134 heimbach 1.1 C parameters for subroutine "blmix"
135    
136     cstar = 10.
137    
138     C-----------------------------------------------------------------------
139    
140     C-- Read settings from model parameter file "data.kpp".
141     READ(UNIT=iUnit,NML=KPP_PARM01,IOSTAT=errIO)
142     IF ( errIO .LT. 0 ) THEN
143     WRITE(msgBuf,'(A)')
144     & 'S/R INI_PARMS'
145     CALL PRINT_ERROR( msgBuf , 1)
146     WRITE(msgBuf,'(A)')
147     & 'Error reading numerical model '
148     CALL PRINT_ERROR( msgBuf , 1)
149     WRITE(msgBuf,'(A)')
150     & 'parameter file "data.kpp"'
151     CALL PRINT_ERROR( msgBuf , 1)
152     WRITE(msgBuf,'(A)')
153     & 'Problem in namelist KPP_PARM01'
154     CALL PRINT_ERROR( msgBuf , 1)
155     C CALL MODELDATA_EXAMPLE( myThid )
156     STOP 'ABNORMAL END: S/R KPP_INIT'
157     ENDIF
158    
159     CLOSE(iUnit)
160    
161     WRITE(msgBuf,'(A)') ' KPP_INIT: finished reading data.kpp'
162     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
163     & SQUEEZE_RIGHT , 1)
164    
165     _END_MASTER(myThid)
166    
167     C-- Everyone else must wait for the parameters to be loaded
168     _BARRIER
169    
170     #endif /* ALLOW_KPP */
171    
172 jmc 1.16 RETURN
173     END
  ViewVC Help
Powered by ViewVC 1.1.22