Changeset 1449

Timestamp:

Nov 8, 2010, 5:19:04 PM (14 years ago)

Author:

jghattas

Message:

Bug correction : avoid small fractions to reappear at second time step.

/M-A Foujols, JG

File:

• LMDZ5/branches/LMDZ5V1.0-dev/libf/phylmd/change_srf_frac_mod.F90 (modified) (4 diffs)

LMDZ5/branches/LMDZ5V1.0-dev/libf/phylmd/change_srf_frac_mod.F90

@@ -9 +9 @@
 ! 
 ! Change Surface Fractions
-!
+! Author J Ghattas 2008
+
   SUBROUTINE change_srf_frac(itime, dtime, jour, &
        pctsrf, alb1, alb2, tsurf, u10m, v10m, pbl_tke)
@@ -76 +77 @@
     END SELECT
 
-    IF (is_modified) THEN
+
 !****************************************************************************************
 ! 2) 
@@ -84 +85 @@
 !
 !****************************************************************************************
+    IF (is_modified) THEN
   
 ! Test and exit if a fraction is negative
@@ -150 +152 @@
        CALL pbl_surface_newfrac(itime, pctsrf, pctsrf_old, tsurf, alb1, alb2, u10m, v10m, pbl_tke)
 
+    ELSE
+       ! No modifcation should be done
+       pctsrf(:,:) = pctsrf_old(:,:)
+
     END IF ! is_modified