Changeset 833 for trunk/NEMO/OPA_SRC/phycst.F90

Timestamp:

2008-03-07T14:51:35+01:00 (16 years ago)

Author:

rblod

Message:

Merge branche dev_002_LIM back to trunk ticket #70 and #71

File:

• trunk/NEMO/OPA_SRC/phycst.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/phycst.F90

@@ -46 +46 @@
       rtt      = 273.16_wp  ,  &  !: triple point of temperature (Kelvin)
       rt0      = 273.15_wp  ,  &  !: freezing point of water (Kelvin)
+#if defined key_lim3
+      rt0_snow = 273.16_wp  ,  &  !: melting point of snow  (Kelvin)
+      rt0_ice  = 273.16_wp  ,  &  !: melting point of ice   (Kelvin)
+#else
       rt0_snow = 273.15_wp  ,  &  !: melting point of snow  (Kelvin)
       rt0_ice  = 273.05_wp  ,  &  !: melting point of ice   (Kelvin)
+#endif
       rau0     = 1020._wp   ,  &  !: volumic mass of reference (kg/m3)
       rauw     = 1000._wp   ,  &  !: density of pure water (kg/m3)
@@ -54 +59 @@
 
    REAL(wp), PUBLIC ::            &  !:
+#if defined key_lim3
+      rcdsn   =   0.31_wp     ,   &  !: thermal conductivity of snow
+      rcdic   =   2.034396_wp ,   &  !: thermal conductivity of fresh ice
+      cpic    = 2067.0        ,   &
+      ! add the following lines
+      lsub    = 2.834e+6      ,   &  !: pure ice latent heat of sublimation (J.kg-1)
+      lfus    = 0.334e+6      ,   &  !: latent heat of fusion of fresh ice   (J.kg-1)
+      rhoic   = 917._wp       ,   &  !: density of sea ice (kg/m3)
+      tmut    =   0.054       ,   &  !: decrease of seawater meltpoint with salinity
+#else
       rcdsn   =   0.22_wp     ,   &  !: conductivity of the snow
       rcdic   =   2.034396_wp ,   &  !: conductivity of the ice
@@ -62 +77 @@
       xsn     =   2.8e+6      ,   &  !: latent heat of sublimation of snow
       rhoic   = 900._wp       ,   &  !: density of sea ice (kg/m3)
+#endif
       rhosn   = 330._wp       ,   &  !: density of snow (kg/m3)
       emic    =   0.97_wp     ,   &  !: emissivity of snow or ice
@@ -169 +185 @@
          WRITE(numout,*) '          thermal conductivity of the snow          = ', rcdsn   , ' J/s/m/K'
          WRITE(numout,*) '          thermal conductivity of the ice           = ', rcdic   , ' J/s/m/K'
+#if defined key_lim3
+         WRITE(numout,*) '          fresh ice specific heat                   = ', cpic    , ' J/kg/K'
+         WRITE(numout,*) '          latent heat of fusion of fresh ice / snow = ', lfus    , ' J/kg'
+         WRITE(numout,*) '          latent heat of subl.  of fresh ice / snow = ', lsub    , ' J/kg'
+#else
          WRITE(numout,*) '          density times specific heat for snow      = ', rcpsn   , ' J/m^3/K' 
          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' 
          WRITE(numout,*) '          latent heat of sublimation of snow        = ', xsn     , ' J/kg' 
+#endif
          WRITE(numout,*) '          density of sea ice                        = ', rhoic   , ' kg/m^3'
          WRITE(numout,*) '          density of snow                           = ', rhosn   , ' kg/m^3'