Sources/dyn3d/initdynav.f

module initdynav_m

  ! This module is clean: no C preprocessor directive, no include line

  implicit none

contains

  subroutine initdynav(day0, anne0, tstep, nq, fileid, t_ops, t_wrt)

    ! From initdynav.F, version 1.1.1.1, 2004/05/19 12:53:05
    !   L. Fairhead, LMD

    USE histcom, ONLY: histbeg_totreg, histdef, histend, histvert
    use calendar, ONLY: ymds2ju
    USE dimens_m, ONLY : llm
    USE paramet_m, ONLY : iip1, jjp1
    USE comconst, ONLY : pi
    USE comvert, ONLY : nivsigs
    USE comgeom, ONLY : rlatu, rlonv
    USE temps, ONLY : itau_dyn
    USE iniadvtrac_m, ONLY : ttext

    !   Routine d'initialisation des ecritures des fichiers histoires LMDZ
    !   au format IOIPSL. Initialisation du fichier histoire moyenne.

    !   Appels succesifs des routines: histbeg
    !                                  histhori
    !                                  histver
    !                                  histdef
    !                                  histend

    !   Entree:
    !      day0,anne0: date de reference
    !      tstep : frequence d'ecriture
    !      t_ops: frequence de l'operation pour IOIPSL
    !      t_wrt: frequence d'ecriture sur le fichier
    !      nq: nombre de traceurs

    !   Sortie:
    !      fileid: ID du fichier netcdf cree

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

    !   Variables locales

    real zjulian
    integer iq
    real rlong(iip1,jjp1), rlat(iip1,jjp1)
    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) = rlonv(ii) * 180. / pi
          rlat(ii,jj)  = rlatu(jj) * 180. / pi
       enddo
    enddo

    call histbeg_totreg('dyn_hist_ave.nc', rlong(:,1), rlat(1,:), &
         1, iip1, 1, jjp1, &
         itau_dyn, zjulian, tstep, thoriid, fileid)


    !  Appel a histvert pour la grille verticale

    call histvert(fileid, 'sigss', 'Niveaux sigma','Pa', &
         llm, nivsigs, zvertiid)

    !  Appels a histdef pour la definition des variables a sauvegarder

    !  Vents U

    write(6,*)'inithistave',tstep
    call histdef(fileid, 'u', 'vents u scalaires moyennes', &
         'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)


    !  Vents V

    call histdef(fileid, 'v', 'vents v scalaires moyennes', &
         'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)


    !  Temperature

    call histdef(fileid, 'temp', 'temperature moyennee', 'K', &
         iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)

    !  Temperature potentielle

    call histdef(fileid, 'theta', 'temperature potentielle', 'K', &
         iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)


    !  Geopotentiel

    call histdef(fileid, 'phi', 'geopotentiel moyenne', '-', &
         iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)

    !  Traceurs

    DO iq=1,nq
       call histdef(fileid, ttext(iq), ttext(iq), '-', &
            iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
            'ave(X)', t_ops, t_wrt)
    enddo

    !  Masse

    call histdef(fileid, 'masse', 'masse', 'kg', &
         iip1, jjp1, thoriid, 1, 1, 1, -99, &
         'ave(X)', t_ops, t_wrt)

    !  Pression au sol

    call histdef(fileid, 'ps', 'pression naturelle au sol', 'Pa', &
         iip1, jjp1, thoriid, 1, 1, 1, -99, &
         'ave(X)', t_ops, t_wrt)

    !  Pression au sol

    call histdef(fileid, 'phis', 'geopotentiel au sol', '-', &
         iip1, jjp1, thoriid, 1, 1, 1, -99, &
         'ave(X)', t_ops, t_wrt)

    call histend(fileid)

  end subroutine initdynav

end module initdynav_m
1	guez	3	module initdynav_m
2
3			! This module is clean: no C preprocessor directive, no include line
4
5			implicit none
6
7			contains
8
9	guez	27	subroutine initdynav(day0, anne0, tstep, nq, fileid, t_ops, t_wrt)
10	guez	3
11	guez	27	! From initdynav.F, version 1.1.1.1, 2004/05/19 12:53:05
12			! L. Fairhead, LMD
13	guez	3
14	guez	30	USE histcom, ONLY: histbeg_totreg, histdef, histend, histvert
15			use calendar, ONLY: ymds2ju
16	guez	27	USE dimens_m, ONLY : llm
17			USE paramet_m, ONLY : iip1, jjp1
18			USE comconst, ONLY : pi
19			USE comvert, ONLY : nivsigs
20			USE comgeom, ONLY : rlatu, rlonv
21			USE temps, ONLY : itau_dyn
22			USE iniadvtrac_m, ONLY : ttext
23	guez	3
24			! Routine d'initialisation des ecritures des fichiers histoires LMDZ
25			! au format IOIPSL. Initialisation du fichier histoire moyenne.
26
27			! Appels succesifs des routines: histbeg
28			! histhori
29			! histver
30			! histdef
31			! histend
32
33			! Entree:
34			! day0,anne0: date de reference
35			! tstep : frequence d'ecriture
36			! t_ops: frequence de l'operation pour IOIPSL
37			! t_wrt: frequence d'ecriture sur le fichier
38			! nq: nombre de traceurs
39
40			! Sortie:
41			! fileid: ID du fichier netcdf cree
42
43			! Arguments
44			integer day0, anne0
45			real, intent(in):: tstep, t_ops, t_wrt
46			integer fileid
47			integer nq
48			integer thoriid, zvertiid
49
50			! Variables locales
51
52			real zjulian
53			integer iq
54			real rlong(iip1,jjp1), rlat(iip1,jjp1)
55			integer ii,jj
56			integer zan, dayref
57
58			!----------------------------------------------------
59
60			! Appel a histbeg: creation du fichier netcdf et initialisations diverses
61
62			zan = anne0
63			dayref = day0
64			CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
65
66			do jj = 1, jjp1
67			do ii = 1, iip1
68			rlong(ii,jj) = rlonv(ii) * 180. / pi
69			rlat(ii,jj) = rlatu(jj) * 180. / pi
70			enddo
71			enddo
72
73	guez	27	call histbeg_totreg('dyn_hist_ave.nc', rlong(:,1), rlat(1,:), &
74	guez	3	1, iip1, 1, jjp1, &
75	guez	27	itau_dyn, zjulian, tstep, thoriid, fileid)
76	guez	3
77	guez	27
78	guez	3	! Appel a histvert pour la grille verticale
79	guez	27
80	guez	3	call histvert(fileid, 'sigss', 'Niveaux sigma','Pa', &
81			llm, nivsigs, zvertiid)
82	guez	27
83	guez	3	! Appels a histdef pour la definition des variables a sauvegarder
84	guez	27
85	guez	3	! Vents U
86	guez	27
87	guez	3	write(6,*)'inithistave',tstep
88			call histdef(fileid, 'u', 'vents u scalaires moyennes', &
89			'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
90	guez	15	'ave(X)', t_ops, t_wrt)
91	guez	3
92	guez	27
93	guez	3	! Vents V
94	guez	27
95	guez	3	call histdef(fileid, 'v', 'vents v scalaires moyennes', &
96			'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
97	guez	15	'ave(X)', t_ops, t_wrt)
98	guez	3
99	guez	27
100	guez	3	! Temperature
101	guez	27
102	guez	3	call histdef(fileid, 'temp', 'temperature moyennee', 'K', &
103			iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
104	guez	15	'ave(X)', t_ops, t_wrt)
105	guez	27
106	guez	3	! Temperature potentielle
107	guez	27
108	guez	3	call histdef(fileid, 'theta', 'temperature potentielle', 'K', &
109			iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
110	guez	15	'ave(X)', t_ops, t_wrt)
111	guez	3
112
113	guez	27
114	guez	3	! Geopotentiel
115	guez	27
116	guez	3	call histdef(fileid, 'phi', 'geopotentiel moyenne', '-', &
117			iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
118	guez	15	'ave(X)', t_ops, t_wrt)
119	guez	27
120	guez	3	! Traceurs
121	guez	27
122	guez	3	DO iq=1,nq
123			call histdef(fileid, ttext(iq), ttext(iq), '-', &
124			iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
125	guez	15	'ave(X)', t_ops, t_wrt)
126	guez	3	enddo
127	guez	27
128	guez	3	! Masse
129	guez	27
130	guez	3	call histdef(fileid, 'masse', 'masse', 'kg', &
131			iip1, jjp1, thoriid, 1, 1, 1, -99, &
132	guez	15	'ave(X)', t_ops, t_wrt)
133	guez	27
134	guez	3	! Pression au sol
135	guez	27
136	guez	3	call histdef(fileid, 'ps', 'pression naturelle au sol', 'Pa', &
137			iip1, jjp1, thoriid, 1, 1, 1, -99, &
138	guez	15	'ave(X)', t_ops, t_wrt)
139	guez	27
140	guez	3	! Pression au sol
141	guez	27
142	guez	3	call histdef(fileid, 'phis', 'geopotentiel au sol', '-', &
143			iip1, jjp1, thoriid, 1, 1, 1, -99, &
144	guez	15	'ave(X)', t_ops, t_wrt)
145	guez	27
146	guez	3	call histend(fileid)
147
148			end subroutine initdynav
149
150			end module initdynav_m