libf/bibio/inithist.f90

module inithist_m

  implicit none

contains

  subroutine inithist(day0, anne0, tstep, nq, t_ops, t_wrt)

    ! From inithist.F, version 1.1.1.1 2004/05/19 12:53:05

    ! Routine d'initialisation des écritures des fichiers histoires LMDZ
    ! au format IOIPSL
    ! Appels successifs des routines : histbeg, histhori, histver,
    ! histdef, histend

    ! Entrées :
    ! day0, anne0: date de référence
    ! tstep : durée du pas de temps en secondes
    ! t_ops : fréquence de l'opération pour IOIPSL
    ! t_wrt : fréquence d'écriture sur le fichier
    ! nq : nombre de traceurs

    ! L. Fairhead, LMD, 03/99

    USE calendar
    use com_io_dyn, only: histid, histvid, histuid
    use histcom
    use dimens_m
    use paramet_m
    use comconst
    use comvert
    use comgeom
    use serre
    use temps
    use ener
    use iniadvtrac_m
    use nr_util, only: pi

    ! Arguments
    integer day0, anne0
    real, intent(in):: tstep, t_ops, t_wrt
    integer nq

    ! Variables locales
    real zjulian
    integer iq
    real rlong(iip1, jjp1), rlat(iip1, jjp1)
    integer uhoriid, vhoriid, thoriid, zvertiid
    integer ii, jj
    integer zan, dayref

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

    ! Appel a histbeg: creation du fichier netcdf et initialisations diverses

    zan = anne0
    dayref = day0
    CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)

    do jj = 1, jjp1
       do ii = 1, iip1
          rlong(ii, jj) = rlonu(ii) * 180. / pi
          rlat(ii, jj) = rlatu(jj) * 180. / pi
       enddo
    enddo

    call histbeg_totreg("dyn_histu.nc", rlong(:,1), rlat(1,:), 1, iip1, 1, &
         jjp1, itau_dyn, zjulian, tstep, uhoriid, histuid)

    do jj = 1, jjm
       do ii = 1, iip1
          rlong(ii, jj) = rlonv(ii) * 180. / pi
          rlat(ii, jj) = rlatv(jj) * 180. / pi
       enddo
    enddo

    ! Creation du fichier histoire pour la grille en V (oblige pour l'instant,
    ! IOIPSL ne permet pas de grilles avec des nombres de point differents dans
    ! un meme fichier)

    call histbeg_totreg('dyn_histv.nc', rlong(:, 1), rlat(1, :jjm), &
         1, iip1, 1, jjm, &
         itau_dyn, zjulian, tstep, vhoriid, histvid)
    !
    ! Appel a histhori pour rajouter les autres grilles horizontales
    !
    do jj = 1, jjp1
       do ii = 1, iip1
          rlong(ii, jj) = rlonv(ii) * 180. / pi
          rlat(ii, jj) = rlatu(jj) * 180. / pi
       enddo
    enddo

    call histbeg_totreg("dyn_hist.nc", rlong(:, 1), rlat(1, :), 1, iip1, 1, &
         jjp1, itau_dyn, zjulian, tstep, thoriid, histid)
    !
    ! Appel a histvert pour la grille verticale
    !
      call histvert(histid, 'presnivs', 'Niveaux pression','mb', llm, &
           presnivs/100., zvertiid,'down')
      call histvert(histvid, 'presnivs', 'Niveaux pression','mb', llm, &
           presnivs/100., zvertiid,'down')
      call histvert(histuid, 'presnivs', 'Niveaux pression','mb', llm, &
           presnivs/100., zvertiid,'down')
    !
    ! Appels a histdef pour la definition des variables a sauvegarder
    !
    ! Vents U
    !
    call histdef(histuid, 'u', 'vent u', 'm/s', &
         iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, &
         'inst(X)', t_ops, t_wrt)
    !
    ! Vents V
    !
    call histdef(histvid, 'v', 'vent v', 'm/s', &
         iip1, jjm, vhoriid, llm, 1, llm, zvertiid, &
         'inst(X)', t_ops, t_wrt)

    !
    ! Temperature potentielle
    !
    call histdef(histid, 'teta', 'temperature potentielle', '-', &
         iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'inst(X)', t_ops, t_wrt)
    !
    ! Geopotentiel
    !
    call histdef(histid, 'phi', 'geopotentiel', '-', &
         iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'inst(X)', t_ops, t_wrt)
    !
    ! Traceurs
    !
    DO iq=1, nq
       call histdef(histid, ttext(iq), ttext(iq), '-', &
            iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
            'inst(X)', t_ops, t_wrt)
    enddo
    !
    ! Masse
    !
    call histdef(histid, 'masse', 'masse', 'kg', &
         iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'inst(X)', t_ops, t_wrt)
    !
    ! Pression au sol
    !
    call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', &
         iip1, jjp1, thoriid, 1, 1, 1, -99, &
         'inst(X)', t_ops, t_wrt)
    !
    ! Geopotentiel au sol
    !
    call histdef(histid, 'phis', 'geopotentiel au sol', '-', &
         iip1, jjp1, thoriid, 1, 1, 1, -99, &
         'inst(X)', t_ops, t_wrt)
    !
    ! Fin
    !
    call histend(histid)
    call histend(histuid)
    call histend(histvid)

  end subroutine inithist

