trunk/dyn3d/dynredem0.f

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 : ju2ymds, ymds2ju
    USE dimens_m, ONLY : iim, jjm, llm, nqmx
    USE paramet_m, ONLY : iip1, jjp1, llmp1
    USE comconst, ONLY : cpp, daysec, dtvr, g, kappa, omeg, rad
    USE comvert, ONLY : ap, bp, nivsig, nivsigs, pa, preff, presnivs
    USE logic, ONLY : fxyhypb, ysinus
    USE comgeom, ONLY : aire_2d, cu_2d, cv_2d, rlatu, rlatv, rlonu, rlonv
    USE serre, ONLY : clat, clon, dzoomx, dzoomy, grossismx, grossismy, &
         taux, tauy
    USE temps, ONLY : annee_ref, day_ref, itaufin, itau_dyn
    USE ener, ONLY : ang0, etot0, ptot0, stot0, ztot0
    USE iniadvtrac_m, ONLY : tname, ttext
    USE netcdf95, ONLY : nf95_close, nf95_create, nf95_def_dim, &
         nf95_def_var, nf95_enddef, nf95_inq_varid, nf95_put_att, &
         nf95_put_var
    USE netcdf, ONLY : nf90_clobber, nf90_float, nf90_global, &
         nf90_unlimited

    CHARACTER (len=*), INTENT (IN) :: fichnom
    INTEGER, INTENT (IN) :: iday_end
    REAL, INTENT (IN) :: phis(:, :)

    !   Local:

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

    !   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.
    END DO
    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.
    END IF

    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_def_var(nid, 'rlonu', nf90_float, idim_rlonu, nvarid)
    CALL nf95_put_att(nid, nvarid, 'title', 'Longitudes des points U')

    CALL nf95_def_var(nid, 'rlatu', nf90_float, idim_rlatu, nvarid)
    CALL nf95_put_att(nid, nvarid, 'title', 'Latitudes des points U')

    CALL nf95_def_var(nid, 'rlonv', nf90_float, idim_rlonv, nvarid)
    CALL nf95_put_att(nid, nvarid, 'title', 'Longitudes des points V')

    CALL nf95_def_var(nid, 'rlatv', nf90_float, idim_rlatv, nvarid)
    CALL nf95_put_att(nid, nvarid, 'title', 'Latitudes des points V')

    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_def_var(nid, 'nivsig', nf90_float, idim_sig, nvarid)
    CALL nf95_put_att(nid, nvarid, 'title', &
         'Numero naturel des couches sigma')

    CALL nf95_def_var(nid, 'ap', nf90_float, idim_sig, nvarid)
    CALL nf95_put_att(nid, nvarid, 'title', 'Coefficient A pour hybride')

    CALL nf95_def_var(nid, 'bp', nf90_float, idim_sig, nvarid)
    CALL nf95_put_att(nid, nvarid, 'title', 'Coefficient B pour hybride')

    CALL nf95_def_var(nid, 'presnivs', nf90_float, idim_s, nvarid)

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

    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')

    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')

    ! Aire de chaque maille:

    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')

    ! Geopentiel au sol:

    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')

    ! Definir les variables pour pouvoir les enregistrer plus tard:

    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))
    END DO

    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')

    CALL nf95_enddef(nid)

    CALL nf95_inq_varid(nid, 'controle', nvarid)
    CALL nf95_put_var(nid, nvarid, tab_cntrl)

    CALL nf95_inq_varid(nid, 'rlonu', nvarid)
    CALL nf95_put_var(nid, nvarid, rlonu)

    CALL nf95_inq_varid(nid, 'rlatu', nvarid)
    CALL nf95_put_var(nid, nvarid, rlatu)

    CALL nf95_inq_varid(nid, 'rlonv', nvarid)
    CALL nf95_put_var(nid, nvarid, rlonv)

    CALL nf95_inq_varid(nid, 'rlatv', nvarid)
    CALL nf95_put_var(nid, nvarid, rlatv)

    CALL nf95_inq_varid(nid, 'nivsigs', nvarid)
    CALL nf95_put_var(nid, nvarid, nivsigs)

    CALL nf95_inq_varid(nid, 'nivsig', nvarid)
    CALL nf95_put_var(nid, nvarid, nivsig)

    CALL nf95_inq_varid(nid, 'ap', nvarid)
    CALL nf95_put_var(nid, nvarid, ap)

    CALL nf95_inq_varid(nid, 'bp', nvarid)
    CALL nf95_put_var(nid, nvarid, bp)

    CALL nf95_inq_varid(nid, 'presnivs', nvarid)
    CALL nf95_put_var(nid, nvarid, presnivs)

    CALL nf95_inq_varid(nid, 'cu', nvarid)
    CALL nf95_put_var(nid, nvarid, cu_2d)

    CALL nf95_inq_varid(nid, 'cv', nvarid)
    CALL nf95_put_var(nid, nvarid, cv_2d)

    CALL nf95_inq_varid(nid, 'aire', nvarid)
    CALL nf95_put_var(nid, nvarid, aire_2d)

    CALL nf95_inq_varid(nid, 'phisinit', nvarid)
    CALL nf95_put_var(nid, nvarid, phis)

    CALL nf95_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	MODULE dynredem0_m
