[Thread Prev][Thread Next][Index]

Re: [ferret_users] Reading ascii file



Hi Pierre,

Here is a solution.
Haven't found simpler.

!####################################################################################################
def symb fileInput=($01%nbp.txt%)

def sym nblines=`{SPAWN:"cat ($fileInput) | wc -l"}`
def sym nbcolumns=`{SPAWN:"head -n 1 ($fileInput) | awk '{ print NF}'"}`

!====================================================
! Read file

def axis/x=1:($nbcolumns):1 xaxis
def axis/t=1:($nblines):1 taxis
def grid/x=xaxis/t=taxis mygrid
file/g=mygrid/column=($nbcolumns) "($fileInput)"

list v1

!====================================================
! First colum is a time axis

def axis/x=1:`($nbcolumns)-1`:1 regions
def axis/t/units=years tyears=v1[i=1]
def grid/t=tyears/x=regions myvargrid

let/title="nbp"/units="" nbp=reshape(v1[i=2:($nbcolumns)], x[g=myvargrid]*t[g=myvargrid])

list nbp
show grid nbp

shade/transpose nbp
     currently SET data sets:
    1> ./nbp.txt  (default)
 name     title                             I         J         K         L
 NBP      nbp                              1:18      ...       ...       1:60

 

yes? list nbp[i=1]
             VARIABLE : nbp
             FILENAME : nbp.txt
 18   / 18: -0.0313



Thanks for you help.
Pierre




Attachment: readAscii.jnl
Description: Binary data


[Thread Prev][Thread Next][Index]