Stage IV for James

From: Li, Yongzuo <yongzuo.li_at_nyahnyahspammersnyahnyah> Date: Tue Jun 15 2010 - 12:50:30 MDT · > Hi,

--
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================
------------------------------
Message: 6
Date: Mon, 14 Jun 2010 16:17:40 -0700
From: "Correia, James" <james.correia@pnl.gov>
Subject: Stage IV precip data: conversion error
To: ncl-talk <ncl-talk@ucar.edu>
Message-ID: <C83C06A4.69F9%james.correia@pnl.gov>
Content-Type: text/plain; charset="iso-8859-1"
Hi Everybody-
I am using NCEP Stage 4 data. I used ncl's .Grb to .nc converter.
I can view the data just fine using ncview, and it appears that I can print
the array values correctly in NCL, BUT I cannot plot the data. There is no
error message so I am struggling to understand what is wrong in order to fix
it.
Perhaps an NCL developer can take a look at the data and assist? I uploaded
the file via ftp: ST2ml2002092008.Grb
Any help is greatly appreciated.
Below is the code.
;************************************************
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/WRF_contributed.ncl"
begin
;************************************************
; open file and read in data
;************************************************
mon = (/"04","05","06","07","08","09"/)
hr =
(/"00","01","02","03","04","05","06","07","08","09","10","11","12","13","14"
,"15","16","17","18","19","20","21","22","23"/)
day =
(/"01","02","03","04","05","06","07","08","09","10","11","12","13","14","15"
,"16","17","18","19","20","21","22","23","24","25","26","27","28","29","30",
"31"/)
suma = new((/881,1121/),"float")
suma = 0.
suma@_FillValue = 1e20
sumh = new((/24,881,1121/),"float")
sumh = 0.
sumh@_FillValue = 1e20
  wks = gsn_open_wks("ps" ,"test")           ; ps,pdf,x11,ncgm,eps
y1 = 2002
m1 = 5
if(m1 .eq. 5)then
nt = 30*24
end if
do yr=y1,y1
do month=m1,m1
if(month .eq. 0)then
nd = 29
end if
if(month .eq. 1)then
nd = 30
end if
if(month .eq. 2)then
nd = 29
end if
if(month .eq. 3)then
nd = 30
end if
if(month .eq. 4)then
nd = 30
end if
if(month .eq. 5)then
nd = 29
end if
nt = -1
do d=19,19
do h=8,8
nt = nt + 1
files = "ST2ml"+yr+mon(month)+day(d)+hr(h)+".nc"
  g     = addfile(files, "r")
print(files+" ")
if(nt .eq. 0)then
xlat = g->gridlat_240
xlon = g->gridlon_240
end if
;need to do qc?
var = g->A_PCP_240_SFC_acc1h
;print(var)
;print(files+" "+max(var)+" "+min(var))
;printVarSummary(var)
suma = suma + var(:,:)
;print(suma)
sumh(h,:,:) = sumh(h,:,:) + var(:,:)
end do
end do
end do
end do
print("max: "+max(suma))
print("min: "+min(suma))
do h=0,23
print(h+" "+"max: "+max(sumh(h,:,:)))
print(h+" "+"min: "+min(sumh(h,:,:)))
end do
;************************************************
; create plots
;************************************************
  gsn_define_colormap(wks ,"BlAqGrYeOrReVi200"); choose colormap
  res                       = True             ; plot mods desired
  res@gsnMaximize           = True             ; uncomment to maximize size
  res@gsnSpreadColors       = True             ; use full range of colormap
  res@cnFillOn              = True             ; color plot desired
  res@cnLinesOn             = False            ; turn off contour lines
  res@cnLineLabelsOn        = False            ; turn off contour labels
  res@lbLabelAutoStride     = True             ; let NCL figure lb stride
  res@cnLevelSelectionMode  = "ExplicitLevels" ; explicit [unequal] cn
levels
  res@mpOutlineBoundarySets = "GeophysicalAndUSStates"
  res@cnLevels              = (/0,.254,1,8,16,32,48,64,80,96,128/)
;  WRF_map_c(gg,res,0)                           ; set map resources
;************************************************
; set True for native mapping (faster plotting)
; set to False othewise
;************************************************
   res@tfDoNDCOverlay       = False
