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

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

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/

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