[ncl-talk] Sub: Array
Adv
advita6 at gmail.com
Sat Oct 31 18:59:50 MDT 2015
Hi NCL user,
I would like to write final average corresponding to two different stations
to single array.
Could someone help me to accomplish this?
Thank you,
Adv
fkey = "minnesota"
strs=asciiread("minnesota50.csv" , -1, "string")
;***********************************************
; user specified YYYYMM and VARIABLE
;***********************************************
yyyymm = 197301
fldstr = "PRCP" ; case sensitive
con = 1.00
dat=new((/1/),float)
;***********************************************
; match all data lines with specified YYYYYY
;***********************************************
seldata = str_match_ic_regex(strs, tostring(yyyymm)) ; 6.3.0
; seldata = str_match_ic_regex(strs, 1973) ; 6.3.0
; print(seldata)
print("=====")
n=(/"THIEF","PARK"/)
print(n)
do i =0, dimsizes(n)-1
strkey = n(i)
print(strkey)
;strkey= "THIEF"
;print(strkey)
;return
; find all lines (rows) which contain the key string
seldata1 = str_match_ic_regex(seldata, strkey) ; 6.3.0 onward
print(seldata1)
print(dimsizes(seldata1))
;***********************************************
; write selected data to ascii file
;***********************************************
seldir = "./"
selfil = fkey+"."+yyyymm+".csv"
system("/bin/rm -f "+seldir+selfil)
asciiwrite(seldir+selfil, seldata1) ; all data for selected yyyymm
;***********************************************
; Which field matches
;***********************************************
FIELDS = (/"STATION","STATION_NAME","ELEVATION","LATITUDE","LONGITUDE" \
,"DATE","MDPR","MDSF","DAPR","DASF","PRCP","SNWD","SNOW" \
,"TMAX","TMIN","TOBS"
/)
nfield = ind(FIELDS.eq.fldstr) + 1 ; field to be extracted
print("fldstr="+fldstr+" corresponds to field number "+nfield)
print("=====")
;***********************************************
; extract the values
;***********************************************
sdat = str_get_field(seldata1 , 2, ",")
;*********Monthly average********************************
data = tofloat( str_get_field(seldata1 , nfield, ",") )*con
data at _FillValue = -9999*con
print(sdat+" : "+data)
print("=====")
;print(data)
;return
tavg=avg(data)
end do
print(tavg)
return
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151031/363122da/attachment.html
More information about the ncl-talk
mailing list