a bit more about dimensions

Eduardo Agosta Scarel eduardo.agosta at GMAIL.COM
Thu Apr 7 11:12:05 EDT 2005


Dear Pesquero, thanks for your help, It was really great.

But I wrote the script bellow following your suggestions, I do not
know why TE'level' and TES'level do not give different results, i mean
de variables saved only last calculation I think.

Could anybody help? What do I do wrong in the script?
Thanks again (and I appreciate a lot this forum!!)

*******************************
'set lat -60 10'
'set lon -120 -30'
'set t 1 610'

'define HR=rhum.1/100'

'define Tair=air.2+273.15'

'define rm=shum.3'

***********************
i = 1
ultimonivel=8
str="1000 925 850 700 500 600 400 300"

while (i <= ultimonivel)

level = subwrd(str,i)
'set lev 'level

'define airlift=1/(1/(Tair-55)-log(HR)/2840)+55'
'define A=0.2854*(1-0.28*10e-3*rm)'
'define B=3.376/airlift-0.00254'
'define C=1+0.81*10e-3*rm'
'define TE1=Tair*pow(1000/lev, A)'
'define TE2=exp(B*rm*C)'
'define TE'level'=TE1*TE2'

'define es=exp(-2937.4/Tair -4.29283*log(Tair) + 22.5518)'

'define presion=lev/10'
'define epsi=0.622'
'define rs=epsi*es/presion'

'define A=0.2854*(1-0.28*10e-3*rs)'

'define B=3.376/Tair-0.00254'
'define C=1+0.81*10e-3*rs'
'define TES1=Tair*pow(1000/lev, A)'
'define TES2=exp(B*rs*C)'
'define TES'level'=TES1*TES2'

i = i+1
endwhile



More information about the gradsusr mailing list