dyn3d/Guide/guide.f

MODULE guide_m

  ! From dyn3d/guide.F, version 1.3 2005/05/25 13:10:09
  ! and dyn3d/guide.h, version 1.1.1.1 2004/05/19 12:53:06

  IMPLICIT NONE

  REAL aire_min, aire_max

CONTAINS

  SUBROUTINE guide(itau, ucov, vcov, teta, q, masse, ps)

    ! Author: F.Hourdin

    USE comconst, ONLY: cpp, daysec, dtvr, kappa
    USE comgeom, ONLY: aire, rlatu, rlonv
    USE conf_gcm_m, ONLY: day_step, iperiod
    use conf_guide_m, only: conf_guide, guide_u, guide_v, guide_t, guide_q, &
         ncep, ini_anal, tau_min_u, tau_max_u, tau_min_v, tau_max_v, &
         tau_min_t, tau_max_t, tau_min_q, tau_max_q, tau_min_p, tau_max_p, &
         online
    USE dimens_m, ONLY: iim, jjm, llm
    USE disvert_m, ONLY: ap, bp, preff, presnivs
    USE exner_hyb_m, ONLY: exner_hyb
    USE inigrads_m, ONLY: inigrads
    use massdair_m, only: massdair
    use netcdf, only: nf90_nowrite, nf90_open, nf90_close, nf90_inq_dimid, &
         nf90_inquire_dimension
    use nr_util, only: pi
    USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1, jjp1, llmp1
    USE q_sat_m, ONLY: q_sat
    use read_reanalyse_m, only: read_reanalyse
    USE serre, ONLY: clat, clon
    use tau2alpha_m, only: tau2alpha, dxdys

    INTEGER, INTENT(IN):: itau

    ! variables dynamiques
    REAL ucov(ip1jmp1, llm), vcov(ip1jm, llm) ! vents covariants
    REAL, intent(inout):: teta(ip1jmp1, llm) ! temperature potentielle 
    REAL q(ip1jmp1, llm) ! temperature potentielle 
    REAL, intent(out):: masse(ip1jmp1, llm) ! masse d'air
    REAL, intent(in):: ps(:, :) ! (iim + 1, jjm + 1) pression au sol

    ! Local:

    ! variables dynamiques pour les reanalyses.
    REAL, save:: ucovrea1(ip1jmp1, llm), vcovrea1(ip1jm, llm) !vts cov reas
    REAL, save:: tetarea1(ip1jmp1, llm) ! temp pot reales
    REAL, save:: qrea1(ip1jmp1, llm) ! temp pot reales
    REAL, save:: ucovrea2(ip1jmp1, llm), vcovrea2(ip1jm, llm) !vts cov reas
    REAL, save:: tetarea2(ip1jmp1, llm) ! temp pot reales
    REAL, save:: qrea2(ip1jmp1, llm) ! temp pot reales
    REAL, save:: masserea2(ip1jmp1, llm) ! masse

    REAL, save:: alpha_q(ip1jmp1)
    REAL, save:: alpha_t(ip1jmp1), alpha_p(ip1jmp1)
    REAL, save:: alpha_u(ip1jmp1), alpha_v(ip1jm)
    REAL dday_step, toto, reste
    real, save:: itau_test
    INTEGER, save:: step_rea, count_no_rea

    INTEGER ilon, ilat
    REAL factt, ztau(ip1jmp1)

    INTEGER ij, l
    INTEGER ncidpl, varidpl, status
    INTEGER rcod, rid
    REAL ditau, tau, a
    INTEGER, SAVE:: nlev

    ! TEST SUR QSAT
    REAL p(iim + 1, jjm + 1, llmp1), pk(ip1jmp1, llm), pks(ip1jmp1)
    REAL pkf(ip1jmp1, llm)
    REAL pres(ip1jmp1, llm)

    REAL qsat(ip1jmp1, llm)
    REAL unskap
    REAL tnat(ip1jmp1, llm)

    LOGICAL:: first = .TRUE.
    CHARACTER(len=10) file
    INTEGER:: igrads = 2
    REAL:: dtgrads = 100.

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

    PRINT *, 'Call sequence information: guide'

    ! calcul de l'humidite saturante

    forall (l = 1: llm + 1) p(:, :, l) = ap(l) + bp(l) * ps
    CALL massdair(p, masse)
    CALL exner_hyb(ps, p, pks, pk, pkf)
    tnat(:, :) = pk(:, :)*teta(:, :)/cpp
    unskap = 1./kappa
    pres(:, :) = preff*(pk(:, :)/cpp)**unskap
    qsat = q_sat(tnat, pres)

    ! initialisations pour la lecture des reanalyses.
    ! alpha determine la part des injections de donnees a chaque etape
    ! alpha=1 signifie pas d'injection
    ! alpha=0 signifie injection totale

    IF (online==-1) THEN
       RETURN
    END IF

    IF (first) THEN
       CALL conf_guide
       file = 'guide'
       CALL inigrads(igrads, rlonv, 180./pi, -180., 180., rlatu, -90., 90., &
            180./pi, presnivs, 1., dtgrads, file, 'dyn_zon ')
       PRINT *, '1: en-ligne, 0: hors-ligne (x=x_rea), -1: climat (x=x_gcm)'
       IF (online==-1) RETURN

       IF (online==1) THEN
          ! Constantes de temps de rappel en jour
          ! 0.1 c'est en gros 2h30. 
          ! 1e10 est une constante infinie donc en gros pas de guidage

          ! coordonnees du centre du zoom
          CALL coordij(clon, clat, ilon, ilat)
          ! aire de la maille au centre du zoom
          aire_min = aire(ilon+(ilat-1)*iip1)
          ! aire maximale de la maille
          aire_max = 0.
          DO ij = 1, ip1jmp1
             aire_max = max(aire_max, aire(ij))
          END DO
          ! factt = pas de temps en fraction de jour
          factt = dtvr*iperiod/daysec

          CALL tau2alpha(3, iip1, jjm, factt, tau_min_v, tau_max_v, alpha_v)
          CALL tau2alpha(2, iip1, jjp1, factt, tau_min_u, tau_max_u, alpha_u)
          CALL tau2alpha(1, iip1, jjp1, factt, tau_min_t, tau_max_t, alpha_t)
          CALL tau2alpha(1, iip1, jjp1, factt, tau_min_p, tau_max_p, alpha_p)
          CALL tau2alpha(1, iip1, jjp1, factt, tau_min_q, tau_max_q, alpha_q)

          CALL dump2d(iip1, jjp1, aire, 'AIRE MAILLe ')
          CALL dump2d(iip1, jjp1, alpha_u, 'COEFF U ')
          CALL dump2d(iip1, jjp1, alpha_t, 'COEFF T ')

          ! Cas ou on force exactement par les variables analysees
       ELSE
          alpha_t = 0.
          alpha_u = 0.
          alpha_v = 0.
          alpha_p = 0.
          ! physic=.false.
       END IF

       itau_test = 1001
       step_rea = 1
       count_no_rea = 0
       ncidpl = -99

       ! itau_test montre si l'importation a deja ete faite au rang itau
       ! lecture d'un fichier netcdf pour determiner le nombre de niveaux
       if (guide_u) then
          if (ncidpl.eq.-99) rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
       endif

       if (guide_v) then
          if (ncidpl.eq.-99) rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
       endif

       if (guide_T) then
          if (ncidpl.eq.-99) rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
       endif

       if (guide_Q) then
          if (ncidpl.eq.-99) rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl)
       endif

       IF (ncep) THEN
          status = nf90_inq_dimid(ncidpl, 'LEVEL', rid)
       ELSE
          status = nf90_inq_dimid(ncidpl, 'PRESSURE', rid)
       END IF
       status = nf90_inquire_dimension(ncidpl, rid, len=nlev)
       PRINT *, 'nlev', nlev
       rcod = nf90_close(ncidpl)
       ! Lecture du premier etat des reanalyses.
       CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, &
            masserea2, nlev)
       qrea2(:, :) = max(qrea2(:, :), 0.1)

       ! Debut de l'integration temporelle:
    END IF ! first

    ! IMPORTATION DES VENTS, PRESSION ET TEMPERATURE REELS:

    ditau = real(itau)
    dday_step = real(day_step)
    WRITE (*, *) 'ditau, dday_step'
    WRITE (*, *) ditau, dday_step
    toto = 4*ditau/dday_step
    reste = toto - aint(toto)

    IF (reste==0.) THEN
       IF (itau_test==itau) THEN
          WRITE (*, *) 'deuxieme passage de advreel a itau=', itau
          STOP
       ELSE
          vcovrea1(:, :) = vcovrea2(:, :)
          ucovrea1(:, :) = ucovrea2(:, :)
          tetarea1(:, :) = tetarea2(:, :)
          qrea1(:, :) = qrea2(:, :)

          PRINT *, 'LECTURE REANALYSES, pas ', step_rea, 'apres ', &
               count_no_rea, ' non lectures'
          step_rea = step_rea + 1
          itau_test = itau
          CALL read_reanalyse(step_rea, ps, ucovrea2, vcovrea2, tetarea2, &
               qrea2, masserea2, nlev)
          qrea2(:, :) = max(qrea2(:, :), 0.1)
          factt = dtvr*iperiod/daysec
          ztau(:) = factt/max(alpha_t(:), 1.E-10)
          CALL wrgrads(igrads, 1, aire, 'aire ', 'aire ')
          CALL wrgrads(igrads, 1, dxdys, 'dxdy ', 'dxdy ')
          CALL wrgrads(igrads, 1, alpha_u, 'au ', 'au ')
          CALL wrgrads(igrads, 1, alpha_t, 'at ', 'at ')
          CALL wrgrads(igrads, 1, ztau, 'taut ', 'taut ')
          CALL wrgrads(igrads, llm, ucov, 'u ', 'u ')
          CALL wrgrads(igrads, llm, ucovrea2, 'ua ', 'ua ')
          CALL wrgrads(igrads, llm, teta, 'T ', 'T ')
          CALL wrgrads(igrads, llm, tetarea2, 'Ta ', 'Ta ')
          CALL wrgrads(igrads, llm, qrea2, 'Qa ', 'Qa ')
          CALL wrgrads(igrads, llm, q, 'Q ', 'Q ')

          CALL wrgrads(igrads, llm, qsat, 'QSAT ', 'QSAT ')

       END IF
    ELSE
       count_no_rea = count_no_rea + 1
    END IF

    ! Guidage
    ! x_gcm = a * x_gcm + (1-a) * x_reanalyses

    IF (ini_anal) PRINT *, 'ATTENTION !!! ON PART DU GUIDAGE'

    ditau = real(itau)
    dday_step = real(day_step)

    tau = 4*ditau/dday_step
    tau = tau - aint(tau)

    ! ucov
    IF (guide_u) THEN
       DO l = 1, llm
          DO ij = 1, ip1jmp1
             a = (1.-tau)*ucovrea1(ij, l) + tau*ucovrea2(ij, l)
             ucov(ij, l) = (1.-alpha_u(ij))*ucov(ij, l) + alpha_u(ij)*a
             IF (first .AND. ini_anal) ucov(ij, l) = a
          END DO
       END DO
    END IF

    IF (guide_t) THEN
       DO l = 1, llm
          DO ij = 1, ip1jmp1
             a = (1.-tau)*tetarea1(ij, l) + tau*tetarea2(ij, l)
             teta(ij, l) = (1.-alpha_t(ij))*teta(ij, l) + alpha_t(ij)*a
             IF (first .AND. ini_anal) teta(ij, l) = a
          END DO
       END DO
    END IF

    IF (guide_q) THEN
       DO l = 1, llm
          DO ij = 1, ip1jmp1
             a = (1.-tau)*qrea1(ij, l) + tau*qrea2(ij, l)
             ! hum relative en % -> hum specif
             a = qsat(ij, l)*a*0.01
             q(ij, l) = (1.-alpha_q(ij))*q(ij, l) + alpha_q(ij)*a
             IF (first .AND. ini_anal) q(ij, l) = a
          END DO
       END DO
    END IF

    ! vcov
    IF (guide_v) THEN
       DO l = 1, llm
          DO ij = 1, ip1jm
             a = (1.-tau)*vcovrea1(ij, l) + tau*vcovrea2(ij, l)
             vcov(ij, l) = (1.-alpha_v(ij))*vcov(ij, l) + alpha_v(ij)*a
             IF (first .AND. ini_anal) vcov(ij, l) = a
          END DO
          IF (first .AND. ini_anal) vcov(ij, l) = a
       END DO
    END IF

    first = .FALSE.

  END SUBROUTINE guide

