put real time axis

I will rearrange this later~!!

let kk = 1
repeat/name=mm/range=1958:2014 (;\
if `mm eq 2014` then;\
use lh_oaflux_`mm`rt.nc;\
else;\
use lh_oaflux_`mm`.nc;\
endif;\
DEFINE AXIS/T=”15-Jan-1958″:”1-Sep-2014″:1/UNIT=month tannual;\
DEFINE GRID/T=tannual gannual;\
LET timestamp = T[G=gannual] * 0;\
LET lh = lhtfl[L=1:12] + timestamp;\
if `mm eq 2014` then;\
SAVE/append/FILE=test.nc lh[L=`kk`:`kk+8`];\
else;\
SAVE/append/FILE=test.nc lh[L=`kk`:`kk+11`];\
endif;\
let kk = `kk + 12`;\
can d/a;\
)

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s