Changeset 1860 for trunk/NEMO/OPA_SRC/DOM
- Timestamp:
- 2010-05-08T11:23:48+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DOM/phycst.F90
r1739 r1860 4 4 !! Definition of of both ocean and ice parameters used in the code 5 5 !!===================================================================== 6 !! History : ! 90-10 (C. Levy - G. Madec) Original code 7 !! ! 91-11 (G. Madec) 8 !! ! 91-12 (M. Imbard) 9 !! 8.5 ! 02-08 (G. Madec, C. Ethe) F90, add ice constants 10 !! 9.0 ! 06-08 (G. Madec) style 6 !! History : OPA ! 1990-10 (C. Levy, G. Madec) Original code 7 !! ! 1991-11 (G. Madec, M. Imbard) 8 !! NEMO 1.0 ! 2002-08 (G. Madec, C. Ethe) F90, add ice constants 9 !! - ! 2006-08 (G. Madec) style 11 10 !!---------------------------------------------------------------------- 12 11 … … 59 58 rcdic = 2.034396_wp , & !: thermal conductivity of fresh ice 60 59 cpic = 2067.0 , & 61 ! add the following lines62 60 lsub = 2.834e+6 , & !: pure ice latent heat of sublimation (J.kg-1) 63 61 lfus = 0.334e+6 , & !: latent heat of fusion of fresh ice (J.kg-1) … … 69 67 rcpsn = 6.9069e+5_wp, & !: density times specific heat for snow 70 68 rcpic = 1.8837e+6_wp, & !: volumetric latent heat fusion of sea ice 71 xlsn = 110.121e+6_wp , & !: volumetric latent heat fusion of snow 72 xlic = 300.33e+6_wp , & !: volumetric latent heat fusion of ice 69 lfus = 0.3337e+6 , & !: latent heat of fusion of fresh ice (J.kg-1) 70 xlsn = 110.121e+6_wp , & !: = lfus * rhosn, volumetric latent heat fusion of snow 71 xlic = 300.33e+6_wp , & !: = lfus * rhosn, volumetric latent heat fusion of ice 73 72 xsn = 2.8e+6 , & !: latent heat of sublimation of snow 74 73 rhoic = 900._wp , & !: volumic mass of sea ice (kg/m3) … … 83 82 stefan = 5.67e-8_wp !: Stefan-Boltzmann constant 84 83 !!---------------------------------------------------------------------- 85 !! OPA 9.0 , LOCEAN-IPSL (2005)84 !! NEMO/OPA 3.2 , LOCEAN-IPSL (2010) 86 85 !! $Id$ 87 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt86 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 88 87 !!---------------------------------------------------------------------- 89 88 … … 160 159 IF(lwp) WRITE(numout,*) ' 1. / ( rau0 * rcp ) = ro0cpr = ', ro0cpr 161 160 161 #if defined key_lim3 162 xlsn = lfus * rhosn ! volumetric latent heat fusion of snow [J/m3] 163 #else 164 lfus = xlsn / rhosn ! latent heat of fusion of fresh ice 165 #endif 166 162 167 IF(lwp) THEN 163 168 WRITE(numout,*)
Note: See TracChangeset
for help on using the changeset viewer.