include file: C_PERMA 4
      REAL HCAPW,HCAPI,RHO_ICE,DPSIDT                                      C_PERMA.2      
      PARAMETER(                                                           C_PERMA.3      
     + HCAPW=4180      ! Specific heat capacity of water (J/kg/K)          C_PERMA.4      
     +,HCAPI=2100      ! Specific heat capacity of ice (J/kg/K)            C_PERMA.5      
     +,RHO_ICE=917     ! Density of ice (kg/m3)                            C_PERMA.6      
     +,DPSIDT=114.3    ! Rate of change of ice potential with temperatur   C_PERMA.7      
C                      ! RHO_ICE*LF/ZERODEGC*1/(RHO_WATER*G) (m/K)         C_PERMA.8      
     +)                                                                    C_PERMA.9      
C-----------------------------------------------------------------------   C_PERMA.10