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, infile, t_ops, t_wrt)

    ! From initdynav.F,v 1.1.1.1 2004/05/19 12:53:05

    USE IOIPSL, only: ymds2ju, histbeg_totreg, histvert, histdef, histend
    use dimens_m
    use paramet_m
    use comconst, only: pi
    use comvert, only: nivsigs
    use logic
    use comgeom
    use serre
    use temps
    use ener
    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:
    !      infile: nom du fichier histoire a creer
    !      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

    !   L. Fairhead, LMD, 03/99

    !   Arguments
    character(len=*) infile
    integer day0, anne0
    real, intent(in):: tstep, t_ops, t_wrt
    integer fileid
    integer nq
    integer thoriid, zvertiid

    !   Variables locales

    integer tau0
    real zjulian
    integer iq
    real rlong(iip1,jjp1), rlat(iip1,jjp1)
    integer ii,jj
    integer zan, dayref

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

    !  Initialisations

    pi = 4. * atan (1.)
    !
    !  Appel a histbeg: creation du fichier netcdf et initialisations diverses
    !         

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

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