- Timestamp:
- 2015-02-09T14:39:07+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r5067 r5070 631 631 rswitch = 1._wp - MAX( 0._wp , SIGN( 1._wp , - a_i_1d(ji) + epsi20 ) ) 632 632 ht_i_1d(ji) = rswitch * zv / MAX( a_i_1d(ji), epsi20 ) 633 ! adjust thickness633 ! retrieve total concentration 634 634 at_i_1d(ji) = a_i_1d(ji) 635 635 END IF … … 651 651 !!------------------------------------------------------------------- 652 652 INTEGER :: ios ! Local integer output status for namelist read 653 NAMELIST/namicethd/ rn_hnewice, nn_frazil, rn_maxfrazb, rn_vfrazb, rn_Cfrazb, & 654 & rn_himin, parsub, rn_betas, & 655 & rn_kappa_i, nn_conv_dif, rn_terr_dif, nn_ice_thcon, & 653 NAMELIST/namicethd/ rn_hnewice, ln_frazil, rn_maxfrazb, rn_vfrazb, rn_Cfrazb, & 654 & rn_himin, parsub, rn_betas, rn_kappa_i, nn_conv_dif, rn_terr_dif, nn_ice_thcon, & 656 655 & nn_monocat 657 656 !!------------------------------------------------------------------- … … 682 681 WRITE(numout,*) 683 682 WRITE(numout,*)' Namelist of ice parameters for ice thermodynamic computation ' 684 WRITE(numout,*)' ice thick. for lateral accretion rn_hnewice 685 WRITE(numout,*)' Frazil ice thickness as a function of wind or not nn_frazil = ', nn_frazil686 WRITE(numout,*)' Maximum proportion of frazil ice collecting at bottom rn_maxfrazb 687 WRITE(numout,*)' Thresold relative drift speed for collection of frazil rn_vfrazb 688 WRITE(numout,*)' Squeezing coefficient for collection of frazil rn_Cfrazb 689 WRITE(numout,*)' minimum ice thickness rn_himin 683 WRITE(numout,*)' ice thick. for lateral accretion rn_hnewice = ', rn_hnewice 684 WRITE(numout,*)' Frazil ice thickness as a function of wind or not ln_frazil = ', ln_frazil 685 WRITE(numout,*)' Maximum proportion of frazil ice collecting at bottom rn_maxfrazb = ', rn_maxfrazb 686 WRITE(numout,*)' Thresold relative drift speed for collection of frazil rn_vfrazb = ', rn_vfrazb 687 WRITE(numout,*)' Squeezing coefficient for collection of frazil rn_Cfrazb = ', rn_Cfrazb 688 WRITE(numout,*)' minimum ice thickness rn_himin = ', rn_himin 690 689 WRITE(numout,*)' numerical carac. of the scheme for diffusion in ice ' 691 690 WRITE(numout,*)' switch for snow sublimation (=1) or not (=0) parsub = ', parsub 692 WRITE(numout,*)' coefficient for ice-lead partition of snowfall rn_betas 693 WRITE(numout,*)' extinction radiation parameter in sea ice rn_kappa_i 691 WRITE(numout,*)' coefficient for ice-lead partition of snowfall rn_betas = ', rn_betas 692 WRITE(numout,*)' extinction radiation parameter in sea ice rn_kappa_i = ', rn_kappa_i 694 693 WRITE(numout,*)' maximal n. of iter. for heat diffusion computation nn_conv_dif = ', nn_conv_dif 695 694 WRITE(numout,*)' maximal err. on T for heat diffusion computation rn_terr_dif = ', rn_terr_dif 696 WRITE(numout,*)' switch for comp. of thermal conductivity in the ice nn_ice_thcon 695 WRITE(numout,*)' switch for comp. of thermal conductivity in the ice nn_ice_thcon = ', nn_ice_thcon 697 696 WRITE(numout,*)' check heat conservation in the ice/snow con_i = ', con_i 698 697 WRITE(numout,*)' virtual ITD mono-category parameterizations (1) or not nn_monocat = ', nn_monocat
Note: See TracChangeset
for help on using the changeset viewer.