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