[gradsusr] Convert TXT to NetCDF

[Guillaume Dueymes]

Dear Dr Gamal,

i'm trying something through IDL. can you give me the dimensions of your
grid you should expect.
ie in Netcdf header :

*dimensions:
       ** x =  ;
        y =  ;
        time = UNLIMITED ; // (1 currently)
variables:
        double lon(y, x) ;
                lon:long_name = "longitude" ;
                lon:units = "degrees_east" ;
                lon:standard_name = "grid_longitude" ;
        double lat(y, x) ;
                lat:long_name = "latitude" ;
                lat:units = "degrees_north" ;
                lat:standard_name = "grid_latitude" ;
        double time(time) ;
                time:units = "day as %Y%m%d.%f" ;
        float preacc(time, y, x) ;
                preacc:long_name = "Precipitation" ;
                preacc:standard_name = "precipitation_flux" ;

*
Cheers

Guillaume D.
*
*
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
>
>
> 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
>
>


-- 
Guillaume Dueymes M.Sc.
Agent de Recherche
Centre ESCER - Université du Québec à Montréal
Environnement Canada - Service Météorologique du Canada
Université du Québec à Montréal
201 avenue Président-Kennedy, local PK-2315
Montréal (Québec) Canada
Tel: (514) 987-3000 poste 1407*
*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120131/1cbcf2fc/attachment-0003.html