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

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

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

revision 1.15 by heimbach, Mon Oct 23 16:28:35 2006 UTC revision 1.16 by jmc, Fri Oct 19 19:11:17 2007 UTC
# Line 42  C--   KPP vertical mixing parameters Line 42  C--   KPP vertical mixing parameters
42       & conc1, conam, concm, conc2, zetam,       & conc1, conam, concm, conc2, zetam,
43       & conas, concs, conc3, zetas,       & conas, concs, conc3, zetas,
44       & Ricr, cekman, cmonob, concv, hbf,       & Ricr, cekman, cmonob, concv, hbf,
45       & zmin, zmax, umin, umax,       & zmin, zmax, umin, umax,
46       & num_v_smooth_Ri, num_v_smooth_BV,       & num_v_smooth_Ri, num_v_smooth_BV,
47       & num_z_smooth_sh, num_m_smooth_sh,       & num_z_smooth_sh, num_m_smooth_sh,
48       & Riinfty, BVSQcon, difm0, difs0, dift0,       & Riinfty, BVSQcon, difm0, difs0, dift0,
# Line 78  C--------------------------------------- Line 78  C---------------------------------------
78  C define some non-dimensional constants and  C define some non-dimensional constants and
79  C the vertical mixing coefficients in m-k-s units  C the vertical mixing coefficients in m-k-s units
80  C-----------------------------------------------------------------------  C-----------------------------------------------------------------------
   
       epsln   = 1.e-20  
       phepsi  = 1.e-10  
       epsilon = 0.1  
       vonk    = 0.40  
       dB_dz   = 5.2e-5  
       conc1   = 5.0  
       conam   = 1.257  
       concm   = 8.380  
       conc2   = 16.0  
       zetam   = -0.2  
       conas   = -28.86  
       concs   = 98.96  
       conc3   = 16.0  
       zetas   = -1.0  
   
 c     parameters for subroutine "bldepth"  
   
       Ricr    = 0.30  
       cekman  = 0.7  
       cmonob  = 1.0  
       concv   = 1.8  
       hbf     = 1.0  
   
 c     parameters and common arrays for subroutines  
 c     "kmixinit" and "wscale"  
   
       zmin    = -4.e-7  
       zmax    = 0.0  
       umin    = 0.0  
       umax    = .04  
