--- trunk/dyn3d/Guide/conf_guide.f	2014/09/19 11:41:35	114
+++ trunk/Sources/dyn3d/Guide/conf_guide.f	2015/09/29 19:48:59	171
@@ -2,32 +2,34 @@
 
   IMPLICIT NONE
 
-  !   Constantes de rappel. Unite : fraction de jour
-  REAL:: tau_min_u = 0.02
+  ! Constantes de rappel, en jours :
+  REAL:: tau_min_u = 0.03
   REAL:: tau_max_u = 10.
-  REAL:: tau_min_v = 0.02
+  REAL:: tau_min_v = 0.03
   REAL:: tau_max_v = 10.
-  REAL:: tau_min_t = 0.02
+  REAL:: tau_min_t = 0.03
   REAL:: tau_max_t = 10.
-  REAL:: tau_min_q = 0.02
+  REAL:: tau_min_q = 0.03
   REAL:: tau_max_q = 10.
-  REAL:: tau_min_p = 0.02
+  REAL:: tau_min_p = 0.03
   REAL:: tau_max_p = 10.
 
   LOGICAL:: ncep = .false. ! Coordonnee vert NCEP ou ECMWF
   LOGICAL:: ini_anal = .false. ! Initial = analyse
-  LOGICAL:: guide_u = .true. ! guidage de u
-  LOGICAL:: guide_v = .true. ! gvidage de v
+  LOGICAL:: guide_u = .false. ! guidage de u
+  LOGICAL:: guide_v = .false. ! gvidage de v
   LOGICAL:: guide_t = .false. ! guidage de T
   LOGICAL:: guide_q = .false. ! guidage de q
 
   logical:: online = .true. ! controle du guide
   ! hors-ligne: x=x_rea
 
-  ! Latitude min et max pour le rappel dans le cas ou on 'a les
-  ! analyses que sur une bande de latitudes.
-  REAL:: lat_min_guide = -90. ! Latitude minimum pour le guidage
-  real:: lat_max_guide = 90. ! Latitude maximum pour le guidage
+  ! Dans le cas où on n'a les analyses que sur une bande de latitudes :
+  REAL, save:: lat_min_guide ! minimum latitude for nudging, in rad
+  real, save:: lat_max_guide ! maximum latitude for nudging, in rad
+
+  logical, save:: ok_guide ! guidage
+  REAL, save:: factt ! pas de temps entre deux appels au guidage, en jours
 
 contains
 
@@ -36,12 +38,22 @@
     ! From LMDZ4/libf/dyn3d/conf_guide.F, version 1.1.1.1 2004/05/19 12:53:07
     !  Parametres de controle du run:
 
+    use abort_gcm_m, only: abort_gcm
+    use comconst, only: daysec, dtvr
+    use conf_gcm_m, only: day_step, iperiod
+    use dynetat0_m, only: grossismx, grossismy
+    use nr_util, only: assert, pi
     use unit_nml_m, only: unit_nml
 
+    ! Local:
+
+    REAL:: lat_min_guide_deg = -90. ! in degrees
+    real:: lat_max_guide_deg = 90. ! in degrees
+
     namelist /conf_guide_nml/ ncep, ini_anal, guide_u, guide_v, guide_t, &
-         online, 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, &
-         lat_min_guide, lat_max_guide
+         guide_q, online, 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, &
+         lat_min_guide_deg, lat_max_guide_deg
 
     !-----------------------------------------------------------------------
 
@@ -51,6 +63,28 @@
     read(unit=*, nml=conf_guide_nml)
     write(unit_nml, nml=conf_guide_nml)
 
+    lat_min_guide = lat_min_guide_deg / 180. * pi
+    lat_max_guide = lat_max_guide_deg / 180. * pi
+
+    ok_guide = any((/guide_u, guide_v, guide_t, guide_q/))
+    if (ok_guide .and. mod(day_step, 4 * iperiod) /= 0) call &
+         abort_gcm("conf_guide", 'ok_guide day_step iperiod')
+
+    if (ok_guide .and. online) then
+       factt = dtvr * iperiod / daysec
+       print *, "factt = ", factt
+       if (abs(grossismx - 1.) >= 0.1 .and. abs(grossismy - 1.) >= 0.1) then
+          if (guide_u) call assert(factt / tau_min_u < 1, &
+               "conf_guide tau_min_u")
+          if (guide_v) call assert(factt / tau_min_v < 1, &
+               "conf_guide tau_min_v")
+          if (guide_t) call assert(factt / tau_min_t < 1, &
+               "conf_guide tau_min_t")
+          if (guide_q) call assert(factt / tau_min_q < 1, &
+               "conf_guide tau_min_q")
+       end if
+    end if
+
   end SUBROUTINE conf_guide
 
 end module conf_guide_m