[ncl-talk] write variable attribute
Luo, Chao
chao.luo at eas.gatech.edu
Fri Aug 28 16:51:58 MDT 2015
Sorry, I made the mistake. It works. This write out metadata, including lat/lon. Is it possible to just write out variable value and attribute, not lat/lon data?
Thanks much,
Chao
----- Original Message -----
From: "Dennis Shea" <shea at ucar.edu>
To: "Chao Luo" <chao.luo at eas.gatech.edu>
Cc: "Alan Brammer" <abrammer at albany.edu>, "ncl-talk Talk" <ncl-talk at ucar.edu>
Sent: Friday, August 28, 2015 1:41:32 PM
Subject: Re: [ncl-talk] write variable attribute
Alan's script should work.
system("/bin/rm -f test.nc")
ncout = addfile("test.nc","c")
x = 0
x!0 = "test"
x&test = 23
x at Some_att = "First Att"
x at Some_other_att = "Second Att"
x at date = systemfunc("date")
ncout->x = x
quit
If the above does not work ... something is wrong
diro = "/data11/cluo/cmaq4.6/data/post/smk31_post/"
filo = "out.nc" ; input file name
fl_nm_out = diro+filo
system ("/bin/rm -f "+fl_nm_out) ; Remove any pre-existing file
ncdf = addfile(fl_nm_out,"c")
[snip]
temp2 = ESMF_regrid(temp2_36,Opt) ; Do the regridding
printVarSummary(temp2) ; this has attributes and possiblt time/level
; coordinate information
; not sure why you are doing this
temp2!0 = "time"
temp2!1 = "lay"
temp2!2 = "lat"
temp2!3 = "lon"
temp2 at long_name = "TEMP2"
temp2 at units = "k"
temp2 at var_desc = "temperature at 2 m"
printVarSummary (temp2)
printMinMax(temp2,True)
ncdf->TFLAG =TFLAG
ncdf->temp2 = temp2 ;***** NOT (/temp2/)
; the (/.../) strip meta data
On Fri, Aug 28, 2015 at 2:21 PM, Luo, Chao <chao.luo at eas.gatech.edu> wrote:
> Thanks very much for your help and time. I still can't write out variable
> attribute even I followed your sample, and I don't what caused the problem.
> Appreciate any suggestions.
>
> Chao
>
> ________________________________
> From: "Alan Brammer" <abrammer at albany.edu>
> To: "Chao Luo" <chao.luo at eas.gatech.edu>
> Cc: "ncl-talk Talk" <ncl-talk at ucar.edu>
> Sent: Friday, August 28, 2015 12:39:28 PM
>
> Subject: Re: [ncl-talk] write variable attribute
>
> The below example is how it should work. A ncl_filedump of your output
> netcdf might be useful for people to see what is missing compared to the
> data in the script.
> Not sure I follow what the problem is. If the metadata is written I don't
> see how the attributes weren't also.
>
>
> $ ncl
>
>
> ncout = addfile("test.nc","c")
> x = 0
> x!0 = "test"
> x&test = 23
> x at Some_att = "First Att"
> x at Some_other_att = "Second Att"
> x at date = systemfunc("date")
> ncout->x = x
> exit
>
> $ ncl_filedump test.nc
> Copyright (C) 1995-2015 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 6.3.0
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
>
> Variable: f
> Type: file
> filename: test
> path: test.nc
> file global attributes:
> dimensions:
> test = 1
> variables:
> integer x ( test )
> date : Fri Aug 28 19:32:35 UTC 2015
> Some_other_att : Second Att
> Some_att : First Att
>
> integer test ( test )
>
>
> On Fri, Aug 28, 2015 at 3:13 PM, Luo, Chao <chao.luo at eas.gatech.edu> wrote:
>>
>> Yes, it write out metadata of TEMP2, but still no TEMP2 attribute write
>> out when I use ncdf->tem2 = temp2. What I want to write out is TEMP2 values
>> and attribute. Thanks!
>>
>> Chao
>>
>> ________________________________
>> From: "Alan Brammer" <abrammer at albany.edu>
>> To: "Chao Luo" <chao.luo at eas.gatech.edu>
>> Cc: "ncl-talk Talk" <ncl-talk at ucar.edu>
>> Sent: Friday, August 28, 2015 11:06:55 AM
>> Subject: Re: [ncl-talk] write variable attribute
>>
>>
>> You've commented out the line that explicitly writes the attributes to the
>> file. Likely because you haven't explicitly defined the variable in the
>> file.
>> But by also surrounding temp2 with (/ /) you're telling ncl to not copy
>> any metadata.
>> That notation is only useful when you have already defined the variable in
>> the file.
>> So
>> ncdf->temp2 = temp2
>> Will write out the contents and the metadata/attributes with the variable.
>>
>>
>>
>> ; filevarattdef(ncdf,"temp2",temp2)
>>
>> ncdf->TFLAG =TFLAG
>> ncdf->temp2 = (/temp2/)
>>
>> On Fri, Aug 28, 2015 at 1:33 PM, Luo, Chao <chao.luo at eas.gatech.edu>
>> wrote:
>>>
>>> Dear All,
>>>
>>> I am working on temperature interpolation from 36km to 4km by
>>> ESMF_regrid, and want to write out the TEMP2 at new resolution,
>>> unfortunately the temp2 attribute didn't write out in netcdf file, I have
>>> tried to define TEMP2 attribute by:
>>>
>>> temp2!0 = "time"
>>> temp2!1 = "lay"
>>> temp2!2 = "lat"
>>> temp2!3 = "lon"
>>> temp2 at long_name = "TEMP2"
>>> temp2 at units = "k"
>>> temp2 at var_desc = "temperature at 2 m"
>>>
>>> and by function of filevarattdef. Both didn't work.
>>>
>>> Thanks very much for and help or suggestions!
>>>
>>> Here is the script I used:
>>>
>>>
>>> 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"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
>>>
>>> begin
>>>
>>> diri_36 = "/data11/cluo/CMAQv5.0.1/scripts/mcip4.1/36km/201107/data/"
>>> diri_4 = "/data12/cluo/smoke3.6.5/data/ge_dat/"
>>> namei_36 = "GRIDCRO2D_148x112_20110701"
>>> namei_4 = "geo_em.d01_1332x1008.nc"
>>>
>>> f1 = addfile(diri_36+namei_36,"r")
>>>
>>> ;output lat/lon
>>> lat2d_36km = f1->LAT(0,0,:,:) ; (TSTEP, LAY, ROW, COL) , 36km
>>> lon2d_36km = f1->LON(0,0,:,:) ; (TSTEP, LAY, ROW, COL) , 36km
>>>
>>> printVarSummary (lon2d_36km)
>>> printMinMax(lon2d_36km,True)
>>>
>>>
>>> f2 = addfile(diri_4+namei_4,"r")
>>> ;input lat/lon
>>>
>>> lat2d_4km = f2->XLAT_M(0,:,:) ; 4km
>>> lon2d_4km = f2->XLONG_M(0,:,:) ; 4km
>>>
>>> printVarSummary (lon2d_4km)
>>> printMinMax(lon2d_4km,True)
>>>
>>> name_in = "METCRO2D_148x112_20110701"
>>>
>>> nlay = 1
>>> ntime = 25
>>>
>>> time = new((/ntime/),integer)
>>> time = ispan(0,ntime-1,1)
>>>
>>> f3 = addfile(diri_36+name_in,"r") temp2_36 = f3->TEMP2(:,:,:,:)
>>> ;(TSTEP, LAY, ROW, COL)
>>> TFLAG = f3->TFLAG(:,0,:)
>>> temp2_36 at _FillValue = 9.999e+36
>>>
>>> printVarSummary (temp2_36)
>>> printMinMax(temp2_36,True)
>>>
>>> diro = "/data11/cluo/cmaq4.6/data/post/smk31_post/"
>>> filo = "out.nc" ; input file name
>>> fl_nm_out = diro+filo
>>> system ("/bin/rm -f "+fl_nm_out) ; Remove any pre-existing file
>>> ncdf = addfile(fl_nm_out,"c")
>>>
>>> ;---Set up options for regridding
>>> Opt = True
>>>
>>> Opt at SrcGridLat = lat2d_36km ; source grid
>>> Opt at SrcGridLon = lon2d_36km
>>> Opt at DstGridLat = lat2d_4km ; destination grid
>>> Opt at DstGridLon = lon2d_4km
>>>
>>> Opt at SrcRegional = True ; Necessary if grids
>>> Opt at DstRegional = True ; are regional
>>>
>>> Opt at InterpMethod = "bilinear" ; "patch", "conserve"
>>>
>>> Opt at ForceOverwrite = True ; Optional, but recommended.
>>> Opt at PrintTimings = True ; Optional.
>>> Opt at Debug = True ; Optional
>>>
>>> temp2 = ESMF_regrid(temp2_36,Opt) ; Do the regridding
>>>
>>> temp2!0 = "time"
>>> temp2!1 = "lay"
>>> temp2!2 = "lat"
>>> temp2!3 = "lon"
>>> temp2 at long_name = "TEMP2"
>>> temp2 at units = "k"
>>> temp2 at var_desc = "temperature at 2 m"
>>>
>>> printVarSummary (temp2)
>>> printMinMax(temp2,True)
>>>
>>> ; filevarattdef(ncdf,"temp2",temp2)
>>>
>>> ncdf->TFLAG =TFLAG
>>> ncdf->temp2 = (/temp2/)
>>>
>>> end
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
More information about the ncl-talk
mailing list