;************************************************
; associate the 2-dimensional coordinates to variables for plotting
; only if res@tfDoNDCOverlay=False
;************************************************
  if (.not.res@tfDoNDCOverlay) then
      lat2d    =  xlat
      lon2d    =  xlon
  end if
lat2d@units = "degrees_north"
lon2d@units = "degrees_east"
;  lata = gg->XLAT(0,:,:)
;  lona = gg->XLONG(0,:,:)
;x2 = 143
;y2 = 118
;x1 = 10
;y1 = 10
;print("wrf")
;  res@mpRightCornerLonF = lona(y2,x2)
;  res@mpRightCornerLatF = lata(y2,x2)
;  res@mpLeftCornerLonF = lona(y1,x1)
;  res@mpLeftCornerLatF = lata(y1,x1)
;lets try the stereographic suggestion from the archives
res@mpProjection = "Stereographic"
res@mpLimitMode = "Corners"
res@mpCenterLonF = 255.
res@mpCenterLatF = 90.
res@cnInfoLabelOn = False
res@mpLeftCornerLatF = lat2d(0,0)
res@mpLeftCornerLonF = lon2d(0,0)
res@mpRightCornerLatF = lat2d(880,1120)
res@mpRightCornerLonF = lon2d(880,1120)
res@gsnAddCyclic = False
;res@mpGridSpacingF = 300.
;res@mpGridLatSpacingF = 3.0
;res@mpGridLonSpacingF = 3.0
suma!0 = "lat"
suma!1 = "lon"
;suma&lat2d = lat2d
;suma&lon2d = lon2d
;sumh&lat2d = lat2d
;sumh&lon2d = lon2d
;map it to the wrfgrid perhaps? Wow that takes entirely too long.
;nlat = fspan(20.,55.,360)
;nlon = fspan(-120.,-60.,611)
;print("here")
;suma2 = rcm2rgrid(xlat,xlon,suma,nlat,nlon,1)
;print("here2")
;suma1 = rgrid2rcm(nlat,nlon,suma2,lata,lona,1)
plot = gsn_csm_contour_map(wks,suma,res)
;delete(suma2)
;delete(suma1)
;suma2 = rcm2rgrid(lat2d,lon2d,sumh,nlat,nlon,1)
;suma1 = rgrid2rcm(nlat,nlon,suma2,lata,lona,1)
;do h=0,23
;plot = gsn_csm_contour_map(wks,sumh(h,:,:),res)
;end do
end
------------------------------
Message: 7
Date: Mon, 14 Jun 2010 17:40:17 -0600
From: Mary Haley <haley@ucar.edu>
Subject: Re: Argument 0 of... was coerced to the
        appropriate type        and thus will not change...
