Changeset 1414 for LMDZ4/trunk/libf/cosp/llnl_stats.F90

Timestamp:

Jul 15, 2010, 5:21:22 PM (14 years ago)

Author:

idelkadi

Message:

Passage a la version cosp.v1.3 pour le Lidar et ISCCP
Corrections de bugs pour ISCCP et optimisation pour le Lidar

File:

• LMDZ4/trunk/libf/cosp/llnl_stats.F90 (modified) (4 diffs)

LMDZ4/trunk/libf/cosp/llnl_stats.F90

@@ -24 +24 @@
 
 MODULE MOD_LLNL_STATS
+  USE MOD_COSP_CONSTANTS
   IMPLICIT NONE
 
@@ -62 +63 @@
    do j = 1, Nlevels, 1
       do k = 1, Ncolumns, 1
-         do i = 1, Npoints, 1 
-            if ((x(i,k,j) >= xmin) .and. (x(i,k,j) <= xmax)) then 
+         do i = 1, Npoints, 1
+            if (x(i,k,j) == R_GROUND) then
+               cosp_cfad(i,:,j) = R_UNDEF
+            elseif ((x(i,k,j) >= xmin) .and. (x(i,k,j) <= xmax)) then 
                ibin = ceiling((x(i,k,j) - bmin)/bwidth)
                if (ibin > Nbins) ibin = Nbins
@@ -72 +75 @@
       enddo  !k
    enddo  !j
-   cosp_cfad = cosp_cfad / Ncolumns
+   where ((cosp_cfad /= R_UNDEF).and.(cosp_cfad /= 0.0)) cosp_cfad = cosp_cfad / Ncolumns
 END FUNCTION COSP_CFAD
 
@@ -98 +101 @@
    integer :: pr,i,j
    
-!    lidar_only_freq_cloud = 0.0
-!    tcc = 0.0
+   lidar_only_freq_cloud = 0.0
+   tcc = 0.0
    do pr=1,Npoints
      do i=1,Ncolumns