Sources/dyn3d/initdynav.f

module initdynav_m

  implicit none

  integer histaveid

contains

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

    ! From initdynav.F, version 1.1.1.1, 2004/05/19 12:53:05
    ! L. Fairhead, LMD
    ! Routine d'initialisation des écritures des fichiers histoires LMDZ
    ! au format IOIPSL. Initialisation du fichier histoire moyenne.

    use calendar, ONLY: ymds2ju
    USE comvert, ONLY : nivsigs
    USE comgeom, ONLY : rlatu, rlonv
    USE dimens_m, ONLY : llm
    USE histbeg_totreg_m, ONLY : histbeg_totreg
    USE histdef_m, ONLY : histdef
    USE histend_m, ONLY : histend
    USE histvert_m, ONLY : histvert
    USE iniadvtrac_m, ONLY : ttext
    USE nr_util, ONLY : pi
    USE paramet_m, ONLY : iip1, jjp1
    USE temps, ONLY : itau_dyn

    integer, intent(in):: day0, anne0 ! date de reference
    real, intent(in):: tstep ! frequence d'ecriture
    integer, intent(in):: nq ! nombre de traceurs
    real, intent(in):: t_ops ! frequence de l'operation pour IOIPSL
    real, intent(in):: t_wrt ! frequence d'ecriture sur le fichier

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

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

    print *, "Call sequence information: initdynav"

    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, histaveid)
    call histvert(histaveid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &
         zvertiid)

    call histdef(histaveid, 'u', 'vents u scalaires moyennes', &
         'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)
    call histdef(histaveid, 'v', 'vents v scalaires moyennes', &
         'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)
    call histdef(histaveid, 'temp', 'temperature moyennee', 'K', &
         iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)
    call histdef(histaveid, 'theta', 'temperature potentielle', 'K', &
         iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)
    call histdef(histaveid, 'phi', 'geopotentiel moyenne', '-', &
         iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', t_ops, t_wrt)

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

    call histdef(histaveid, 'masse', 'masse', 'kg', &
         iip1, jjp1, thoriid, 1, 1, 1, -99, &
         'ave(X)', t_ops, t_wrt)
    call histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', &
         iip1, jjp1, thoriid, 1, 1, 1, -99, &
         'ave(X)', t_ops, t_wrt)
    call histdef(histaveid, 'phis', 'geopotentiel au sol', '-', &
         iip1, jjp1, thoriid, 1, 1, 1, -99, &
         'ave(X)', t_ops, t_wrt)

    call histend(histaveid)

  end subroutine initdynav

end module initdynav_m
1	guez	3	module initdynav_m
2
3			implicit none
4
5	guez	56	integer histaveid
6
7	guez	3	contains
8
9	guez	56	subroutine initdynav(day0, anne0, tstep, nq, 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	guez	40	! L. Fairhead, LMD
13			! Routine d'initialisation des écritures des fichiers histoires LMDZ
14			! au format IOIPSL. Initialisation du fichier histoire moyenne.
15	guez	3
16	guez	30	use calendar, ONLY: ymds2ju
17	guez	39	USE comvert, ONLY : nivsigs
18			USE comgeom, ONLY : rlatu, rlonv
19	guez	27	USE dimens_m, ONLY : llm
20	guez	61	USE histbeg_totreg_m, ONLY : histbeg_totreg
21			USE histdef_m, ONLY : histdef
22			USE histend_m, ONLY : histend
23			USE histvert_m, ONLY : histvert
24	guez	39	USE iniadvtrac_m, ONLY : ttext
25			USE nr_util, ONLY : pi
26	guez	27	USE paramet_m, ONLY : iip1, jjp1
27			USE temps, ONLY : itau_dyn
28	guez	3
29	guez	40	integer, intent(in):: day0, anne0 ! date de reference
30			real, intent(in):: tstep ! frequence d'ecriture
31	guez	56	integer, intent(in):: nq ! nombre de traceurs
32	guez	40	real, intent(in):: t_ops ! frequence de l'operation pour IOIPSL
33			real, intent(in):: t_wrt ! frequence d'ecriture sur le fichier
34	guez	3
35	guez	40	! Variables locales
36	guez	3	integer thoriid, zvertiid
37			real zjulian
38			integer iq
39	guez	40	real rlong(iip1, jjp1), rlat(iip1, jjp1)
40			integer ii, jj
41	guez	3	integer zan, dayref
42
43			!----------------------------------------------------
44
45	guez	40	print *, "Call sequence information: initdynav"
46	guez	3
47			zan = anne0
48			dayref = day0
49			CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
50
51			do jj = 1, jjp1
52			do ii = 1, iip1
53	guez	40	rlong(ii, jj) = rlonv(ii) * 180. / pi
54			rlat(ii, jj) = rlatu(jj) * 180. / pi
55	guez	3	enddo
56			enddo
57
58	guez	40	call histbeg_totreg('dyn_hist_ave.nc', rlong(:, 1), rlat(1, :), 1, iip1, &
59	guez	56	1, jjp1, itau_dyn, zjulian, tstep, thoriid, histaveid)
60			call histvert(histaveid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &
61	guez	40	zvertiid)
62	guez	3
63	guez	56	call histdef(histaveid, 'u', 'vents u scalaires moyennes', &
64	guez	3	'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
65	guez	15	'ave(X)', t_ops, t_wrt)
66	guez	56	call histdef(histaveid, 'v', 'vents v scalaires moyennes', &
67	guez	3	'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
68	guez	15	'ave(X)', t_ops, t_wrt)
69	guez	56	call histdef(histaveid, 'temp', 'temperature moyennee', 'K', &
70	guez	3	iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
71	guez	15	'ave(X)', t_ops, t_wrt)
72	guez	56	call histdef(histaveid, 'theta', 'temperature potentielle', 'K', &
73	guez	3	iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
74	guez	15	'ave(X)', t_ops, t_wrt)
75	guez	56	call histdef(histaveid, 'phi', 'geopotentiel moyenne', '-', &
76	guez	3	iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
77	guez	15	'ave(X)', t_ops, t_wrt)
78	guez	27
79	guez	40	! Traceurs
80			DO iq = 1, nq
81	guez	56	call histdef(histaveid, ttext(iq), ttext(iq), '-', &
82	guez	3	iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
83	guez	15	'ave(X)', t_ops, t_wrt)
84	guez	3	enddo
85	guez	27
86	guez	56	call histdef(histaveid, 'masse', 'masse', 'kg', &
87	guez	3	iip1, jjp1, thoriid, 1, 1, 1, -99, &
88	guez	15	'ave(X)', t_ops, t_wrt)
89	guez	56	call histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', &
90	guez	3	iip1, jjp1, thoriid, 1, 1, 1, -99, &
91	guez	15	'ave(X)', t_ops, t_wrt)
92	guez	56	call histdef(histaveid, 'phis', 'geopotentiel au sol', '-', &
93	guez	3	iip1, jjp1, thoriid, 1, 1, 1, -99, &
94	guez	15	'ave(X)', t_ops, t_wrt)
95	guez	27
96	guez	56	call histend(histaveid)
97	guez	3
98			end subroutine initdynav
99
100			end module initdynav_m