trunk/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	module read_serre_m
2
3	implicit none
4
5	contains
6
7	subroutine read_serre
8
9	use dynetat0_m, only: clon, clat, grossismx, grossismy, dzoomx, dzoomy, &
10	taux, tauy
11	use unit_nml_m, only: unit_nml
12	use nr_util, only: assert, pi
13
14	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
17	namelist /serre_nml/ clon_deg, clat_deg, grossismx, grossismy, dzoomx, &
18	dzoomy, taux, tauy
19
20	!-------------------------------------------------
21
22	! 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