Hello,
I am having problems with the ability to neglect values at a particular
index. I have a 1-D array containing 20 values. If any index (i) of that
array corresponds to a missing/fill value (-999), I do not want to use that
index in furture calculations. Any thoughts on how this may be accomplished?
Below is a subsection of my script. Thanks!
Sincerely,
Karen A. Kosiba
angle = new((/20/),float)
azAngle = new((/20/),float)
vel_circle = new((/20/),float)
alpha = new((/20/),float)
aVec = new((/20/),float)
bVec = new((/20/),float)
cVec = new((/20/),float)
aVecSq = new((/20/),float)
bVecSq = new((/20/),float)
vel_circle@_FillValue = -999
do i=0,19
angle(i) = (2*3.14*(i))/20
rr(i) = dim_rMid+RR*sin(angle(i))
az(i) = dim_azMid+RAZ*cos(angle(i))
azAngle(i) = az(floattoint(az(i)))
azAngle(i) = azAngle(i)*(3.14/180)
vel_circle(i) = vel(floattoint(az(i)),floattoint(rr(i)))*.01
if (vel_circle(i).neq.-999) then
aVec(i) = sin(angle(i)-azAngle(i))
bVec(i) = cos(angle(i)-azAngle(i))
cVec(i) = vel_circle(i)-15*cos(1.4-azAngle(i))
aVecSq(i) = (aVec(i))^2
bVecSq(i) = (bVec(i))^2
end if
end do
********************************
Karen A. Kosiba
Graduate Student
Purdue University, Dept. of EAS
********************************
_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Fri Feb 11 2005 - 09:18:10 MST