Changeset 1450 for LMDZ4/branches


Ignore:
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:
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.