Wednesday, April 14, 2010

some new idea for journal

i've got several new idea for journal writing...but , yet, i have little time to spend investigating 'em (at least until this July-after conference and also when my 'real' proposal done...)..well, since i'm under scholarship now, i feel more motivated- i hope this will last longer..heeehoooo~... don't get me wrong, the motivation can be towards a lot of thing, not only toward my research, which, my research is still the main reason i'm doing this...

if they were rejected, i will post some of the ideas here..that would be a great start and here is the best platform, i guess..

me: still under 'not enough sleep' influence.... still counting time to write something heavy in here..;p

Sunday, April 11, 2010

**update**

woah..this blog is 'spose to be a platform to write sth about earth science, my idea , etc etc...but i am too busy with everything..everything..haaaa...maybe later, after july perhaps? aor in between, i'll try to update, while updating my other blog..heh..

me: time to eat..lapar~

Saturday, April 10, 2010

WRF-Teleconnection-ObservationVsModelling

three things in my mindd....
not to forget the scholarship thingy...
and possibility of one friend come to see what happening with PRECIS
but the WRF is running and each running gonna take 17 hours (give and take) - slow processor, por favor.
and there's 11 runs to go...so, 11x17=187 hrs...yea baby!

me: tired eyes, staring at the LCD long enuf!

Wednesday, March 31, 2010

Innovation and Creativity Expo

I'll be at this expo for 3 afternoon...another time stealer of the month. anyhow, i'll try to enjoy the exhibition, since exhibition is not something new to me...and who knows, if there's a time where less people visiting, i can read something or do something...

Still, the progress, of the research....
since i cannot really relate the siberian high index ( pressure average) with Arctic Oscillation Index - correlation coefficient lower than i hope it would be (unless i made a mistake) - i've decided to look at different way..just see the average pressure at negative and positive AO. generally it gave some agreement with hypothesis, but ENSO index! gosh, i still cannot relate the SHI to ENSO with popular beliefs..argh

me: setting up booth now..

Wednesday, March 24, 2010

Focus!

Lack of them lately..urgh..
Now, 2 things in my mind...
1. teleconnection project
2. my dissertation

got my dissertation specific topic, but then, i need to finish this lagging teleconnection research. such a slooowwwwwwwwwwww...slow slow slow slow slow....sigh~

me: pls, motivation!

Tuesday, February 2, 2010

Another GrADS probs..

huh...i've came across this prob before but since it was a few weeks ago (and i was in rushing, ) i did not found what alteration i did. but then, i just found the conflict. felt so silly..hahaha..

after finish doing gs file for 5 day average, and output the DAT file and read it with CTL file, this GS file created to print the single value in text file..

*************************************************************************************
* untuk print index into text file


'reinit'
'open f:\ncep\scripts\work\slppentad\pentad_shi.ctl'


pen=1

while(pen<=2847)
'set t 'pen (forgot to put this here,before!!! argh!!)
'set x 1 144'
'define av'pen '=aave(slpp,lon=80,lon=120,lat=40,lat=60)'
pen=pen+1

endwhile


pen=1
while(pen<=2847)

'd av'pen
res=subwrd(result,4)
ret=write('f:\ncep\scripts\work\siberianhighindex.txt',res)
pen=pen+1
endwhile

*************************************************************************************

me: too much things to think, i guess??

Saturday, January 16, 2010

GrADS in LINUX for WRF

Got one problem today while scripting in GrADS, in hand-on training. we are animating and printing streamline and rainfall information.
the script is like this (shorter version):

********************************************************************************
'reinit'
'open test1.ctl'
tt=2

while (tt<74)
'c'
'set t'tt
'set gxout shaded'
'define rain = rainc+rainnc'
'set lev 1000'
'd rain'tt
'set lev 850'
'set gxout stream'
'd u;v'
'printim test'tt'.png white'
tt=tt+1
endwhile
******************************************************************
note that rain and rainnc data are available in level 1000 only.
and u;v data available in several atmospheric level (19 to be exact)

The output data works for t=2, but then the rain data became 'undefined' (for the rest of t)...
After a few time of figuring out how to solve the problem, we came with moving 'set lev 1000'
to above 'define rain=rainc+rainnc' ,

My theory is the 'define rain' data need also to be set to 'set lev 1000' for the 2nd looping, because the second execution bringing the value of level=850 (that being used for wind output) into the 'define rain' calculation..so, the defined 'rain' is now in value at 850 hPa, instead of level 1000...

so, we adjusted >>

********************************************************************************
'reinit'
'open test1.ctl'
tt=2

while (tt<74)
'c'
'set t'tt
'set lev 1000'
'set gxout shaded'
'define rain = rainc+rainnc'
!****!
'd rain'tt
'set lev 850'
'set gxout stream'
'd u;v'
'printim test'tt'.png white'
tt=tt+1
endwhile
******************************************************************


me: GrADS help is so limited especially in simple-easy problem, i'll try to write something more complicated whenever i'm in intermediate-advanced level..hihihi