Ignore:
Timestamp:
2017-09-07T20:08:11+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part5 - reaching the end

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd.F90

    r8506 r8512  
    2626   USE sbc_oce , ONLY : sss_m, sst_m, e3t_m, utau, vtau, ssu_m, ssv_m, frq_m, qns_tot, qsr_tot, sprecip, ln_cpl 
    2727   USE sbc_ice , ONLY : qsr_oce, qns_oce, qemp_oce, qsr_ice, qns_ice, dqns_ice, evap_ice, qprec_ice, qevap_ice, & 
    28       &                 fr1_i0, fr2_i0, nn_limflx 
     28      &                 fr1_i0, fr2_i0 
    2929   USE ice1D          ! thermodynamic sea-ice variables 
    3030   USE icethd_dif     ! vertical diffusion 
     
    9292      IF( kt == nit000 .AND. lwp ) THEN 
    9393         WRITE(numout,*) 
    94          WRITE(numout,*)' icethd : sea-ice thermodynamics' 
    95          WRITE(numout,*)' ~~~~~~~~' 
     94         WRITE(numout,*) 'ice_thd: sea-ice thermodynamics' 
     95         WRITE(numout,*) '~~~~~~~' 
    9696      ENDIF 
    9797       
     
    104104      ! computation of friction velocity at T points 
    105105      !---------------------------------------------! 
    106       IF( ln_limdyn ) THEN 
     106      IF( ln_icedyn ) THEN 
    107107         zu_io(:,:) = u_ice(:,:) - ssu_m(:,:) 
    108108         zv_io(:,:) = v_ice(:,:) - ssv_m(:,:) 
     
    532532      !!              thermodynamics 
    533533      !! 
    534       !! ** Method  :   Read the namicethd namelist and check the ice-thermo 
     534      !! ** Method  :   Read the namice_thd namelist and check the ice-thermo 
    535535      !!              parameter values called at the first timestep (nit000) 
    536536      !! 
     
    539539      INTEGER  ::   ios   ! Local integer output status for namelist read 
    540540      !! 
    541       NAMELIST/namicethd/ rn_kappa_i, nn_ice_thcon, ln_dqnsice, rn_cdsn,                                  & 
     541      NAMELIST/namice_thd/ ln_limthd, rn_kappa_i, nn_ice_thcon, ln_dqnsice, rn_cdsn,                       & 
    542542         &                ln_limdH, rn_betas,                                                             & 
    543543         &                ln_limdA, rn_beta, rn_dmin,                                                     & 
     
    546546      !!------------------------------------------------------------------- 
    547547      ! 
    548       REWIND( numnam_ice_ref )              ! Namelist namicethd in reference namelist : Ice thermodynamics 
    549       READ  ( numnam_ice_ref, namicethd, IOSTAT = ios, ERR = 901) 
    550 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicethd in reference namelist', lwp ) 
    551  
    552       REWIND( numnam_ice_cfg )              ! Namelist namicethd in configuration namelist : Ice thermodynamics 
    553       READ  ( numnam_ice_cfg, namicethd, IOSTAT = ios, ERR = 902 ) 
    554 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicethd in configuration namelist', lwp ) 
    555       IF(lwm) WRITE ( numoni, namicethd ) 
     548      REWIND( numnam_ice_ref )              ! Namelist namice_thd in reference namelist : Ice thermodynamics 
     549      READ  ( numnam_ice_ref, namice_thd, IOSTAT = ios, ERR = 901) 
     550901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_thd in reference namelist', lwp ) 
     551 
     552      REWIND( numnam_ice_cfg )              ! Namelist namice_thd in configuration namelist : Ice thermodynamics 
     553      READ  ( numnam_ice_cfg, namice_thd, IOSTAT = ios, ERR = 902 ) 
     554902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_thd in configuration namelist', lwp ) 
     555      IF(lwm) WRITE ( numoni, namice_thd ) 
    556556      ! 
    557557      ! 
     
    559559         WRITE(numout,*) 'ice_thd_init : Ice Thermodynamics' 
    560560         WRITE(numout,*) '~~~~~~~~~~~~~' 
    561          WRITE(numout,*)'   Namelist namicethd' 
    562          WRITE(numout,*)'   -- icethd_dif --' 
    563          WRITE(numout,*)'      extinction radiation parameter in sea ice               rn_kappa_i   = ', rn_kappa_i 
    564          WRITE(numout,*)'      switch for comp. of thermal conductivity in the ice     nn_ice_thcon = ', nn_ice_thcon 
    565          WRITE(numout,*)'      change the surface non-solar flux with Tsu or not       ln_dqnsice   = ', ln_dqnsice 
    566          WRITE(numout,*)'      thermal conductivity of the snow                        rn_cdsn      = ', rn_cdsn 
    567          WRITE(numout,*)'   -- icethd_dh --' 
    568          WRITE(numout,*)'      activate ice thick change from top/bot (T) or not (F)   ln_limdH     = ', ln_limdH 
    569          WRITE(numout,*)'      coefficient for ice-lead partition of snowfall          rn_betas     = ', rn_betas 
    570          WRITE(numout,*)'   -- icethd_da --' 
    571          WRITE(numout,*)'      activate lateral melting (T) or not (F)                 ln_limdA     = ', ln_limdA 
    572          WRITE(numout,*)'      Coef. beta for lateral melting param.                   rn_beta      = ', rn_beta 
    573          WRITE(numout,*)'      Minimum floe diameter for lateral melting param.        rn_dmin      = ', rn_dmin 
    574          WRITE(numout,*)'   -- icethd_lac --' 
    575          WRITE(numout,*)'      activate ice growth in open-water (T) or not (F)        ln_limdO     = ', ln_limdO 
    576          WRITE(numout,*)'      ice thickness for lateral accretion                     rn_hnewice   = ', rn_hnewice 
    577          WRITE(numout,*)'      Frazil ice thickness as a function of wind or not       ln_frazil    = ', ln_frazil 
    578          WRITE(numout,*)'      Maximum proportion of frazil ice collecting at bottom   rn_maxfrazb  = ', rn_maxfrazb 
    579          WRITE(numout,*)'      Threshold relative drift speed for collection of frazil rn_vfrazb    = ', rn_vfrazb 
    580          WRITE(numout,*)'      Squeezing coefficient for collection of frazil          rn_Cfrazb    = ', rn_Cfrazb 
    581          WRITE(numout,*)'   -- iceitd --' 
    582          WRITE(numout,*)'      minimum ice thickness                                   rn_himin     = ', rn_himin  
    583          WRITE(numout,*)'   -- icestp --' 
    584          WRITE(numout,*)'      Multicategory heat flux formulation                     nn_limflx    = ', nn_limflx 
     561         WRITE(numout,*) '   Namelist namice_thd' 
     562         WRITE(numout,*) '      Ice thermodynamics (T) or not (F)                       ln_limthd    = ', ln_limthd 
     563         WRITE(numout,*) '   -- icethd_dif --' 
     564         WRITE(numout,*) '      extinction radiation parameter in sea ice               rn_kappa_i   = ', rn_kappa_i 
     565         WRITE(numout,*) '      switch for comp. of thermal conductivity in the ice     nn_ice_thcon = ', nn_ice_thcon 
     566         WRITE(numout,*) '      change the surface non-solar flux with Tsu or not       ln_dqnsice   = ', ln_dqnsice 
     567         WRITE(numout,*) '      thermal conductivity of the snow                        rn_cdsn      = ', rn_cdsn 
     568         WRITE(numout,*) '   -- icethd_dh --' 
     569         WRITE(numout,*) '      activate ice thick change from top/bot (T) or not (F)   ln_limdH     = ', ln_limdH 
     570         WRITE(numout,*) '      coefficient for ice-lead partition of snowfall          rn_betas     = ', rn_betas 
     571         WRITE(numout,*) '   -- icethd_da --' 
     572         WRITE(numout,*) '      activate lateral melting (T) or not (F)                 ln_limdA     = ', ln_limdA 
     573         WRITE(numout,*) '      Coef. beta for lateral melting param.                   rn_beta      = ', rn_beta 
     574         WRITE(numout,*) '      Minimum floe diameter for lateral melting param.        rn_dmin      = ', rn_dmin 
     575         WRITE(numout,*) '   -- icethd_lac --' 
     576         WRITE(numout,*) '      activate ice growth in open-water (T) or not (F)        ln_limdO     = ', ln_limdO 
     577         WRITE(numout,*) '      ice thickness for lateral accretion                     rn_hnewice   = ', rn_hnewice 
     578         WRITE(numout,*) '      Frazil ice thickness as a function of wind or not       ln_frazil    = ', ln_frazil 
     579         WRITE(numout,*) '      Maximum proportion of frazil ice collecting at bottom   rn_maxfrazb  = ', rn_maxfrazb 
     580         WRITE(numout,*) '      Threshold relative drift speed for collection of frazil rn_vfrazb    = ', rn_vfrazb 
     581         WRITE(numout,*) '      Squeezing coefficient for collection of frazil          rn_Cfrazb    = ', rn_Cfrazb 
     582         WRITE(numout,*) '   -- iceitd --' 
     583         WRITE(numout,*) '      minimum ice thickness                                   rn_himin     = ', rn_himin  
     584         WRITE(numout,*) '   -- icestp --' 
     585         WRITE(numout,*) '      Multicategory heat flux formulation                     nn_limflx    = ', nn_limflx 
    585586      ENDIF 
    586587      ! 
Note: See TracChangeset for help on using the changeset viewer.