Changeset 2293 for branches/nemo_v3_3_beta/NEMOGCM/NEMO
- Timestamp:
- 2010-10-19T14:35:49+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90
r2287 r2293 44 44 45 45 ! !!! ** Namelist namzdf_gls ** 46 LOGICAL :: ln_crban = .FALSE.! =T use Craig and Banner scheme46 LOGICAL :: ln_crban = .FALSE. ! =T use Craig and Banner scheme 47 47 LOGICAL :: ln_length_lim = .FALSE. ! use limit on the dissipation rate understable stratification (Galperin et al., 1988) 48 LOGICAL :: ln_sigpsi = .FALSE.! Activate Burchard (2003) modification for k-eps closure AND wave breaking mixing49 REAL(wp) :: rn_epsmin = 1.e-12_wp! minimum value of dissipation (m2/s3)50 REAL(wp) :: rn_emin = 1.e-6_wp! minimum value of TKE (m2/s2)48 LOGICAL :: ln_sigpsi = .FALSE. ! Activate Burchard (2003) modification for k-eps closure AND wave breaking mixing 49 REAL(wp) :: rn_epsmin = 1.e-12_wp ! minimum value of dissipation (m2/s3) 50 REAL(wp) :: rn_emin = 1.e-6_wp ! minimum value of TKE (m2/s2) 51 51 INTEGER :: nn_tkebc_surf = 0 ! TKE surface boundary condition (=0/1) 52 52 INTEGER :: nn_tkebc_bot = 0 ! TKE bottom boundary condition (=0/1) … … 54 54 INTEGER :: nn_psibc_bot = 0 ! PSI bottom boundary condition (=0/1) 55 55 INTEGER :: nn_stab_func = 0 ! stability functions G88, KC or Canuto (=0/1/2) 56 INTEGER :: nn_clos = 0! closure 0/1/2/3 MY82/k-eps/k-w/gen57 REAL(wp) :: clim_galp = 0.53_wp! Holt 2008 value for k-eps: 0.26758 REAL(wp) :: hsro = 0.003_wp! Minimum surface roughness59 REAL(wp) :: rn_charn = 2.e+5_wp! Charnock constant for surface breaking waves mixing : 1400. (standard) or 2.e5 (Stacey value)60 REAL(wp) :: rn_crban = 100._wp! Craig and Banner constant for surface breaking waves mixing61 REAL(wp) :: rn_cm_sf = 0.73_wp! Shear free turbulence parameters62 REAL(wp) :: rn_a_sf = -2.0_wp! Must be negative -2 < rn_a_sf < -163 REAL(wp) :: rn_l_sf = 0.2_wp! 0 <rn_l_sf<vkarmn64 REAL(wp) :: rn_ghmin = -0.2865 REAL(wp) :: rn_gh0 = 0.032966 REAL(wp) :: rn_ghcri = 0.0367 REAL(wp) :: rn_a1 = 0.92_wp68 REAL(wp) :: rn_a2 = 0.74_wp69 REAL(wp) :: rn_b1 = 16.60_wp70 REAL(wp) :: rn_b2 = 10.10_wp71 REAL(wp) :: rn_e1 = 1.80_wp72 REAL(wp) :: rn_e2 = 1.33_wp73 REAL(wp) :: rn_l1 = 0.107_wp74 REAL(wp) :: rn_l2 = 0.0032_wp75 REAL(wp) :: rn_l3 = 0.0864_wp76 REAL(wp) :: rn_l4 = 0.12_wp77 REAL(wp) :: rn_l5 = 11.9_wp78 REAL(wp) :: rn_l6 = 0.4_wp79 REAL(wp) :: rn_l7 = 0.0_wp80 REAL(wp) :: rn_l8 = 0.48_wp81 REAL(wp) :: rn_m1 = 0.127_wp82 REAL(wp) :: rn_m2 = 0.00336_wp83 REAL(wp) :: rn_m3 = 0.0906_wp84 REAL(wp) :: rn_m4 = 0.101_wp85 REAL(wp) :: rn_m5 = 11.2_wp86 REAL(wp) :: rn_m6 = 0.4_wp87 REAL(wp) :: rn_m7 = 0.0_wp88 REAL(wp) :: rn_m8 = 0.318_wp56 INTEGER :: nn_clos = 0 ! closure 0/1/2/3 MY82/k-eps/k-w/gen 57 REAL(wp) :: rn_clim_galp = 0.53_wp ! Holt 2008 value for k-eps: 0.267 58 REAL(wp) :: hsro = 0.003_wp ! Minimum surface roughness 59 REAL(wp) :: rn_charn = 2.e+5_wp ! Charnock constant for surface breaking waves mixing : 1400. (standard) or 2.e5 (Stacey value) 60 REAL(wp) :: rn_crban = 100._wp ! Craig and Banner constant for surface breaking waves mixing 61 REAL(wp) :: rn_cm_sf = 0.73_wp ! Shear free turbulence parameters 62 REAL(wp) :: rn_a_sf = -2.0_wp ! Must be negative -2 < rn_a_sf < -1 63 REAL(wp) :: rn_l_sf = 0.2_wp ! 0 <rn_l_sf<vkarmn 64 REAL(wp) :: rn_ghmin = -0.28 65 REAL(wp) :: rn_gh0 = 0.0329 66 REAL(wp) :: rn_ghcri = 0.03 67 REAL(wp) :: rn_a1 = 0.92_wp 68 REAL(wp) :: rn_a2 = 0.74_wp 69 REAL(wp) :: rn_b1 = 16.60_wp 70 REAL(wp) :: rn_b2 = 10.10_wp 71 REAL(wp) :: rn_e1 = 1.80_wp 72 REAL(wp) :: rn_e2 = 1.33_wp 73 REAL(wp) :: rn_l1 = 0.107_wp 74 REAL(wp) :: rn_l2 = 0.0032_wp 75 REAL(wp) :: rn_l3 = 0.0864_wp 76 REAL(wp) :: rn_l4 = 0.12_wp 77 REAL(wp) :: rn_l5 = 11.9_wp 78 REAL(wp) :: rn_l6 = 0.4_wp 79 REAL(wp) :: rn_l7 = 0.0_wp 80 REAL(wp) :: rn_l8 = 0.48_wp 81 REAL(wp) :: rn_m1 = 0.127_wp 82 REAL(wp) :: rn_m2 = 0.00336_wp 83 REAL(wp) :: rn_m3 = 0.0906_wp 84 REAL(wp) :: rn_m4 = 0.101_wp 85 REAL(wp) :: rn_m5 = 11.2_wp 86 REAL(wp) :: rn_m6 = 0.4_wp 87 REAL(wp) :: rn_m7 = 0.0_wp 88 REAL(wp) :: rn_m8 = 0.318_wp 89 89 REAL(wp) :: rn_c02 90 90 REAL(wp) :: rn_c02r … … 781 781 ! Galperin criterium (NOTE : Not required if the proper value of C3 in stable cases is calculated) 782 782 zrn2 = MAX( rn2(ji,jj,jk), rsmall ) 783 mxln(ji,jj,jk) = MIN( clim_galp*SQRT( 2.*en(ji,jj,jk)/zrn2 ), mxln(ji,jj,jk) )783 mxln(ji,jj,jk) = MIN( rn_clim_galp*SQRT( 2.*en(ji,jj,jk)/zrn2 ), mxln(ji,jj,jk) ) 784 784 END DO 785 785 END DO … … 933 933 !! 934 934 NAMELIST/namzdf_gls/rn_emin, rn_epsmin, ln_length_lim, & 935 & clim_galp, ln_crban, ln_sigpsi,&936 & rn_crban, rn_charn, &937 & nn_tkebc_surf, nn_tkebc_bot, &938 & nn_psibc_surf, nn_psibc_bot, &935 & rn_clim_galp, ln_crban, ln_sigpsi, & 936 & rn_crban, rn_charn, & 937 & nn_tkebc_surf, nn_tkebc_bot, & 938 & nn_psibc_surf, nn_psibc_bot, & 939 939 & nn_stab_func, nn_clos 940 940 !!---------------------------------------------------------- … … 958 958 WRITE(numout,*) ' Bottom roughness (m) rn_hbro = ', rn_hbro 959 959 WRITE(numout,*) ' Limit dissipation rate under stable stratification ln_length_lim = ',ln_length_lim 960 WRITE(numout,*) ' Galperin length scale limitation coef (Standard: 0.53, Holt: 0.26) clim_galp = ',clim_galp960 WRITE(numout,*) ' Galperin length scale limitation coef (Standard: 0.53, Holt: 0.26) rn_clim_galp = ', rn_clim_galp 961 961 WRITE(numout,*) ' TKE Surface boundary condition nn_tkebc_surf = ', nn_tkebc_surf 962 962 WRITE(numout,*) ' TKE Bottom boundary condition nn_tkebc_bot = ', nn_tkebc_bot … … 1206 1206 rn_c03 = rn_c0**3._wp 1207 1207 rn_c04 = rn_c0**4._wp 1208 rn_c03_sqrt2_galp = rn_c03 / SQRT(2._wp) / clim_galp1208 rn_c03_sqrt2_galp = rn_c03 / SQRT(2._wp) / rn_clim_galp 1209 1209 rn_sbc_mb = 0.5 * (15.8*rn_crban)**(2./3.) ! Surf. bound. cond. from Mellor and Blumberg 1210 1210 rn_sbc_std = 3.75 ! Surf. bound. cond. standard (prod=diss)
Note: See TracChangeset
for help on using the changeset viewer.