Changeset 833 for trunk/NEMO/OPA_SRC/phycst.F90
- Timestamp:
- 2008-03-07T14:51:35+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/phycst.F90
r719 r833 46 46 rtt = 273.16_wp , & !: triple point of temperature (Kelvin) 47 47 rt0 = 273.15_wp , & !: freezing point of water (Kelvin) 48 #if defined key_lim3 49 rt0_snow = 273.16_wp , & !: melting point of snow (Kelvin) 50 rt0_ice = 273.16_wp , & !: melting point of ice (Kelvin) 51 #else 48 52 rt0_snow = 273.15_wp , & !: melting point of snow (Kelvin) 49 53 rt0_ice = 273.05_wp , & !: melting point of ice (Kelvin) 54 #endif 50 55 rau0 = 1020._wp , & !: volumic mass of reference (kg/m3) 51 56 rauw = 1000._wp , & !: density of pure water (kg/m3) … … 54 59 55 60 REAL(wp), PUBLIC :: & !: 61 #if defined key_lim3 62 rcdsn = 0.31_wp , & !: thermal conductivity of snow 63 rcdic = 2.034396_wp , & !: thermal conductivity of fresh ice 64 cpic = 2067.0 , & 65 ! add the following lines 66 lsub = 2.834e+6 , & !: pure ice latent heat of sublimation (J.kg-1) 67 lfus = 0.334e+6 , & !: latent heat of fusion of fresh ice (J.kg-1) 68 rhoic = 917._wp , & !: density of sea ice (kg/m3) 69 tmut = 0.054 , & !: decrease of seawater meltpoint with salinity 70 #else 56 71 rcdsn = 0.22_wp , & !: conductivity of the snow 57 72 rcdic = 2.034396_wp , & !: conductivity of the ice … … 62 77 xsn = 2.8e+6 , & !: latent heat of sublimation of snow 63 78 rhoic = 900._wp , & !: density of sea ice (kg/m3) 79 #endif 64 80 rhosn = 330._wp , & !: density of snow (kg/m3) 65 81 emic = 0.97_wp , & !: emissivity of snow or ice … … 169 185 WRITE(numout,*) ' thermal conductivity of the snow = ', rcdsn , ' J/s/m/K' 170 186 WRITE(numout,*) ' thermal conductivity of the ice = ', rcdic , ' J/s/m/K' 187 #if defined key_lim3 188 WRITE(numout,*) ' fresh ice specific heat = ', cpic , ' J/kg/K' 189 WRITE(numout,*) ' latent heat of fusion of fresh ice / snow = ', lfus , ' J/kg' 190 WRITE(numout,*) ' latent heat of subl. of fresh ice / snow = ', lsub , ' J/kg' 191 #else 171 192 WRITE(numout,*) ' density times specific heat for snow = ', rcpsn , ' J/m^3/K' 172 193 WRITE(numout,*) ' density times specific heat for ice = ', rcpic , ' J/m^3/K' … … 174 195 WRITE(numout,*) ' volumetric latent heat fusion of snow = ', xlsn , ' J/m' 175 196 WRITE(numout,*) ' latent heat of sublimation of snow = ', xsn , ' J/kg' 197 #endif 176 198 WRITE(numout,*) ' density of sea ice = ', rhoic , ' kg/m^3' 177 199 WRITE(numout,*) ' density of snow = ', rhosn , ' kg/m^3'
Note: See TracChangeset
for help on using the changeset viewer.