[gradsusr] Convert TXT to NetCDF

Dr.Gamal El Afandi gamalafandy at yahoo.com
Mon Jan 30 22:55:42 EST 2012



Dear all,
Could you please let me know if anyone have idea about converting the TXT or ASCII data into NetCDF?
Please attached is the file required to be transformed into NetCDF.
Please below is the program I used to convert from binary to TXT
 
>PROGRAM convert_rfe_bin2asc
>
>c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
>c  Program: convert_rfe2.f
>c
>c  Description: This program reads in a daily RFE (rainfall estimate) file
>c              (version 2) in .bin format and creates an ascii integer
>c              file.  Each file contains rfe coordinate data ordered by rows.
>c              Rainfall resolution is 0.1 degree with a domain of
>c              20W-55E, 40S-40N.
>c
>c  Inputs: 'input_file' = the 10-day binary precip file to be converted
>c
>c  Output: 'output.txt' = the output ascii converted text file
>c
>c  Modified: July, 2010 by Nick Novella (nicholas.novella at noaa.gov)
>c
>c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
>
>REAL*4 rfe(751,801),arr(3,601551)
>real*4 lat,lon
>INTEGER:: i,j,ind
>character*25:: input_file
>
>c #####  CHANGE FILE NAME BELOW TO INPUT FILE  #####
>
>input_file="all_products.bin.20100621"
>
>c ##################################################
>
>c #####  Initialize the arrays  #####
>do i = 1,751 
>do j = 1,801 
>  rfe (i,j) = -999
>end do
>end do
>
>c #####  Read in the input rainfall data  #####
>open
 (unit=22,file=input_file,access="direct",
>1	convert='big_endian',status="old",recl=751*801*4)
>read  (22, rec=1) rfe 
>
>c #####  Open the output file  #####
>open (88,file="output.txt",access='sequential',
>1	status='unknown',form='formatted')
>
>c #####  The next code converts the input rainfall data. Each row is then written
>c        to the output file.
>
>ind=1 
>  	lat= -40.0
>lon= -20.0  
>do j=1,801
>do i=1,751
>          arr(1,ind)=lat
>          arr(2,ind)=lon
>          arr(3,ind)=rfe(i,j)
>
>write (88,499) arr(1,ind), arr(2,ind), arr(3,ind)
>499	format(1X,f6.2,1X,f6.2,1X,f8.2)
>
>   ind=ind+1
>   lon=lon+0.1
>end do
>        lat=lat+0.1
>        lon=-20.0
>end do
>
>
>
>
>
>      END PROGRAM
>
>
>With my kind regards
>Gamal
>___________________________________________________________
>Gamal Salah El Afandi, Ph.D. 
>College of Agricultural, Environment and Nutrition Sciences  
>College of Engineering  
>Tuskegee University, Tuskegee Alabama USA
>
>http://www.tuskegee.edu/
>http://www.tuskegee.edu/academics/colleges/caens.aspx
>http://www.tuskegee.edu/academics/colleges/ceps.aspx
>http://maindb.unfccc.int/public/roe/expert_list.pl?id_country=66&;
>http://www.itap.purdue.edu/pto/NATO_KEEO/partCoDir_en.html
>http://www.ucar.edu/governance/iap/iapreps.shtml 
>
>
>   
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120130/f170468f/attachment-0003.html 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: output.txt
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20120130/f170468f/attachment-0003.txt 


More information about the gradsusr mailing list