Changeset 2655 for LMDZ5/trunk/libf

Timestamp:

Oct 10, 2016, 9:47:19 AM (8 years ago)

Author:

Ehouarn Millour

Message:

Correction in pbl_surface_mod to better account and set default values for ocean and sea ice when such new fractions are created.
NB: This correction changes test-bench results.
FC

File:

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

LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -3095 +3095 @@
                 endif
                 mfois(nsrf) = mfois(nsrf) + 1
+                ! 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
+                      alb_dir(i,k,nsrf) = 0.06 ! typical Ocean albedo
+                      alb_dif(i,k,nsrf) = 0.06
+                    END DO
+                ELSE IF (nsrf.EQ.is_sic) THEN
+                    tsurf(i,nsrf) = 273.15 ! Melting ice 
+                    DO k=1,nsw
+                      alb_dir(i,k,nsrf) = 0.3 ! thin ice
+                      alb_dif(i,k,nsrf) = 0.3
+                    END DO
+                END IF
+                ! F. Codron
              ELSE
                 ! The continents have changed. The new fraction receives the mean sum of the existent fractions