/[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 265 by guez, Tue Mar 20 09:35:59 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      USE dimensions, ONLY : iim, jjm
43      USE dynetat0_m, ONLY : rlatu, rlatv, xprimu, grossismx      USE dynetat0_m, ONLY : rlatu, rlatv, xprimu, grossismx
44      use inifgn_m, only: inifgn      use inifgn_m, only: inifgn
45      use inifilr_hemisph_m, only: inifilr_hemisph      use inifilr_hemisph_m, only: inifilr_hemisph
# Line 72  contains Line 72  contains
72    
73      rlamda = iim / pi / colat0 * grossismx / sqrt(- eignvl(2: iim))      rlamda = iim / pi / colat0 * grossismx / sqrt(- eignvl(2: iim))
74      print *, "1 / rlamda(iim) = ", 1. / rlamda(iim)      print *, "1 / rlamda(iim) = ", 1. / rlamda(iim)
   
75      ! This is demonstrated in the notes but just to be sure:      ! This is demonstrated in the notes but just to be sure:
76      call assert(rlamda(iim) * colat0 >= 1. - epsilon(0.), &      call assert(rlamda(iim) * colat0 >= 1. - 2. * epsilon(0.), &
77           "inifilr rlamda(iim) * colat0")           "inifilr rlamda(iim) * colat0")
78    
79      call new_unit(unit)      call new_unit(unit)
Legend:
Removed from v.169  
changed lines
  Added in v.265
  ViewVC Help
Powered by ViewVC 1.1.21