Changeset 1450 for LMDZ4

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

Location:

LMDZ4/branches/LMDZ4_AR5/libf

Files:

• dyn3d/grid_noro.F (modified) (1 diff)
• dyn3dpar/grid_noro.F (modified) (1 diff)
• phylmd/change_srf_frac_mod.F90 (modified) (4 diffs)

LMDZ4/branches/LMDZ4_AR5/libf/dyn3d/grid_noro.F

@@ -458 +458 @@
 C MAKE A MOVING AVERAGE OVER 9 GRIDPOINTS OF THE X FIELDS
 
-      PARAMETER (ISMo=400,JSMo=200)
-      REAL X(IMAR,JMAR),XF(ISMo,JSMo)
+      REAL X(IMAR,JMAR),XF(IMAR,JMAR)
       real WEIGHTpb(-1:1,-1:1)
 
-      if(imar.gt.ismo) stop'surdimensionner ismo dans mva9 (grid_noro)'
-      if(jmar.gt.jsmo) stop'surdimensionner jsmo dans mva9 (grid_noro)'
-      
       SUM=0.
       DO IS=-1,1

LMDZ4/branches/LMDZ4_AR5/libf/dyn3dpar/grid_noro.F

@@ -458 +458 @@
 C MAKE A MOVING AVERAGE OVER 9 GRIDPOINTS OF THE X FIELDS
 
-      PARAMETER (ISMo=300,JSMo=200)
-      REAL X(IMAR,JMAR),XF(ISMo,JSMo)
+      REAL X(IMAR,JMAR),XF(IMAR,JMAR)
       real WEIGHTpb(-1:1,-1:1)
 
-      if(imar.gt.ismo) stop'surdimensionner ismo dans mva9 (grid_noro)'
-      if(jmar.gt.jsmo) stop'surdimensionner jsmo dans mva9 (grid_noro)'
-      
       SUM=0.
       DO IS=-1,1

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