libf/dyn3d/dynredem0.f90

module dynredem0_m

  IMPLICIT NONE

contains

  SUBROUTINE dynredem0(fichnom, iday_end, phis)

    ! From dyn3d/dynredem.F, v 1.2 2004/06/22 11:45:30
    ! Ecriture du fichier de redémarrage au format NetCDF (initialisation)

    USE IOIPSL, only: ymds2ju, ju2ymds
    use dimens_m, only: iim, jjm, llm, nqmx
    use paramet_m, only: iip1, jjp1, llmp1
    use comconst, only: rad, cpp, daysec, dtvr, kappa, g, omeg
    use comvert, only: pa, bp, ap, nivsigs, preff, presnivs, nivsig
    use logic
    use comgeom
    use serre
    use temps, only: annee_ref, day_ref, itaufin, itau_dyn
    use ener
    use advtrac_m, only: tname, ttext
    use netcdf95, only: nf95_create, nf95_put_att, nf95_def_dim, &
         nf95_def_var, NF95_ENDDEF, NF95_PUT_VAR
    use netcdf, only: NF90_CLOBBER, NF90_GLOBAL, NF90_UNLIMITED, nf90_float

    CHARACTER(len=*), intent(in):: fichnom
    INTEGER, intent(in):: iday_end
    REAL, intent(in):: phis(:, :)

    !   Local:

    include "netcdf.inc"

    INTEGER iq, l
    INTEGER length
    PARAMETER (length = 100)
    REAL tab_cntrl(length) ! tableau des parametres du run
    INTEGER ierr

    !   Variables locales pour NetCDF:

    INTEGER dims2(2), dims3(3), dims4(4)
    INTEGER idim_index
    INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
    INTEGER idim_s, idim_sig
    INTEGER idim_tim
    INTEGER nid, nvarid

    REAL zjulian, hours
    INTEGER yyears0, jjour0, mmois0
    character(len=30) unites

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

    print *, "Call sequence information: dynredem0"

    call ymds2ju(annee_ref, 1, iday_end, 0.0, zjulian)
    call ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)

    DO l=1, length
       tab_cntrl(l) = 0.
    ENDDO
    tab_cntrl(1)  = REAL(iim)
    tab_cntrl(2)  = REAL(jjm)
    tab_cntrl(3)  = REAL(llm)
    tab_cntrl(4)  = REAL(day_ref)
    tab_cntrl(5)  = REAL(annee_ref)
    tab_cntrl(6)  = rad
    tab_cntrl(7)  = omeg
    tab_cntrl(8)  = g
    tab_cntrl(9)  = cpp
    tab_cntrl(10) = kappa
    tab_cntrl(11) = daysec
    tab_cntrl(12) = dtvr
    tab_cntrl(13) = etot0
    tab_cntrl(14) = ptot0
    tab_cntrl(15) = ztot0
    tab_cntrl(16) = stot0
    tab_cntrl(17) = ang0
    tab_cntrl(18) = pa
    tab_cntrl(19) = preff

    ! Paramètres  pour le zoom :

    tab_cntrl(20)  = clon
    tab_cntrl(21)  = clat
    tab_cntrl(22)  = grossismx
    tab_cntrl(23)  = grossismy

    IF ( fxyhypb )   THEN
       tab_cntrl(24) = 1.
       tab_cntrl(25) = dzoomx
       tab_cntrl(26) = dzoomy
       tab_cntrl(27) = 0.
       tab_cntrl(28) = taux
       tab_cntrl(29) = tauy
    ELSE
       tab_cntrl(24) = 0.
       tab_cntrl(25) = dzoomx
       tab_cntrl(26) = dzoomy
       tab_cntrl(27) = 0.
       tab_cntrl(28) = 0.
       tab_cntrl(29) = 0.
       IF( ysinus )  tab_cntrl(27) = 1.
    ENDIF

    tab_cntrl(30) = REAL(iday_end)
    tab_cntrl(31) = REAL(itau_dyn + itaufin)

    call nf95_create(fichnom, NF90_CLOBBER, nid)
    call nf95_put_att(nid, NF90_GLOBAL, "title", &
         "Fichier de démarrage dynamique")

    ! Definir les dimensions du fichiers:

    call nf95_def_dim(nid, "index", length, idim_index)
    call NF95_DEF_DIM(nid, "rlonu", iip1, idim_rlonu)
    call NF95_DEF_DIM(nid, "rlatu", jjp1, idim_rlatu)
    call NF95_DEF_DIM(nid, "rlonv", iip1, idim_rlonv)
    call NF95_DEF_DIM(nid, "rlatv", jjm, idim_rlatv)
    call NF95_DEF_DIM(nid, "sigs", llm, idim_s)
    call NF95_DEF_DIM(nid, "sig", llmp1, idim_sig)
    call NF95_DEF_DIM(nid, "temps", NF90_UNLIMITED, idim_tim)

    ! Definir et enregistrer certains champs invariants:

    call nf95_def_var(nid, "controle", NF90_FLOAT, idim_index, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Parametres de controle")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, tab_cntrl)

    ierr = NF_REDEF (nid)
    call nf95_def_var(nid, "rlonu", NF90_FLOAT, idim_rlonu, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Longitudes des points U")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, rlonu)

    ierr = NF_REDEF (nid)
    call nf95_def_var(nid, "rlatu", NF90_FLOAT, idim_rlatu, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Latitudes des points U")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, rlatu)

    ierr = NF_REDEF (nid)
    call nf95_def_var(nid, "rlonv", NF90_FLOAT, idim_rlonv, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Longitudes des points V")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, rlonv)

    ierr = NF_REDEF (nid)
    call nf95_def_var(nid, "rlatv", NF90_FLOAT, idim_rlatv, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Latitudes des points V")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, rlatv)

    ierr = NF_REDEF (nid)
    call nf95_def_var(nid, "nivsigs", NF90_FLOAT, idim_s, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches s")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, nivsigs)

    ierr = NF_REDEF (nid)
    call nf95_def_var(nid, "nivsig", NF90_FLOAT, idim_sig, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches sigma")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, nivsig)

    ierr = NF_REDEF (nid)
    call nf95_def_var(nid, "ap", NF90_FLOAT, idim_sig, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Coefficient A pour hybride")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, ap)

    ierr = NF_REDEF (nid)
    call nf95_def_var(nid, "bp", NF90_FLOAT, idim_sig, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Coefficient B pour hybride")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, bp)

    ierr = NF_REDEF (nid)
    call nf95_def_var(nid, "presnivs", NF90_FLOAT, idim_s, nvarid)
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, presnivs)

    ! Coefficients de passage cov. <-> contra. <--> naturel

    ierr = NF_REDEF (nid)
    dims2(1) = idim_rlonu
    dims2(2) = idim_rlatu
    call nf95_def_var(nid, "cu", NF90_FLOAT, dims2, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour U")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, cu_2d)

    ierr = NF_REDEF (nid)
    dims2(1) = idim_rlonv
    dims2(2) = idim_rlatv
    call nf95_def_var(nid, "cv", NF90_FLOAT, dims2, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour V")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, cv_2d)

    ! Aire de chaque maille:

    ierr = NF_REDEF (nid)
    dims2(1) = idim_rlonv
    dims2(2) = idim_rlatu
    call nf95_def_var(nid, "aire", NF90_FLOAT, dims2, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Aires de chaque maille")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, aire_2d)

    ! Geopentiel au sol:

    ierr = NF_REDEF (nid)
    dims2(1) = idim_rlonv
    dims2(2) = idim_rlatu
    call nf95_def_var(nid, "phisinit", NF90_FLOAT, dims2, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Geopotentiel au sol")
    call NF95_ENDDEF(nid)
    call NF95_PUT_VAR(nid, nvarid, phis)

    ! Definir les variables pour pouvoir les enregistrer plus tard:

    ierr = NF_REDEF (nid) ! entrer dans le mode de definition

    call nf95_def_var(nid, "temps", NF90_FLOAT, idim_tim, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Temps de simulation")
    write(unites, 200)yyears0, mmois0, jjour0
200 format('days since ', i4, '-', i2.2, '-', i2.2, ' 00:00:00')
    call nf95_put_att(nid, nvarid, "units", unites)


    dims4(1) = idim_rlonu
    dims4(2) = idim_rlatu
    dims4(3) = idim_s
    dims4(4) = idim_tim
    call nf95_def_var(nid, "ucov", NF90_FLOAT, dims4, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Vitesse U")

    dims4(1) = idim_rlonv
    dims4(2) = idim_rlatv
    dims4(3) = idim_s
    dims4(4) = idim_tim
    call nf95_def_var(nid, "vcov", NF90_FLOAT, dims4, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Vitesse V")

    dims4(1) = idim_rlonv
    dims4(2) = idim_rlatu
    dims4(3) = idim_s
    dims4(4) = idim_tim
    call nf95_def_var(nid, "teta", NF90_FLOAT, dims4, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Temperature")

    dims4(1) = idim_rlonv
    dims4(2) = idim_rlatu
    dims4(3) = idim_s
    dims4(4) = idim_tim
    DO iq=1, nqmx
       call nf95_def_var(nid, tname(iq), NF90_FLOAT, dims4, nvarid)
       call nf95_put_att(nid, nvarid, "title", ttext(iq))
    ENDDO

    dims4(1) = idim_rlonv
    dims4(2) = idim_rlatu
    dims4(3) = idim_s
    dims4(4) = idim_tim
    call nf95_def_var(nid, "masse", NF90_FLOAT, dims4, nvarid)
    call nf95_put_att(nid, nvarid, "title", "C est quoi ?")

    dims3(1) = idim_rlonv
    dims3(2) = idim_rlatu
    dims3(3) = idim_tim
    call nf95_def_var(nid, "ps", NF90_FLOAT, dims3, nvarid)
    call nf95_put_att(nid, nvarid, "title", "Pression au sol")

    ierr = NF_ENDDEF(nid) ! sortir du mode de definition
    ierr = NF_CLOSE(nid) ! fermer le fichier

    PRINT*, 'iim, jjm, llm, iday_end', iim, jjm, llm, iday_end
    PRINT*, 'rad, omeg, g, cpp, kappa', rad, omeg, g, cpp, kappa

  END SUBROUTINE dynredem0

end module dynredem0_m
1	guez	7	module dynredem0_m
2	guez	3
3	guez	7	IMPLICIT NONE
4	guez	3
5	guez	7	contains
6	guez	3
7	guez	7	SUBROUTINE dynredem0(fichnom, iday_end, phis)
8	guez	3
9	guez	7	! From dyn3d/dynredem.F, v 1.2 2004/06/22 11:45:30
10			! Ecriture du fichier de redémarrage au format NetCDF (initialisation)
11	guez	3
12	guez	7	USE IOIPSL, only: ymds2ju, ju2ymds
13			use dimens_m, only: iim, jjm, llm, nqmx
14	guez	10	use paramet_m, only: iip1, jjp1, llmp1
15	guez	7	use comconst, only: rad, cpp, daysec, dtvr, kappa, g, omeg
16			use comvert, only: pa, bp, ap, nivsigs, preff, presnivs, nivsig
17			use logic
18			use comgeom
19			use serre
20			use temps, only: annee_ref, day_ref, itaufin, itau_dyn
21			use ener
22			use advtrac_m, only: tname, ttext
23			use netcdf95, only: nf95_create, nf95_put_att, nf95_def_dim, &
24			nf95_def_var, NF95_ENDDEF, NF95_PUT_VAR
25			use netcdf, only: NF90_CLOBBER, NF90_GLOBAL, NF90_UNLIMITED, nf90_float
26	guez	5
27	guez	7	CHARACTER(len=*), intent(in):: fichnom
28			INTEGER, intent(in):: iday_end
29			REAL, intent(in):: phis(:, :)
30	guez	5
31	guez	7	! Local:
32	guez	3
33	guez	7	include "netcdf.inc"
34	guez	5
35	guez	7	INTEGER iq, l
36			INTEGER length
37			PARAMETER (length = 100)
38			REAL tab_cntrl(length) ! tableau des parametres du run
39			INTEGER ierr
40	guez	3
41	guez	7	! Variables locales pour NetCDF:
42	guez	3
43	guez	7	INTEGER dims2(2), dims3(3), dims4(4)
44			INTEGER idim_index
45			INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
46			INTEGER idim_s, idim_sig
47			INTEGER idim_tim
48			INTEGER nid, nvarid
49	guez	3
50	guez	7	REAL zjulian, hours
51			INTEGER yyears0, jjour0, mmois0
52			character(len=30) unites
53	guez	3
54	guez	7	!-----------------------------------------------------------------------
55	guez	3
56	guez	7	print *, "Call sequence information: dynredem0"
57	guez	3
58	guez	7	call ymds2ju(annee_ref, 1, iday_end, 0.0, zjulian)
59			call ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)
60	guez	5
61	guez	7	DO l=1, length
62			tab_cntrl(l) = 0.
63			ENDDO
64			tab_cntrl(1) = REAL(iim)
65			tab_cntrl(2) = REAL(jjm)
66			tab_cntrl(3) = REAL(llm)
67			tab_cntrl(4) = REAL(day_ref)
68			tab_cntrl(5) = REAL(annee_ref)
69			tab_cntrl(6) = rad
70			tab_cntrl(7) = omeg
71			tab_cntrl(8) = g
72			tab_cntrl(9) = cpp
73			tab_cntrl(10) = kappa
74			tab_cntrl(11) = daysec
75			tab_cntrl(12) = dtvr
76			tab_cntrl(13) = etot0
77			tab_cntrl(14) = ptot0
78			tab_cntrl(15) = ztot0
79			tab_cntrl(16) = stot0
80			tab_cntrl(17) = ang0
81			tab_cntrl(18) = pa
82			tab_cntrl(19) = preff
83	guez	5
84	guez	7	! Paramètres pour le zoom :
85	guez	3
86	guez	7	tab_cntrl(20) = clon
87			tab_cntrl(21) = clat
88			tab_cntrl(22) = grossismx
89			tab_cntrl(23) = grossismy
90	guez	5
91	guez	7	IF ( fxyhypb ) THEN
92			tab_cntrl(24) = 1.
93			tab_cntrl(25) = dzoomx
94			tab_cntrl(26) = dzoomy
95			tab_cntrl(27) = 0.
96			tab_cntrl(28) = taux
97			tab_cntrl(29) = tauy
98			ELSE
99			tab_cntrl(24) = 0.
100			tab_cntrl(25) = dzoomx
101			tab_cntrl(26) = dzoomy
102			tab_cntrl(27) = 0.
103			tab_cntrl(28) = 0.
104			tab_cntrl(29) = 0.
105			IF( ysinus ) tab_cntrl(27) = 1.
106			ENDIF
107	guez	5
108	guez	7	tab_cntrl(30) = REAL(iday_end)
109			tab_cntrl(31) = REAL(itau_dyn + itaufin)
110	guez	5
111	guez	7	call nf95_create(fichnom, NF90_CLOBBER, nid)
112			call nf95_put_att(nid, NF90_GLOBAL, "title", &
113			"Fichier de démarrage dynamique")
114	guez	5
115	guez	7	! Definir les dimensions du fichiers:
116	guez	5
117	guez	7	call nf95_def_dim(nid, "index", length, idim_index)
118			call NF95_DEF_DIM(nid, "rlonu", iip1, idim_rlonu)
119			call NF95_DEF_DIM(nid, "rlatu", jjp1, idim_rlatu)
120			call NF95_DEF_DIM(nid, "rlonv", iip1, idim_rlonv)
121			call NF95_DEF_DIM(nid, "rlatv", jjm, idim_rlatv)
122			call NF95_DEF_DIM(nid, "sigs", llm, idim_s)
123			call NF95_DEF_DIM(nid, "sig", llmp1, idim_sig)
124			call NF95_DEF_DIM(nid, "temps", NF90_UNLIMITED, idim_tim)
125	guez	5
126	guez	7	! Definir et enregistrer certains champs invariants:
127	guez	5
128	guez	7	call nf95_def_var(nid, "controle", NF90_FLOAT, idim_index, nvarid)
129			call nf95_put_att(nid, nvarid, "title", "Parametres de controle")
130			call NF95_ENDDEF(nid)
131			call NF95_PUT_VAR(nid, nvarid, tab_cntrl)
132	guez	5
133	guez	7	ierr = NF_REDEF (nid)
134			call nf95_def_var(nid, "rlonu", NF90_FLOAT, idim_rlonu, nvarid)
135			call nf95_put_att(nid, nvarid, "title", "Longitudes des points U")
136			call NF95_ENDDEF(nid)
137			call NF95_PUT_VAR(nid, nvarid, rlonu)
138	guez	5
139	guez	7	ierr = NF_REDEF (nid)
140			call nf95_def_var(nid, "rlatu", NF90_FLOAT, idim_rlatu, nvarid)
141			call nf95_put_att(nid, nvarid, "title", "Latitudes des points U")
142			call NF95_ENDDEF(nid)
143			call NF95_PUT_VAR(nid, nvarid, rlatu)
144	guez	5
145	guez	7	ierr = NF_REDEF (nid)
146			call nf95_def_var(nid, "rlonv", NF90_FLOAT, idim_rlonv, nvarid)
147			call nf95_put_att(nid, nvarid, "title", "Longitudes des points V")
148			call NF95_ENDDEF(nid)
149			call NF95_PUT_VAR(nid, nvarid, rlonv)
150	guez	5
151	guez	7	ierr = NF_REDEF (nid)
152			call nf95_def_var(nid, "rlatv", NF90_FLOAT, idim_rlatv, nvarid)
153			call nf95_put_att(nid, nvarid, "title", "Latitudes des points V")
154			call NF95_ENDDEF(nid)
155			call NF95_PUT_VAR(nid, nvarid, rlatv)
156	guez	5
157	guez	7	ierr = NF_REDEF (nid)
158			call nf95_def_var(nid, "nivsigs", NF90_FLOAT, idim_s, nvarid)
159			call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches s")
160			call NF95_ENDDEF(nid)
161			call NF95_PUT_VAR(nid, nvarid, nivsigs)
162	guez	5
163	guez	7	ierr = NF_REDEF (nid)
164			call nf95_def_var(nid, "nivsig", NF90_FLOAT, idim_sig, nvarid)
165			call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches sigma")
166			call NF95_ENDDEF(nid)
167			call NF95_PUT_VAR(nid, nvarid, nivsig)
168	guez	5
169	guez	7	ierr = NF_REDEF (nid)
170			call nf95_def_var(nid, "ap", NF90_FLOAT, idim_sig, nvarid)
171			call nf95_put_att(nid, nvarid, "title", "Coefficient A pour hybride")
172			call NF95_ENDDEF(nid)
173			call NF95_PUT_VAR(nid, nvarid, ap)
174	guez	5
175	guez	7	ierr = NF_REDEF (nid)
176			call nf95_def_var(nid, "bp", NF90_FLOAT, idim_sig, nvarid)
177			call nf95_put_att(nid, nvarid, "title", "Coefficient B pour hybride")
178			call NF95_ENDDEF(nid)
179			call NF95_PUT_VAR(nid, nvarid, bp)
180	guez	5
181	guez	7	ierr = NF_REDEF (nid)
182			call nf95_def_var(nid, "presnivs", NF90_FLOAT, idim_s, nvarid)
183			call NF95_ENDDEF(nid)
184			call NF95_PUT_VAR(nid, nvarid, presnivs)
185	guez	5
186	guez	7	! Coefficients de passage cov. <-> contra. <--> naturel
187	guez	5
188	guez	7	ierr = NF_REDEF (nid)
189			dims2(1) = idim_rlonu
190			dims2(2) = idim_rlatu
191			call nf95_def_var(nid, "cu", NF90_FLOAT, dims2, nvarid)
192			call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour U")
193			call NF95_ENDDEF(nid)
194			call NF95_PUT_VAR(nid, nvarid, cu_2d)
195	guez	5
196	guez	7	ierr = NF_REDEF (nid)
197			dims2(1) = idim_rlonv
198			dims2(2) = idim_rlatv
199			call nf95_def_var(nid, "cv", NF90_FLOAT, dims2, nvarid)
200			call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour V")
201			call NF95_ENDDEF(nid)
202			call NF95_PUT_VAR(nid, nvarid, cv_2d)
203	guez	5
204	guez	7	! Aire de chaque maille:
205	guez	5
206	guez	7	ierr = NF_REDEF (nid)
207			dims2(1) = idim_rlonv
208			dims2(2) = idim_rlatu
209			call nf95_def_var(nid, "aire", NF90_FLOAT, dims2, nvarid)
210			call nf95_put_att(nid, nvarid, "title", "Aires de chaque maille")
211			call NF95_ENDDEF(nid)
212			call NF95_PUT_VAR(nid, nvarid, aire_2d)
213	guez	5
214	guez	7	! Geopentiel au sol:
215	guez	5
216	guez	7	ierr = NF_REDEF (nid)
217			dims2(1) = idim_rlonv
218			dims2(2) = idim_rlatu
219			call nf95_def_var(nid, "phisinit", NF90_FLOAT, dims2, nvarid)
220			call nf95_put_att(nid, nvarid, "title", "Geopotentiel au sol")
221			call NF95_ENDDEF(nid)
222			call NF95_PUT_VAR(nid, nvarid, phis)
223	guez	5
224	guez	7	! Definir les variables pour pouvoir les enregistrer plus tard:
225	guez	3
226	guez	7	ierr = NF_REDEF (nid) ! entrer dans le mode de definition
227
228			call nf95_def_var(nid, "temps", NF90_FLOAT, idim_tim, nvarid)
229			call nf95_put_att(nid, nvarid, "title", "Temps de simulation")
230			write(unites, 200)yyears0, mmois0, jjour0
231	guez	5	200 format('days since ', i4, '-', i2.2, '-', i2.2, ' 00:00:00')
232	guez	7	call nf95_put_att(nid, nvarid, "units", unites)
233	guez	5
234
235	guez	7	dims4(1) = idim_rlonu
236			dims4(2) = idim_rlatu
237			dims4(3) = idim_s
238			dims4(4) = idim_tim
239			call nf95_def_var(nid, "ucov", NF90_FLOAT, dims4, nvarid)
240			call nf95_put_att(nid, nvarid, "title", "Vitesse U")
241	guez	5
242	guez	7	dims4(1) = idim_rlonv
243			dims4(2) = idim_rlatv
244			dims4(3) = idim_s
245			dims4(4) = idim_tim
246			call nf95_def_var(nid, "vcov", NF90_FLOAT, dims4, nvarid)
247			call nf95_put_att(nid, nvarid, "title", "Vitesse V")
248	guez	5
249	guez	7	dims4(1) = idim_rlonv
250			dims4(2) = idim_rlatu
251			dims4(3) = idim_s
252			dims4(4) = idim_tim
253			call nf95_def_var(nid, "teta", NF90_FLOAT, dims4, nvarid)
254			call nf95_put_att(nid, nvarid, "title", "Temperature")
255	guez	5
256	guez	7	dims4(1) = idim_rlonv
257			dims4(2) = idim_rlatu
258			dims4(3) = idim_s
259			dims4(4) = idim_tim
260			DO iq=1, nqmx
261			call nf95_def_var(nid, tname(iq), NF90_FLOAT, dims4, nvarid)
262			call nf95_put_att(nid, nvarid, "title", ttext(iq))
263			ENDDO
264	guez	5
265	guez	7	dims4(1) = idim_rlonv
266			dims4(2) = idim_rlatu
267			dims4(3) = idim_s
268			dims4(4) = idim_tim
269			call nf95_def_var(nid, "masse", NF90_FLOAT, dims4, nvarid)
270			call nf95_put_att(nid, nvarid, "title", "C est quoi ?")
271	guez	5
272	guez	7	dims3(1) = idim_rlonv
273			dims3(2) = idim_rlatu
274			dims3(3) = idim_tim
275			call nf95_def_var(nid, "ps", NF90_FLOAT, dims3, nvarid)
276			call nf95_put_att(nid, nvarid, "title", "Pression au sol")
277	guez	5
278	guez	7	ierr = NF_ENDDEF(nid) ! sortir du mode de definition
279			ierr = NF_CLOSE(nid) ! fermer le fichier
280	guez	3
281	guez	7	PRINT*, 'iim, jjm, llm, iday_end', iim, jjm, llm, iday_end
282			PRINT*, 'rad, omeg, g, cpp, kappa', rad, omeg, g, cpp, kappa
283	guez	3
284	guez	7	END SUBROUTINE dynredem0
285
286			end module dynredem0_m