program 

c*************************************************************************************
c  	   to read  training_data SeeBorV4.1 with 2003 BF emissivity Version 2
c
c          example to compile: f77 read_SeeBorV5_TD.f -o read_SeeBorV5_TDS.e         
c
c	    July 2009 Created by Eva Borbas, University of Wisconsin, SSEC/CIMSS
c*************************************************************************************

          integer nl,leng, nrec,nemis
	  parameter(nl=101,nlayer=nl-1,leng=338)
	  parameter(nrec=15704,nemis=10)
	  real prof(leng)
     	  
	  real p(nl),t(nl),w(nl),o(nl)
	  
	  real emis(nemis),freq(nemis)
 
         REAL igbp_type
	 REAL alon,jday
	 integer nday, year,month,day,hour
      
        data p/0.005,.016,.038,.077,.137,.224,.345,.506,.714,
     &  .975,1.297,1.687,2.153,2.701,3.340,4.077,4.920,
     &  5.878,6.957,8.165,9.512,11.004,12.649,14.456,16.432,
     &  18.585,20.922,23.453,26.183,29.121,32.274,35.651,39.257,
     &  43.100,47.188,51.528,56.126,60.989,66.125,71.540,77.240,
     &  83.231,89.520,96.114,103.017,110.237,117.777,125.646,133.846,
     &  142.385,151.266,160.496,170.078,180.018,190.320,200.989,212.028,
     &  223.441,235.234,247.408,259.969,272.919,286.262,300.000,314.137,
     &  328.675,343.618,358.966,374.724,390.893,407.474,424.470,441.882,
     &  459.712,477.961,496.630,515.720,535.232,555.167,575.525,596.306,
     &  617.511,639.140,661.192,683.667,706.565,729.886,753.628,777.790,
     &  802.371,827.371,852.788,878.620,904.866,931.524,958.591,986.067,
     &  1013.948,1042.232,1070.917,1100.000/


	 open(2,file='SeeBorV5.1_Training_data_Emis10inf2004.bin',recl=leng*4,
     &	 access='direct',status='old')
           
	  do i=1,nrec
	   
	  	read(2,rec=i) prof

		alon=prof(3*nl+2)
		year=prof(313)
		day=prof(315)
		hour=prof(316)
 
   		alat=prof(3*nl+1)
		amonth=prof(3*nl+11)
		month=amonth
	  
	  	igbp_type=prof(310)

		if(i.eq.2) then
			write(*,*) prof
		endif	
		
	  
	  	do j=1,nemis
	  		freq(j)=prof(3*nl+14+j)
			emis(j)=prof(3*nl+14+nemis+j)
c	  		write(*,'(2i8,2f8.2,i5)') 
c     &			i,j,freq(j),emis(j),month
	  	enddo

          	sfc_height=prof(311)
	  
	  	DO J=1,NL
             		T(J)=prof(J)
             		W(J)=prof(J+nl)*1000.0
             		O(J)=prof(J+2*nl)
	     		if(O(J).lt.0.) then 
	         		O(J)=0.
	    		endif 

         	ENDDO

          enddo

	 close(2)
	   
	 stop
	 end