Changeset 1450


Ignore:
Timestamp:
11/09/10 11:10:11 (7 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:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4_AR5/libf/dyn3d/grid_noro.F

    r773 r1450  
    458458C MAKE A MOVING AVERAGE OVER 9 GRIDPOINTS OF THE X FIELDS 
    459459 
    460       PARAMETER (ISMo=400,JSMo=200) 
    461       REAL X(IMAR,JMAR),XF(ISMo,JSMo) 
     460      REAL X(IMAR,JMAR),XF(IMAR,JMAR) 
    462461      real WEIGHTpb(-1:1,-1:1) 
    463462 
    464       if(imar.gt.ismo) stop'surdimensionner ismo dans mva9 (grid_noro)' 
    465       if(jmar.gt.jsmo) stop'surdimensionner jsmo dans mva9 (grid_noro)' 
    466        
    467463      SUM=0. 
    468464      DO IS=-1,1 
  • LMDZ4/branches/LMDZ4_AR5/libf/dyn3dpar/grid_noro.F

    r764 r1450  
    458458C MAKE A MOVING AVERAGE OVER 9 GRIDPOINTS OF THE X FIELDS 
    459459 
    460       PARAMETER (ISMo=300,JSMo=200) 
    461       REAL X(IMAR,JMAR),XF(ISMo,JSMo) 
     460      REAL X(IMAR,JMAR),XF(IMAR,JMAR) 
    462461      real WEIGHTpb(-1:1,-1:1) 
    463462 
    464       if(imar.gt.ismo) stop'surdimensionner ismo dans mva9 (grid_noro)' 
    465       if(jmar.gt.jsmo) stop'surdimensionner jsmo dans mva9 (grid_noro)' 
    466        
    467463      SUM=0. 
    468464      DO IS=-1,1 
  • LMDZ4/branches/LMDZ4_AR5/libf/phylmd/change_srf_frac_mod.F90

    r996 r1450  
    99!  
    1010! Change Surface Fractions 
    11 ! 
     11! Author J Ghattas 2008 
     12 
    1213  SUBROUTINE change_srf_frac(itime, dtime, jour, & 
    1314       pctsrf, alb1, alb2, tsurf, u10m, v10m, pbl_tke) 
     
    7677    END SELECT 
    7778 
    78     IF (is_modified) THEN 
     79 
    7980!**************************************************************************************** 
    8081! 2)  
     
    8485! 
    8586!**************************************************************************************** 
     87    IF (is_modified) THEN 
    8688   
    8789! Test and exit if a fraction is negative 
     
    150152       CALL pbl_surface_newfrac(itime, pctsrf, pctsrf_old, tsurf, alb1, alb2, u10m, v10m, pbl_tke) 
    151153 
     154    ELSE 
     155       ! No modifcation should be done 
     156       pctsrf(:,:) = pctsrf_old(:,:) 
     157 
    152158    END IF ! is_modified 
    153159 
Note: See TracChangeset for help on using the changeset viewer.