[gradsusr] Convert TXT to NetCDF

[Arlindo da Silva]

On Mon, Jan 30, 2012 at 10:55 PM, Dr.Gamal El Afandi
<gamalafandy at yahoo.com>wrote:

> Dear all,
> Could you please let me know if anyone have idea about converting the TXTor 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
>

Why did you get ASCII involved at all? It appears that your binary file can
be directly read by grads (with the proper
ctl<http://grads.iges.org/grads/gadoc/descriptorfile.html>).
Once grads can open your file you can convert it to netcdf with
sdfwrite<http://grads.iges.org/grads/gadoc/gradcomdsdfwrite.html>,
or if you have opengrads, with lats4d<http://opengrads.org/doc/scripts/lats4d/>
.

   Arlindo



>
> 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
>
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>


-- 
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120131/7416e4e9/attachment-0003.html