/[lmdze]/trunk/Sources/dyn3d/read_serre.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/read_serre.f

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

revision 112 by guez, Thu Sep 18 13:36:51 2014 UTC revision 119 by guez, Wed Jan 7 14:34:57 2015 UTC
# Line 5  module serre Line 5  module serre
5    REAL:: clon = 0. ! longitude of the center of the zoom, in degrees    REAL:: clon = 0. ! longitude of the center of the zoom, in degrees
6    real:: clat = 0. ! latitude of the center of the zoom, in degrees    real:: clat = 0. ! latitude of the center of the zoom, in degrees
7    
8    real, save:: transx, transy    real:: grossismx = 1., grossismy = 1.
9    real, save:: alphax, alphay ! anciennes formulations des grossissements    ! facteurs de grossissement du zoom, selon la longitude et la latitude
10    real, save:: pxo, pyo    ! = 2 si 2 fois, = 3 si 3 fois, etc.
11    
12    real:: grossismx = 1. ! facteur de grossissement du zoom, selon la longitude    real:: dzoomx = 0., dzoomy = 0.
13    real:: grossismy = 1. ! facteur de grossissement du zoom, selon la latitude    ! extensions en longitude et latitude de la zone du zoom (fractions
14      ! de la zone totale)
15    real:: dzoomx = 0.  
16    ! extension en longitude de la zone du zoom (fraction de la zone totale)    real:: taux = 3., tauy = 3.
17      ! raideur de la transition de l'intérieur à l'extérieur du zoom
18    real:: dzoomy = 0.    
   ! extension en latitude de la zone du zoom (fraction de la zone totale)  
   
   real:: taux = 3. ! raideur du zoom en X  
   real:: tauy = 3. ! raideur du zoom en Y  
   
19  contains  contains
20    
21    subroutine read_serre    subroutine read_serre
22    
23      use unit_nml_m, only: unit_nml      use unit_nml_m, only: unit_nml
24        use nr_util, only: assert
25    
26      namelist /serre_nml/ clon, clat, grossismx, grossismy, dzoomx, dzoomy, &      namelist /serre_nml/ clon, clat, grossismx, grossismy, dzoomx, dzoomy, &
27           taux, tauy           taux, tauy
# Line 36  contains Line 32  contains
32      read(unit=*, nml=serre_nml)      read(unit=*, nml=serre_nml)
33      write(unit_nml, nml=serre_nml)      write(unit_nml, nml=serre_nml)
34    
35      IF (grossismx < 1.) THEN      call assert(grossismx >= 1. .and. grossismy >= 1., "read_serre grossism")
        PRINT *, 'Error: grossismx < 1'  
        STOP 1  
     ELSE  
        alphax = 1. - 1. / grossismx  
     ENDIF  
     IF (grossismy < 1.) THEN  
        PRINT *, 'Error: grossismy < 1'  
        STOP 1  
     ELSE  
        alphay = 1. - 1. / grossismy  
     ENDIF  
     PRINT *, 'alphax = ', alphax  
     PRINT *, 'alphay = ', alphay  
36    
37    end subroutine read_serre    end subroutine read_serre
38    
Legend:
Removed from v.112  
changed lines
  Added in v.119
  ViewVC Help
Powered by ViewVC 1.1.21