Sources/dyn3d/read_serre.f

module serre

  implicit none

  REAL:: clon = 0. ! longitude of the center of the zoom, in rad
  real:: clat = 0. ! latitude of the center of the zoom, in rad

  real:: grossismx = 1., grossismy = 1.
  ! facteurs de grossissement du zoom, selon la longitude et la latitude
  ! = 2 si 2 fois, = 3 si 3 fois, etc.

  real:: dzoomx = 0.2, dzoomy = 0.2
  ! extensions en longitude et latitude de la zone du zoom (fractions
  ! de la zone totale)

  real:: taux = 3., tauy = 3. 
  ! raideur de la transition de l'intérieur à l'extérieur du zoom
  
contains

  subroutine read_serre

    use unit_nml_m, only: unit_nml
    use nr_util, only: assert, pi

    REAL:: clon_deg = 0. ! longitude of the center of the zoom, in degrees
    real:: clat_deg = 0. ! latitude of the center of the zoom, in degrees

    namelist /serre_nml/ clon_deg, clat_deg, grossismx, grossismy, dzoomx, &
         dzoomy, taux, tauy

    !-------------------------------------------------

    print *, "Enter namelist 'serre_nml'."
    read(unit=*, nml=serre_nml)
    write(unit_nml, nml=serre_nml)

    call assert(grossismx >= 1. .and. grossismy >= 1., "read_serre grossism")
    call assert(dzoomx > 0., dzoomx < 1., dzoomy < 1., &
         "read_serre dzoomx dzoomy")
    clon = clon_deg / 180. * pi
    clat = clat_deg / 180. * pi

  end subroutine read_serre

end module serre
1	guez	3	module serre
2
3			implicit none
4
5	guez	127	REAL:: clon = 0. ! longitude of the center of the zoom, in rad
6			real:: clat = 0. ! latitude of the center of the zoom, in rad
7	guez	3
8	guez	119	real:: grossismx = 1., grossismy = 1.
9			! facteurs de grossissement du zoom, selon la longitude et la latitude
10			! = 2 si 2 fois, = 3 si 3 fois, etc.
11	guez	3
12	guez	122	real:: dzoomx = 0.2, dzoomy = 0.2
13	guez	119	! extensions en longitude et latitude de la zone du zoom (fractions
14			! de la zone totale)
15	guez	3
16	guez	119	real:: taux = 3., tauy = 3.
17			! raideur de la transition de l'intérieur à l'extérieur du zoom
18
19	guez	112	contains
20
21			subroutine read_serre
22
23			use unit_nml_m, only: unit_nml
24	guez	127	use nr_util, only: assert, pi
25	guez	112
26	guez	127	REAL:: clon_deg = 0. ! longitude of the center of the zoom, in degrees
27			real:: clat_deg = 0. ! latitude of the center of the zoom, in degrees
28	guez	112
29	guez	127	namelist /serre_nml/ clon_deg, clat_deg, grossismx, grossismy, dzoomx, &
30			dzoomy, taux, tauy
31
32	guez	112	!-------------------------------------------------
33
34			print *, "Enter namelist 'serre_nml'."
35			read(unit=*, nml=serre_nml)
36			write(unit_nml, nml=serre_nml)
37
38	guez	113	call assert(grossismx >= 1. .and. grossismy >= 1., "read_serre grossism")
39	guez	122	call assert(dzoomx > 0., dzoomx < 1., dzoomy < 1., &
40			"read_serre dzoomx dzoomy")
41	guez	127	clon = clon_deg / 180. * pi
42			clat = clat_deg / 180. * pi
43	guez	112
44			end subroutine read_serre
45
46	guez	3	end module serre