Sources/dyn3d/read_serre.f

module read_serre_m

  implicit none

contains

  subroutine read_serre

    use dynetat0_m, only: clon, clat, grossismx, grossismy, dzoomx, dzoomy, &
         taux, tauy
    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

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

    ! Default values:
    grossismx = 1.
    grossismy = 1.
    dzoomx = 0.2
    dzoomy = 0.2
    taux = 3.
    tauy = 3. 

    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 read_serre_m
1	guez	139	module read_serre_m
2	guez	3
3			implicit none
4
5	guez	112	contains
6
7			subroutine read_serre
8
9	guez	139	use dynetat0_m, only: clon, clat, grossismx, grossismy, dzoomx, dzoomy, &
10			taux, tauy
11	guez	112	use unit_nml_m, only: unit_nml
12	guez	127	use nr_util, only: assert, pi
13	guez	112
14	guez	127	REAL:: clon_deg = 0. ! longitude of the center of the zoom, in degrees
15			real:: clat_deg = 0. ! latitude of the center of the zoom, in degrees
16	guez	112
17	guez	127	namelist /serre_nml/ clon_deg, clat_deg, grossismx, grossismy, dzoomx, &
18			dzoomy, taux, tauy
19
20	guez	112	!-------------------------------------------------
21
22	guez	139	! Default values:
23			grossismx = 1.
24			grossismy = 1.
25			dzoomx = 0.2
26			dzoomy = 0.2
27			taux = 3.
28			tauy = 3.
29
30	guez	112	print *, "Enter namelist 'serre_nml'."
31			read(unit=*, nml=serre_nml)
32			write(unit_nml, nml=serre_nml)
33
34	guez	113	call assert(grossismx >= 1. .and. grossismy >= 1., "read_serre grossism")
35	guez	122	call assert(dzoomx > 0., dzoomx < 1., dzoomy < 1., &
36			"read_serre dzoomx dzoomy")
37	guez	127	clon = clon_deg / 180. * pi
38			clat = clat_deg / 180. * pi
39	guez	112
40			end subroutine read_serre
41
42	guez	139	end module read_serre_m