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

Diff of /trunk/dyn3d/read_serre.f

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

trunk/dyn3d/serre.f revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC trunk/dyn3d/read_serre.f revision 254 by guez, Mon Feb 5 10:39:38 2018 UTC
# Line 1  Line 1 
1  module serre  module read_serre_m
2    
3    implicit none    implicit none
4    
5    REAL:: clon = 0. ! longitude of the center of the zoom, in degrees  contains
   real:: clat = 0. ! latitude of the center of the zoom, in degrees  
6    
7    real, save:: transx, transy    subroutine read_serre
   real, save:: alphax, alphay ! anciennes formulations des grossissements  
   real, save:: pxo, pyo  
8    
9    real:: grossismx = 1. ! facteur de grossissement du zoom, selon la longitude      use dynetat0_m, only: clon, clat, grossismx, grossismy, dzoomx, dzoomy, &
10    real:: grossismy = 1. ! facteur de grossissement du zoom, selon la latitude           taux, tauy
11        use unit_nml_m, only: unit_nml
12        use nr_util, only: assert, pi
13    
14    real:: dzoomx = 0.      REAL:: clon_deg = 0. ! longitude of the center of the zoom, in degrees
15    ! extension en longitude de la zone du zoom (fraction de la zone totale)      real:: clat_deg = 0. ! latitude of the center of the zoom, in degrees
16    
17    real:: dzoomy = 0.      namelist /serre_nml/ clon_deg, clat_deg, grossismx, grossismy, dzoomx, &
18    ! extension en latitude de la zone du zoom (fraction de la zone totale)           dzoomy, taux, tauy
19    
20    real:: taux = 3. ! raideur du zoom en X      !-------------------------------------------------
   real:: tauy = 3. ! raideur du zoom en Y  
21    
22  end module serre      ! Default values:
23        grossismx = 1.
24        grossismy = 1.
25        dzoomx = 0.2
26        dzoomy = 0.2
27        taux = 3.
28        tauy = 3.
29    
30        print *, "Enter namelist 'serre_nml'."
31        read(unit=*, nml=serre_nml)
32        write(unit_nml, nml=serre_nml)
33    
34        call assert(grossismx >= 1. .and. grossismy >= 1., "read_serre grossism")
35        call assert(dzoomx > 0., dzoomx < 1., dzoomy < 1., &
36             "read_serre dzoomx dzoomy")
37        clon = clon_deg / 180. * pi
38        clat = clat_deg / 180. * pi
39    
40      end subroutine read_serre
41    
42    end module read_serre_m
Legend:
Removed from v.82  
changed lines
  Added in v.254
  ViewVC Help
Powered by ViewVC 1.1.21