--- trunk/dyn3d/guide.f	2014/07/15 13:43:24	102
+++ trunk/dyn3d/guide.f	2014/09/11 15:09:15	107
@@ -18,18 +18,17 @@
     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
+         tau_min_t, tau_max_t, tau_min_q, tau_max_q, online
     USE dimens_m, ONLY: iim, jjm, llm
     USE disvert_m, ONLY: ap, bp, preff, presnivs
     use dump2d_m, only: dump2d
     USE exner_hyb_m, ONLY: exner_hyb
     USE inigrads_m, ONLY: inigrads
-    use massdair_m, only: massdair
-    use netcdf, only: nf90_nowrite, nf90_close, nf90_inq_dimid
-    use netcdf95, only: nf95_inquire_dimension, nf95_open
+    use netcdf, only: nf90_nowrite
+    use netcdf95, only: nf95_close, nf95_inq_dimid, nf95_inquire_dimension, &
+         nf95_open
     use nr_util, only: pi
-    USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1, jjp1, llmp1
+    USE paramet_m, ONLY: iip1, ip1jmp1, jjp1, llmp1
     USE q_sat_m, ONLY: q_sat
     use read_reanalyse_m, only: read_reanalyse
     USE serre, ONLY: clat, clon
@@ -67,7 +66,7 @@
     ! alpha=1 signifie pas d'injection
     ! alpha=0 signifie injection totale
     REAL, save:: alpha_q(iim + 1, jjm + 1)
-    REAL, save:: alpha_t(iim + 1, jjm + 1), alpha_p(ip1jmp1)
+    REAL, save:: alpha_t(iim + 1, jjm + 1)
     REAL, save:: alpha_u(iim + 1, jjm + 1), alpha_v(iim + 1, jjm)
 
     INTEGER, save:: step_rea, count_no_rea
@@ -77,8 +76,8 @@
     real ztau(iim + 1, jjm + 1)
 
     INTEGER ij, l
-    INTEGER ncidpl, status
-    INTEGER rcod, rid
+    INTEGER ncidpl
+    INTEGER rid
     REAL tau
     INTEGER, SAVE:: nlev
 
@@ -115,21 +114,20 @@
 
           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 tau2alpha(3, factt, tau_min_v, tau_max_v, alpha_v)
+          CALL tau2alpha(2, factt, tau_min_u, tau_max_u, alpha_u)
+          CALL tau2alpha(1, factt, tau_min_t, tau_max_t, alpha_t)
+          CALL tau2alpha(1, 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 ')
        ELSE
-          ! Cas ou on force exactement par les variables analysees
+          ! Cas où on force exactement par les variables analysées
           alpha_t = 0.
           alpha_u = 0.
           alpha_v = 0.
-          alpha_p = 0.
+          alpha_q = 0.
        END IF
 
        step_rea = 1
@@ -143,13 +141,13 @@
        if (guide_Q) call nf95_open('hur.nc',nf90_nowrite, ncidpl)
 
        IF (ncep) THEN
-          status = nf90_inq_dimid(ncidpl, 'LEVEL', rid)
+          call nf95_inq_dimid(ncidpl, 'LEVEL', rid)
        ELSE
-          status = nf90_inq_dimid(ncidpl, 'PRESSURE', rid)
+          call nf95_inq_dimid(ncidpl, 'PRESSURE', rid)
        END IF
        call nf95_inquire_dimension(ncidpl, rid, nclen=nlev)
        PRINT *, 'nlev', nlev
-       rcod = nf90_close(ncidpl)
+       call nf95_close(ncidpl)
        ! Lecture du premier etat des reanalyses.
        CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, &
             masserea2, nlev)