To: Yi Wang <yiwang6@hdfgroup.org>
Cc: "ncl-talk@ucar.edu Talk" <ncl-talk@ucar.edu>
Message-ID: <9A96F7E8-38E5-46C6-9699-01EFE6EE3910@ucar.edu>
Content-Type: text/plain; charset="us-ascii"
On Jun 14, 2010, at 4:46 PM, Mary Haley wrote:
> Hi Yi,
>
> This warning is coming from you calling an NCL function or procedure
> and passing a variable whose type is not what the function is
> expecting. In this case, it's the first variable being passed
> (argument 0).
>
> NCL will try to coerce your variable to the correct type, but it's
> telling you that if this is a variable that is supposed to be
> modified by the function, that's not going to happen.
>
> Usually these errors are harmless, but I highly recommend fixing
> them, as you don't want to get in the habit of having these errors
> in your code.
>
> Do you know which function or procedure is causing this error?
>
> If not, you can try this: comment out the "begin" and "end" for your
> main code (if you have them), and run your script with the "-x"
> option:
>
>    ncl -x   your_script.ncl
>
> This will echo every line as it is being executed, and you should
> then see where the warning is coming from.
>
> If the function/procedure is inside a do loop or an if statement,
> the whole thing will be echoed as one block, unfortunately, but
> hopefully you can still see where the problem is.
>
> --Mary
>
>
> On Jun 14, 2010, at 3:01 PM, Yi Wang wrote:
>
>> Hi friends,
>>
>> When I run SW_S2B01667E.20031000005-20031000147.20061182243.ncl
>> (see attachment) I got a warning:
>> Argument 0 of the current function or procedure was coerced to the
>> appropriate type and thus will not change if the function or
>> procedure modifies its value.
>>
>> What is that? Do I need to care it?
>>
>> The corresponding hdf file is below:
>> ftp://ftp.hdfgroup.uiuc.edu/pub/outgoing/veer/MAP_DATA/files/PODAAC/SW_S2B01667E.20031000005-20031000147.20061182243.hdf
>>
>> Thanks a lot.
>>
>> Best,
>> Yi Wang
>> <SW_S2B01667E.
>> 20031000005
>> -20031000147.20061182243
>> .ncl>_______________________________________________
>> 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/20100614/16dbc0f6/attachment.html
------------------------------
Message: 8
Date: Mon, 14 Jun 2010 17:43:18 -0600
From: David Brown <dbrown@ucar.edu>
Subject: Re: Stage IV precip data: conversion error
To: "Correia, James" <james.correia@pnl.gov>
Cc: ncl-talk <ncl-talk@ucar.edu>
Message-ID: <00556DE8-689F-44A7-9636-2BB8F79EC885@ucar.edu>
Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes
Hi James,
When I tried to run this script, I did get errors as follows:
(0)     check_for_y_lat_coord: Warning: Data either does not contain a
valid latitude coordinate array or doesn't contain one at all.
(0)     A valid latitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0)         'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0)     check_for_lon_coord: Warning: Data either does not contain a valid
longitude coordinate array or doesn't contain one at all.
(0)     A valid longitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0)         'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
This indicates you do not have valid coordinates set.
You have commented out the statements
  ;suma&lat2d = lat2d
  ;suma&lon2d = lon2d
Since your coordinates are 2d you cannot set these arrays as the
coordinate variables of your data variable, which is what the '&'
syntax is trying to do.
Actually "lat2d" and "lon2d" are special *attributes* recognized by
the gsn routines. I was able to get a plot by changing those 2 lines
to these:
suma&lat2d = lat2d
suma&lon2d = lon2d
  -dave
