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