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 comvert, ONLY : ap, bp, nivsig, nivsigs, pa, preff, presnivs
    USE comgeom, ONLY : aire_2d, cu_2d, cv_2d, rlatu, rlatv, rlonu, rlonv
    USE dimens_m, ONLY : iim, jjm, llm, nqmx
    USE ener, ONLY : ang0, etot0, ptot0, stot0, ztot0
    USE calendar, ONLY : ju2ymds, ymds2ju
    USE iniadvtrac_m, ONLY : tname, ttext
    USE logic, ONLY : fxyhypb, ysinus
    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
    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

    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 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 :: 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, 'nivsigs', nf90_float, idim_s, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Numero naturel des couches s')

    CALL nf95_def_var(ncid, 'nivsig', nf90_float, idim_sig, varid)
    CALL nf95_put_att(ncid, varid, 'title', &
         'Numero naturel des couches sigma')

    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

    dims2(1) = idim_rlonu
    dims2(2) = idim_rlatu
    CALL nf95_def_var(ncid, 'cu', nf90_float, dims2, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Coefficient de passage pour U')

    dims2(1) = idim_rlonv
    dims2(2) = idim_rlatv
    CALL nf95_def_var(ncid, 'cv', nf90_float, dims2, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Coefficient de passage pour V')

    ! Aire de chaque maille:

    dims2(1) = idim_rlonv
    dims2(2) = idim_rlatu
    CALL nf95_def_var(ncid, 'aire', nf90_float, dims2, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Aires de chaque maille')

    ! Geopentiel au sol:

    dims2(1) = idim_rlonv
    dims2(2) = idim_rlatu
    CALL nf95_def_var(ncid, 'phisinit', nf90_float, dims2, 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)


    dims4(1) = idim_rlonu
    dims4(2) = idim_rlatu
    dims4(3) = idim_s
    dims4(4) = idim_tim
    CALL nf95_def_var(ncid, 'ucov', nf90_float, dims4, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Vitesse U')

    dims4(1) = idim_rlonv
    dims4(2) = idim_rlatv
    dims4(3) = idim_s
    dims4(4) = idim_tim
    CALL nf95_def_var(ncid, 'vcov', nf90_float, dims4, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'Vitesse V')

    dims4(1) = idim_rlonv
    dims4(2) = idim_rlatu
    dims4(3) = idim_s
    dims4(4) = idim_tim
    CALL nf95_def_var(ncid, 'teta', nf90_float, dims4, varid)
    CALL nf95_put_att(ncid, varid, '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(ncid, tname(iq), nf90_float, dims4, varid)
       CALL nf95_put_att(ncid, varid, '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(ncid, 'masse', nf90_float, dims4, varid)
    CALL nf95_put_att(ncid, varid, 'title', 'C est quoi ?')

    dims3(1) = idim_rlonv
    dims3(2) = idim_rlatu
    dims3(3) = idim_tim
    CALL nf95_def_var(ncid, 'ps', nf90_float, dims3, 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, 'nivsigs', varid)
    CALL nf95_put_var(ncid, varid, nivsigs)

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

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