On Jun 14, 2010, at 5:17 PM, Correia, James wrote:
> Hi Everybody-
> I am using NCEP Stage 4 data. I used ncl's .Grb to .nc converter.
>
> I can view the data just fine using ncview, and it appears that I
> can print
> the array values correctly in NCL, BUT I cannot plot the data. There
> is no
> error message so I am struggling to understand what is wrong in
> order to fix
> it.
>
> Perhaps an NCL developer can take a look at the data and assist? I
> uploaded
> the file via ftp: ST2ml2002092008.Grb
>
> Any help is greatly appreciated.
>
> Below is the code.
>
> ;************************************************
> 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/WRF_contributed.ncl"
> begin
> ;************************************************
> ; open file and read in data
> ;************************************************
>
> mon = (/"04","05","06","07","08","09"/)
> hr =
> (/"00
> ","01
> ","02","03","04","05","06","07","08","09","10","11","12","13","14"
> ,"15","16","17","18","19","20","21","22","23"/)
> day =
> (/"01
> ","02
> ","03","04","05","06","07","08","09","10","11","12","13","14","15"
> ,"16
> ","17
> ","18","19","20","21","22","23","24","25","26","27","28","29","30",
> "31"/)
>
> suma = new((/881,1121/),"float")
> suma = 0.
> suma@_FillValue = 1e20
> sumh = new((/24,881,1121/),"float")
> sumh = 0.
> sumh@_FillValue = 1e20
>
>  wks = gsn_open_wks("ps" ,"test")           ; ps,pdf,x11,ncgm,eps
>
> y1 = 2002
> m1 = 5
> if(m1 .eq. 5)then
> nt = 30*24
> end if
>
> do yr=y1,y1
> do month=m1,m1
> if(month .eq. 0)then
> nd = 29
> end if
> if(month .eq. 1)then
> nd = 30
> end if
> if(month .eq. 2)then
> nd = 29
> end if
> if(month .eq. 3)then
> nd = 30
> end if
> if(month .eq. 4)then
> nd = 30
> end if
> if(month .eq. 5)then
> nd = 29
> end if
> nt = -1
> do d=19,19
> do h=8,8
> nt = nt + 1
>
> files = "ST2ml"+yr+mon(month)+day(d)+hr(h)+".nc"
>  g     = addfile(files, "r")
> print(files+" ")
> if(nt .eq. 0)then
> xlat = g->gridlat_240
> xlon = g->gridlon_240
> end if
>
> ;need to do qc?
>
>
> var = g->A_PCP_240_SFC_acc1h
> ;print(var)
> ;print(files+" "+max(var)+" "+min(var))
> ;printVarSummary(var)
> suma = suma + var(:,:)
> ;print(suma)
>
> sumh(h,:,:) = sumh(h,:,:) + var(:,:)
>
> end do
> end do
> end do
> end do
>
> print("max: "+max(suma))
> print("min: "+min(suma))
>
> do h=0,23
> print(h+" "+"max: "+max(sumh(h,:,:)))
> print(h+" "+"min: "+min(sumh(h,:,:)))
> end do
>
>
> ;************************************************
> ; create plots
> ;************************************************
>  gsn_define_colormap(wks ,"BlAqGrYeOrReVi200"); choose colormap
>
>  res                       = True             ; plot mods desired
>  res@gsnMaximize           = True             ; uncomment to
> maximize size
>  res@gsnSpreadColors       = True             ; use full range of
> colormap
>  res@cnFillOn              = True             ; color plot desired
>  res@cnLinesOn             = False            ; turn off contour lines
>  res@cnLineLabelsOn        = False            ; turn off contour
> labels
>  res@lbLabelAutoStride     = True             ; let NCL figure lb
> stride
>  res@cnLevelSelectionMode  = "ExplicitLevels" ; explicit [unequal] cn
> levels
>  res@mpOutlineBoundarySets = "GeophysicalAndUSStates"
>  res@cnLevels              = (/0,.254,1,8,16,32,48,64,80,96,128/)
>
> ;  WRF_map_c(gg,res,0)                           ; set map resources
> ;************************************************
> ; set True for native mapping (faster plotting)
> ; set to False othewise
> ;************************************************
>   res@tfDoNDCOverlay       = False
>
> ;************************************************
> ; associate the 2-dimensional coordinates to variables for plotting
> ; only if res@tfDoNDCOverlay=False
> ;************************************************
>  if (.not.res@tfDoNDCOverlay) then
>      lat2d    =  xlat
>      lon2d    =  xlon
>  end if
> lat2d@units = "degrees_north"
> lon2d@units = "degrees_east"
>
> ;  lata = gg->XLAT(0,:,:)
> ;  lona = gg->XLONG(0,:,:)
> ;x2 = 143
> ;y2 = 118
> ;x1 = 10
> ;y1 = 10
> ;print("wrf")
> ;  res@mpRightCornerLonF = lona(y2,x2)
> ;  res@mpRightCornerLatF = lata(y2,x2)
> ;  res@mpLeftCornerLonF = lona(y1,x1)
> ;  res@mpLeftCornerLatF = lata(y1,x1)
>
> ;lets try the stereographic suggestion from the archives
> res@mpProjection = "Stereographic"
> res@mpLimitMode = "Corners"
> res@mpCenterLonF = 255.
> res@mpCenterLatF = 90.
> res@cnInfoLabelOn = False
> res@mpLeftCornerLatF = lat2d(0,0)
> res@mpLeftCornerLonF = lon2d(0,0)
> res@mpRightCornerLatF = lat2d(880,1120)
> res@mpRightCornerLonF = lon2d(880,1120)
> res@gsnAddCyclic = False
> ;res@mpGridSpacingF = 300.
> ;res@mpGridLatSpacingF = 3.0
> ;res@mpGridLonSpacingF = 3.0
>
> suma!0 = "lat"
> suma!1 = "lon"
> ;suma&lat2d = lat2d
> ;suma&lon2d = lon2d
> ;sumh&lat2d = lat2d
> ;sumh&lon2d = lon2d
>
> ;map it to the wrfgrid perhaps? Wow that takes entirely too long.
> ;nlat = fspan(20.,55.,360)
> ;nlon = fspan(-120.,-60.,611)
> ;print("here")
> ;suma2 = rcm2rgrid(xlat,xlon,suma,nlat,nlon,1)
> ;print("here2")
> ;suma1 = rgrid2rcm(nlat,nlon,suma2,lata,lona,1)
>
> plot = gsn_csm_contour_map(wks,suma,res)
>
> ;delete(suma2)
> ;delete(suma1)
> ;suma2 = rcm2rgrid(lat2d,lon2d,sumh,nlat,nlon,1)
> ;suma1 = rgrid2rcm(nlat,nlon,suma2,lata,lona,1)
>
> ;do h=0,23
> ;plot = gsn_csm_contour_map(wks,sumh(h,:,:),res)
> ;end do
>
> end
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
------------------------------
Message: 9
Date: Tue, 15 Jun 2010 11:12:59 -0300
From: Mateus Teixeira <mateus.teixeira@gmail.com>
Subject: Re: [ncl-install] Moving NCL installation
        directory
