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

Monday, January 11, 2010

WRF training


Having training on January's weekend. Saturday and Sunday morning. WRF(Weather Research and Forecasting) and MM5 (Mesoscale Modelling).. Nice one. On weekend, because the trainers are from Malaysian Meteorological Dept, and they are working on weekdays.. anyhow, i'm looking forward for the end result of this training.. then i can start doing modelling! wow...can't wait (even though i feel lousy right now..aarrgghh!)

and currently doing training, accessing linux server system using windoes, xmanager. cool~

Tuesday, December 22, 2009

Current Doing


Yeah! Finally i can update this blog with some useful info by me..ha ha
Currently i'm doing pentad-climatology for Sea Level Pressure, and try to concentrate on longitude range 60oE to 140oE. and now doing some time series to see the pattern of the Siberian high. still in data-extraction stage. i hope can go so much faster. please 'distraction', leave me alone!
.......still having errors! AAARGHH!!
me: fast, faster!!