trunk/dyn3d/dynredem0.f

MODULE dynredem0_m

  IMPLICIT NONE

CONTAINS

  SUBROUTINE dynredem0(fichnom, iday_end, phis)

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

    USE comconst, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad
    USE comgeom, ONLY: aire_2d, cu_2d, cv_2d, rlatu, rlatv, rlonu, rlonv
    use conf_gcm_m, ONLY: fxyhypb, ysinus
    USE dimens_m, ONLY: iim, jjm, llm, nqmx
    USE disvert_m, ONLY: ap, bp, pa, preff, presnivs
    USE ener, ONLY: ang0, etot0, ptot0, stot0, ztot0
    USE iniadvtrac_m, ONLY: tname, ttext
    USE ju2ymds_m, ONLY: ju2ymds
    USE netcdf, ONLY: nf90_clobber, nf90_float, nf90_global, nf90_unlimited
    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 paramet_m, ONLY: iip1, jjp1, llmp1
    USE serre, ONLY: clat, clon, dzoomx, dzoomy, grossismx, grossismy, taux, &
         tauy
    USE temps, ONLY: annee_ref, day_ref
    use ymds2ju_m, only: ymds2ju

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

    ! Local:

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

    ! Variables locales pour NetCDF:

    INTEGER idim_index
    INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
    INTEGER idim_s, idim_sig
    INTEGER idim_tim
    INTEGER ncid, varid

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

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

    PRINT *, 'Call sequence information: dynredem0'

    CALL ymds2ju(annee_ref, 1, iday_end, 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)

    CALL nf95_create(fichnom, nf90_clobber, ncid)
    CALL nf95_put_att(ncid, nf90_global, 'title', &
         'Fichier de démarrage dynamique')

    ! Definir les dimensions du fichiers:

    CALL nf95_def_dim(ncid, 'index', length, idim_index)
    CALL nf95_def_dim(ncid, 'rlonu', iip1, idim_rlonu)
    CALL nf95_def_dim(ncid, 'rlatu', jjp1, idim_rlatu)
    CALL nf95_def_dim(ncid, 'rlonv', iip1, idim_rlonv)
    CALL nf95_def_dim(ncid, 'rlatv', jjm, idim_rlatv)
    CALL nf95_def_dim(ncid, 'sigs', llm, idim_s)
    CALL nf95_def_dim(ncid, 'sig', llmp1, idim_sig)
    CALL nf95_def_dim(ncid, 'temps', nf90_unlimited, idim_tim)

    ! Definir et enregistrer certains champs invariants:

    CALL nf95_def_var(ncid, 'controle', nf90_float, idim_index, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Parametres de controle')

    CALL nf95_def_var(ncid, 'rlonu', nf90_float, idim_rlonu, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Longitudes des points U')

    CALL nf95_def_var(ncid, 'rlatu', nf90_float, idim_rlatu, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Latitudes des points U')

    CALL nf95_def_var(ncid, 'rlonv', nf90_float, idim_rlonv, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Longitudes des points V')

    CALL nf95_def_var(ncid, 'rlatv', nf90_float, idim_rlatv, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Latitudes des points V')

    CALL nf95_def_var(ncid, 'ap', nf90_float, idim_sig, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Coefficient A pour hybride')

    CALL nf95_def_var(ncid, 'bp', nf90_float, idim_sig, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Coefficient B pour hybride')

    CALL nf95_def_var(ncid, 'presnivs', nf90_float, idim_s, varid)

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

    CALL nf95_def_var(ncid, 'cu', nf90_float, (/idim_rlonu, idim_rlatu/), varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Coefficient de passage pour U')

    CALL nf95_def_var(ncid, 'cv', nf90_float, (/idim_rlonv, idim_rlatv/), varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Coefficient de passage pour V')

    ! Aire de chaque maille:

    CALL nf95_def_var(ncid, 'aire', nf90_float, (/idim_rlonv, idim_rlatu/), &
         varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Aires de chaque maille')

    ! Geopentiel au sol:

    CALL nf95_def_var(ncid, 'phisinit', nf90_float, &
         (/idim_rlonv, idim_rlatu/), varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Geopotentiel au sol')

    ! Definir les variables pour pouvoir les enregistrer plus tard:

    CALL nf95_def_var(ncid, 'temps', nf90_float, idim_tim, varid)
    CALL nf95_put_att(ncid, varid, '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(ncid, varid, 'units', unites)

    CALL nf95_def_var(ncid, 'ucov', nf90_float, &
         (/idim_rlonu, idim_rlatu, idim_s, idim_tim/), varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Vitesse U')

    CALL nf95_def_var(ncid, 'vcov', nf90_float, &
         (/idim_rlonv, idim_rlatv, idim_s, idim_tim/), varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Vitesse V')

    CALL nf95_def_var(ncid, 'teta', nf90_float, &
         (/idim_rlonv, idim_rlatu, idim_s, idim_tim/), varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Temperature')

    DO iq = 1, nqmx
       CALL nf95_def_var(ncid, tname(iq), nf90_float, &
            (/idim_rlonv, idim_rlatu, idim_s, idim_tim/), varid)
       CALL nf95_put_att(ncid, varid, 'title', ttext(iq))
    END DO

    CALL nf95_def_var(ncid, 'masse', nf90_float, &
         (/idim_rlonv, idim_rlatu, idim_s, idim_tim/), varid)
    CALL nf95_put_att(ncid, varid, 'title', 'C est quoi ?')

    CALL nf95_def_var(ncid, 'ps', nf90_float, &
         (/idim_rlonv, idim_rlatu, idim_tim/), varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Pression au sol')

    CALL nf95_enddef(ncid)

    CALL nf95_inq_varid(ncid, 'controle', varid)
    CALL nf95_put_var(ncid, varid, tab_cntrl)

    CALL nf95_inq_varid(ncid, 'rlonu', varid)
    CALL nf95_put_var(ncid, varid, rlonu)

    CALL nf95_inq_varid(ncid, 'rlatu', varid)
    CALL nf95_put_var(ncid, varid, rlatu)

    CALL nf95_inq_varid(ncid, 'rlonv', varid)
    CALL nf95_put_var(ncid, varid, rlonv)

    CALL nf95_inq_varid(ncid, 'rlatv', varid)
    CALL nf95_put_var(ncid, varid, rlatv)

    CALL nf95_inq_varid(ncid, 'ap', varid)
    CALL nf95_put_var(ncid, varid, ap)

    CALL nf95_inq_varid(ncid, 'bp', varid)
    CALL nf95_put_var(ncid, varid, bp)

    CALL nf95_inq_varid(ncid, 'presnivs', varid)
    CALL nf95_put_var(ncid, varid, presnivs)

    CALL nf95_inq_varid(ncid, 'cu', varid)
    CALL nf95_put_var(ncid, varid, cu_2d)

    CALL nf95_inq_varid(ncid, 'cv', varid)
    CALL nf95_put_var(ncid, varid, cv_2d)

    CALL nf95_inq_varid(ncid, 'aire', varid)
    CALL nf95_put_var(ncid, varid, aire_2d)

    CALL nf95_inq_varid(ncid, 'phisinit', varid)
    CALL nf95_put_var(ncid, varid, phis)

    CALL nf95_close(ncid) ! 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	18	MODULE dynredem0_m
2	guez	3
3	guez	7	IMPLICIT NONE
4	guez	3
5	guez	18	CONTAINS
6	guez	3
7	guez	7	SUBROUTINE dynredem0(fichnom, iday_end, phis)
8	guez	3
9	guez	23	! From dyn3d/dynredem.F, version 1.2 2004/06/22 11:45:30
10	guez	7	! Ecriture du fichier de redémarrage au format NetCDF (initialisation)
11	guez	3
12	guez	62	USE comconst, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad
13			USE comgeom, ONLY: aire_2d, cu_2d, cv_2d, rlatu, rlatv, rlonu, rlonv
14	guez	92	use conf_gcm_m, ONLY: fxyhypb, ysinus
15	guez	62	USE dimens_m, ONLY: iim, jjm, llm, nqmx
16	guez	92	USE disvert_m, ONLY: ap, bp, pa, preff, presnivs
17	guez	62	USE ener, ONLY: ang0, etot0, ptot0, stot0, ztot0
18			USE iniadvtrac_m, ONLY: tname, ttext
19	guez	92	USE ju2ymds_m, ONLY: ju2ymds
20			USE netcdf, ONLY: nf90_clobber, nf90_float, nf90_global, nf90_unlimited
21	guez	62	USE netcdf95, ONLY: nf95_close, nf95_create, nf95_def_dim, nf95_def_var, &
22			nf95_enddef, nf95_inq_varid, nf95_put_att, nf95_put_var
23			USE paramet_m, ONLY: iip1, jjp1, llmp1
24			USE serre, ONLY: clat, clon, dzoomx, dzoomy, grossismx, grossismy, taux, &
25			tauy
26			USE temps, ONLY: annee_ref, day_ref
27	guez	92	use ymds2ju_m, only: ymds2ju
28	guez	5
29	guez	62	CHARACTER(len=*), INTENT(IN):: fichnom
30			INTEGER, INTENT(IN):: iday_end
31			REAL, INTENT(IN):: phis(:, :)
32	guez	5
33	guez	62	! Local:
34	guez	3
35	guez	62	INTEGER iq, l
36	guez	22	INTEGER, PARAMETER:: length = 100
37	guez	62	REAL tab_cntrl(length) ! tableau des paramètres du run
38	guez	5
39	guez	62	! Variables locales pour NetCDF:
40	guez	3
41	guez	62	INTEGER idim_index
42			INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
43			INTEGER idim_s, idim_sig
44			INTEGER idim_tim
45			INTEGER ncid, varid
46	guez	3
47	guez	62	REAL zjulian, hours
48			INTEGER yyears0, jjour0, mmois0
49			CHARACTER(len=30) unites
50	guez	3
51	guez	7	!-----------------------------------------------------------------------
52	guez	3
53	guez	18	PRINT *, 'Call sequence information: dynredem0'
54	guez	3
55	guez	23	CALL ymds2ju(annee_ref, 1, iday_end, 0., zjulian)
56	guez	18	CALL ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)
57	guez	5
58	guez	18	DO l = 1, length
59	guez	7	tab_cntrl(l) = 0.
60	guez	18	END DO
61			tab_cntrl(1) = real(iim)
62			tab_cntrl(2) = real(jjm)
63			tab_cntrl(3) = real(llm)
64			tab_cntrl(4) = real(day_ref)
65			tab_cntrl(5) = real(annee_ref)
66			tab_cntrl(6) = rad
67			tab_cntrl(7) = omeg
68			tab_cntrl(8) = g
69			tab_cntrl(9) = cpp
70	guez	7	tab_cntrl(10) = kappa
71			tab_cntrl(11) = daysec
72			tab_cntrl(12) = dtvr
73			tab_cntrl(13) = etot0
74			tab_cntrl(14) = ptot0
75			tab_cntrl(15) = ztot0
76			tab_cntrl(16) = stot0
77			tab_cntrl(17) = ang0
78			tab_cntrl(18) = pa
79			tab_cntrl(19) = preff
80	guez	5
81	guez	62	! Paramètres pour le zoom :
82	guez	3
83	guez	18	tab_cntrl(20) = clon
84			tab_cntrl(21) = clat
85			tab_cntrl(22) = grossismx
86			tab_cntrl(23) = grossismy
87	guez	5
88	guez	18	IF (fxyhypb) THEN
89	guez	7	tab_cntrl(24) = 1.
90			tab_cntrl(25) = dzoomx
91			tab_cntrl(26) = dzoomy
92			tab_cntrl(27) = 0.
93			tab_cntrl(28) = taux
94			tab_cntrl(29) = tauy
95			ELSE
96			tab_cntrl(24) = 0.
97			tab_cntrl(25) = dzoomx
98			tab_cntrl(26) = dzoomy
99			tab_cntrl(27) = 0.
100			tab_cntrl(28) = 0.
101			tab_cntrl(29) = 0.
102	guez	18	IF (ysinus) tab_cntrl(27) = 1.
103			END IF
104	guez	5
105	guez	18	tab_cntrl(30) = real(iday_end)
106	guez	5
107	guez	39	CALL nf95_create(fichnom, nf90_clobber, ncid)
108			CALL nf95_put_att(ncid, nf90_global, 'title', &
109	guez	18	'Fichier de démarrage dynamique')
110	guez	5
111	guez	7	! Definir les dimensions du fichiers:
112	guez	5
113	guez	39	CALL nf95_def_dim(ncid, 'index', length, idim_index)
114			CALL nf95_def_dim(ncid, 'rlonu', iip1, idim_rlonu)
115			CALL nf95_def_dim(ncid, 'rlatu', jjp1, idim_rlatu)
116			CALL nf95_def_dim(ncid, 'rlonv', iip1, idim_rlonv)
117			CALL nf95_def_dim(ncid, 'rlatv', jjm, idim_rlatv)
118			CALL nf95_def_dim(ncid, 'sigs', llm, idim_s)
119			CALL nf95_def_dim(ncid, 'sig', llmp1, idim_sig)
120			CALL nf95_def_dim(ncid, 'temps', nf90_unlimited, idim_tim)
121	guez	5
122	guez	7	! Definir et enregistrer certains champs invariants:
123	guez	5
124	guez	39	CALL nf95_def_var(ncid, 'controle', nf90_float, idim_index, varid)
125			CALL nf95_put_att(ncid, varid, 'title', 'Parametres de controle')
126	guez	5
127	guez	39	CALL nf95_def_var(ncid, 'rlonu', nf90_float, idim_rlonu, varid)
128			CALL nf95_put_att(ncid, varid, 'title', 'Longitudes des points U')
129	guez	5
130	guez	39	CALL nf95_def_var(ncid, 'rlatu', nf90_float, idim_rlatu, varid)
131			CALL nf95_put_att(ncid, varid, 'title', 'Latitudes des points U')
132	guez	5
133	guez	39	CALL nf95_def_var(ncid, 'rlonv', nf90_float, idim_rlonv, varid)
134			CALL nf95_put_att(ncid, varid, 'title', 'Longitudes des points V')
135	guez	5
136	guez	39	CALL nf95_def_var(ncid, 'rlatv', nf90_float, idim_rlatv, varid)
137			CALL nf95_put_att(ncid, varid, 'title', 'Latitudes des points V')
138	guez	5
139	guez	39	CALL nf95_def_var(ncid, 'ap', nf90_float, idim_sig, varid)
140			CALL nf95_put_att(ncid, varid, 'title', 'Coefficient A pour hybride')
141	guez	5
142	guez	39	CALL nf95_def_var(ncid, 'bp', nf90_float, idim_sig, varid)
143			CALL nf95_put_att(ncid, varid, 'title', 'Coefficient B pour hybride')
144	guez	5
145	guez	39	CALL nf95_def_var(ncid, 'presnivs', nf90_float, idim_s, varid)
146	guez	5
147	guez	7	! Coefficients de passage cov. <-> contra. <--> naturel
148	guez	5
149	guez	62	CALL nf95_def_var(ncid, 'cu', nf90_float, (/idim_rlonu, idim_rlatu/), varid)
150	guez	39	CALL nf95_put_att(ncid, varid, 'title', 'Coefficient de passage pour U')
151	guez	5
152	guez	62	CALL nf95_def_var(ncid, 'cv', nf90_float, (/idim_rlonv, idim_rlatv/), varid)
153	guez	39	CALL nf95_put_att(ncid, varid, 'title', 'Coefficient de passage pour V')
154	guez	5
155	guez	7	! Aire de chaque maille:
156	guez	5
157	guez	62	CALL nf95_def_var(ncid, 'aire', nf90_float, (/idim_rlonv, idim_rlatu/), &
158			varid)
159	guez	39	CALL nf95_put_att(ncid, varid, 'title', 'Aires de chaque maille')
160	guez	5
161	guez	7	! Geopentiel au sol:
162	guez	5
163	guez	62	CALL nf95_def_var(ncid, 'phisinit', nf90_float, &
164			(/idim_rlonv, idim_rlatu/), varid)
165	guez	39	CALL nf95_put_att(ncid, varid, 'title', 'Geopotentiel au sol')
166	guez	5
167	guez	7	! Definir les variables pour pouvoir les enregistrer plus tard:
168	guez	3
169	guez	39	CALL nf95_def_var(ncid, 'temps', nf90_float, idim_tim, varid)
170			CALL nf95_put_att(ncid, varid, 'title', 'Temps de simulation')
171	guez	62	WRITE(unites, 200) yyears0, mmois0, jjour0
172	guez	18	200 FORMAT ('days since ', I4, '-', I2.2, '-', I2.2, ' 00:00:00')
173	guez	39	CALL nf95_put_att(ncid, varid, 'units', unites)
174	guez	7
175	guez	62	CALL nf95_def_var(ncid, 'ucov', nf90_float, &
176			(/idim_rlonu, idim_rlatu, idim_s, idim_tim/), varid)
177	guez	39	CALL nf95_put_att(ncid, varid, 'title', 'Vitesse U')
178	guez	5
179	guez	62	CALL nf95_def_var(ncid, 'vcov', nf90_float, &
180			(/idim_rlonv, idim_rlatv, idim_s, idim_tim/), varid)
181	guez	39	CALL nf95_put_att(ncid, varid, 'title', 'Vitesse V')
182	guez	5
183	guez	62	CALL nf95_def_var(ncid, 'teta', nf90_float, &
184			(/idim_rlonv, idim_rlatu, idim_s, idim_tim/), varid)
185	guez	39	CALL nf95_put_att(ncid, varid, 'title', 'Temperature')
186	guez	5
187	guez	18	DO iq = 1, nqmx
188	guez	62	CALL nf95_def_var(ncid, tname(iq), nf90_float, &
189			(/idim_rlonv, idim_rlatu, idim_s, idim_tim/), varid)
190	guez	39	CALL nf95_put_att(ncid, varid, 'title', ttext(iq))
191	guez	18	END DO
192	guez	5
193	guez	62	CALL nf95_def_var(ncid, 'masse', nf90_float, &
194			(/idim_rlonv, idim_rlatu, idim_s, idim_tim/), varid)
195	guez	39	CALL nf95_put_att(ncid, varid, 'title', 'C est quoi ?')
196	guez	5
197	guez	62	CALL nf95_def_var(ncid, 'ps', nf90_float, &
198			(/idim_rlonv, idim_rlatu, idim_tim/), varid)
199	guez	39	CALL nf95_put_att(ncid, varid, 'title', 'Pression au sol')
200	guez	5
201	guez	39	CALL nf95_enddef(ncid)
202	guez	3
203	guez	39	CALL nf95_inq_varid(ncid, 'controle', varid)
204			CALL nf95_put_var(ncid, varid, tab_cntrl)
205	guez	3
206	guez	39	CALL nf95_inq_varid(ncid, 'rlonu', varid)
207			CALL nf95_put_var(ncid, varid, rlonu)
208	guez	18
209	guez	39	CALL nf95_inq_varid(ncid, 'rlatu', varid)
210			CALL nf95_put_var(ncid, varid, rlatu)
211	guez	18
212	guez	39	CALL nf95_inq_varid(ncid, 'rlonv', varid)
213			CALL nf95_put_var(ncid, varid, rlonv)
214	guez	18
215	guez	39	CALL nf95_inq_varid(ncid, 'rlatv', varid)
216			CALL nf95_put_var(ncid, varid, rlatv)
217	guez	18
218	guez	39	CALL nf95_inq_varid(ncid, 'ap', varid)
219			CALL nf95_put_var(ncid, varid, ap)
220	guez	18
221	guez	39	CALL nf95_inq_varid(ncid, 'bp', varid)
222			CALL nf95_put_var(ncid, varid, bp)
223	guez	18
224	guez	39	CALL nf95_inq_varid(ncid, 'presnivs', varid)
225			CALL nf95_put_var(ncid, varid, presnivs)
226	guez	18
227	guez	39	CALL nf95_inq_varid(ncid, 'cu', varid)
228			CALL nf95_put_var(ncid, varid, cu_2d)
229	guez	18
230	guez	39	CALL nf95_inq_varid(ncid, 'cv', varid)
231			CALL nf95_put_var(ncid, varid, cv_2d)
232	guez	18
233	guez	39	CALL nf95_inq_varid(ncid, 'aire', varid)
234			CALL nf95_put_var(ncid, varid, aire_2d)
235	guez	18
236	guez	39	CALL nf95_inq_varid(ncid, 'phisinit', varid)
237			CALL nf95_put_var(ncid, varid, phis)
238	guez	18
239	guez	39	CALL nf95_close(ncid) ! fermer le fichier
240	guez	18
241			PRINT *, 'iim, jjm, llm, iday_end', iim, jjm, llm, iday_end
242			PRINT *, 'rad, omeg, g, cpp, kappa', rad, omeg, g, cpp, kappa
243
244	guez	7	END SUBROUTINE dynredem0
245
246	guez	18	END MODULE dynredem0_m