To: Helen Parish <hparish@ess.ucla.edu>, ncl-talk@ucar.edu
Message-ID:
        <AANLkTil033r1wdnhgvpXjjcZ18tF4mTJylQ2-3EX3pry@mail.gmail.com>
Content-Type: text/plain; charset="iso-8859-1"
Hi Helen,
If I have understood correctly, you want to move to another directory
because of memory. I think that the memory that you refer is about disk
space. As soon as I know, the directory of the ncl files does not affect it.
If you are experiencing problems with your disk space executing your scripts
in your home directory, try to move only your data and script files to
another directory, keeping the ncl files in the initial location, i.e., the
locations in your old .cshrc.
If you are experiencing problems with memory, I think that you should talk
with your system administrator, since the location of a program does not
interfere in the memory usage.
I hope have helped you.
Best regards,
Mateus
2010/6/15 Helen Parish <hparish@ess.ucla.edu>
> Mateus,
>
> Here are the lines from the old .cshrc file which work:
>
> setenv NCARG_ROOT /Users/helenparish/nc
> set path=(/Users/helenparish/nc/bin $path)
> set path=(/Users/helenparish/bin $path)
>
> Here are the lines from the new .cshrc:
>
> setenv NCARG_ROOT /Volumes/"HP Mac Disk 2"/nc
> set path=(/Volumes/"HP Mac Disk 2"/nc/bin $path)
> set path=(/Volumes/"HP Mac Disk 2"/bin $path)
>
> Note that the directory /bin simply contains the following:
>
> g95 -> /Users/helenparish/g95-install/bin/powerpc-apple-darwin6.8-g95
>
> I was not sure exactly how to alter that line when moving directories, even
> though I had copied the
> g95 directory to the new disk (as well as leaving a copy in the old home
> directory) so I also tried the following:
>
> setenv NCARG_ROOT /Volumes/"HP Mac Disk 2"/nc
> set path=(/Volumes/"HP Mac Disk 2"/nc/bin $path)
> set path=(/Users/helenparish/bin $path)
>
> keeping the /bin in its original location.
>
> I need to move ncl so that it will run from another directory, because it
> does not have sufficient memory
> when running from the home directory (especially given the sizes of the
> files that i need to process using
> ncl). I am not sure what happens even if I am still attempting to run g95
> from the home directory given that I cannot
> alter the g95 -> statement above, or whether this matters in terms of how
> much memory is used in the home
> directory. I am not sure if this has anything to do with the reason that
> the ncl command is not understood though. The
> command ncl is not recognised as soon as I type it in.
>
> Any help would be appreciated.
>
> Thanks,
> Helen.
>
>
> On Jun 12, 2010, at 3:07 PM, Mateus Teixeira wrote:
>
> Helen,
>
> Can you provide the lines in both new and old .cshrc files?
>
> One alternative is to set NCARG_ROOT and PATH in the /etc/profile or other
> general configuration file available in your system. To avoid problems with
> different users I usually set these environmental variables in the
> /etc/profile of the system.
>
> Mateus
>
>
>
> 2010/6/12 Helen Parish <hparish@ess.ucla.edu>
>
>> I have been using NCL successfully for a while, with it installed in
>> my home (login) directory. I now want to be able to run the code from
>> another directory on the same computer.
>>
>> I have moved the /nc directory to the disk where I want to use the
>> code from. In my home directory I have changed the statements in
>> the .cshrc file beginning "setenv NCARG_ROOT"  and "set path =" , to
>> contain the new directory paths.
>> However, the command "ncl" is no longer recognised. If I move the /nc
>> directory back to its original location and replace the new .cshrc
>> with the old one then ncl works again.
>>
>> Does anyone know if there is something else I need to do in order to
>> get ncl working from a new directory ?.
>>
>> Thanks,
>> Helen.
>> _______________________________________________
>> ncl-install mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-install
>>
>
>
>
> --
> Mateus da Silva Teixeira
>
> Registered Linux User #466740 (http://counter.li.org/)
>
> --
> This message has been scanned for viruses and
> dangerous content by *MailScanner* <http://www.mailscanner.info/>, and is
> believed to be clean.
>
>
>
--
Mateus da Silva Teixeira
Registered Linux User #466740 (http://counter.li.org/)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100615/a686797b/attachment.html
------------------------------
Message: 10
Date: Tue, 15 Jun 2010 10:13:37 -0500
From: Yi Wang <yiwang6@hdfgroup.org>
Subject: Re: Argument 0 of... was coerced to the
        appropriate type        and thus will not change...
