Hi Otacilio,
It's difficult to debug the problem when it's such a complex script, without having your data.
I tried a simple script doing similar things, and the use of viewports itself doesn't seem to cause what you are seeing. Based on the example from color_vector_demo.jnl, my script sets a region, sets viewports, calls the scripts margins.jnl, color_vector.jnl, and land.jnl, and writes some labels. My test is below.
I would suggest you move the frame command and the main label outside the repeat loops. Does that change the problem? If not then gradually simplify your script to see if you can see which part of it is causing the vector arrows to disappear. What is in the script 'go brasil' ? If you comment out that command does it change anything? If you don't use viewports, but just put a pause statement in the inner loop so you see each plot as it's made, do they have vector arrows?
Here is my script.
can mode logo
set win/siz=0.5
set data monthly_navy_winds
let uu = uwnd[d=1]
let vv = vwnd[d=1]
set reg/x=160:360/y=-30:70/t=15-JUL-1989
!--define color_vector fields as in the demo
define symbol vec_cmd = VECTOR/hlimits=160:360:20/ vlimits=-30:70:10/nolab
define symbol vec_u = uu ! u var name, no [],(),*,/,+,-, ^ stuff
define symbol vec_v = vv ! v var name, no [],(),*,/,+,-, ^ stuff
define symbol vec_len = 10 ! value for /LENGTH qualifier
define symbol vec_xskip = 2 ! value for /XSKIP qualifier
define symbol vec_yskip = 2 ! value for /YSKIP qualifier
define symbol vec_arrowkey = YES ! YES or NO, for scale/key vector
define symbol vec_veckey = 325,78,,(f4.1, " m/s")! arguments to "PPL VECKEY" command
define symbol vec_lo = 0 ! low value for color levels
define symbol vec_hi = 9 ! high value for color levels
define symbol vec_delta = 0.5 ! delta value for color levels
define symbol vec_pal = rainbow ! palette file name
define symbol vec_colorkey = YES ! YES or NO, for color bar/key
define symbol vec_shakey = 1,1,0.12,2 ! arguments to "PPL SHAKEY" command
define symbol vec_otherppl = none !other ppl commands (should start with "ppl"),
! for example
! ppl xfor (f5.1, ''lonE'')
set view upper
go margins 1 0.9 0.9 1.2
go color_vector
go land
label/nouser `($PPL$XLEN)/2`, -.8, 0, 0, .2, @ASlabel1
set view lower
go margins 1 0.9 0.9 1.2
go color_vector
go land
label/nouser `($PPL$XLEN)/2`, -.8, 0, 0, .2, @ASlabel2
set view full
label/nou `($PPL$XLEN)/2`, 7.3, 0, 0, .18, @C07 Wind Speed xx mbar
frame/file=color_vec_2.gif
On 3/4/2013 4:34 AM, Otacilio Leandro wrote:
Dear ferret users,
I am trying to analise a six month wind climatlogy over South America.
For this i made a loop to make easy to make the figures for several experiments that i'm runing.
The problem is:
When combining viewport with the script color_vector, Ferret is clearing the first viewport when the second is filled and so on and the result is attached.
Thank you!
==================== My script =============
! Script to analise the wind speed at pressure levels
!go startup.jnl
set window/aspect=.8
let modo = "$1"
IF `$1 eq 200` THENlet alt = 8ELIF `$1 eq 500` THENlet alt = 12ELIF `$1 eq 700` THENlet alt = 14ELIF `$1 eq 850` THENlet alt = 15ELIF `$1 eq 925` THENlet alt = 16ELIF `$1 eq 1000` THENlet alt = 17elsesay "Inform the correct level.... quitting"exitENDIF
use uv_aabnra.pm.nc ! HPR 50%use uv_aabnxa.pm.nc ! HPR 100%use uv_aabnva.pm.nc ! HPR OFFuse uv_aabnsa.pm.nc ! CONTROLEuse uv_aabwua.pm.nc ! BIOMASS OFF 50%
let imes = {"JAN", "FEB", "MAR", "APR", "MAY", "JUN", "JUL", "AUG", "SEP", "OCT", "NOV", "DEC"}
define region/x=85W:32W/y=15N:55S/k=`alt` SAset region/@SA
def axis/calendar=360_day/t="16-JAN-1979":"16-NOV-2005":1/unit=month tax
! Several experiments
let exp_name = {"hpr_50","hpr_100","hpr_off","control","biomass_off"}
repeat/range=1:5/name=a (\let var = exp_name[i=`a`] ;\let u_`var` = u_1[d=`a`,gt=tax] ;\let v_`var` = v_1[d=`a`,gt=tax,gxy=u_1] ;\)
! Climatologies
use climatological_axes
repeat/range=1:5/name=a (\let var = exp_name[i=`a`] ;\let u_`var`_clim = u_`var`[gt=month_360_day@mod] ;\let v_`var`_clim = v_`var`[gt=month_360_day@mod] ;\)
go landscape3x2.jnllet posicao = {"ul6", "um6", "ur6","ll6", "lm6", "lr6"}
ppl color,6,80,80,80ppl axlsze .16 .16
let marg = "1 1.2 1.2 .4"
let variavel1 = {"HPR_50","HPR_100","HPR_OFF","CONTROL","BIOMASS_OFF"}let nome_var1 = {"HPR 50","HPR 100","HPR OFF","CONTROL","BIOMASS OFF"}
! COLOR VECTOR
define symbol vec_cmd = VECTOR/nolab/setdefine symbol vec_u = uudefine symbol vec_v = vvdefine symbol vec_len = 10define symbol vec_xskip = 1define symbol vec_yskip = 1define symbol vec_arrowkey = YES!define symbol vec_veckey = 325,78,(f4.1, " m/s")define symbol vec_veckey = 310,19,,(f4.1, " m s^-^1")define symbol vec_lo = 0define symbol vec_hi = 20define symbol vec_delta = 1define symbol vec_pal = padraodefine symbol vec_colorkey = YESdefine symbol vec_shakey = 1,1,0.12,2define symbol vec_otherppl = nonePPL AXNMTC 3,3PPL AXLINT 2,2
! Test for the 1st experiment ===> "m" loop on experiments and "a" loop for monthsrepeat/range=1:1/name=m ( \repeat/range=6:11/name=a ( \let var1 = variavel1[i=`m`] ;\let nome1 = nome_var1[i=`m`] ;\;\let mes = imes[i=`a`] ;\let pos = posicao[i=`a-5`] ;\;\ ! MONTH ;\set view `pos` ; go margins 1 .9 .9 1.2 ;\;\let uu=U_`var1`_CLIM[l=`a`] ; let vv=V_`var1`_CLIM[l=`a`] ;\go color_vector ; go brasil ; go land 1 " " 1 ;\label/nou `($PPL$XLEN)/2`, -.8, 0, 0, .2, @AS`mes` ;\;\set view full ;\label/nou `($PPL$XLEN)/2`, 7.3, 0, 0, .18, @C07 Wind Speed `modo` mbar ;\;\frame/file=color-vec-exemple.gif ;\))
exit
--
Otacilio Leandro de Menezes Neto
Instituto Nacional de Pesquisas Espaciais
Doutorando em Meteorologia
Contato: (12) 3186-9537 / (12) 8194-1455
Email alternativo: otacilio.neto@xxxxxxxxxxxxx