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

Diff of /trunk/filtrez/inifilr.f

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

trunk/Sources/filtrez/inifilr.f revision 169 by guez, Mon Sep 14 17:13:16 2015 UTC trunk/filtrez/inifilr.f revision 313 by guez, Mon Dec 10 15:54:30 2018 UTC
# Line 39  contains Line 39  contains
39      ! (stretching coefficient) is set equal to zero for the regular      ! (stretching coefficient) is set equal to zero for the regular
40      ! grid case.      ! grid case.
41    
42      USE dimens_m, ONLY : iim, jjm      ! Libraries:
     USE dynetat0_m, ONLY : rlatu, rlatv, xprimu, grossismx  
     use inifgn_m, only: inifgn  
     use inifilr_hemisph_m, only: inifilr_hemisph  
43      use jumble, only: new_unit      use jumble, only: new_unit
44      use nr_util, only: pi, ifirstloc, assert      use nr_util, only: pi, ifirstloc, assert
45    
46        USE dimensions, ONLY: iim, jjm
47        USE dynetat0_m, ONLY: rlatu, rlatv, xprimu
48        USE dynetat0_chosen_m, ONLY: grossismx
49        use inifgn_m, only: inifgn
50        use inifilr_hemisph_m, only: inifilr_hemisph
51    
52      ! Local:      ! Local:
53    
54      REAL dlatu(jjm)      REAL dlatu(jjm)
# Line 72  contains Line 75  contains
75    
76      rlamda = iim / pi / colat0 * grossismx / sqrt(- eignvl(2: iim))      rlamda = iim / pi / colat0 * grossismx / sqrt(- eignvl(2: iim))
77      print *, "1 / rlamda(iim) = ", 1. / rlamda(iim)      print *, "1 / rlamda(iim) = ", 1. / rlamda(iim)
   
78      ! This is demonstrated in the notes but just to be sure:      ! This is demonstrated in the notes but just to be sure:
79      call assert(rlamda(iim) * colat0 >= 1. - epsilon(0.), &      call assert(rlamda(iim) * colat0 >= 1. - 2. * epsilon(0.), &
80           "inifilr rlamda(iim) * colat0")           "inifilr rlamda(iim) * colat0")
81    
82      call new_unit(unit)      call new_unit(unit)
Legend:
Removed from v.169  
changed lines
  Added in v.313
  ViewVC Help
Powered by ViewVC 1.1.21