Subversion Repositories lagranto.um

Rev

Go to most recent revision | Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
16 michaesp 1 
#!/bin/csh
2 
#
3 
# Script to create secondary data files (S-files) from P-files
4 
# H.Wernli	1. April 96
5 
 
6 
unalias cd
7 
 
8 
echo "file in p2s.csh" $1
9 
 
10 
if ( $#argv < 2 ) then
11 
  echo "Usage: p2s p_filename mode"
12 
  echo "       if mode=2: put PS, TH, PV, RH, THW, DHR, DPVR on S-files"
13 
  echo "       if mode=5: put PS, TH, PV, RH, THW, DTHWDP on S-files"
14 
  echo "       if mode=6: put PS, TH, RH on S-files"
15 
  echo "       if mode=1: put PS, TH, RH, THE, THW on S-files"
16 
  echo "       if mode=4: put PS, TH, PV, RH, DTHDP, DTHWDP, THW, VO, RI on S-files"
17 
  exit 1
18 
endif
19 
 
20 
set file=$1
21 
set mode=$2
22 
 
23 
set gva=$LAGRANTO/goodies/getvars
24 
 
25 
set inputfile=fort.9
26 
 
27 
if ( ! -f $file ) then
28 
  echo "*** error: P-file $file does not exist ***"
29 
  exit 1
30 
endif
31 
 
32 
# check whether Z is already defined on P-file (important for mode=10)
33 
 
34 
set z=`$gva $file | grep Z | wc -c`
35 
if ( $z == 0 ) then
36 
  set zdef=0
37 
else
38 
  set zdef=1
39 
endif
40 
echo zdef $zdef
41 
 
42 
# determine the qmode
43 
 
44 
set qd=`$gva $file | grep QD | wc -c`
45 
if ( $qd == 0 ) then
46 
  set qmode=1
47 
else
48 
  set qmode=2
49 
endif
50 
 
51 
# prepare inputfile
52 
 
53 
\rm -f $inputfile
54 
touch $inputfile
55 
 
56 
#set date=`echo $file | cut -c 2-12`
57 
set date=`echo $file | cut -c 2-14`
58 
 
59 
echo $date >> $inputfile
60 
echo $mode >> $inputfile
61 
echo $qmode >> $inputfile
62 
if ( $mode == 10 ) echo $zdef >> $inputfile
63 
 
64 
# execute program ptos
65 
 
66 
$LAGRANTO/ppconvert/ptos
67 
 
68 
#\rm -f $inputfile