Script carta superficie

[waldenio almeida]

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
  TITLE 
r

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
==============================================================================