To: Mary Haley <haley@ucar.edu>, Shea Dennis <shea@ucar.edu>
Cc: "ncl-talk@ucar.edu Talk" <ncl-talk@ucar.edu>
Message-ID:
        <AANLkTilmhov550HaAvkyOxar7-Cz8qSut0uZBX2Q7O_d@mail.gmail.com>
Content-Type: text/plain; charset="iso-8859-1"
Hi Mary and Dennis,
Great! Now the warning message does not show up again. It looks that
"printVarSummary(variable) " is really useful. Thanks a lot. :)
Best,
Yi
On Mon, Jun 14, 2010 at 6:40 PM, Mary Haley <haley@ucar.edu> wrote:
>
>  On Jun 14, 2010, at 4:46 PM, Mary Haley wrote:
>
>  Hi Yi,
>
> This warning is coming from you calling an NCL function or procedure and
> passing a variable whose type is not what the function is expecting. In this
> case, it's the first variable being passed (argument 0).
>
> NCL will try to coerce your variable to the correct type, but it's telling
> you that if this is a variable that is supposed to be modified by the
> function, that's not going to happen.
>
> Usually these errors are harmless, but I highly recommend fixing them, as
> you don't want to get in the habit of having these errors in your code.
>
> Do you know which function or procedure is causing this error?
>
> If not, you can try this: comment out the "begin" and "end" for your main
> code (if you have them), and run your script with the "-x" option:
>
>    ncl -x   your_script.ncl
>
> This will echo every line as it is being executed, and you should then see
> where the warning is coming from.
>
> If the function/procedure is inside a do loop or an if statement, the whole
> thing will be echoed as one block, unfortunately, but hopefully you can
> still see where the problem is.
>
> --Mary
>
>
>  On Jun 14, 2010, at 3:01 PM, Yi Wang wrote:
>
>  Hi friends,
>
> When I run SW_S2B01667E.20031000005-20031000147.20061182243.ncl (see
> attachment) I got a warning:
> Argument 0 of the current function or procedure was coerced to the
> appropriate type and thus will not change if the function or procedure
> modifies its value.
>
> What is that? Do I need to care it?
>
> The corresponding hdf file is below:
>
> ftp://ftp.hdfgroup.uiuc.edu/pub/outgoing/veer/MAP_DATA/files/PODAAC/SW_S2B01667E.20031000005-20031000147.20061182243.hdf
>
> Thanks a lot.
>
> Best,
> Yi Wang
> <SW_S2B01667E.20031000005-20031000147.20061182243.ncl>
> _______________________________________________
> 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/20100615/c8ad38a0/attachment.html
------------------------------
_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
End of ncl-talk Digest, Vol 79, Issue 31
****************************************
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk