trunk/dyn3d/dynetat0_chosen.f90

module dynetat0_chosen_m

  IMPLICIT NONE

  integer, protected, save:: day_ref ! jour de l'ann\'ee de l'\'etat initial
  ! (= 350 si 20 d\'ecembre par exemple)

  integer, protected, save:: annee_ref
  ! Annee de l'etat initial (avec 4 chiffres)

  REAL, protected, save:: clon ! longitude of the center of the zoom, in rad
  real, protected, save:: clat ! latitude of the center of the zoom, in rad

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

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

  real, protected, save:: taux, tauy
  ! raideur de la transition de l'int\'erieur \`a l'ext\'erieur du zoom
  
  real, protected, save:: pa ! in Pa

contains

  SUBROUTINE dynetat0_chosen

    ! This procedure reads the initial state of the atmosphere. Values
    ! that were chosen in ce0l.

    ! Libraries:
    use netcdf, only: NF90_NOWRITE
    use netcdf95, only: nf95_open, nf95_inq_varid, NF95_CLOSE, NF95_Gw_VAR
    use nr_util, only: assert

    use comconst, only: dtvr
    use conf_gcm_m, only: raz_date
    use dimensions, only: iim, jjm, llm
    use unit_nml_m, only: unit_nml

    ! Local: 
    REAL, allocatable:: tab_cntrl(:) ! tableau des param\`etres du run
    INTEGER ncid, varid

    namelist /dynetat0_nml/ day_ref, annee_ref

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

    print *, "Call sequence information: dynetat0_chosen"

    ! Fichier \'etat initial :
    call nf95_open("start.nc", NF90_NOWRITE, ncid)

    call nf95_inq_varid(ncid, "controle", varid)
    call NF95_Gw_VAR(ncid, varid, tab_cntrl)

    call assert(int(tab_cntrl(1)) == iim, "dynetat0_chosen tab_cntrl iim") 
    call assert(int(tab_cntrl(2)) == jjm, "dynetat0_chosen tab_cntrl jjm") 
    call assert(int(tab_cntrl(3)) == llm, "dynetat0_chosen tab_cntrl llm") 

    IF (dtvr /= tab_cntrl(12)) THEN
       print *, 'Warning: the time steps from day_step and "start.nc" ' // &
            'are different.'
       print *, 'dtvr from day_step: ', dtvr
       print *, 'dtvr from "start.nc": ', tab_cntrl(12)
       print *, 'Using the value from day_step.'
    ENDIF

    pa = tab_cntrl(18)

    clon = tab_cntrl(20)
    clat = tab_cntrl(21)
    grossismx = tab_cntrl(22)
    grossismy = tab_cntrl(23)
    dzoomx = tab_cntrl(25)
    dzoomy = tab_cntrl(26)
    taux = tab_cntrl(28)
    tauy = tab_cntrl(29)

    if (raz_date) then
       ! Default values:
       day_ref = 1
       annee_ref = 1998
       
       print *, "Enter namelist 'dynetat0_nml'."
       read(unit = *, nml = dynetat0_nml)
       write(unit_nml, nml = dynetat0_nml)
    else
       day_ref = tab_cntrl(4)
       annee_ref = tab_cntrl(5)
    end if

    call NF95_CLOSE(ncid)

  END SUBROUTINE dynetat0_chosen

  !********************************************************************

  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
    namelist /dynetat0_nml/ day_ref, annee_ref

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

    ! 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

    ! Default values:
    day_ref = 1
    annee_ref = 1998

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

    pa = 5e4

  end subroutine read_serre

end module dynetat0_chosen_m
1	guez	313	module dynetat0_chosen_m
2
3			IMPLICIT NONE
4
5			integer, protected, save:: day_ref ! jour de l'ann\'ee de l'\'etat initial
6			! (= 350 si 20 d\'ecembre par exemple)
7
8			integer, protected, save:: annee_ref
9			! Annee de l'etat initial (avec 4 chiffres)
10
11			REAL, protected, save:: clon ! longitude of the center of the zoom, in rad
12			real, protected, save:: clat ! latitude of the center of the zoom, in rad
13
14			real, protected, save:: grossismx, grossismy
15			! facteurs de grossissement du zoom, selon la longitude et la latitude
16			! = 2 si 2 fois, = 3 si 3 fois, etc.
17
18			real, protected, save:: dzoomx, dzoomy
19			! extensions en longitude et latitude de la zone du zoom (fractions
20			! de la zone totale)
21
22			real, protected, save:: taux, tauy
23			! raideur de la transition de l'int\'erieur \`a l'ext\'erieur du zoom
24
25			real, protected, save:: pa ! in Pa
26
27			contains
28
29			SUBROUTINE dynetat0_chosen
30
31			! This procedure reads the initial state of the atmosphere. Values
32			! that were chosen in ce0l.
33
34			! Libraries:
35			use netcdf, only: NF90_NOWRITE
36			use netcdf95, only: nf95_open, nf95_inq_varid, NF95_CLOSE, NF95_Gw_VAR
37			use nr_util, only: assert
38
39			use comconst, only: dtvr
40			use conf_gcm_m, only: raz_date
41			use dimensions, only: iim, jjm, llm
42			use unit_nml_m, only: unit_nml
43
44			! Local:
45			REAL, allocatable:: tab_cntrl(:) ! tableau des param\`etres du run
46			INTEGER ncid, varid
47
48			namelist /dynetat0_nml/ day_ref, annee_ref
49
50			!-----------------------------------------------------------------------
51
52			print *, "Call sequence information: dynetat0_chosen"
53
54			! Fichier \'etat initial :
55			call nf95_open("start.nc", NF90_NOWRITE, ncid)
56
57			call nf95_inq_varid(ncid, "controle", varid)
58			call NF95_Gw_VAR(ncid, varid, tab_cntrl)
59
60			call assert(int(tab_cntrl(1)) == iim, "dynetat0_chosen tab_cntrl iim")
61			call assert(int(tab_cntrl(2)) == jjm, "dynetat0_chosen tab_cntrl jjm")
62			call assert(int(tab_cntrl(3)) == llm, "dynetat0_chosen tab_cntrl llm")
63
64			IF (dtvr /= tab_cntrl(12)) THEN
65			print *, 'Warning: the time steps from day_step and "start.nc" ' // &
66			'are different.'
67			print *, 'dtvr from day_step: ', dtvr
68			print *, 'dtvr from "start.nc": ', tab_cntrl(12)
69			print *, 'Using the value from day_step.'
70			ENDIF
71
72			pa = tab_cntrl(18)
73
74			clon = tab_cntrl(20)
75			clat = tab_cntrl(21)
76			grossismx = tab_cntrl(22)
77			grossismy = tab_cntrl(23)
78			dzoomx = tab_cntrl(25)
79			dzoomy = tab_cntrl(26)
80			taux = tab_cntrl(28)
81			tauy = tab_cntrl(29)
82
83			if (raz_date) then
84			! Default values:
85			day_ref = 1
86			annee_ref = 1998
87
88			print *, "Enter namelist 'dynetat0_nml'."
89			read(unit = *, nml = dynetat0_nml)
90			write(unit_nml, nml = dynetat0_nml)
91			else
92			day_ref = tab_cntrl(4)
93			annee_ref = tab_cntrl(5)
94			end if
95
96			call NF95_CLOSE(ncid)
97
98			END SUBROUTINE dynetat0_chosen
99
100			!********************************************************************
101
102			subroutine read_serre
103
104			use unit_nml_m, only: unit_nml
105			use nr_util, only: assert, pi
106
107			REAL:: clon_deg = 0. ! longitude of the center of the zoom, in degrees
108			real:: clat_deg = 0. ! latitude of the center of the zoom, in degrees
109
110			namelist /serre_nml/ clon_deg, clat_deg, grossismx, grossismy, dzoomx, &
111			dzoomy, taux, tauy
112			namelist /dynetat0_nml/ day_ref, annee_ref
113
114			!-------------------------------------------------
115
116			! Default values:
117			grossismx = 1.
118			grossismy = 1.
119			dzoomx = 0.2
120			dzoomy = 0.2
121			taux = 3.
122			tauy = 3.
123
124			print *, "Enter namelist 'serre_nml'."
125			read(unit=*, nml=serre_nml)
126			write(unit_nml, nml=serre_nml)
127
128			call assert(grossismx >= 1. .and. grossismy >= 1., "read_serre grossism")
129			call assert(dzoomx > 0., dzoomx < 1., dzoomy < 1., &
130			"read_serre dzoomx dzoomy")
131			clon = clon_deg / 180. * pi
132			clat = clat_deg / 180. * pi
133
134			! Default values:
135			day_ref = 1
136			annee_ref = 1998
137
138			print *, "Enter namelist 'dynetat0_nml'."
139			read(unit = *, nml = dynetat0_nml)
140			write(unit_nml, nml = dynetat0_nml)
141
142			pa = 5e4
143
144			end subroutine read_serre
145
146			end module dynetat0_chosen_m