[ncl-talk] query ( CONVERSION ASCII TO NetCDF)
Kunal Bali
kunal.bali9 at gmail.com
Wed Sep 24 06:33:38 MDT 2014
Dear NCL user
i tried these command to save my plot in pdf file. but its not working. pdf
creating but without map. Here i attached my scripts for plotting map.
please look at this and let me know what am i doing wrong.
wks_type = "pdf"
wks_type at wkPaperSize = "A4"
wks = gsn_open_wks(wks_type,"sample1")
My aim to just save my results in pdf format.
Thaank YOu
Kunal Bali
Research Scholar
Radio & Atmospheric Science Division
CSIR - National Physical Laboratory
New Delhi - 110012
On Wed, Sep 24, 2014 at 5:14 PM, Karin Meier-Fleischer <
meier-fleischer at dkrz.de> wrote:
> You set wks to x11 and create a plot. The second setting of wks do not
> have any influence.
> Change the the first wks setting to the pdf settings and it should work.
>
> ;---Open workstation and change color map
> wks_type = "pdf"
> wks_type at wkPaperSize = "A4"
> wks = gsn_open_wks(wks_type,"sample1")
>
>
> ;---Set some resources
> res = True
>
> res at gsnMaximize = True ; maximize plot in frame
>
> res at cnFillOn = True ; turn on contour fill
> res at cnLinesOn = False ; turn off contour lines
> res at cnLineLabelsOn = False ; turn off line labels
>
> res at tiMainString = filename
>
> res at gsnAddCyclic = False
>
> plot = gsn_csm_contour_map(wks,totc,res)
>
>
> Bye,
> Karin
>
> Am 24.09.14 13:16, schrieb Kunal Bali:
>
> please find the attached file
>
> Thank You
>
> Kunal Bali
> Research Scholar
> Radio & Atmospheric Science Division
> CSIR - National Physical Laboratory
> New Delhi - 110012
>
>
>
>
>
>
> On Wed, Sep 24, 2014 at 4:22 PM, Karin Meier-Fleischer <
> meier-fleischer at dkrz.de> wrote:
>
>> 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> 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 <%2B49%20%280%2940%20460094%20126>
>>> Fax: +49 (0)40 460094 270 <%2B49%20%280%2940%20460094%20270>
>>> 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
>>>
>>>
>>>
>>>
>>>
>>
>>
>> _______________________________________________
>> 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
>>
>>
>> _______________________________________________
>> 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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140924/c40e7e33/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: TotCh4_plot1.ncl
Type: application/octet-stream
Size: 943 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140924/c40e7e33/attachment.obj
More information about the ncl-talk
mailing list