Changeset 1450 for LMDZ4/branches/LMDZ4_AR5/libf/phylmd/change_srf_frac_mod.F90

Timestamp:

Nov 9, 2010, 11:10:11 AM (14 years ago)

Author:

jghattas

Message:

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

grid_noro : Deleted dimesion using local parameters. Instead using dimesions already existing as input arguments. This modification avoids need for editing the code before running a higer resolution (above 300x200).

/M-A Foujols, JG

File:

• LMDZ4/branches/LMDZ4_AR5/libf/phylmd/change_srf_frac_mod.F90 (modified) (4 diffs)

LMDZ4/branches/LMDZ4_AR5/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