Dear Anley,
thank you for the reply, but it doesn't work yet. I have tried with /ASPECT and changing the relation between x and y axis. Heigth is in meters. This is the whole script: Thank you. Carlos use onda3.nc let /title="full geopotential at layer boundaries" /unit="m" geopot_model = (phb+ph)/9.81 let staggered_z=z[gz=bottom_top_stag]-0.5 define axis /z my_staggered_z=staggered_z let geopot_stag=geopot_model[gz=my_staggered_z@asn] let /title=”Heigth” /unit=”m” geopotential_f=geopot_stag[gz=bottom_top] !------------------------------------------------------------ PPL axlsze 0.16,0.16 PPL labset 0.16,0.16,0.16,0.16 PPL shakey 0.16 PPL title temperature PPL SHASET RESET !-------------------------------------------------------------- let /title="viento" viento=((u^2+v^2)^0.5) let direccion=atan(u/v) let my_t='T' let theta=my_t+300 let ptot=(PB+P)/100 let temp=(theta*((ptot/1000)^(2/7)))-273.15 cancel mode logo set region /j=20:70/i=60/l=600 REPEAT /l=600:900 (vector /ASPECT /vlim=400:4000 u,v,xlat,geopotential_f); On 01/10/2013 08:31 PM, Ansley Manke wrote: Carlos, -- Carlos Román Cascón Departamento de Geofísica y Meteorología Facultad de Ciencias Físicas Universidad Complutense de Madrid Madrid, Spain. |