81    
82  c     parameters for subroutine "Ri_iwmix"        epsln   =   1. _d -20
83          phepsi  =   1. _d -10
84          epsilon =   1. _d -1
85          vonk    =   0.4 _d 0
86          dB_dz   =   5.2 _d -5
87          conc1   =   5. _d 0
88          conam   =   1.257 _d 0
89          concm   =   8.380 _d 0
90          conc2   =  16. _d 0
91          zetam   =  -0.2 _d 0
92          conas   = -28.86 _d 0
93          concs   =  98.96 _d 0
94          conc3   =  16. _d 0
95          zetas   =  -1. _d 0
96    
97    C     parameters for subroutine "bldepth"
98    
99          Ricr    = 0.3 _d 0
100          cekman  = 0.7 _d 0
101          cmonob  = 1.  _d 0
102          concv   = 1.8 _d 0
103          hbf     = 1.  _d 0
104    
105    C     parameters and common arrays for subroutines
106    C     "kmixinit" and "wscale"
107    
108          zmin    = -4. _d -7
109          zmax    =  0. _d 0
110          umin    =  0. _d 0
111          umax    =  4. _d -2
112    
113    C     parameters for subroutine "Ri_iwmix"
114    
115        num_v_smooth_Ri = 0        num_v_smooth_Ri = 0
116        num_v_smooth_BV = 0        num_v_smooth_BV = 0
117        num_z_smooth_sh = 0        num_z_smooth_sh = 0
118        num_m_smooth_sh = 0        num_m_smooth_sh = 0
119        Riinfty = 0.7        Riinfty =  0.7 _d 0
120        BVSQcon = -0.2e-4        BVSQcon = -0.2 _d -4
121    
122        difm0   = 0.005        difm0   = 5. _d -3
123        difs0   = 0.005        difs0   = 5. _d -3
124        dift0   = 0.005        dift0   = 5. _d -3
125    
126        difmcon = 0.1        difmcon = 0.1 _d 0
127        difscon = 0.1        difscon = 0.1 _d 0
128        diftcon = 0.1        diftcon = 0.1 _d 0
129    
130  C     parameters for subroutine "blmix"  C     parameters for subroutine "blmix"
131    
# Line 165  C      CALL MODELDATA_EXAMPLE( myThid ) Line 165  C      CALL MODELDATA_EXAMPLE( myThid )
165  C       Define grid types for KPP variables  C       Define grid types for KPP variables
166          CALL MNC_CW_ADD_VNAME('KPPviscAz', 'Cen_xy_Hn__C__t',          CALL MNC_CW_ADD_VNAME('KPPviscAz', 'Cen_xy_Hn__C__t',
167       &       4,5, myThid)       &       4,5, myThid)
168          CALL MNC_CW_ADD_VATTR_TEXT('KPPviscAz','units','m^2/s',          CALL MNC_CW_ADD_VATTR_TEXT('KPPviscAz','units','m^2/s',
169       &       myThid)       &       myThid)
170          CALL MNC_CW_ADD_VATTR_TEXT('KPPviscAz','long_name',          CALL MNC_CW_ADD_VATTR_TEXT('KPPviscAz','long_name',
171       &       'KPP_vertical_eddy_viscosity_coefficient', myThid)       &       'KPP_vertical_eddy_viscosity_coefficient', myThid)
172          CALL MNC_CW_ADD_VATTR_TEXT('KPPviscAz',          CALL MNC_CW_ADD_VATTR_TEXT('KPPviscAz',
173       &       'coordinates','XC YC RC iter', myThid)       &       'coordinates','XC YC RC iter', myThid)
174            
175          CALL MNC_CW_ADD_VNAME('KPPdiffKzS', 'Cen_xy_Hn__C__t',          CALL MNC_CW_ADD_VNAME('KPPdiffKzS', 'Cen_xy_Hn__C__t',
176       &       4,5, myThid)       &       4,5, myThid)
177          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzS','units','m^2/s',          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzS','units','m^2/s',
178       &       myThid)       &       myThid)
179          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzS','long_name',          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzS','long_name',
180       &       'KPP_salt-tracer_vertical_diffusion_coefficient',       &       'KPP_salt-tracer_vertical_diffusion_coefficient',
181       &       myThid)       &       myThid)
182          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzS',          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzS',
183       &       'coordinates','XC YC RC iter', myThid)       &       'coordinates','XC YC RC iter', myThid)
184            
185          CALL MNC_CW_ADD_VNAME('KPPdiffKzT', 'Cen_xy_Hn__C__t',          CALL MNC_CW_ADD_VNAME('KPPdiffKzT', 'Cen_xy_Hn__C__t',
186       &       4,5, myThid)       &       4,5, myThid)
187          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzT','units','m^2/s',          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzT','units','m^2/s',
188       &       myThid)       &       myThid)
189          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzT','long_name',          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzT','long_name',
190       &       'KPP_vertical_heat_diffusion_coefficient', myThid)       &       'KPP_vertical_heat_diffusion_coefficient', myThid)
191          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzT',          CALL MNC_CW_ADD_VATTR_TEXT('KPPdiffKzT',
192       &       'coordinates','XC YC RC iter', myThid)       &       'coordinates','XC YC RC iter', myThid)
193            
194          CALL MNC_CW_ADD_VNAME('KPPGHAT', 'Cen_xy_Hn__C__t',          CALL MNC_CW_ADD_VNAME('KPPGHAT', 'Cen_xy_Hn__C__t',
195       &       4,5, myThid)       &       4,5, myThid)
196          CALL MNC_CW_ADD_VATTR_TEXT('KPPGHAT','units','s/m^2',          CALL MNC_CW_ADD_VATTR_TEXT('KPPGHAT','units','s/m^2',
197       &       myThid)       &       myThid)
198          CALL MNC_CW_ADD_VATTR_TEXT('KPPGHAT','long_name',          CALL MNC_CW_ADD_VATTR_TEXT('KPPGHAT','long_name',
199       &       'KPP_nonlocal_transport_coefficient', myThid)       &       'KPP_nonlocal_transport_coefficient', myThid)
200          CALL MNC_CW_ADD_VATTR_TEXT('KPPGHAT',          CALL MNC_CW_ADD_VATTR_TEXT('KPPGHAT',
201       &       'coordinates','XC YC RC iter', myThid)       &       'coordinates','XC YC RC iter', myThid)
202            
203          CALL MNC_CW_ADD_VNAME('KPPHBL', 'Cen_xy_Hn__-__t',          CALL MNC_CW_ADD_VNAME('KPPHBL', 'Cen_xy_Hn__-__t',
204       &       3,4, myThid)       &       3,4, myThid)
205          CALL MNC_CW_ADD_VATTR_TEXT('KPPHBL','units','m',          CALL MNC_CW_ADD_VATTR_TEXT('KPPHBL','units','m',
206       &       myThid)       &       myThid)
207          CALL MNC_CW_ADD_VATTR_TEXT('KPPHBL','long_name',          CALL MNC_CW_ADD_VATTR_TEXT('KPPHBL','long_name',
208       &       'KPP_boundary_layer_depth', myThid)       &       'KPP_boundary_layer_depth', myThid)
209          CALL MNC_CW_ADD_VATTR_TEXT('KPPHBL',          CALL MNC_CW_ADD_VATTR_TEXT('KPPHBL',
210       &       'coordinates','XC YC iter', myThid)       &       'coordinates','XC YC iter', myThid)
211            
212          CALL MNC_CW_ADD_VNAME('KPPMLD', 'Cen_xy_Hn__-__t',          CALL MNC_CW_ADD_VNAME('KPPMLD', 'Cen_xy_Hn__-__t',
213       &       3,4, myThid)       &       3,4, myThid)
214          CALL MNC_CW_ADD_VATTR_TEXT('KPPMLD','units','m',          CALL MNC_CW_ADD_VATTR_TEXT('KPPMLD','units','m',
215       &       myThid)       &       myThid)
216          CALL MNC_CW_ADD_VATTR_TEXT('KPPMLD','long_name',          CALL MNC_CW_ADD_VATTR_TEXT('KPPMLD','long_name',
217       &       'KPP_mixed_layer_depth', myThid)       &       'KPP_mixed_layer_depth', myThid)
218          CALL MNC_CW_ADD_VATTR_TEXT('KPPMLD',          CALL MNC_CW_ADD_VATTR_TEXT('KPPMLD',
219       &       'coordinates','XC YC iter', myThid)       &       'coordinates','XC YC iter', myThid)
220            
221          CALL MNC_CW_ADD_VNAME('KPPFRAC', 'Cen_xy_Hn__-__t',          CALL MNC_CW_ADD_VNAME('KPPFRAC', 'Cen_xy_Hn__-__t',
222       &       3,4, myThid)       &       3,4, myThid)
223          CALL MNC_CW_ADD_VATTR_TEXT('KPPFRAC','units','dimless',          CALL MNC_CW_ADD_VATTR_TEXT('KPPFRAC','units','dimless',
224       &       myThid)       &       myThid)
225          CALL MNC_CW_ADD_VATTR_TEXT('KPPFRAC','long_name',          CALL MNC_CW_ADD_VATTR_TEXT('KPPFRAC','long_name',
226       &       'KPP_short-wave_fraction_penetrating_mixing_layer',       &       'KPP_short-wave_fraction_penetrating_mixing_layer',
227       &       myThid)       &       myThid)
228          CALL MNC_CW_ADD_VATTR_TEXT('KPPFRAC',          CALL MNC_CW_ADD_VATTR_TEXT('KPPFRAC',
229       &       'coordinates','XC YC iter', myThid)       &       'coordinates','XC YC iter', myThid)
# Line 235  C--   Everyone else must wait for the pa Line 235  C--   Everyone else must wait for the pa
235    
236  #endif /* ALLOW_KPP */  #endif /* ALLOW_KPP */
237    
238        return        RETURN
239        end        END
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22