end module inithist_m
1	guez	3	module inithist_m
2
3			implicit none
4
5			contains
6
7	guez	56	subroutine inithist(day0, anne0, tstep, nq, t_ops, t_wrt)
8	guez	3
9	guez	39	! From inithist.F, version 1.1.1.1 2004/05/19 12:53:05
10	guez	3
11	guez	39	! Routine d'initialisation des écritures des fichiers histoires LMDZ
12			! au format IOIPSL
13			! Appels successifs des routines : histbeg, histhori, histver,
14			! histdef, histend
15	guez	3
16	guez	39	! Entrées :
17			! day0, anne0: date de référence
18			! tstep : durée du pas de temps en secondes
19			! t_ops : fréquence de l'opération pour IOIPSL
20			! t_wrt : fréquence d'écriture sur le fichier
21			! nq : nombre de traceurs
22	guez	3
23	guez	39	! L. Fairhead, LMD, 03/99
24	guez	3
25	guez	30	USE calendar
26	guez	56	use com_io_dyn, only: histid, histvid, histuid
27	guez	30	use histcom
28	guez	3	use dimens_m
29			use paramet_m
30			use comconst
31			use comvert
32			use comgeom
33			use serre
34			use temps
35			use ener
36	guez	18	use iniadvtrac_m
37	guez	39	use nr_util, only: pi
38	guez	3
39	guez	39	! Arguments
40	guez	3	integer day0, anne0
41			real, intent(in):: tstep, t_ops, t_wrt
42			integer nq
43
44	guez	39	! Variables locales
45	guez	3	real zjulian
46			integer iq
47	guez	39	real rlong(iip1, jjp1), rlat(iip1, jjp1)
48	guez	3	integer uhoriid, vhoriid, thoriid, zvertiid
49	guez	39	integer ii, jj
50	guez	3	integer zan, dayref
51
52			!-----------------------------------------------------------------------
53
54	guez	39	! Appel a histbeg: creation du fichier netcdf et initialisations diverses
55	guez	3
56			zan = anne0
57			dayref = day0
58			CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
59
60			do jj = 1, jjp1
61			do ii = 1, iip1
62	guez	39	rlong(ii, jj) = rlonu(ii) * 180. / pi
63			rlat(ii, jj) = rlatu(jj) * 180. / pi
64	guez	3	enddo
65			enddo
66
67	guez	56	call histbeg_totreg("dyn_histu.nc", rlong(:,1), rlat(1,:), 1, iip1, 1, &
68			jjp1, itau_dyn, zjulian, tstep, uhoriid, histuid)
69	guez	3
70			do jj = 1, jjm
71			do ii = 1, iip1
72	guez	39	rlong(ii, jj) = rlonv(ii) * 180. / pi
73			rlat(ii, jj) = rlatv(jj) * 180. / pi
74	guez	3	enddo
75			enddo
76
77	guez	56	! Creation du fichier histoire pour la grille en V (oblige pour l'instant,
78			! IOIPSL ne permet pas de grilles avec des nombres de point differents dans
79			! un meme fichier)
80
81	guez	39	call histbeg_totreg('dyn_histv.nc', rlong(:, 1), rlat(1, :jjm), &
82	guez	3	1, iip1, 1, jjm, &
83	guez	56	itau_dyn, zjulian, tstep, vhoriid, histvid)
84	guez	3	!
85	guez	39	! Appel a histhori pour rajouter les autres grilles horizontales
86	guez	3	!
87			do jj = 1, jjp1
88			do ii = 1, iip1
89	guez	39	rlong(ii, jj) = rlonv(ii) * 180. / pi
90			rlat(ii, jj) = rlatu(jj) * 180. / pi
91	guez	3	enddo
92			enddo
93
94	guez	56	call histbeg_totreg("dyn_hist.nc", rlong(:, 1), rlat(1, :), 1, iip1, 1, &
95			jjp1, itau_dyn, zjulian, tstep, thoriid, histid)
96	guez	3	!
97	guez	39	! Appel a histvert pour la grille verticale
98	guez	3	!
99	guez	56	call histvert(histid, 'presnivs', 'Niveaux pression','mb', llm, &
100			presnivs/100., zvertiid,'down')
101			call histvert(histvid, 'presnivs', 'Niveaux pression','mb', llm, &
102			presnivs/100., zvertiid,'down')
103			call histvert(histuid, 'presnivs', 'Niveaux pression','mb', llm, &
104			presnivs/100., zvertiid,'down')
105	guez	3	!
106	guez	39	! Appels a histdef pour la definition des variables a sauvegarder
107	guez	3	!
108	guez	39	! Vents U
109	guez	3	!
110	guez	56	call histdef(histuid, 'u', 'vent u', 'm/s', &
111	guez	3	iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, &
112	guez	15	'inst(X)', t_ops, t_wrt)
113	guez	3	!
114	guez	39	! Vents V
115	guez	3	!
116	guez	56	call histdef(histvid, 'v', 'vent v', 'm/s', &
117	guez	3	iip1, jjm, vhoriid, llm, 1, llm, zvertiid, &
118	guez	15	'inst(X)', t_ops, t_wrt)
119	guez	3
120			!
121	guez	39	! Temperature potentielle
122	guez	3	!
123	guez	56	call histdef(histid, 'teta', 'temperature potentielle', '-', &
124	guez	3	iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
125	guez	15	'inst(X)', t_ops, t_wrt)
126	guez	3	!
127	guez	39	! Geopotentiel
128	guez	3	!
129	guez	56	call histdef(histid, 'phi', 'geopotentiel', '-', &
130	guez	3	iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
131	guez	15	'inst(X)', t_ops, t_wrt)
132	guez	3	!
133	guez	39	! Traceurs
134	guez	3	!
135	guez	39	DO iq=1, nq
136	guez	56	call histdef(histid, ttext(iq), ttext(iq), '-', &
137	guez	3	iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
138	guez	15	'inst(X)', t_ops, t_wrt)
139	guez	3	enddo
140			!
141	guez	39	! Masse
142	guez	3	!
143	guez	56	call histdef(histid, 'masse', 'masse', 'kg', &
144			iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
145	guez	15	'inst(X)', t_ops, t_wrt)
146	guez	3	!
147	guez	39	! Pression au sol
148	guez	3	!
149	guez	56	call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', &
150	guez	3	iip1, jjp1, thoriid, 1, 1, 1, -99, &
151	guez	15	'inst(X)', t_ops, t_wrt)
152	guez	3	!
153	guez	56	! Geopotentiel au sol
154	guez	3	!
155	guez	56	call histdef(histid, 'phis', 'geopotentiel au sol', '-', &
156	guez	3	iip1, jjp1, thoriid, 1, 1, 1, -99, &
157	guez	15	'inst(X)', t_ops, t_wrt)
158	guez	3	!
159	guez	39	! Fin
160	guez	3	!
161	guez	56	call histend(histid)
162			call histend(histuid)
163			call histend(histvid)
164	guez	3
165			end subroutine inithist
166
167			end module inithist_m