/[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 167 by guez, Mon Aug 24 16:30:33 2015 UTC revision 168 by guez, Wed Sep 9 10:41:47 2015 UTC
# Line 57  contains Line 57  contains
57    
58      REAL dlatu(jjm)      REAL dlatu(jjm)
59      REAL rlamda(2:iim) ! > 0, in descending order      REAL rlamda(2:iim) ! > 0, in descending order
60      real eignvl(iim) ! eigenvalues sorted in descending order (<= 0)      real eignvl(iim) ! eigenvalues (<= 0) sorted in descending order
61      INTEGER j, unit      INTEGER j, unit
62      REAL colat0 ! > 0      REAL colat0 ! > 0
63      integer j1 ! index of smallest positive latitude      integer j1 ! index of smallest positive latitude
# Line 77  contains Line 77  contains
77      PRINT *, 'colat0 = ', colat0      PRINT *, 'colat0 = ', colat0
78    
79      rlamda = iim / (pi * colat0 / grossismx) / sqrt(- eignvl(2: iim))      rlamda = iim / (pi * colat0 / grossismx) / sqrt(- eignvl(2: iim))
80        print *, "1 / rlamda(iim) = ", 1. / rlamda(iim)
81      call new_unit(unit)      call new_unit(unit)
82      open(unit, file = "modfrst.csv", status = "replace", action = "write")      open(unit, file = "modfrst.csv", status = "replace", action = "write")
83      write(unit, fmt = *) '"rlat (degrees)" modfrst' ! title line      write(unit, fmt = *) '"rlat (degrees)" modfrst' ! title line
84    
    ! D\'etermination de jfilt[ns][uv] :  
   
85      j1 = ifirstloc(rlatu <= 0.)      j1 = ifirstloc(rlatu <= 0.)
86    
87      call inifilr_hemisph(rlatu(j1 - 1:2:- 1), colat0, rlamda, unit, eignfnv, &      call inifilr_hemisph(rlatu(j1 - 1:2:- 1), colat0, rlamda, unit, eignfnv, &
Legend:
Removed from v.167  
changed lines
  Added in v.168
  ViewVC Help
Powered by ViewVC 1.1.21