Se coloco proj=SAT no gpmap e gdplot, o dado do modelo não aparece.
Alguém sabe porque isso ocorre?
Obrigado!
2007/9/19, Rafael C.Borges <rafael.ciram@xxxxxxxxx>:
>
> Olá pessoal ...
> Alguém pode me ajudar!
> Ao gerar uma figura "imagem de satélite+PMSL(AVN)", tenho problema com as
> projeções.
>
> Segue abaixo o script:
>
> ##################################################
> #!/bin/csh -f
> source ~gempak/Gemenviron
> #
> # GOES-12 variable set in Gemenviron
> # MODEL - GFS(AVN) - PMSL
> #
> ## Datas
> set YMD=`date -u '+%Y%m%d'`
> set YY=`date -u '+%Y'`
> set MM=`date -u '+%m'`
> set DD=`date -u '+%d'`
> set HH=`date -u '+%H'`
> set min=`date -u '+%M'`
> ##
> #
> set FONTE=/data/ldm/gempak/images/sat/GOES-12/4km/IR
> set DESTINO=/home/meteop/i/meteorologia/satelite_metar
> mkdir $DESTINO
> cd $FONTE
> #
> set CURRENT_SAT=`ls IR_* | tail -1`
> set STATION=dset
> set SAIDA=${DESTINO}/porraaaaaaaa${CURRENT_SAT}'Z'.gif
> set ARQSAIDA = 'gif|'${SAIDA}'|1500;1300'
> set GAREA='-42;-75;-18;-38'
>
> if (-e ${SAIDA}) then
> rm -f ${SAIDA}
> endif
> #if(-e sat.gif ) then
> # rm -f sat.gif*
> #endif
>
> # generate gifs in /meteorologia/SatMetar
> cd $DESTINO
>
> foreach FILE ($CURRENT_SAT)
> ############ IMAGEM SATELITE GOES12#####################
> gpmap << EOF
> MAP = 1
> DEVICE = ${ARQSAIDA}
> SATFIL = /data/ldm/gempak/images/sat/GOES-12/4km/IR/$FILE
> TEXT = 1.5/2/2/hw
> PROJ = SAT
> GAREA = ${GAREA}
> CLEAR = y
> LATLON = 0
> TITLE = 1/-0.1/GOES12-${CURRENT_SAT}Z \ EPAGRI/CIRAM
> \$mapfil = hipowo.cia+histbr.uni
> r
>
> e
> EOF
> echo " "
> ############### MODEL GFS(AVN) ########################
> gdplot << EOF
> GDFILE = /data/ldm/gempak/model/gfs/2007091906_thin.gem
> GDATTIM = LAST
> GLEVEL = 0
> GVCORD = NONE
> PANEL = 0
> SKIP = 0
> SCALE = 999
> GFUNC = PMSL
> CTYPE = C
> CONTUR = 0
> CINT = 0
> LINE = 20/1/2/2
> FINT = 0
> FLINE = 10-20
> HILO = 26;3/A#;B#/1018-1070;900-1012//30;30/y
> HLSYM =
> CLRBAR =
> GVECT =
> WIND =
> REFVEC =
> TITLE = 1
> TEXT = 1
> CLEAR = NO
> GAREA = ${GAREA}
> IJSKIP =
> PROJ = SAT
> MAP = 1
> MSCALE = 0
> LATLON = 0
> DEVICE = ${ARQSAIDA}
> STNPLT =
> r
>
> e
> EOF
> echo " "
> #### Apaga arquivos nts ####
> rm -f last.nts gemglb.nts
> rm -f LAST.NTS GEMGLB.NTS
> gpend
> end
> cleanup
> exit
> #####################################################
>
> Obrigado pela atenção.
>
> Att.
> Rafael
> --
> Rafael Censi Borges
> Epagri/Ciram
> +55 (48)3239-8064/8066
> http://ciram.epagri.rct-sc.br/
