trunk/dyn3d/init_dynzon.f

module init_dynzon_m

  IMPLICIT NONE

  integer, parameter:: ntr = 5
  integer, parameter:: nQ = 7
  integer ncum, fileid
  character(len=10) znom(ntr, nQ)
  character(len=4), parameter:: nom(nQ) = (/'T   ', 'gz  ', 'K   ', 'ang ', &
       'u   ', 'ovap', 'un  '/)

contains

  SUBROUTINE init_dynzon(dt_app)

    ! From LMDZ4/libf/dyn3d/bilan_dyn.F, version 1.5 2005/03/16 10:12:17

    USE conf_gcm_m, ONLY : day_step, iperiod, periodav
    USE histcom, ONLY : histbeg_totreg, histdef, histend, histvert
    USE calendar, ONLY : ymds2ju
    USE dimens_m, ONLY : jjm, llm
    USE comvert, ONLY : presnivs
    USE comgeom, ONLY : rlatv
    USE temps, ONLY : annee_ref, day_ref, itau_dyn
    USE nr_util, ONLY : pi

    ! Arguments:

    real, intent(in):: dt_app

    ! Local:

    real dt_cum
    character(len=5), parameter:: unites(nQ) = (/'K    ', 'm2/s2', 'm2/s2', &
         'ang  ', 'm/s  ', 'kg/kg', 'un   '/)

    ! champs de tansport en moyenne zonale
    integer itr

    character(len=26) znoml(ntr, nQ)
    character(len=12) zunites(ntr, nQ)

    character(len=3), parameter:: ctrs(ntr) = (/'   ', 'TOT', 'MMC', 'TRS', &
         'STN'/)

    integer iQ

    ! Initialisation du fichier contenant les moyennes zonales.

    integer thoriid, zvertiid

    real zjulian
    integer zan, dayref

    real rlong(jjm), rlatg(jjm)

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

    print *, "Call sequence information: init_dynzon"

    ! initialisation des fichiers
    ! ncum est la frequence de stokage en pas de temps
    ncum = day_step / iperiod * periodav
    dt_cum = ncum * dt_app

    ! Initialisation du fichier contenant les moyennes zonales

    zan = annee_ref
    dayref = day_ref
    CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)

    rlong = 0.
    rlatg = rlatv*180./pi

    call histbeg_totreg('dynzon', rlong(:1), rlatg, 1, 1, 1, jjm, itau_dyn, &
         zjulian, dt_cum, thoriid, fileid)

    ! Appel à histvert pour la grille verticale

    call histvert(fileid, 'presnivs', 'Niveaux sigma', 'mb', llm, presnivs, &
         zvertiid)

    ! Appels à histdef pour la définition des variables à sauvegarder
    do iQ = 1, nQ
       do itr = 1, ntr
          if (itr == 1) then
             znom(itr, iQ) = nom(iQ)
             znoml(itr, iQ) = nom(iQ)
             zunites(itr, iQ) = unites(iQ)
          else
             znom(itr, iQ) = ctrs(itr)//'v'//nom(iQ)
             znoml(itr, iQ) = 'transport : v * '//nom(iQ)//' '//ctrs(itr)
             zunites(itr, iQ) = 'm/s * '//unites(iQ)
          endif
       enddo
    enddo

    ! Déclarations des champs avec dimension verticale
    do iQ = 1, nQ
       do itr = 1, ntr
          call histdef(fileid, znom(itr, iQ), znoml(itr, iQ), &
               zunites(itr, iQ), 1, jjm, thoriid, llm, 1, llm, zvertiid, &
               'ave(X)', dt_cum, dt_cum)
       enddo
       ! Déclarations pour les fonctions de courant
       call histdef(fileid, 'psi'//nom(iQ), 'stream fn. '//znoml(2, iQ), &
            zunites(2, iQ), 1, jjm, thoriid, llm, 1, llm, zvertiid, &
            'ave(X)', dt_cum, dt_cum)
    enddo

    ! Déclarations pour les champs de transport d'air
    call histdef(fileid, 'masse', 'masse', &
         'kg', 1, jjm, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', dt_cum, dt_cum)
    call histdef(fileid, 'v', 'v', &
         'm/s', 1, jjm, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', dt_cum, dt_cum)
    ! Déclarations pour les fonctions de courant
    call histdef(fileid, 'psi', 'stream fn. MMC ', 'mega t/s', &
         1, jjm, thoriid, llm, 1, llm, zvertiid, &
         'ave(X)', dt_cum, dt_cum)

    ! Déclaration des champs 1D de transport en latitude
    do iQ = 1, nQ
       do itr = 2, ntr
          call histdef(fileid, 'a'//znom(itr, iQ), znoml(itr, iQ), &
               zunites(itr, iQ), 1, jjm, thoriid, 1, 1, 1, -99, &
               'ave(X)', dt_cum, dt_cum)
       enddo
    enddo

    CALL histend(fileid)

  end SUBROUTINE init_dynzon

end module init_dynzon_m
1	module init_dynzon_m
2
3	IMPLICIT NONE
4
5	integer, parameter:: ntr = 5
6	integer, parameter:: nQ = 7
7	integer ncum, fileid
8	character(len=10) znom(ntr, nQ)
9	character(len=4), parameter:: nom(nQ) = (/'T ', 'gz ', 'K ', 'ang ', &
10	'u ', 'ovap', 'un '/)
11
12	contains
13
14	SUBROUTINE init_dynzon(dt_app)
15
16	! From LMDZ4/libf/dyn3d/bilan_dyn.F, version 1.5 2005/03/16 10:12:17
17
18	USE conf_gcm_m, ONLY : day_step, iperiod, periodav
19	USE histcom, ONLY : histbeg_totreg, histdef, histend, histvert
20	USE calendar, ONLY : ymds2ju
21	USE dimens_m, ONLY : jjm, llm
22	USE comvert, ONLY : presnivs
23	USE comgeom, ONLY : rlatv
24	USE temps, ONLY : annee_ref, day_ref, itau_dyn
25	USE nr_util, ONLY : pi
26
27	! Arguments:
28
29	real, intent(in):: dt_app
30
31	! Local:
32
33	real dt_cum
34	character(len=5), parameter:: unites(nQ) = (/'K ', 'm2/s2', 'm2/s2', &
35	'ang ', 'm/s ', 'kg/kg', 'un '/)
36
37	! champs de tansport en moyenne zonale
38	integer itr
39
40	character(len=26) znoml(ntr, nQ)
41	character(len=12) zunites(ntr, nQ)
42
43	character(len=3), parameter:: ctrs(ntr) = (/' ', 'TOT', 'MMC', 'TRS', &
44	'STN'/)
45
46	integer iQ
47
48	! Initialisation du fichier contenant les moyennes zonales.
49
50	integer thoriid, zvertiid
51
52	real zjulian
53	integer zan, dayref
54
55	real rlong(jjm), rlatg(jjm)
56
57	!-----------------------------------------------------------------
58
59	print *, "Call sequence information: init_dynzon"
60
61	! initialisation des fichiers
62	! ncum est la frequence de stokage en pas de temps
63	ncum = day_step / iperiod * periodav
64	dt_cum = ncum * dt_app
65
66	! Initialisation du fichier contenant les moyennes zonales
67
68	zan = annee_ref
69	dayref = day_ref
70	CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
71
72	rlong = 0.
73	rlatg = rlatv*180./pi
74
75	call histbeg_totreg('dynzon', rlong(:1), rlatg, 1, 1, 1, jjm, itau_dyn, &
76	zjulian, dt_cum, thoriid, fileid)
77
78	! Appel à histvert pour la grille verticale
79
80	call histvert(fileid, 'presnivs', 'Niveaux sigma', 'mb', llm, presnivs, &
81	zvertiid)
82
83	! Appels à histdef pour la définition des variables à sauvegarder
84	do iQ = 1, nQ
85	do itr = 1, ntr
86	if (itr == 1) then
87	znom(itr, iQ) = nom(iQ)
88	znoml(itr, iQ) = nom(iQ)
89	zunites(itr, iQ) = unites(iQ)
90	else
91	znom(itr, iQ) = ctrs(itr)//'v'//nom(iQ)
92	znoml(itr, iQ) = 'transport : v * '//nom(iQ)//' '//ctrs(itr)
93	zunites(itr, iQ) = 'm/s * '//unites(iQ)
94	endif
95	enddo
96	enddo
97
98	! Déclarations des champs avec dimension verticale
99	do iQ = 1, nQ
100	do itr = 1, ntr
101	call histdef(fileid, znom(itr, iQ), znoml(itr, iQ), &
102	zunites(itr, iQ), 1, jjm, thoriid, llm, 1, llm, zvertiid, &
103	'ave(X)', dt_cum, dt_cum)
104	enddo
105	! Déclarations pour les fonctions de courant
106	call histdef(fileid, 'psi'//nom(iQ), 'stream fn. '//znoml(2, iQ), &
107	zunites(2, iQ), 1, jjm, thoriid, llm, 1, llm, zvertiid, &
108	'ave(X)', dt_cum, dt_cum)
109	enddo
110
111	! Déclarations pour les champs de transport d'air
112	call histdef(fileid, 'masse', 'masse', &
113	'kg', 1, jjm, thoriid, llm, 1, llm, zvertiid, &
114	'ave(X)', dt_cum, dt_cum)
115	call histdef(fileid, 'v', 'v', &
116	'm/s', 1, jjm, thoriid, llm, 1, llm, zvertiid, &
117	'ave(X)', dt_cum, dt_cum)
118	! Déclarations pour les fonctions de courant
119	call histdef(fileid, 'psi', 'stream fn. MMC ', 'mega t/s', &
120	1, jjm, thoriid, llm, 1, llm, zvertiid, &
121	'ave(X)', dt_cum, dt_cum)
122
123	! Déclaration des champs 1D de transport en latitude
124	do iQ = 1, nQ
125	do itr = 2, ntr
126	call histdef(fileid, 'a'//znom(itr, iQ), znoml(itr, iQ), &
127	zunites(itr, iQ), 1, jjm, thoriid, 1, 1, 1, -99, &
128	'ave(X)', dt_cum, dt_cum)
129	enddo
130	enddo
131
132	CALL histend(fileid)
133
134	end SUBROUTINE init_dynzon
135
136	end module init_dynzon_m