/[lmdze]/trunk/filtrez/inifilr.f
ViewVC logotype

Diff of /trunk/filtrez/inifilr.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 131 by guez, Thu Sep 18 19:56:46 2014 UTC revision 132 by guez, Fri Mar 20 16:31:06 2015 UTC
# Line 43  contains Line 43  contains
43      USE serre, ONLY : grossismx      USE serre, ONLY : grossismx
44    
45      ! Local:      ! Local:
46      REAL dlonu(iim), dlatu(jjm)      REAL dlatu(jjm)
47      REAL rlamda(2: iim), eignvl(iim)      REAL rlamda(2: iim), eignvl(iim)
48    
49      REAL lamdamax, cof      REAL lamdamax, cof
50      INTEGER i, j, modemax, imx, k, kf      INTEGER i, j, modemax, imx, k, kf
51      REAL dymin, dxmin, colat0      REAL dymin, colat0
52      REAL eignft(iim, iim), coff      REAL eignft(iim, iim), coff
53    
54      !-----------------------------------------------------------      !-----------------------------------------------------------
55    
56      print *, "Call sequence information: inifilr"      print *, "Call sequence information: inifilr"
57    
     DO i = 1, iim  
        dlonu(i) = xprimu(i)  
     END DO  
   
58      CALL inifgn(eignvl)      CALL inifgn(eignvl)
59    
60      PRINT *, 'EIGNVL '      PRINT *, 'EIGNVL '
# Line 78  contains Line 74  contains
74         dlatu(j) = rlatu(j) - rlatu(j+1)         dlatu(j) = rlatu(j) - rlatu(j+1)
75      END DO      END DO
76    
     dxmin = dlonu(1)  
     DO i = 2, iim  
        dxmin = min(dxmin, dlonu(i))  
     END DO  
77      dymin = dlatu(1)      dymin = dlatu(1)
78      DO j = 2, jjm      DO j = 2, jjm
79         dymin = min(dymin, dlatu(j))         dymin = min(dymin, dlatu(j))
80      END DO      END DO
81    
82      colat0 = min(0.5, dymin/dxmin)      colat0 = min(0.5, dymin / minval(xprimu(:iim)))
83    
84      PRINT *, 'colat0 = ', colat0      PRINT *, 'colat0 = ', colat0
85    
Legend:
Removed from v.131  
changed lines
  Added in v.132
  ViewVC Help
Powered by ViewVC 1.1.21