Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90

Timestamp:

2012-01-28T17:44:18+01:00 (12 years ago)

Author:

rblod

Message:

Merge of 3.4beta into the trunk

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90

@@ -9 +9 @@
    !!              -   !  2006-08  (G. Madec)  style 
    !!             3.2  !  2006-08  (S. Masson, G. Madec)  suppress useless variables + style 
+   !!             3.4  !  2011-11  (C. Harris)  minor changes for CICE constants 
    !!----------------------------------------------------------------------
 
@@ -48 +49 @@
 #endif
 
+#if defined key_cice
+   REAL(wp), PUBLIC ::   rau0     = 1026._wp      !: reference volumic mass (density)  (kg/m3)
+#else
    REAL(wp), PUBLIC ::   rau0     = 1035._wp      !: reference volumic mass (density)  (kg/m3)
+#endif
    REAL(wp), PUBLIC ::   rau0r                    !: reference specific volume         (m3/kg)
    REAL(wp), PUBLIC ::   rcp      =    4.e+3_wp   !: ocean specific heat
    REAL(wp), PUBLIC ::   ro0cpr                   !: = 1. / ( rau0 * rcp )
 
-#if defined key_lim3
+#if defined key_lim3 || defined key_cice
    REAL(wp), PUBLIC ::   rcdsn   =   0.31_wp      !: thermal conductivity of snow
    REAL(wp), PUBLIC ::   rcdic   =   2.034396_wp  !: thermal conductivity of fresh ice
@@ -100 +105 @@
       rsiyea = 365.25 * rday * 2. * rpi / 6.283076
       rsiday = rday / ( 1. + rday / rsiyea )
+#if defined key_cice
+      omega =  7.292116e-05
+#else
       omega  = 2. * rpi / rsiday 
+#endif
 
       rau0r  = 1. /   rau0  
@@ -155 +164 @@
          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'
+#elif defined key_cice
+         WRITE(numout,*) '          latent heat of fusion of fresh ice / snow = ', lfus    , ' J/kg'
 #else
          WRITE(numout,*) '          density times specific heat for snow      = ', rcpsn   , ' J/m^3/K'