Index: BOL/trunk/Class_Reg/geo2reg.F90
===================================================================
--- BOL/trunk/Class_Reg/geo2reg.F90	(revision 611)
+++ BOL/trunk/Class_Reg/geo2reg.F90	(revision 612)
@@ -45,5 +45,5 @@
   real, allocatable, dimension(:)   :: xpdfmean
 
-  integer nb_bin
+  integer nb_bin, nbr_inf, nbr_sup
   real min_bin, max_bin, step_bin, w_mult, x_mult
   real xpdftot(mmax)
@@ -467,5 +467,4 @@
 ! initialisations:
 
-
   allocate(xpdfmean(nb_bin))
   allocate(xpdf(nb_bin, itime))
@@ -474,8 +473,4 @@
   allocate(w_binw(nb_bin, itime))
   allocate(x_binw(nb_bin, itime))
-
-   
-   
-
 
 !-- temporal loop:
@@ -538,4 +533,5 @@
       enddo
     enddo
+    nbr_inf = 0; nbr_sup = 0 
 ! tropical belt:
     do j = 1, jm
@@ -549,5 +545,13 @@
             x1 = x(i,j)*x_mult
 ! bin w500:
-            ir = INT((w1-min_bin)/step_bin) + 1
+            if (w1 < min_bin) then
+              ir = 1
+              nbr_inf = nbr_inf + 1
+            else if (w1 > max_bin) then
+              ir = nb_bin
+              nbr_sup = nbr_sup + 1
+            else
+              ir = INT((w1-min_bin)/step_bin) + 1
+            endif
 ! monthly PDF :
             xpdfmeantot = xpdfmeantot + surf(i,j)
@@ -565,4 +569,6 @@
       endif ! lat
     enddo ! j
+    write(lunout,*)'nbre de points ou w500 < ',min_bin,' = ', nbr_inf
+    write(lunout,*)'nbre de points ou w500 > ',max_bin,' = ', nbr_sup 
   enddo ! m
 
@@ -626,4 +632,5 @@
     stop
   endif
+  deallocate(var_dim)