- Timestamp:
- 2014-06-04T19:18:57+02:00 (10 years ago)
- Location:
- branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM
- Files:
-
- 8 edited
-
CONFIG/ORCA2_LIM3/EXP00/namelist_ice (modified) (3 diffs)
-
CONFIG/ORCA2_LIM3/EXP00/namelist_ice_lim3 (modified) (3 diffs)
-
NEMO/LIM_SRC_3/ice.F90 (modified) (4 diffs)
-
NEMO/LIM_SRC_3/limdyn.F90 (modified) (2 diffs)
-
NEMO/LIM_SRC_3/limitd_me.F90 (modified) (7 diffs)
-
NEMO/LIM_SRC_3/limthd.F90 (modified) (3 diffs)
-
NEMO/LIM_SRC_3/limthd_lac.F90 (modified) (2 diffs)
-
NEMO/LIM_SRC_3/thd_ice.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_ice
r4654 r4659 44 44 !----------------------------------------------------------------------- 45 45 epsd = 1.0e-20 ! tolerance parameter 46 alpha = 0.5 ! coefficient for semi-implicit coriolis47 dm = 0.6e+03 ! diffusion constant for dynamics48 nbiter = 1 ! number of sub-time steps for relaxation49 nbitdr = 500 ! maximum number of iterations for relaxation50 46 om = 0.5 ! relaxation constant 51 resl = 5.0e-05 ! maximum value for the residual of relaxation52 47 cw = 5.0e-03 ! drag coefficient for oceanic stress 53 48 angvg = 0.0 ! turning angle for oceanic stress 54 49 pstar = 2.0e+04 ! 1st bulk-rheology parameter 55 50 c_rhg = 20.0 ! 2nd bulk-rhelogy parameter 56 etamn = 0.0e+07 ! minimun value for viscosity57 51 creepl = 1.0e-12 ! creep limit 58 52 ecc = 2.0 ! eccentricity of the elliptical yield curve … … 64 58 / 65 59 !----------------------------------------------------------------------- 66 &namicetrp ! ice transport67 !-----------------------------------------------------------------------68 bound = 0. ! boundary conditions (=0.0 no-slip, =1.0 free-slip)69 /70 !-----------------------------------------------------------------------71 60 &namicethd ! ice thermodynamic 72 61 !----------------------------------------------------------------------- 73 62 hmelt = -0.15 ! maximum melting at the bottom 74 hiccrit = 0.1 , 0.1 ! ice thickness for lateral accretion in the Northern (Southern) Hemisphere63 hiccrit = 0.1 ! ice thickness for lateral accretion 75 64 ! caution 1.0, 1.0 best value to be used!!! (gilles G.) ???? 76 fraz_swi = 0 .0! use of frazil ice collection thickness in function of wind (1.0) or not (0.0)65 fraz_swi = 0 ! use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 77 66 maxfrazb = 0.0 ! maximum portion of frazil ice collecting at the ice bottom 78 67 vfrazb = 0.4166667 ! thresold drift speed for frazil ice collecting at the ice bottom 79 68 Cfrazb = 5.0 ! squeezing coefficient for frazil ice collecting at the ice bottom 80 hicmin = 0.2 ! ice thickness corr. to max. energy stored in brine pocket81 69 hiclim = 0.10 ! minimum ice thickness 82 sbeta = 1. ! numerical caracteritic of the scheme for diffusion in ice83 ! Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)84 parlat = 0.0 ! percentage of energy used for lateral ablation85 hakspl = 0.25 ! slope of distr. for Hakkinen-Mellor's lateral melting86 hibspl = 0.5 ! slope of distribution for Hibler's lateral melting87 exld = 2.0 ! exponent for leads-closure rate88 hakdif = 1.0 ! coefficient for diffusions of ice and snow89 thth = 0.2 ! threshold thickness for comp. of eq. thermal conductivity90 70 hnzst = 0.1 ! thickness of the surf. layer in temp. computation 91 71 parsub = 1.0 ! switch for snow sublimation or not 92 alphs = 1.0 ! coefficient for snow density when snow ice formation93 72 betas = 0.6 ! exponent in lead-ice fractionation of snow precipitation 0.66 94 73 ! betas = 1 -> equipartition, betas < 1 -> more on leads … … 127 106 astar = 0.05 ! equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice) 128 107 Hstar = 100.0 ! parameter determining the maximum thickness of ridged ice 129 raft swi= 1 ! rafting or not108 raft_swi = 1 ! rafting or not 130 109 hparmeter = 0.75 ! threshold thickness for rafting or not 131 110 Craft = 5.0 ! coefficient used in the rafting function 132 111 ridge_por = 0.3 ! initial porosity of the ridged ice (typically 0.30) 133 sal_max_ridge = 15.0 ! maximum ridged ice salinity134 112 partfun_swi = 1 ! participation function linear, TH75 (0) or exponential Letal07 (1) 135 transfun_swi = 0 ! transfer function uniform of H80 (0) or exponential Letal07 (1)136 113 brinstren_swi = 0 ! (1) use brine volume to diminish ice strength 137 114 / -
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_ice_lim3
r4654 r4659 44 44 !----------------------------------------------------------------------- 45 45 epsd = 1.0e-20 ! tolerance parameter 46 alpha = 0.5 ! coefficient for semi-implicit coriolis47 dm = 0.6e+03 ! diffusion constant for dynamics48 nbiter = 1 ! number of sub-time steps for relaxation49 nbitdr = 500 ! maximum number of iterations for relaxation50 46 om = 0.5 ! relaxation constant 51 resl = 5.0e-05 ! maximum value for the residual of relaxation52 47 cw = 5.0e-03 ! drag coefficient for oceanic stress 53 48 angvg = 0.0 ! turning angle for oceanic stress 54 49 pstar = 2.0e+04 ! 1st bulk-rheology parameter 55 50 c_rhg = 20.0 ! 2nd bulk-rhelogy parameter 56 etamn = 0.0e+07 ! minimun value for viscosity57 51 creepl = 1.0e-12 ! creep limit 58 52 ecc = 2.0 ! eccentricity of the elliptical yield curve … … 64 58 / 65 59 !----------------------------------------------------------------------- 66 &namicetrp ! ice transport67 !-----------------------------------------------------------------------68 bound = 0. ! boundary conditions (=0.0 no-slip, =1.0 free-slip)69 /70 !-----------------------------------------------------------------------71 60 &namicethd ! ice thermodynamic 72 61 !----------------------------------------------------------------------- 73 62 hmelt = -0.15 ! maximum melting at the bottom 74 hiccrit = 0.1 , 0.1 ! ice thickness for lateral accretion in the Northern (Southern) Hemisphere63 hiccrit = 0.1 ! ice thickness for lateral accretion 75 64 ! caution 1.0, 1.0 best value to be used!!! (gilles G.) ???? 76 fraz_swi = 0 .0! use of frazil ice collection thickness in function of wind (1.0) or not (0.0)65 fraz_swi = 0 ! use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 77 66 maxfrazb = 0.0 ! maximum portion of frazil ice collecting at the ice bottom 78 67 vfrazb = 0.4166667 ! thresold drift speed for frazil ice collecting at the ice bottom 79 68 Cfrazb = 5.0 ! squeezing coefficient for frazil ice collecting at the ice bottom 80 hicmin = 0.2 ! ice thickness corr. to max. energy stored in brine pocket81 69 hiclim = 0.10 ! minimum ice thickness 82 sbeta = 1. ! numerical caracteritic of the scheme for diffusion in ice83 ! Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)84 parlat = 0.0 ! percentage of energy used for lateral ablation85 hakspl = 0.25 ! slope of distr. for Hakkinen-Mellor's lateral melting86 hibspl = 0.5 ! slope of distribution for Hibler's lateral melting87 exld = 2.0 ! exponent for leads-closure rate88 hakdif = 1.0 ! coefficient for diffusions of ice and snow89 thth = 0.2 ! threshold thickness for comp. of eq. thermal conductivity90 70 hnzst = 0.1 ! thickness of the surf. layer in temp. computation 91 71 parsub = 1.0 ! switch for snow sublimation or not 92 alphs = 1.0 ! coefficient for snow density when snow ice formation93 72 betas = 0.6 ! exponent in lead-ice fractionation of snow precipitation 0.66 94 73 ! betas = 1 -> equipartition, betas < 1 -> more on leads … … 127 106 astar = 0.05 ! equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice) 128 107 Hstar = 100.0 ! parameter determining the maximum thickness of ridged ice 129 raft swi= 1 ! rafting or not108 raft_swi = 1 ! rafting or not 130 109 hparmeter = 0.75 ! threshold thickness for rafting or not 131 110 Craft = 5.0 ! coefficient used in the rafting function 132 111 ridge_por = 0.3 ! initial porosity of the ridged ice (typically 0.30) 133 sal_max_ridge = 15.0 ! maximum ridged ice salinity134 112 partfun_swi = 1 ! participation function linear, TH75 (0) or exponential Letal07 (1) 135 transfun_swi = 0 ! transfer function uniform of H80 (0) or exponential Letal07 (1)136 113 brinstren_swi = 0 ! (1) use brine volume to diminish ice strength 137 114 / -
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/ice.F90
r4649 r4659 166 166 167 167 ! !!** ice-dynamic namelist (namicedyn) ** 168 INTEGER , PUBLIC :: nbiter = 1 !: number of sub-time steps for relaxation169 INTEGER , PUBLIC :: nbitdr = 250 !: maximum number of iterations for relaxation170 168 INTEGER , PUBLIC :: nevp = 400 !: number of iterations for subcycling 171 169 INTEGER , PUBLIC :: nlay_i = 5 !: number of layers in the ice … … 173 171 ! !!** ice-dynamic namelist (namicedyn) ** 174 172 REAL(wp), PUBLIC :: epsd = 1.0e-20_wp !: tolerance parameter for dynamic 175 REAL(wp), PUBLIC :: alpha = 0.5_wp !: coefficient for semi-implicit coriolis176 REAL(wp), PUBLIC :: dm = 0.6e+03_wp !: diffusion constant for dynamics177 173 REAL(wp), PUBLIC :: om = 0.5_wp !: relaxation constant 178 REAL(wp), PUBLIC :: resl = 5.0e-05_wp !: maximum value for the residual of relaxation179 174 REAL(wp), PUBLIC :: cw = 5.0e-03_wp !: drag coefficient for oceanic stress 180 175 REAL(wp), PUBLIC :: angvg = 0._wp !: turning angle for oceanic stress 181 176 REAL(wp), PUBLIC :: pstar = 1.0e+04_wp !: determines ice strength (N/M), Hibler JPO79 182 177 REAL(wp), PUBLIC :: c_rhg = 20._wp !: determines changes in ice strength 183 REAL(wp), PUBLIC :: etamn = 0.0e+07_wp !: minimun value for viscosity : has to be 0184 178 REAL(wp), PUBLIC :: creepl = 2.0e-08_wp !: creep limit : has to be under 1.0e-9 185 179 REAL(wp), PUBLIC :: ecc = 2._wp !: eccentricity of the elliptical yield curve … … 220 214 REAL(wp), PUBLIC :: Craft = 5.0_wp !: coefficient for smoothness of the hyperbolic tangent in rafting 221 215 REAL(wp), PUBLIC :: ridge_por = 0.0_wp !: initial porosity of ridges (0.3 regular value) 222 REAL(wp), PUBLIC :: sal_max_ridge = 15.0_wp !: maximum ridged ice salinity (ppt) 223 REAL(wp), PUBLIC :: betas = 1.0_wp !: coef. for partitioning of snowfall between leads and sea ice 224 REAL(wp), PUBLIC :: kappa_i = 1.0_wp !: coef. for the extinction of radiation Grenfell et al. (2006) [1/m] 216 REAL(wp), PUBLIC :: betas = 1.0_wp !: coef. for partitioning of snowfall between leads and sea ice 217 REAL(wp), PUBLIC :: kappa_i = 1.0_wp !: coef. for the extinction of radiation Grenfell et al. (2006) [1/m] 225 218 REAL(wp), PUBLIC :: nconv_i_thd = 50_wp !: maximal number of iterations for heat diffusion 226 219 REAL(wp), PUBLIC :: maxer_i_thd = 1.0e-4_wp !: maximal tolerated error (C) for heat diffusion … … 228 221 ! !!** ice-mechanical redistribution namelist (namiceitdme) 229 222 INTEGER , PUBLIC :: ridge_scheme_swi = 0 !: scheme used for ice ridging 230 INTEGER , PUBLIC :: raft swi= 1 !: rafting of ice or not223 INTEGER , PUBLIC :: raft_swi = 1 !: rafting of ice or not 231 224 INTEGER , PUBLIC :: partfun_swi = 1 !: participation function: =0 Thorndike et al. (1975), =1 Lipscomb et al. (2007) 232 INTEGER , PUBLIC :: transfun_swi = 0 !: transfer function: =0 Hibler 1980, =1 Lipscomb et al. 2007233 225 INTEGER , PUBLIC :: brinstren_swi = 0 !: use brine volume to diminish ice strength 234 226 -
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90
r4649 r4659 240 240 !! ** input : Namelist namicedyn 241 241 !!------------------------------------------------------------------- 242 NAMELIST/namicedyn/ epsd, alpha, & 243 & dm, nbiter, nbitdr, om, resl, cw, angvg, pstar, & 244 & c_rhg, etamn, creepl, ecc, ahi0, & 242 NAMELIST/namicedyn/ epsd, om, cw, angvg, pstar, & 243 & c_rhg, creepl, ecc, ahi0, & 245 244 & nevp, telast, alphaevp, hminrhg 246 245 !!------------------------------------------------------------------- … … 254 253 WRITE(numout,*) '~~~~~~~~~~~~' 255 254 WRITE(numout,*) ' tolerance parameter epsd = ', epsd 256 WRITE(numout,*) ' coefficient for semi-implicit coriolis alpha = ', alpha257 WRITE(numout,*) ' diffusion constant for dynamics dm = ', dm258 WRITE(numout,*) ' number of sub-time steps for relaxation nbiter = ', nbiter259 WRITE(numout,*) ' maximum number of iterations for relaxation nbitdr = ', nbitdr260 255 WRITE(numout,*) ' relaxation constant om = ', om 261 WRITE(numout,*) ' maximum value for the residual of relaxation resl = ', resl262 256 WRITE(numout,*) ' drag coefficient for oceanic stress cw = ', cw 263 257 WRITE(numout,*) ' turning angle for oceanic stress angvg = ', angvg 264 258 WRITE(numout,*) ' first bulk-rheology parameter pstar = ', pstar 265 259 WRITE(numout,*) ' second bulk-rhelogy parameter c_rhg = ', c_rhg 266 WRITE(numout,*) ' minimun value for viscosity etamn = ', etamn267 260 WRITE(numout,*) ' creep limit creepl = ', creepl 268 261 WRITE(numout,*) ' eccentricity of the elliptical yield curve ecc = ', ecc -
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90
r4649 r4659 618 618 !!---------------------------------------------------------------------! 619 619 INTEGER :: ji,jj, jl ! dummy loop indices 620 INTEGER :: krdg_index !621 620 REAL(wp) :: Gstari, astari, hi, hrmean, zdummy ! local scalar 622 621 REAL(wp), POINTER, DIMENSION(:,:) :: zworka ! temporary array used here … … 694 693 !----------------------------------------------------------------- 695 694 696 krdg_index = 1 697 698 IF( krdg_index == 0 ) THEN !--- Linear formulation (Thorndike et al., 1975) 695 IF( partfun_swi == 0 ) THEN !--- Linear formulation (Thorndike et al., 1975) 699 696 DO jl = 0, ice_cat_bounds(1,2) ! only undeformed ice participates 700 697 DO jj = 1, jpj … … 724 721 END DO 725 722 ! 726 ENDIF ! krdg_index727 728 IF( raft swi == 1 ) THEN ! Ridging and rafting ice participation functions723 ENDIF ! partfun_swi 724 725 IF( raft_swi == 1 ) THEN ! Ridging and rafting ice participation functions 729 726 ! 730 727 DO jl = 1, jpl … … 742 739 END DO ! jl 743 740 744 ELSE ! raft swi = 0741 ELSE ! raft_swi = 0 745 742 ! 746 743 DO jl = 1, jpl … … 750 747 ENDIF 751 748 752 IF ( raft swi == 1 ) THEN749 IF ( raft_swi == 1 ) THEN 753 750 754 751 IF( MAXVAL(aridge + araft - athorn(:,:,1:jpl)) .GT. epsi10 ) THEN … … 1348 1345 !! ** input : Namelist namiceitdme 1349 1346 !!------------------------------------------------------------------- 1350 NAMELIST/namiceitdme/ ridge_scheme_swi, Cs, Cf, fsnowrdg, fsnowrft,& 1351 Gstar, astar, & 1352 Hstar, raftswi, hparmeter, Craft, ridge_por, & 1353 sal_max_ridge, partfun_swi, transfun_swi, & 1354 brinstren_swi 1347 NAMELIST/namiceitdme/ ridge_scheme_swi, Cs, Cf, fsnowrdg, fsnowrft, & 1348 & Gstar, astar, Hstar, raft_swi, hparmeter, Craft, ridge_por, & 1349 & partfun_swi, brinstren_swi 1355 1350 !!------------------------------------------------------------------- 1356 1351 ! … … 1370 1365 WRITE(numout,*)' Equivalent to G* for an exponential part function astar ', astar 1371 1366 WRITE(numout,*)' Quantity playing a role in max ridged ice thickness Hstar ', Hstar 1372 WRITE(numout,*)' Rafting of ice sheets or not raft swi ', raftswi1367 WRITE(numout,*)' Rafting of ice sheets or not raft_swi ', raft_swi 1373 1368 WRITE(numout,*)' Parmeter thickness (threshold between ridge-raft) hparmeter ', hparmeter 1374 1369 WRITE(numout,*)' Rafting hyperbolic tangent coefficient Craft ', Craft 1375 1370 WRITE(numout,*)' Initial porosity of ridges ridge_por ', ridge_por 1376 WRITE(numout,*)' Maximum salinity of ridging ice sal_max_ridge ', sal_max_ridge1377 1371 WRITE(numout,*)' Switch for part. function (0) linear (1) exponential partfun_swi ', partfun_swi 1378 WRITE(numout,*)' Switch for tran. function (0) linear (1) exponential transfun_swi ', transfun_swi1379 1372 WRITE(numout,*)' Switch for including brine volume in ice strength comp. brinstren_swi ', brinstren_swi 1380 1373 ENDIF -
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r4649 r4659 610 610 !!------------------------------------------------------------------- 611 611 NAMELIST/namicethd/ hmelt , hiccrit, fraz_swi, maxfrazb, vfrazb, Cfrazb, & 612 & hicmin, hiclim, & 613 & sbeta , parlat, hakspl, hibspl, exld, & 614 & hakdif, hnzst , thth , parsub, alphs, betas, & 612 & hiclim, hnzst, parsub, betas, & 615 613 & kappa_i, nconv_i_thd, maxer_i_thd, thcon_i_swi 616 614 !!------------------------------------------------------------------- … … 629 627 WRITE(numout,*)' Namelist of ice parameters for ice thermodynamic computation ' 630 628 WRITE(numout,*)' maximum melting at the bottom hmelt = ', hmelt 631 WRITE(numout,*)' ice thick. for lateral accretion in NH (SH) hiccrit(1/2)= ', hiccrit629 WRITE(numout,*)' ice thick. for lateral accretion hiccrit = ', hiccrit 632 630 WRITE(numout,*)' Frazil ice thickness as a function of wind or not fraz_swi = ', fraz_swi 633 631 WRITE(numout,*)' Maximum proportion of frazil ice collecting at bottom maxfrazb = ', maxfrazb 634 632 WRITE(numout,*)' Thresold relative drift speed for collection of frazil vfrazb = ', vfrazb 635 633 WRITE(numout,*)' Squeezing coefficient for collection of frazil Cfrazb = ', Cfrazb 636 WRITE(numout,*)' ice thick. corr. to max. energy stored in brine pocket hicmin = ', hicmin637 634 WRITE(numout,*)' minimum ice thickness hiclim = ', hiclim 638 635 WRITE(numout,*)' numerical carac. of the scheme for diffusion in ice ' 639 WRITE(numout,*)' Cranck-Nicholson (=0.5), implicit (=1), explicit (=0) sbeta = ', sbeta640 WRITE(numout,*)' percentage of energy used for lateral ablation parlat = ', parlat641 WRITE(numout,*)' slope of distr. for Hakkinen-Mellor lateral melting hakspl = ', hakspl642 WRITE(numout,*)' slope of distribution for Hibler lateral melting hibspl = ', hibspl643 WRITE(numout,*)' exponent for leads-closure rate exld = ', exld644 WRITE(numout,*)' coefficient for diffusions of ice and snow hakdif = ', hakdif645 WRITE(numout,*)' threshold thick. for comp. of eq. thermal conductivity zhth = ', thth646 636 WRITE(numout,*)' thickness of the surf. layer in temp. computation hnzst = ', hnzst 647 637 WRITE(numout,*)' switch for snow sublimation (=1) or not (=0) parsub = ', parsub 648 WRITE(numout,*)' coefficient for snow density when snow ice formation alphs = ', alphs649 638 WRITE(numout,*)' coefficient for ice-lead partition of snowfall betas = ', betas 650 639 WRITE(numout,*)' extinction radiation parameter in sea ice (1.0) kappa_i = ', kappa_i … … 655 644 ENDIF 656 645 ! 657 rcdsn = hakdif * rcdsn658 rcdic = hakdif * rcdic659 !660 646 END SUBROUTINE lim_thd_init 661 647 -
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_lac.F90
r4649 r4659 157 157 158 158 ! Default new ice thickness 159 hicol(:,:) = hiccrit (1)160 161 IF( fraz_swi == 1 ._wp) THEN159 hicol(:,:) = hiccrit 160 161 IF( fraz_swi == 1 ) THEN 162 162 163 163 !-------------------- … … 322 322 !---------------------- 323 323 DO ji = 1, nbpac 324 zh_newice(ji) = hiccrit (1)325 END DO 326 IF( fraz_swi == 1 .0) zh_newice(:) = hicol_b(:)324 zh_newice(ji) = hiccrit 325 END DO 326 IF( fraz_swi == 1 ) zh_newice(:) = hicol_b(:) 327 327 328 328 !---------------------- -
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/thd_ice.F90
r4649 r4659 20 20 ! !!! ** ice-thermo namelist (namicethd) ** 21 21 REAL(wp), PUBLIC :: hmelt = -0.15 !: maximum melting at the bottom; active only for one category 22 REAL(wp), PUBLIC :: hicmin = 0.2 !: (REMOVE)23 22 REAL(wp), PUBLIC :: hiclim = 0.05 !: minimum ice thickness 24 REAL(wp), PUBLIC :: sbeta = 1.0 !: numerical scheme for diffusion in ice (REMOVE)25 REAL(wp), PUBLIC :: parlat = 0.0 !: (REMOVE)26 REAL(wp), PUBLIC :: hakspl = 0.5 !: (REMOVE)27 REAL(wp), PUBLIC :: hibspl = 0.5 !: (REMOVE)28 REAL(wp), PUBLIC :: exld = 2.0 !: (REMOVE)29 REAL(wp), PUBLIC :: hakdif = 1.0 !: (REMOVE)30 REAL(wp), PUBLIC :: thth = 0.2 !: (REMOVE)31 23 REAL(wp), PUBLIC :: hnzst = 0.1 !: thick. of the surf. layer in temp. comp. 32 24 REAL(wp), PUBLIC :: parsub = 1.0 !: switch for snow sublimation or not 33 REAL(wp), PUBLIC :: alphs = 1.0 !: coef. for snow density when snow-ice formation34 REAL(wp), PUBLIC :: fraz_swi= 1.0 !: use of frazil ice collection in function of wind (1.0) or not (0.0)35 25 REAL(wp), PUBLIC :: maxfrazb= 0.7 !: maximum portion of frazil ice collecting at the ice bottom 36 26 REAL(wp), PUBLIC :: vfrazb = 0.41667 !: threshold drift speed for collection of bottom frazil ice 37 27 REAL(wp), PUBLIC :: Cfrazb = 5.0 !: squeezing coefficient for collection of bottom frazil ice 28 REAL(wp), PUBLIC :: hiccrit = 0.3 !: ice th. for lateral accretion in the NH (SH) (m) 38 29 39 REAL(wp), PUBLIC, DIMENSION(2) :: hiccrit = (/0.3,0.3/) !: ice th. for lateral accretion in the NH (SH) (m)30 INTEGER , PUBLIC :: fraz_swi= 1 !: use of frazil ice collection in function of wind (1) or not (0) 40 31 41 32 !!-----------------------------
Note: See TracChangeset
for help on using the changeset viewer.