END MODULE guide_m
1	guez	20	MODULE guide_m
2	guez	3
3	guez	29	! From dyn3d/guide.F, version 1.3 2005/05/25 13:10:09
4			! and dyn3d/guide.h, version 1.1.1.1 2004/05/19 12:53:06
5	guez	3
6	guez	37	IMPLICIT NONE
7
8	guez	36	REAL aire_min, aire_max
9	guez	3
10	guez	20	CONTAINS
11	guez	3
12	guez	29	SUBROUTINE guide(itau, ucov, vcov, teta, q, masse, ps)
13	guez	3
14	guez	29	! Author: F.Hourdin
15	guez	3
16	guez	83	USE comconst, ONLY: cpp, daysec, dtvr, kappa
17			USE comgeom, ONLY: aire, rlatu, rlonv
18			USE conf_gcm_m, ONLY: day_step, iperiod
19	guez	44	use conf_guide_m, only: conf_guide, guide_u, guide_v, guide_t, guide_q, &
20	guez	85	ncep, ini_anal, tau_min_u, tau_max_u, tau_min_v, tau_max_v, &
21	guez	44	tau_min_t, tau_max_t, tau_min_q, tau_max_q, tau_min_p, tau_max_p, &
22			online
23	guez	88	USE dimens_m, ONLY: iim, jjm, llm
24	guez	83	USE disvert_m, ONLY: ap, bp, preff, presnivs
25			USE exner_hyb_m, ONLY: exner_hyb
26			USE inigrads_m, ONLY: inigrads
27	guez	67	use massdair_m, only: massdair
28	guez	44	use netcdf, only: nf90_nowrite, nf90_open, nf90_close, nf90_inq_dimid, &
29			nf90_inquire_dimension
30	guez	39	use nr_util, only: pi
31	guez	83	USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1, jjp1, llmp1
32			USE q_sat_m, ONLY: q_sat
33	guez	88	use read_reanalyse_m, only: read_reanalyse
34	guez	83	USE serre, ONLY: clat, clon
35	guez	44	use tau2alpha_m, only: tau2alpha, dxdys
36	guez	3
37	guez	83	INTEGER, INTENT(IN):: itau
38
39	guez	44	! variables dynamiques
40	guez	83	REAL ucov(ip1jmp1, llm), vcov(ip1jm, llm) ! vents covariants
41	guez	29	REAL, intent(inout):: teta(ip1jmp1, llm) ! temperature potentielle
42	guez	36	REAL q(ip1jmp1, llm) ! temperature potentielle
43	guez	70	REAL, intent(out):: masse(ip1jmp1, llm) ! masse d'air
44	guez	88	REAL, intent(in):: ps(:, :) ! (iim + 1, jjm + 1) pression au sol
45	guez	3
46	guez	83	! Local:
47
48	guez	44	! variables dynamiques pour les reanalyses.
49			REAL, save:: ucovrea1(ip1jmp1, llm), vcovrea1(ip1jm, llm) !vts cov reas
50			REAL, save:: tetarea1(ip1jmp1, llm) ! temp pot reales
51			REAL, save:: qrea1(ip1jmp1, llm) ! temp pot reales
52			REAL, save:: ucovrea2(ip1jmp1, llm), vcovrea2(ip1jm, llm) !vts cov reas
53			REAL, save:: tetarea2(ip1jmp1, llm) ! temp pot reales
54			REAL, save:: qrea2(ip1jmp1, llm) ! temp pot reales
55			REAL, save:: masserea2(ip1jmp1, llm) ! masse
56	guez	3
57	guez	44	REAL, save:: alpha_q(ip1jmp1)
58			REAL, save:: alpha_t(ip1jmp1), alpha_p(ip1jmp1)
59			REAL, save:: alpha_u(ip1jmp1), alpha_v(ip1jm)
60			REAL dday_step, toto, reste
61			real, save:: itau_test
62			INTEGER, save:: step_rea, count_no_rea
63	guez	3
64	guez	36	INTEGER ilon, ilat
65			REAL factt, ztau(ip1jmp1)
66	guez	3
67	guez	36	INTEGER ij, l
68	guez	44	INTEGER ncidpl, varidpl, status
69	guez	36	INTEGER rcod, rid
70			REAL ditau, tau, a
71	guez	44	INTEGER, SAVE:: nlev
72	guez	3
73	guez	44	! TEST SUR QSAT
74	guez	88	REAL p(iim + 1, jjm + 1, llmp1), pk(ip1jmp1, llm), pks(ip1jmp1)
75	guez	36	REAL pkf(ip1jmp1, llm)
76			REAL pres(ip1jmp1, llm)
77	guez	3
78	guez	36	REAL qsat(ip1jmp1, llm)
79			REAL unskap
80			REAL tnat(ip1jmp1, llm)
81	guez	3
82	guez	37	LOGICAL:: first = .TRUE.
83	guez	44	CHARACTER(len=10) file
84			INTEGER:: igrads = 2
85			REAL:: dtgrads = 100.
86	guez	3
87	guez	29	!-----------------------------------------------------------------------
88	guez	3
89	guez	29	PRINT *, 'Call sequence information: guide'
90	guez	3
91	guez	29	! calcul de l'humidite saturante
92	guez	3
93	guez	88	forall (l = 1: llm + 1) p(:, :, l) = ap(l) + bp(l) * ps
94	guez	29	CALL massdair(p, masse)
95			CALL exner_hyb(ps, p, pks, pk, pkf)
96			tnat(:, :) = pk(:, :)*teta(:, :)/cpp
97			unskap = 1./kappa
98			pres(:, :) = preff(pk(:, :)/cpp)*unskap
99			qsat = q_sat(tnat, pres)
100	guez	3
101	guez	44	! initialisations pour la lecture des reanalyses.
102			! alpha determine la part des injections de donnees a chaque etape
103			! alpha=1 signifie pas d'injection
104			! alpha=0 signifie injection totale
105	guez	3
106	guez	29	IF (online==-1) THEN
107			RETURN
108			END IF
109	guez	3
110	guez	29	IF (first) THEN
111			CALL conf_guide
112			file = 'guide'
113			CALL inigrads(igrads, rlonv, 180./pi, -180., 180., rlatu, -90., 90., &
114			180./pi, presnivs, 1., dtgrads, file, 'dyn_zon ')
115	guez	44	PRINT *, '1: en-ligne, 0: hors-ligne (x=x_rea), -1: climat (x=x_gcm)'
116			IF (online==-1) RETURN
117	guez	3
118	guez	29	IF (online==1) THEN
119	guez	44	! Constantes de temps de rappel en jour
120			! 0.1 c'est en gros 2h30.
121			! 1e10 est une constante infinie donc en gros pas de guidage
122	guez	3
123	guez	44	! coordonnees du centre du zoom
124	guez	29	CALL coordij(clon, clat, ilon, ilat)
125	guez	44	! aire de la maille au centre du zoom
126	guez	29	aire_min = aire(ilon+(ilat-1)*iip1)
127	guez	44	! aire maximale de la maille
128	guez	29	aire_max = 0.
129			DO ij = 1, ip1jmp1
130			aire_max = max(aire_max, aire(ij))
131			END DO
132	guez	44	! factt = pas de temps en fraction de jour
133	guez	29	factt = dtvr*iperiod/daysec
134	guez	3
135	guez	29	CALL tau2alpha(3, iip1, jjm, factt, tau_min_v, tau_max_v, alpha_v)
136			CALL tau2alpha(2, iip1, jjp1, factt, tau_min_u, tau_max_u, alpha_u)
137			CALL tau2alpha(1, iip1, jjp1, factt, tau_min_t, tau_max_t, alpha_t)
138			CALL tau2alpha(1, iip1, jjp1, factt, tau_min_p, tau_max_p, alpha_p)
139			CALL tau2alpha(1, iip1, jjp1, factt, tau_min_q, tau_max_q, alpha_q)
140	guez	3
141	guez	29	CALL dump2d(iip1, jjp1, aire, 'AIRE MAILLe ')
142	guez	44	CALL dump2d(iip1, jjp1, alpha_u, 'COEFF U ')
143			CALL dump2d(iip1, jjp1, alpha_t, 'COEFF T ')
144	guez	3
145	guez	44	! Cas ou on force exactement par les variables analysees
146	guez	29	ELSE
147			alpha_t = 0.
148			alpha_u = 0.
149			alpha_v = 0.
150			alpha_p = 0.
151	guez	44	! physic=.false.
152	guez	29	END IF
153	guez	3
154	guez	29	itau_test = 1001
155			step_rea = 1
156			count_no_rea = 0
157			ncidpl = -99
158	guez	3
159	guez	44	! itau_test montre si l'importation a deja ete faite au rang itau
160	guez	29	! lecture d'un fichier netcdf pour determiner le nombre de niveaux
161			if (guide_u) then
162			if (ncidpl.eq.-99) rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
163			endif
164	guez	3
165	guez	29	if (guide_v) then
166			if (ncidpl.eq.-99) rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
167			endif
168	guez	3
169	guez	29	if (guide_T) then
170			if (ncidpl.eq.-99) rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
171			endif
172	guez	3
173	guez	29	if (guide_Q) then
174			if (ncidpl.eq.-99) rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl)
175			endif
176	guez	3
177	guez	29	IF (ncep) THEN
178	guez	44	status = nf90_inq_dimid(ncidpl, 'LEVEL', rid)
179	guez	29	ELSE
180	guez	44	status = nf90_inq_dimid(ncidpl, 'PRESSURE', rid)
181	guez	29	END IF
182	guez	44	status = nf90_inquire_dimension(ncidpl, rid, len=nlev)
183	guez	29	PRINT *, 'nlev', nlev
184			rcod = nf90_close(ncidpl)
185	guez	44	! Lecture du premier etat des reanalyses.
186	guez	29	CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, &
187	guez	88	masserea2, nlev)
188	guez	29	qrea2(:, :) = max(qrea2(:, :), 0.1)
189	guez	3
190	guez	44	! Debut de l'integration temporelle:
191	guez	29	END IF ! first
192	guez	3
193	guez	29	! IMPORTATION DES VENTS, PRESSION ET TEMPERATURE REELS:
194	guez	3
195	guez	29	ditau = real(itau)
196			dday_step = real(day_step)
197			WRITE (, ) 'ditau, dday_step'
198			WRITE (, ) ditau, dday_step
199			toto = 4*ditau/dday_step
200			reste = toto - aint(toto)
201	guez	3
202	guez	29	IF (reste==0.) THEN
203			IF (itau_test==itau) THEN
204			WRITE (, ) 'deuxieme passage de advreel a itau=', itau
205			STOP
206			ELSE
207			vcovrea1(:, :) = vcovrea2(:, :)
208			ucovrea1(:, :) = ucovrea2(:, :)
209			tetarea1(:, :) = tetarea2(:, :)
210			qrea1(:, :) = qrea2(:, :)
211	guez	3
212	guez	29	PRINT *, 'LECTURE REANALYSES, pas ', step_rea, 'apres ', &
213			count_no_rea, ' non lectures'
214			step_rea = step_rea + 1
215			itau_test = itau
216			CALL read_reanalyse(step_rea, ps, ucovrea2, vcovrea2, tetarea2, &
217	guez	88	qrea2, masserea2, nlev)
218	guez	29	qrea2(:, :) = max(qrea2(:, :), 0.1)
219			factt = dtvr*iperiod/daysec
220			ztau(:) = factt/max(alpha_t(:), 1.E-10)
221	guez	44	CALL wrgrads(igrads, 1, aire, 'aire ', 'aire ')
222			CALL wrgrads(igrads, 1, dxdys, 'dxdy ', 'dxdy ')
223			CALL wrgrads(igrads, 1, alpha_u, 'au ', 'au ')
224			CALL wrgrads(igrads, 1, alpha_t, 'at ', 'at ')
225			CALL wrgrads(igrads, 1, ztau, 'taut ', 'taut ')
226			CALL wrgrads(igrads, llm, ucov, 'u ', 'u ')
227			CALL wrgrads(igrads, llm, ucovrea2, 'ua ', 'ua ')
228			CALL wrgrads(igrads, llm, teta, 'T ', 'T ')
229			CALL wrgrads(igrads, llm, tetarea2, 'Ta ', 'Ta ')
230			CALL wrgrads(igrads, llm, qrea2, 'Qa ', 'Qa ')
231			CALL wrgrads(igrads, llm, q, 'Q ', 'Q ')
232	guez	3
233	guez	44	CALL wrgrads(igrads, llm, qsat, 'QSAT ', 'QSAT ')
234	guez	3
235	guez	29	END IF
236			ELSE
237			count_no_rea = count_no_rea + 1
238			END IF
239	guez	3
240	guez	44	! Guidage
241			! x_gcm = a * x_gcm + (1-a) * x_reanalyses
242	guez	3
243	guez	29	IF (ini_anal) PRINT *, 'ATTENTION !!! ON PART DU GUIDAGE'
244	guez	3
245	guez	29	ditau = real(itau)
246			dday_step = real(day_step)
247	guez	3
248	guez	29	tau = 4*ditau/dday_step
249			tau = tau - aint(tau)
250	guez	3
251	guez	44	! ucov
252	guez	29	IF (guide_u) THEN
253			DO l = 1, llm
254			DO ij = 1, ip1jmp1
255			a = (1.-tau)ucovrea1(ij, l) + tauucovrea2(ij, l)
256			ucov(ij, l) = (1.-alpha_u(ij))ucov(ij, l) + alpha_u(ij)a
257			IF (first .AND. ini_anal) ucov(ij, l) = a
258			END DO
259			END DO
260			END IF
261	guez	3
262	guez	29	IF (guide_t) THEN
263			DO l = 1, llm
264			DO ij = 1, ip1jmp1
265			a = (1.-tau)tetarea1(ij, l) + tautetarea2(ij, l)
266			teta(ij, l) = (1.-alpha_t(ij))teta(ij, l) + alpha_t(ij)a
267			IF (first .AND. ini_anal) teta(ij, l) = a
268			END DO
269			END DO
270			END IF
271	guez	3
272	guez	29	IF (guide_q) THEN
273			DO l = 1, llm
274			DO ij = 1, ip1jmp1
275			a = (1.-tau)qrea1(ij, l) + tauqrea2(ij, l)
276	guez	44	! hum relative en % -> hum specif
277	guez	29	a = qsat(ij, l)a0.01
278			q(ij, l) = (1.-alpha_q(ij))q(ij, l) + alpha_q(ij)a
279			IF (first .AND. ini_anal) q(ij, l) = a
280			END DO
281			END DO
282			END IF
283	guez	3
284	guez	29	! vcov
285			IF (guide_v) THEN
286			DO l = 1, llm
287			DO ij = 1, ip1jm
288			a = (1.-tau)vcovrea1(ij, l) + tauvcovrea2(ij, l)
289			vcov(ij, l) = (1.-alpha_v(ij))vcov(ij, l) + alpha_v(ij)a
290			IF (first .AND. ini_anal) vcov(ij, l) = a
291			END DO
292			IF (first .AND. ini_anal) vcov(ij, l) = a
293			END DO
294			END IF
295	guez	3
296	guez	29	first = .FALSE.
297	guez	3
298	guez	29	END SUBROUTINE guide
299	guez	20
300			END MODULE guide_m