Re: read a variable from another NCL script

x seems to be defined twice in your script. You do not delete it before defining it the second time, it seems. I would fix that first and re-try.
Cheers,
Jatin.

On 04/12/2013, at 1:29 PM, Yi-Chih Huang wrote:

> Hello,
>
> In fact, helper_libs.ncl includes many NCL scripts. I replaced helper_libs.ncl with those in $NCARG_ROOT. Some errors showed up as follows. The command about the error on line 27 is "x = f->$varName(0)$" that was right when executing separately. I wonder why x should be defined when the script called a procedure?
>
> Thanks much,
>
> Yi-Chih
>
> ######
> ;load "$SysE/lib/ncl/helper_libs.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>
> ######
> yhuang@arc-env:/fs3/yhuang/scripts$ ncl t.ncl
> Copyright (C) 1995-2013 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 6.1.2
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> fatal:syntax error: line 91 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl before or near undef
> undef
> ----^
>
> fatal:syntax error: line 92 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl before or near size_array
> function size_array
> ------------------^
>
> fatal:Function identifier is defined
> fatal:Variable (x) is undefined
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 94 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
>
> warning:Dimension (0) has not been defined
>
> fatal:Dimension sizes of left hand side and right hand side of assignment do not match
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 27 in file t.ncl
>
>
>
> On Wed, Dec 4, 2013 at 10:54 AM, Rick Brownrigg <brownrig@ucar.edu> wrote:
> Hi,
>
> Its hard to pinpoint the problem exactly from what you sent, but the error messages seem pretty explicit. You might open in an editor your file " $SysE/lib/ncl/helpers/analysis_helper.ncl" and look at the line numbers mentioned in the errors to find the problems. Looks like at least two variables (dirName and aString) are referenced before they are ever defined.
>
> Hope that helps…
> Rick
>
> On Dec 3, 2013, at 6:25 PM, Yi-Chih Huang <dscpln@gmail.com> wrote:
>
>> Hello,
>>
>> I used t_nmm(ty_nmm), procedure and load to read the variable because another NCL script is in a different directory as follows. These two scripts are almost the same. I am trying to draw two lines on the same plot to see the detail. Originally these two scripts separately worked just fine without any error message. But after I changed one script to procedure, there were numerous error message showing up that never happened before. What do you think the problem is in the scripts?
>>
>> Thanks much,
>>
>> Yi-Chih
>>
>> .....
>> fatal:error in statement
>> fatal:Variable (dirName) is undefined
>> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 99 in file $SysE/lib/ncl/helpers/analysis_helper.ncl
>>
>> fatal:syntax error: line 101 in file $SysE/lib/ncl/helpers/analysis_helper.ncl before or near pp
>> procedure pp
>> -----------^
>>
>> fatal:error in statement
>> fatal:Variable (aString) is undefined
>> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 103 in file $SysE/lib/ncl/helpers/analysis_helper.ncl
>>
>> fatal:syntax error: line 105 in file $SysE/lib/ncl/helpers/analysis_helper.ncl before or near TerminateProgram
>> procedure TerminateProgram
>> -------------------------^
>>
>> fatal:error in statement
>> (0) Program Terminated
>>
>> #####
>> load "$SysE/lib/ncl/helper_libs.ncl"
>> load "/fs3/yhuang/nmm/t.ncl"
>>
>> begin
>> latS = -30.0
>> latN = 30.0
>> lonE = 30.0
>> lonW = 160.0
>>
>> dir1 = "/fs3/SysE_DB/nmm/CFSR/Monthly/Mean/"
>> yrStrt = 1981
>> yrLast = 2006
>> year = ispan(yrStrt, yrLast,1) ; integer
>> yrDir = year + "/" ; string
>> nyr = dimsizes(year)
>>
>> lvl = (/"200","500"/)
>> nlvl = dimsizes(lvl)
>> varName= (/"t","lat"/)
>>
>> dir = str_concat(dir1+yrDir(0))
>> fil = systemfunc("cd "+dir+" ; ls t_200.nc")
>> f = addfile (dir+fil(0), "r")
>> x = f->$varName(0)$
>> dims = getfilevardims(f,varName(0))
>> xs = x(5:8,{latS:latN},{lonE:lonW}) ; Indian Ocean
>> xy = dim_avg_n_Wrap( xs, 0 )
>> ds = dimsizes(xy)
>> xx = new((/nyr,nlvl,ds(0),ds(1)/),float,"No_FillValue")
>>
>> delete (fil)
>> do ny= 0,nyr-1
>> dir = str_concat(dir1+yrDir(ny))
>> fil = systemfunc("cd "+dir+" ; ls t_200.nc t_500.nc")
>> nlvl= dimsizes(fil)
>>
>> do nl= 0,nlvl-1
>> f = addfile (dir+fil(nl), "r")
>>
>> x = f->$varName(0)$
>> dims = getfilevardims(f,varName(0))
>> xs = x(5:8,{latS:latN},{lonE:lonW}) ; Indian Ocean
>> xy = dim_avg_n_Wrap( xs, 0 )
>> xx(ny,nl,:,:)= xy
>> end do
>> end do
>>
>> tl= dim_avg_n_Wrap( xx, 0 )
>> tl!0 = "level"
>> t= dim_avg_n_Wrap( tl, 0 )
>>
>> ty = dim_avg_n_Wrap( t, 1 )
>> f = addfile (dir+fil(0), "r")
>> la = f->$varName(1)$
>> dims = getfilevardims(f,varName(1))
>> lat = la({$dims(1)$|30:-30})
>>
>> t_nmm(ty_nmm)
>>
>> wks= gsn_open_wks("pdf","t_CFSR")
>> gsn_define_colormap(wks,"BlueRed") ; use the BlueRed colormap (instead of default colormap)
>>
>> res = True
>> res@gsnAddCyclic = False
>> res@mpMinLatF = -30
>> res@mpMaxLatF = 30
>> res@mpMinLonF = 30
>> res@mpMaxLonF = 160
>> res@cnFillOn = True
>> res@gsnSpreadColors = True ; use full colormap
>> res@gsnSpreadColorStart = 10 ; for BlueRed
>> res@gsnSpreadColorEnd = 250 ; for BlueRed
>> res@cnSpanFillPalette = True
>> res@cnLinesOn = False ; turn off the contour lines
>> res@cnLevelSelectionMode= "ManualLevels"
>> res@cnMinLevelValF = 238
>> res@cnMaxLevelValF = 249
>> res@cnLevelSpacingF = 0.5
>>
>> res@lbLabelBarOn = True
>> res@pmLabelBarWidthF= 0.8
>> res@lbLabelStride = 4
>> res@mpFillOn = False
>> res@mpOutlineBoundarySets = "National" ; turns on country boundaries.
>> res@gsnLeftString = "CFSR: Average Monthly Temperature between 200-500 hPa during JJAS"
>> res@cnFillDrawOrder = "Predraw"
>>
>> plot = gsn_csm_contour_map(wks,t,res)
>>
>> wks = gsn_open_wks ("pdf","tXY_CFSR")
>> resL = True
>> resL@tiMainString = "CFSR: Average Monthly Temperature during JJAS"
>> resL@xyLineColor = "Brown"
>> resR@xyLineColor = "NavyBlue"
>> plot = gsn_csm_xy2(wks,ty&lat,ty({lat|-30:30}),ty_nmm({lat|-30:30}),resL,resR)
>> end
>>
>> #####
>> procedure t_nmm (ty)
>>
>> load "$SysE/lib/ncl/helper_libs.ncl"
>>
>> begin
>> latS = -30.0
>> latN = 30.0
>> lonE = 30.0
>> lonW = 160.0
>>
>> dir1 = "/fs3/SysE_DB/nmm/NMM/Monthly/Mean/"
>> yrStrt = 1981
>> yrLast = 2006
>> year = ispan(yrStrt, yrLast,1) ; integer
>> yrDir = year + "/" ; string
>> nyr = dimsizes(year)
>>
>> lvl = (/"200","500"/)
>> nlvl= dimsizes(lvl)
>> varName= (/"temp","lat"/)
>>
>> dir = str_concat(dir1+yrDir(0))
>> fil = systemfunc("cd "+dir+" ; ls temp_200.nc")
>> f = addfile (dir+fil(0), "r")
>> x = f->$varName(0)$
>> dims = getfilevardims(f,varName(0))
>> xs = x(:,{latS:latN},{lonE:lonW}) ; Indian Ocean
>> xy = dim_avg_n_Wrap( xs, 0 )
>> ds = dimsizes(xy)
>> xx = new((/nyr,nlvl,ds(0),ds(1)/),float,"No_FillValue")
>>
>> delete (fil)
>> do ny= 0,nyr-1
>> dir = str_concat(dir1+yrDir(ny))
>> fil = systemfunc("cd "+dir+" ; ls temp_2??.nc temp_5??.nc")
>> nlvl = dimsizes(fil)
>>
>> do nl= 0,nlvl-1
>> f = addfile (dir+fil(nl), "r")
>>
>> x = f->$varName(0)$
>> dims = getfilevardims(f,varName(0))
>> xs = x(:,{latS:latN},{lonE:lonW}) ; Indian Ocean
>> xy = dim_avg_n_Wrap( xs, 0 )
>> xx(ny,nl,:,:)= xy
>> end do
>> end do
>>
>> tl= dim_avg_n_Wrap( xx, 0 )
>> tl!0 = "level"
>> t = dim_avg_n_Wrap( tl, 0 )
>>
>> ty = dim_avg_n_Wrap( t, 1 )
>> f = addfile (dir+fil(0), "r")
>> la = f->$varName(1)$
>> dims = getfilevardims(f,varName(1))
>> lat = la({$dims(1)$|30:-30}) ; get latitudes in the domain
>> printVarSummary(t)
>> printVarSummary(ty)
>> printVarSummary(lat)
>>
>> wks= gsn_open_wks("pdf","t_nmm")
>> gsn_define_colormap(wks,"BlueRed")
>>
>> res = True
>> res@gsnAddCyclic = False
>> res@mpMinLatF = -30
>> res@mpMaxLatF = 30
>> res@mpMinLonF = 30
>> res@mpMaxLonF = 160
>> res@cnFillOn = True
>> res@gsnSpreadColors = True ; use full colormap
>> res@gsnSpreadColorStart = 10 ; for BlueRed
>> res@gsnSpreadColorEnd = 250 ; for BlueRed
>> res@cnSpanFillPalette = True
>> res@cnLinesOn = False ; turn off the contour lines
>> res@cnLevelSelectionMode= "ManualLevels"
>> res@cnMinLevelValF = 238
>> res@cnMaxLevelValF = 249
>> res@cnLevelSpacingF = 0.5
>>
>> res@lbLabelBarOn = True
>> res@pmLabelBarWidthF= 0.8
>> res@lbLabelStride = 4
>> res@mpFillOn = False
>> res@mpOutlineBoundarySets = "National" ; turns on country boundaries.
>> res@gsnLeftString = "NMM: Average Monthly Temperature between 200-500 hPa during JJAS"
>> res@cnFillDrawOrder = "Predraw"
>>
>> plot = gsn_csm_contour_map(wks,t,res)
>>
>> wks = gsn_open_wks ("pdf","tXY_nmm")
>> xyres = True
>> xyres@tiMainString = "CFSR: Average Monthly Temperature during JJAS"
>> plot = gsn_csm_xy (wks,ty&lat,ty({lat|-30:30}),xyres)
>> end
>>
>>
>> On Wed, Dec 4, 2013 at 1:48 AM, Dennis Shea <shea@ucar.edu> wrote:
>>
>> A counterparts to a fortran subroutine and function
>> are an NCL procedure and function. These are both
>> discussed in the Mini-Language Manual and the NCL
>> Reference Manual.
>>
>> Many functions and procedures are in
>>
>> $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl
>> $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl
>> $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
>> $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
>>
>> %> less $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
>>
>> ===
>> WRAPIT is for interlanguage communication:
>> NCL which is written in C and Fortran
>>
>>
>>
>> On 12/3/13 12:24 AM, Yi-Chih Huang wrote:
>> Hello,
>>
>> I am trying to read a variable from another NCL script. WRAPIT does
>> not seem to support this function from the Mini-Language Manual. What is a
>> simple way to read variable from another NCL script? Is there a
>> counterpart of Fortran subroutine on NCL?
>>
>> Thanks much,
>>
>> Yi-Chih
>>
>>
>>
>> _______________________________________________
>> 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
Received on Tue Dec 3 19:38:54 2013