Script carta superficie

• To: idd-brasil@xxxxxxxxxxxxxxxx
• Subject: Script carta superficie
• From: waldenio almeida <gambi@xxxxxxxxxxxxx>
• Date: Mon, 16 Jul 2007 17:03:30 -0300

Marcelo Romero de Moraes escreveu:

Oi pessoal,

Tudo bem com vocês ? Alguém poderia mandar os scripts do curso para a lista.

Abraços à todos.

Marcelo R. de Moraes
EPAGRI/CIRAM
Florianópolis-SC
48 3239-8053

#! /bin/csh
#################################################################
#
#                  CARTAS SINOTICAS com o GEMPAK
#
# Este script gera as cartas sinoticas do Brasil em GIF
# a partir dos dados decodificados do GEMPAK
#
#
#################################################################
#
if ($1 == 'normal') then
 set myself = $0
 set mytype = $1
 set YMD=`date -u '+%Y%m%d'`
 set YY=`date -u '+%Y'`
 set MM=`date -u '+%m'`
 set DD=`date -u '+%d'`
 set HORA = $2
else if ($1 == 'manual') then
 set myself = $0
 set mytype = $1
 set YMD  = $2$3$4
 set YY = $2
 set MM = $3
 set DD = $4
 set HORA = $5
else
 echo " Usage: "
 echo "       $0 normal hh (manual yyyy mm dd hh)"
 exit
endif
#
# Seta as variaveis necessarias para localizar o GEMPAK
#
source ~gempak/Gemenviron

# Diretorio com os dados decodificados no formato GEMPAK
# mais o arquivo fonte para todas as figuras

set FONTE=/data/ldm/gempak
set SFFILE=${FONTE}/syn/${YMD}_syn.gem

# Diretorio de destino para todas as figuras

set DESTINO=$HOME
set DIRRUN=$HOME

cd $DIRRUN

#################################################################
#
#                Cria mapa com dados synop do Brasil
#
#################################################################

# Define o arquivo de saida
set SAIDA=${DESTINO}/synop${YMD}${HORA}.gif

if (-e ${SAIDA}) then
 rm -f ${SAIDA}
endif

# Define area de dados e grafico
set STATION=dset
set GAREA='-60;-90;10;-20'
set ARQSAIDA = 'gif|'${SAIDA}'|2000;2000'
#set ARQSAIDA = xw
#

gpcolor <<EOF
 colors = 101=white;1=black
 DEVICE = ${ARQSAIDA}
r

e
EOF

echo " "

gpmap << EOF
 CLEAR  = YES
 GAREA  = ${GAREA}
 PROJ   = CED
 DEVICE = ${ARQSAIDA}
 LATLON = 32/10/1/1/20;20
 \$MAPFIL=histbr.uni+base
 MAP    = 1/10/1 + 1//2
 TEXT   = 2/22/sw
 TITLE  = 1/-3/SYNOP ${HORA}Z ${DD}/${MM}/${YY}
 r

 e
EOF
echo " "

sfmap << EOF
 CLEAR  = NO
 colors = 1
 SFFILE = ${SFFILE}
 AREA   = ${STATION}
 GAREA  = ${GAREA}
 DATTIM = ${HORA}
 SFPARM = SKYC;TMPC;WSYM;PMSL;;DWPC;BRBK
 OUTPUT = xw
 PROJ   = CED
 DEVICE = ${ARQSAIDA}
 MAP    = 0
 LATLON = 0/10/1/1/20;20
 FILTER = 0.8
 TEXT   = 1.0/1/1/hw

e
EOF
echo " "

set SFFILE=${FONTE}/ship6hr/${YMD}${HORA}_ship.gem

sfmap << EOF
 CLEAR  = NO
 SFFILE = ${SFFILE}
 AREA   = ${STATION}
 GAREA  = ${GAREA}
 DATTIM = ${HORA}
 SFPARM = SKYC;TMPC;WSYM;PMSL;;DWPC;BRBK
 OUTPUT = xw
 PROJ   = CED
 DEVICE = ${ARQSAIDA}
 MAP    = 0
 LATLON = 0/10/1/1/20;20
 FILTER = 0.8
r

e
EOF

echo " "

gpend

rm -f gemglb.nts last.nts

exit(0)

# end of shell....

==============================================================================
To unsubscribe idd-brasil, visit:
http://www.unidata.ucar.edu/mailing-list-delete-form.html
==============================================================================

• References:
  • scripts
    • From: Marcelo Romero de Moraes