[ncl-talk] query ( CONVERSION ASCII TO NetCDF)
Kunal Bali
kunal.bali9 at gmail.com
Wed Sep 24 04:34:42 MDT 2014
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> 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>:
>
> 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> 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>:
>>
>> 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> 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>:
>>>
>>> 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> 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>:
>>>>
>>>> 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> 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","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> 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] *On
>>>>>> Behalf Of *Kunal Bali
>>>>>> *Sent:* Monday, September 22, 2014 8:53 AM
>>>>>> *To:* Karin Meier-Fleischer
>>>>>> *Cc:* 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","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> 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") ; remove any pre-existing file
>>>>>>
>>>>>> ncdf = addfile("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>:
>>>>>>
>>>>>> 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
>>>> Fax: +49 (0)40 460094 270
>>>> E-Mail: 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
> Fax: +49 (0)40 460094 270
> E-Mail: 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
>
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140924/7ed59b61/attachment.html
More information about the ncl-talk
mailing list