!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  2013-01-14 by YZhang 
!  read_GL11_centroids.f
! -- sample to read the mean-centroid ASCII dataset for official ISCCP 
!    Global Clusters based on 1983 to 2009 D1 data product
!
! Import:
!  The mean-centroid ASCII dataset for a total of 12 clusters (11 cloudy
! and 1 clear-sky scenes). Each of 11 cloudy clusters consists of 42 vectors
! of cloud fraction (in %) in the same order as ISCCP D1's parameters 
! #30d to #71d (six cloud optical thickness ranges over seven cloud top 
! pressure ranges from TOA to surface) with the last, 12th cluster for 
! clear sky (all = 0.):  
! 
!   ~/ctd_GL_8309_43_o2.asc  
!
! Export:
! -- to be customized by the users
!-------------------------------------------------
      parameter (mclus=12, nele=42)
      REAL*4    ctd_in(mclus,nele)
      character*80, dimension(1), parameter ::
     &   ufin = (/'/Users/yzhang11/WS/km_dat/ctd_GL_8309_43_o2.asc'/)
      Logical*4 Lexist
!-----------------------------------------------------------
       inquire(file=ufin, exist=Lexist)
       if(Lexist .NEQV. .true.)then
          print *,'??? Non-existent ufin=',ufin
          stop '11111'
       end if
!
      print *,'input #22=',ufin
      OPEN(22,FILE=ufin,status='old',iostat=irc)
!
      DO i=1,mclus
        read(22,*) (ctd_in(i,k),k=1,nele)
      print *,ctd_in(i,:)
      pause 
      ENDDO

      CLOSE(22) 
!
      stop
      end