Changeset 3402 for LMDZ6/trunk/libf

Timestamp:

Oct 16, 2018, 2:35:21 PM (6 years ago)

Author:

lguez

Message:

Set temperature of appearing sea ice to -1.8 Celsius degrees instead
of 0 Celsius degrees, following advice of F. Codron.

File:

• LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -3227 +3227 @@
                 ! F. Codron sensible default values for ocean and sea ice
                 IF (nsrf.EQ.is_oce) THEN
-                    tsurf(i,nsrf) = 271.35 ! Freezing sea water
-                    DO k=1,nsw
+                   tsurf(i,nsrf) = 271.35
+                   ! (temperature of sea water under sea ice, so that
+                   ! is also the temperature of appearing sea water)
+                   DO k=1,nsw
                       alb_dir(i,k,nsrf) = 0.06 ! typical Ocean albedo
                       alb_dif(i,k,nsrf) = 0.06
-                    ENDDO
+                   ENDDO
                 ELSE IF (nsrf.EQ.is_sic) THEN
-                    tsurf(i,nsrf) = 273.15 ! Melting ice 
-                    DO k=1,nsw
+                   tsurf(i,nsrf) = 271.35
+                   ! (Temperature at base of sea ice. Surface
+                   ! temperature could be higher, up to 0 Celsius
+                   ! degrees. We set it to -1.8 Celsius degrees for
+                   ! consistency with the ocean slab model.)
+                   DO k=1,nsw
                       alb_dir(i,k,nsrf) = 0.3 ! thin ice
                       alb_dif(i,k,nsrf) = 0.3
-                    ENDDO
+                   ENDDO
                 ENDIF
-                ! F. Codron
              ELSE
                 ! The continents have changed. The new fraction receives the mean sum of the existent fractions