2
3	IMPLICIT NONE
4
5	CONTAINS
6
7	SUBROUTINE dynredem0(fichnom, iday_end, phis)
8
9	! 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
12	USE ioipsl, ONLY : ju2ymds, ymds2ju
13	USE dimens_m, ONLY : iim, jjm, llm, nqmx
14	USE paramet_m, ONLY : iip1, jjp1, llmp1
15	USE comconst, ONLY : cpp, daysec, dtvr, g, kappa, omeg, rad
16	USE comvert, ONLY : ap, bp, nivsig, nivsigs, pa, preff, presnivs
17	USE logic, ONLY : fxyhypb, ysinus
18	USE comgeom, ONLY : aire_2d, cu_2d, cv_2d, rlatu, rlatv, rlonu, rlonv
19	USE serre, ONLY : clat, clon, dzoomx, dzoomy, grossismx, grossismy, &
20	taux, tauy
21	USE temps, ONLY : annee_ref, day_ref, itaufin, itau_dyn
22	USE ener, ONLY : ang0, etot0, ptot0, stot0, ztot0
23	USE iniadvtrac_m, ONLY : tname, ttext
24	USE netcdf95, ONLY : nf95_close, nf95_create, nf95_def_dim, &
25	nf95_def_var, nf95_enddef, nf95_inq_varid, nf95_put_att, &
26	nf95_put_var
27	USE netcdf, ONLY : nf90_clobber, nf90_float, nf90_global, &
28	nf90_unlimited
29
30	CHARACTER (len=*), INTENT (IN) :: fichnom
31	INTEGER, INTENT (IN) :: iday_end
32	REAL, INTENT (IN) :: phis(:, :)
33
34	! Local:
35
36	INTEGER :: iq, l
37	INTEGER :: length
38	PARAMETER (length=100)
39	REAL :: tab_cntrl(length) ! tableau des parametres du run
40
41	! Variables locales pour NetCDF:
42
43	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
50	REAL :: zjulian, hours
51	INTEGER :: yyears0, jjour0, mmois0
52	CHARACTER (len=30) :: unites
53
54	!-----------------------------------------------------------------------
55
56	PRINT *, 'Call sequence information: dynredem0'
57
58	CALL ymds2ju(annee_ref, 1, iday_end, 0.0, zjulian)
59	CALL ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)
60
61	DO l = 1, length
62	tab_cntrl(l) = 0.
63	END DO
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
84	! Paramètres pour le zoom :
85
86	tab_cntrl(20) = clon
87	tab_cntrl(21) = clat
88	tab_cntrl(22) = grossismx
89	tab_cntrl(23) = grossismy
90
91	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	END IF
107
108	tab_cntrl(30) = real(iday_end)
109	tab_cntrl(31) = real(itau_dyn+itaufin)
110
111	CALL nf95_create(fichnom, nf90_clobber, nid)
112	CALL nf95_put_att(nid, nf90_global, 'title', &
113	'Fichier de démarrage dynamique')
114
115	! Definir les dimensions du fichiers:
116
117	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
126	! Definir et enregistrer certains champs invariants:
127
128	CALL nf95_def_var(nid, 'controle', nf90_float, idim_index, nvarid)
129	CALL nf95_put_att(nid, nvarid, 'title', 'Parametres de controle')
130
131	CALL nf95_def_var(nid, 'rlonu', nf90_float, idim_rlonu, nvarid)
132	CALL nf95_put_att(nid, nvarid, 'title', 'Longitudes des points U')
133
134	CALL nf95_def_var(nid, 'rlatu', nf90_float, idim_rlatu, nvarid)
135	CALL nf95_put_att(nid, nvarid, 'title', 'Latitudes des points U')
136
137	CALL nf95_def_var(nid, 'rlonv', nf90_float, idim_rlonv, nvarid)
138	CALL nf95_put_att(nid, nvarid, 'title', 'Longitudes des points V')
139
140	CALL nf95_def_var(nid, 'rlatv', nf90_float, idim_rlatv, nvarid)
141	CALL nf95_put_att(nid, nvarid, 'title', 'Latitudes des points V')
142
143	CALL nf95_def_var(nid, 'nivsigs', nf90_float, idim_s, nvarid)
144	CALL nf95_put_att(nid, nvarid, 'title', 'Numero naturel des couches s')
145
146	CALL nf95_def_var(nid, 'nivsig', nf90_float, idim_sig, nvarid)
147	CALL nf95_put_att(nid, nvarid, 'title', &
148	'Numero naturel des couches sigma')
149
150	CALL nf95_def_var(nid, 'ap', nf90_float, idim_sig, nvarid)
151	CALL nf95_put_att(nid, nvarid, 'title', 'Coefficient A pour hybride')
152
153	CALL nf95_def_var(nid, 'bp', nf90_float, idim_sig, nvarid)
154	CALL nf95_put_att(nid, nvarid, 'title', 'Coefficient B pour hybride')
155
156	CALL nf95_def_var(nid, 'presnivs', nf90_float, idim_s, nvarid)
157
158	! Coefficients de passage cov. <-> contra. <--> naturel
159
160	dims2(1) = idim_rlonu
161	dims2(2) = idim_rlatu
162	CALL nf95_def_var(nid, 'cu', nf90_float, dims2, nvarid)
163	CALL nf95_put_att(nid, nvarid, 'title', 'Coefficient de passage pour U')
164
165	dims2(1) = idim_rlonv
166	dims2(2) = idim_rlatv
167	CALL nf95_def_var(nid, 'cv', nf90_float, dims2, nvarid)
168	CALL nf95_put_att(nid, nvarid, 'title', 'Coefficient de passage pour V')
169
170	! Aire de chaque maille:
171
172	dims2(1) = idim_rlonv
173	dims2(2) = idim_rlatu
174	CALL nf95_def_var(nid, 'aire', nf90_float, dims2, nvarid)
175	CALL nf95_put_att(nid, nvarid, 'title', 'Aires de chaque maille')
176
177	! Geopentiel au sol:
178
179	dims2(1) = idim_rlonv
180	dims2(2) = idim_rlatu
181	CALL nf95_def_var(nid, 'phisinit', nf90_float, dims2, nvarid)
182	CALL nf95_put_att(nid, nvarid, 'title', 'Geopotentiel au sol')
183
184	! Definir les variables pour pouvoir les enregistrer plus tard:
185
186	CALL nf95_def_var(nid, 'temps', nf90_float, idim_tim, nvarid)
187	CALL nf95_put_att(nid, nvarid, 'title', 'Temps de simulation')
188	WRITE (unites, 200) yyears0, mmois0, jjour0
189	200 FORMAT ('days since ', I4, '-', I2.2, '-', I2.2, ' 00:00:00')
190	CALL nf95_put_att(nid, nvarid, 'units', unites)
191
192
193	dims4(1) = idim_rlonu
194	dims4(2) = idim_rlatu
195	dims4(3) = idim_s
196	dims4(4) = idim_tim
197	CALL nf95_def_var(nid, 'ucov', nf90_float, dims4, nvarid)
198	CALL nf95_put_att(nid, nvarid, 'title', 'Vitesse U')
199
200	dims4(1) = idim_rlonv
201	dims4(2) = idim_rlatv
202	dims4(3) = idim_s
203	dims4(4) = idim_tim
204	CALL nf95_def_var(nid, 'vcov', nf90_float, dims4, nvarid)
205	CALL nf95_put_att(nid, nvarid, 'title', 'Vitesse V')
206
207	dims4(1) = idim_rlonv
208	dims4(2) = idim_rlatu
209	dims4(3) = idim_s
210	dims4(4) = idim_tim
211	CALL nf95_def_var(nid, 'teta', nf90_float, dims4, nvarid)
212	CALL nf95_put_att(nid, nvarid, 'title', 'Temperature')
213
214	dims4(1) = idim_rlonv
215	dims4(2) = idim_rlatu
216	dims4(3) = idim_s
217	dims4(4) = idim_tim
218	DO iq = 1, nqmx
219	CALL nf95_def_var(nid, tname(iq), nf90_float, dims4, nvarid)
220	CALL nf95_put_att(nid, nvarid, 'title', ttext(iq))
221	END DO
222
223	dims4(1) = idim_rlonv
224	dims4(2) = idim_rlatu
225	dims4(3) = idim_s
226	dims4(4) = idim_tim
227	CALL nf95_def_var(nid, 'masse', nf90_float, dims4, nvarid)
228	CALL nf95_put_att(nid, nvarid, 'title', 'C est quoi ?')
229
230	dims3(1) = idim_rlonv
231	dims3(2) = idim_rlatu
232	dims3(3) = idim_tim
233	CALL nf95_def_var(nid, 'ps', nf90_float, dims3, nvarid)
234	CALL nf95_put_att(nid, nvarid, 'title', 'Pression au sol')
235
236	CALL nf95_enddef(nid)
237
238	CALL nf95_inq_varid(nid, 'controle', nvarid)
239	CALL nf95_put_var(nid, nvarid, tab_cntrl)
240
241	CALL nf95_inq_varid(nid, 'rlonu', nvarid)
242	CALL nf95_put_var(nid, nvarid, rlonu)
243
244	CALL nf95_inq_varid(nid, 'rlatu', nvarid)
245	CALL nf95_put_var(nid, nvarid, rlatu)
246
247	CALL nf95_inq_varid(nid, 'rlonv', nvarid)
248	CALL nf95_put_var(nid, nvarid, rlonv)
249
250	CALL nf95_inq_varid(nid, 'rlatv', nvarid)
251	CALL nf95_put_var(nid, nvarid, rlatv)
252
253	CALL nf95_inq_varid(nid, 'nivsigs', nvarid)
254	CALL nf95_put_var(nid, nvarid, nivsigs)
255
256	CALL nf95_inq_varid(nid, 'nivsig', nvarid)
257	CALL nf95_put_var(nid, nvarid, nivsig)
258
259	CALL nf95_inq_varid(nid, 'ap', nvarid)
260	CALL nf95_put_var(nid, nvarid, ap)
261
262	CALL nf95_inq_varid(nid, 'bp', nvarid)
263	CALL nf95_put_var(nid, nvarid, bp)
264
265	CALL nf95_inq_varid(nid, 'presnivs', nvarid)
266	CALL nf95_put_var(nid, nvarid, presnivs)
267
268	CALL nf95_inq_varid(nid, 'cu', nvarid)
269	CALL nf95_put_var(nid, nvarid, cu_2d)
270
271	CALL nf95_inq_varid(nid, 'cv', nvarid)
272	CALL nf95_put_var(nid, nvarid, cv_2d)
273
274	CALL nf95_inq_varid(nid, 'aire', nvarid)
275	CALL nf95_put_var(nid, nvarid, aire_2d)
276
277	CALL nf95_inq_varid(nid, 'phisinit', nvarid)
278	CALL nf95_put_var(nid, nvarid, phis)
279
280	CALL nf95_close(nid) ! fermer le fichier
281
282	PRINT *, 'iim, jjm, llm, iday_end', iim, jjm, llm, iday_end
283	PRINT *, 'rad, omeg, g, cpp, kappa', rad, omeg, g, cpp, kappa
284
285	END SUBROUTINE dynredem0
286
287	END MODULE dynredem0_m