Changeset 612 for BOL/trunk

Timestamp:

Mar 25, 2005, 4:04:06 PM (20 years ago)

Author:

Laurent Fairhead

Message:

Probleme sur les valeurs de w500 qui peuvent etre hors de l'intervalle
defini pour la creation des bins
LF

File:

• BOL/trunk/Class_Reg/geo2reg.F90 (modified) (7 diffs)

BOL/trunk/Class_Reg/geo2reg.F90

@@ -45 +45 @@
   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 +467 @@
 ! initialisations:
 
-
   allocate(xpdfmean(nb_bin))
   allocate(xpdf(nb_bin, itime))
@@ -474 +473 @@
   allocate(w_binw(nb_bin, itime))
   allocate(x_binw(nb_bin, itime))
-
-   
-   
-
 
 !-- temporal loop:
@@ -538 +533 @@
       enddo
     enddo
+    nbr_inf = 0; nbr_sup = 0 
 ! tropical belt:
     do j = 1, jm
@@ -549 +545 @@
             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 +569 @@
       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 +632 @@
     stop
   endif
+  deallocate(var_dim)