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

changes in style - part6 - more clarity (still not finished)

File:
1 edited

Legend:

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

    r8498 r8531  
    1717   USE ice1D          ! sea-ice: thermodynamic 1D variables 
    1818   ! 
     19   USE in_out_manager ! I/O manager 
    1920   USE lib_mpp        ! MPP library 
    2021   USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)   
     
    2324   PRIVATE 
    2425 
    25    PUBLIC   ice_thd_da   ! called by icethd.F90 
     26   PUBLIC   ice_thd_da        ! called by icethd.F90 
     27   PUBLIC   ice_thd_da_init   ! called by icestp.F90 
     28 
     29   ! ** namelist (namthd_da) ** 
     30   REAL(wp) ::   rn_beta     ! coef. beta for lateral melting param. 
     31   REAL(wp) ::   rn_dmin     ! minimum floe diameter for lateral melting param. 
    2632 
    2733   !!---------------------------------------------------------------------- 
     
    149155      ! 
    150156   END SUBROUTINE ice_thd_da 
    151     
     157 
     158   SUBROUTINE ice_thd_da_init 
     159      !!----------------------------------------------------------------------- 
     160      !!                   ***  ROUTINE ice_thd_da_init ***  
     161      !!                  
     162      !! ** Purpose :   Physical constants and parameters associated with 
     163      !!                ice thermodynamics 
     164      !! 
     165      !! ** Method  :   Read the namthd_da namelist and check the parameters 
     166      !!                called at the first timestep (nit000) 
     167      !! 
     168      !! ** input   :   Namelist namthd_da 
     169      !!------------------------------------------------------------------- 
     170      INTEGER  ::   ios   ! Local integer output status for namelist read 
     171      !! 
     172      NAMELIST/namthd_da/ rn_beta, rn_dmin 
     173      !!------------------------------------------------------------------- 
     174      ! 
     175      REWIND( numnam_ice_ref )              ! Namelist namthd_da in reference namelist : Ice thermodynamics 
     176      READ  ( numnam_ice_ref, namthd_da, IOSTAT = ios, ERR = 901) 
     177901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namthd_da in reference namelist', lwp ) 
     178 
     179      REWIND( numnam_ice_cfg )              ! Namelist namthd_da in configuration namelist : Ice thermodynamics 
     180      READ  ( numnam_ice_cfg, namthd_da, IOSTAT = ios, ERR = 902 ) 
     181902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namthd_da in configuration namelist', lwp ) 
     182      IF(lwm) WRITE ( numoni, namthd_da ) 
     183      ! 
     184      ! 
     185      IF(lwp) THEN                          ! control print 
     186         WRITE(numout,*) 'ice_thd_da_init: Ice lateral melting' 
     187         WRITE(numout,*) '~~~~~~~~~~~~~~~' 
     188         WRITE(numout,*) '   Namelist namthd_da:' 
     189         WRITE(numout,*) '      Coef. beta for lateral melting param.               rn_beta = ', rn_beta 
     190         WRITE(numout,*) '      Minimum floe diameter for lateral melting param.    rn_dmin = ', rn_dmin 
     191      ENDIF 
     192      ! 
     193   END SUBROUTINE ice_thd_da_init 
     194   
    152195#else 
    153196   !!---------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.