[idd-brasil] Satelite+pmsl - erro na proj.

• To: idd-brasil@xxxxxxxxxxxxxxxx
• Subject: [idd-brasil] Satelite+pmsl - erro na proj.
• From: "Rafael C.Borges" <rafael.ciram@xxxxxxxxx>
• Date: Wed, 19 Sep 2007 09:30:31 -0300

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/

• Follow-Ups:
  • Re: [idd-brasil] Satelite+pmsl - erro na proj.
    • From: Rafael C.Borges