trunk/dyn3d/serre.f

module serre

  implicit none

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

  real:: grossismx = 1. ! facteur de grossissement du zoom, selon la longitude
  real:: grossismy = 1. ! facteur de grossissement du zoom, selon la latitude

  real:: dzoomx = 0.
  ! extension en longitude de la zone du zoom (fraction de la zone totale)

  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

contains

  subroutine read_serre

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

    namelist /serre_nml/ clon, clat, 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")

  end subroutine read_serre

end module serre
1	guez	3	module serre
2
3			implicit none
4
5	guez	39	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
7	guez	3
8	guez	39	real:: grossismx = 1. ! facteur de grossissement du zoom, selon la longitude
9			real:: grossismy = 1. ! facteur de grossissement du zoom, selon la latitude
10	guez	3
11	guez	39	real:: dzoomx = 0.
12			! extension en longitude de la zone du zoom (fraction de la zone totale)
13	guez	3
14	guez	39	real:: dzoomy = 0.
15			! extension en latitude de la zone du zoom (fraction de la zone totale)
16	guez	3
17	guez	39	real:: taux = 3. ! raideur du zoom en X
18			real:: tauy = 3. ! raideur du zoom en Y
19	guez	3
20	guez	112	contains
21
22			subroutine read_serre
23
24			use unit_nml_m, only: unit_nml
25	guez	113	use nr_util, only: assert
26	guez	112
27			namelist /serre_nml/ clon, clat, grossismx, grossismy, dzoomx, dzoomy, &
28			taux, tauy
29
30			!-------------------------------------------------
31
32			print *, "Enter namelist 'serre_nml'."
33			read(unit=*, nml=serre_nml)
34			write(unit_nml, nml=serre_nml)
35
36	guez	113	call assert(grossismx >= 1. .and. grossismy >= 1., "read_serre grossism")
37	guez	112
38			end subroutine read_serre
39
40	guez	3	end module serre