[ncl-talk] query ( CONVERSION ASCII TO NetCDF)
Karin Meier-Fleischer
meier-fleischer at dkrz.de
Wed Sep 24 04:52:01 MDT 2014
Please provide the script you are using!
Am 24.09.14 12:34, schrieb Kunal Bali:
> Dear Karin
>
> Thank You once again. Now its coming properly.
> One more problem i am facing right now is that i am trying to save the
> output result in pdf.
> I am using the commands before END command
> wks_type = "pdf"
> wks_type at wkPaperSize = "A4"
> wks = gsn_open_wks(wks_type,"test2")
>
> These command are generating the pdf file in directory but blank
> without any graphs (0 bytes).
>
> So how is the correct way to save the result in pdf or png format or
> how to use above commands correctly.
>
> Thank You
>
>
>
> Kunal Bali
> Research Scholar
> Radio & Atmospheric Science Division
> CSIR - National Physical Laboratory
> New Delhi - 110012
>
>
>
>
>
>
> On Wed, Sep 24, 2014 at 3:30 PM, Karin Meier-Fleischer
> <meier-fleischer at dkrz.de <mailto:meier-fleischer at dkrz.de>> wrote:
>
> Hi Kunal,
>
> didn't you use the isnan_ieee section in the script? You have NaNs
> in your data and they should be set to _FillValue.
>
> if (any(isnan_ieee(data))) then
> value = 1.e20
> replace_ieeenan (data, value, 0)
> data at _FillValue = value
> end if
>
> Again, without your script we could not help properly!
>
> Bye,
> Karin
>
> Am 24.09.2014 um 11:40 schrieb Kunal Bali <kunal.bali9 at gmail.com
> <mailto:kunal.bali9 at gmail.com>>:
>
>> Dear NCL users
>>
>> With the help of NCL users Now i am able to convert ASCII to
>> netcdf format. But when i tried to plot this .nc file the contour
>> map is not coming out in proper way.
>>
>> But when i use the original .nc file (without any conversion)
>> then the plot is coming nicely. Here is attached the file with
>> figure and all commands which i am using to plot.
>>
>> Thank You
>>
>> Kunal Bali
>> Research Scholar
>> Radio & Atmospheric Science Division
>> CSIR - National Physical Laboratory
>> New Delhi - 110012
>>
>>
>>
>>
>>
>>
>> On Wed, Sep 24, 2014 at 1:54 PM, Karin Meier-Fleischer
>> <meier-fleischer at dkrz.de <mailto:meier-fleischer at dkrz.de>> wrote:
>>
>> Hi Kunal,
>>
>> you read the data into an array called data with nlines=64800
>> and ncolumns=3 ---> data(nlines,ncolumns)
>>
>> The 1st column, all lines: data(:,0)
>> The 2nd column, all lines: data(:,1)
>> The 3rd column, all lines: data(:,2)
>>
>> nLongitude = 360
>> Longitude = data(0:nLongitude-1,1) ---> read 2nd column,
>> first 360 lines
>>
>> Please, take a look at
>>
>> http://www.ncl.ucar.edu/Document/Manuals/Getting_Started/basics.shtml#ArraySubscripting
>>
>> Bye,
>> Karin
>>
>> Am 24.09.2014 um 07:10 schrieb Kunal Bali
>> <kunal.bali9 at gmail.com <mailto:kunal.bali9 at gmail.com>>:
>>
>>> Dear*Karin Meier-Fleischer
>>> *
>>>
>>> Thank You so much for helping me. Your run your scripts and
>>> working well.
>>> I just want know what is the means of this line >
>>> *Longitude = data(*0:nLongitude-1*,1)
>>>
>>> *
>>> *Thank You
>>> *
>>>
>>> Kunal Bali
>>> Research Scholar
>>> Radio & Atmospheric Science Division
>>> CSIR - National Physical Laboratory
>>> New Delhi - 110012
>>>
>>>
>>>
>>>
>>>
>>>
>>> On Wed, Sep 24, 2014 at 1:01 AM, Karin Meier-Fleischer
>>> <meier-fleischer at dkrz.de <mailto:meier-fleischer at dkrz.de>>
>>> wrote:
>>>
>>> Hi Kunal,
>>>
>>> one problems occurred while reading the datacaused by
>>> the variable name TotCH4_A in the
>>> header line of your data file 14-01.txt. NCL reads the
>>> value 4 in the name and writes it as first
>>> value into the data array. After renaming to TotCH_A in
>>> 14-01.txt file everything works well.
>>>
>>> You missed to convert the one dimensional variable
>>> TotCH4_A (in your NCL script) to a two
>>> dimensional array. Also as described in the mails before
>>> you should assign latitude and
>>> longitude to the named variables of TotCH4_A2D. See the
>>> script below.
>>>
>>> It is easier to write/save the NCL commands in a text
>>> file (e.g. script.ncl) and run it
>>> using 'ncl script.ncl'
>>>
>>> Bye, Karin
>>>
>>> ;-----------------------------------------------------------------------------
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>
>>> begin
>>>
>>> nLatitude = 180
>>> nLongitude = 360
>>> data =
>>> asciiread("14-01.txt",(/nLatitude*nLongitude,3/),"float")
>>>
>>> Latitude = data(::nLongitude,0)
>>> Longitude = data(0:nLongitude-1,1)
>>> TotCH4_A = data(:,2)
>>> TotCH4_A2D = onedtond(TotCH4_A,(/nLatitude,nLongitude/))
>>>
>>> if (any(isnan_ieee(TotCH4_A2D))) then
>>> value = 1.e20
>>> replace_ieeenan (TotCH4_A2D, value, 0)
>>> TotCH4_A2D at _FillValue = value
>>> end if
>>>
>>> printMinMax(Latitude,0)
>>> printMinMax(Longitude,0)
>>> printMinMax(TotCH4_A,0)
>>>
>>> ; Assign named dimensions
>>>
>>> TotCH4_A2D!0 = "lat"
>>> TotCH4_A2D!1 = "lon"
>>>
>>> ; Assign coordinate variables
>>>
>>> TotCH4_A2D&lat = Latitude
>>> TotCH4_A2D&lon = Longitude
>>>
>>> TotCH4_A2D&lat at units = "degrees_north"
>>> TotCH4_A2D&lon at units = "degrees_east"
>>> printVarSummary(TotCH4_A2D)
>>>
>>> ; delete existing netcdf file to create a new one
>>> system("rm -rf Sample.nc")
>>> ncdf = addfile("Sample.nc","c")
>>>
>>> ; write data to netcdf file
>>>
>>> ncdf->TotCH4_A = TotCH4_A2D
>>>
>>> ; plot data
>>>
>>> wks = gsn_open_wks("x11","test")
>>> plot = gsn_csm_contour_map(wks,TotCH4_A2D,False)
>>>
>>> end
>>>
>>> ;-----------------------------------------------------------------------------
>>>
>>>
>>>
>>> Am 23.09.2014 um 19:45 schrieb Kunal Bali
>>> <kunal.bali9 at gmail.com <mailto:kunal.bali9 at gmail.com>>:
>>>
>>>> Here i attached scripts which i am trying to get the
>>>> ascii to netcdf conversion.
>>>>
>>>> Kunal Bali
>>>> Research Scholar
>>>> Radio & Atmospheric Science Division
>>>> CSIR - National Physical Laboratory
>>>> New Delhi - 110012
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> On Tue, Sep 23, 2014 at 11:03 PM, Karin Meier-Fleischer
>>>> <meier-fleischer at dkrz.de
>>>> <mailto:meier-fleischer at dkrz.de>> wrote:
>>>>
>>>> Hi Kunal,
>>>>
>>>> could you provide the data file and your script.
>>>>
>>>> Bye,
>>>> Karin
>>>>
>>>>
>>>> Am 23.09.2014 um 19:05 schrieb Kunal Bali
>>>> <kunal.bali9 at gmail.com <mailto:kunal.bali9 at gmail.com>>:
>>>>
>>>>> Dear NCL users
>>>>>
>>>>> I tried previous mail scripts for converting ASCII
>>>>> to netcdf format , the file is generating but its
>>>>> without data points.Only variable names generating.
>>>>>
>>>>> Kunal Bali
>>>>> Research Scholar
>>>>> Radio & Atmospheric Science Division
>>>>> CSIR - National Physical Laboratory
>>>>> New Delhi - 110012
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> On Tue, Sep 23, 2014 at 1:28 PM, Karin
>>>>> Meier-Fleischer <meier-fleischer at dkrz.de
>>>>> <mailto:meier-fleischer at dkrz.de>> wrote:
>>>>>
>>>>> Hi Kunal,
>>>>>
>>>>> the answer was given in an earlier mail:
>>>>>
>>>>> 1. read the ASCII data variable, lat and lon
>>>>> 2. create a output file: netcdf =
>>>>> addfile("simple.nc <http://simple.nc/>","c")
>>>>> 3. write your data variable to the output
>>>>> file: ncdf->values = values
>>>>>
>>>>>
>>>>> Please read the 'output to netcdf' examples
>>>>> for more details:
>>>>> http://ncl.ucar.edu/Applications/o-netcdf.shtml
>>>>>
>>>>> Bye,
>>>>> Karin
>>>>>
>>>>> Am 23.09.14 07:04, schrieb Kunal Bali:
>>>>>> Dear NCL users
>>>>>>
>>>>>> I am still not getting how to do the
>>>>>> conversion from ASCII to netcdf ??
>>>>>>
>>>>>> Kunal Bali
>>>>>> Research Scholar
>>>>>> Radio & Atmospheric Science Division
>>>>>> CSIR - National Physical Laboratory
>>>>>> New Delhi - 110012
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Mon, Sep 22, 2014 at 7:35 PM, Yuqiang
>>>>>> Zhang <yuqiangzhang.thu at gmail.com
>>>>>> <mailto:yuqiangzhang.thu at gmail.com>> wrote:
>>>>>>
>>>>>> Hi Kunal,
>>>>>>
>>>>>>
>>>>>> In line 30, "*ncl 30> TotCH4_A =
>>>>>> data(:,2)* ". You variable "TotCH4_A" is
>>>>>> only assigned one dimension's value, so
>>>>>> you can't assign two dimension through
>>>>>> line 33 & 34.
>>>>>>
>>>>>>
>>>>>> Regards,
>>>>>>
>>>>>> Yuqiang
>>>>>>
>>>>>>
>>>>>> *From:*ncl-talk-bounces at ucar.edu
>>>>>> <mailto:ncl-talk-bounces at ucar.edu>
>>>>>> [mailto:ncl-talk-bounces at ucar.edu
>>>>>> <mailto:ncl-talk-bounces at ucar.edu>] *On
>>>>>> Behalf Of *Kunal Bali
>>>>>> *Sent:* Monday, September 22, 2014 8:53 AM
>>>>>> *To:* Karin Meier-Fleischer
>>>>>> *Cc:* ncl-talk at ucar.edu
>>>>>> <mailto:ncl-talk at ucar.edu>
>>>>>> *Subject:* Re: [ncl-talk] query (
>>>>>> CONVERSION ASCII TO NetCDF)
>>>>>>
>>>>>>
>>>>>> Dear NCL user
>>>>>>
>>>>>> I tried this script from NCL but still
>>>>>> getting error
>>>>>>
>>>>>>
>>>>>> ncl 28> Longitude = data(:,1)
>>>>>> ncl 29> Latitude = data(:,0)
>>>>>> ncl 30> TotCH4_A = data(:,2)
>>>>>> ncl 31> ncdf =
>>>>>> addfile("/home/kunal/simple1.nc
>>>>>> <http://simple1.nc/>","c")
>>>>>> ncl 32> filedimdef(ncdf,"time",-1,True)
>>>>>> ncl 33> TotCH4_A!0 = "Latitude"
>>>>>> ncl 34> TotCH4_A!1 = "Longitude"
>>>>>> *fatal:Variable (TotCH4_A) has (1)
>>>>>> dimensions can not write to dimension (1)
>>>>>> fatal:["Execute.c":8565]:Execute: Error
>>>>>> occurred at or near line 34
>>>>>> *
>>>>>>
>>>>>>
>>>>>> Kunal Bali
>>>>>>
>>>>>> Research Scholar
>>>>>>
>>>>>> Radio & Atmospheric Science Division
>>>>>>
>>>>>> CSIR - National Physical Laboratory
>>>>>>
>>>>>> New Delhi - 110012
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Mon, Sep 22, 2014 at 5:32 PM, Karin
>>>>>> Meier-Fleischer <meier-fleischer at dkrz.de
>>>>>> <mailto:meier-fleischer at dkrz.de>> wrote:
>>>>>>
>>>>>> Hi Kunal,
>>>>>>
>>>>>>
>>>>>> Please, take a look at the 'output
>>>>>> netcdf' examples page, method_1
>>>>>>
>>>>>>
>>>>>> http://ncl.ucar.edu/Applications/method_1.shtml
>>>>>>
>>>>>> Try this:
>>>>>>
>>>>>>
>>>>>>
>>>>>> lat = data(:,0)
>>>>>>
>>>>>> lat = data(:,1)
>>>>>>
>>>>>> values= data(:,2)
>>>>>>
>>>>>> system
>>>>>> <http://ncl.ucar.edu/Document/Functions/Built-in/system.shtml>("/bin/rm
>>>>>> -f simple.nc <http://simple.nc/>") ;
>>>>>> remove any pre-existing file
>>>>>>
>>>>>> ncdf = addfile("simple.nc <http://simple.nc/>" ,"c") ; output netCDF file
>>>>>>
>>>>>> filedimdef <http://ncl.ucar.edu/Document/Functions/Built-in/filedimdef.shtml>(ncdf,"time",-1,True) ; recommended
>>>>>>
>>>>>> values!0 = "lat"
>>>>>>
>>>>>> values!1 = "lon"
>>>>>>
>>>>>> values&lat = lat
>>>>>>
>>>>>> values&lon = lon
>>>>>>
>>>>>> ncdf->values = values
>>>>>>
>>>>>>
>>>>>> Bye, Karin
>>>>>>
>>>>>>
>>>>>>
>>>>>> Am 22.09.2014 um 11:40 schrieb Kunal Bali
>>>>>> <kunal.bali9 at gmail.com
>>>>>> <mailto:kunal.bali9 at gmail.com>>:
>>>>>>
>>>>>> Dear NCL users
>>>>>>
>>>>>> I want to convert ASCII to .nc file
>>>>>> format. Could any one please provides
>>>>>> me the link for this query or give
>>>>>> any detail about this ??
>>>>>>
>>>>>> I have data in 3 column ( Latitude
>>>>>> Longitude Values)
>>>>>>
>>>>>> More Detail
>>>>>>
>>>>>> *Variable: data
>>>>>> Type: float
>>>>>> Total Size: 777600 bytes
>>>>>> 194400 values
>>>>>> Number of Dimensions: 2
>>>>>> Dimensions and sizes: [64800] x [3]
>>>>>> Coordinates:
>>>>>> Number Of Attributes: 1
>>>>>> _FillValue : 9.96921e+36
>>>>>>
>>>>>> *
>>>>>>
>>>>>> *So is there any ways to convert
>>>>>> ASCII to netcdf format ???*
>>>>>>
>>>>>> *Thank You*
>>>>>>
>>>>>>
>>>>>> Kunal Bali
>>>>>>
>>>>>> Research Scholar
>>>>>>
>>>>>> Radio & Atmospheric Science Division
>>>>>>
>>>>>> CSIR - National Physical Laboratory
>>>>>>
>>>>>> New Delhi - 110012
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> ncl-talk mailing list
>>>>>> List instructions, subscriber
>>>>>> options, unsubscribe:
>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> ncl-talk mailing list
>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options,
>>>>> unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>> ____________________________________
>>>> Dipl. Geophys. Karin Meier-Fleischer
>>>> Visualisierung
>>>> Abteilung Anwendungen
>>>>
>>>> Deutsches Klimarechenzentrum GmbH (DKRZ)
>>>> Bundesstrasse 45a * D-20146 Hamburg
>>>> Germany
>>>>
>>>> Phone:+49 (0)40 460094 126
>>>> <tel:%2B49%20%280%2940%20460094%20126>
>>>> Fax:+49 (0)40 460094 270
>>>> <tel:%2B49%20%280%2940%20460094%20270>
>>>> E-Mail:meier-fleischer at dkrz.de
>>>> <mailto:meier-fleischer at dkrz.de>
>>>> URL:http://www.dkrz.de/
>>>>
>>>> Geschäftsführer: Prof. Dr. Thomas Ludwig
>>>> Sitz der Gesellschaft: Hamburg
>>>> Amtsgericht Hamburg HRB 39784
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> <14-01.txt><ASCII to
>>>> netcdf.ncl>_______________________________________________
>>>> ncl-talk mailing list
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>>
>> <contour_error.png>_______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> ____________________________________
> Dipl. Geophys. Karin Meier-Fleischer
> Visualisierung
> Abteilung Anwendungen
>
> Deutsches Klimarechenzentrum GmbH (DKRZ)
> Bundesstrasse 45a * D-20146 Hamburg
> Germany
>
> Phone:+49 (0)40 460094 126 <tel:%2B49%20%280%2940%20460094%20126>
> Fax:+49 (0)40 460094 270 <tel:%2B49%20%280%2940%20460094%20270>
> E-Mail:meier-fleischer at dkrz.de <mailto:meier-fleischer at dkrz.de>
> URL:http://www.dkrz.de/
>
> Geschäftsführer: Prof. Dr. Thomas Ludwig
> Sitz der Gesellschaft: Hamburg
> Amtsgericht Hamburg HRB 39784
>
>
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
--
Dipl. Geophys. Karin Meier-Fleischer
Visualization
Application Support
Deutsches Klimarechenzentrum GmbH (DKRZ)
Bundesstrasse 45a - D20146 Hamburg - Germany
Phone: +49 (0)40 460094 126
Fax: +49 (0)40 460094 270
E-Mail: meier-fleischer at dkrz.de
URL: www.dkrz.de
Geschäftsführer: Prof. Dr. Thomas Ludwig
Sitz der Gesellschaft: Hamburg
Amtsgericht Hamburg HRB 39784
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140924/708be210/attachment.html
More information about the ncl-talk
mailing list