#!/usr/bin/ksh
#
rm *.key namelist 
ulimit -c 0
set -kx
#
MLEVDIR_1=~opdata/smhi/hirlamh44l31/model          ; export MLEVDIR_1
MLEVKEY_1=$PWD                                    ; export MLEVKEY_1
PLEVDIR_1=~opdata/smhi/hirlamh44l31/postpr         ; export PLEVDIR_1
PLEVKEY_1=$PWD                                    ; export PLEVKEY_1

#
YEAR1=`date +"%Y"`
YR1=`date +"%y"`
MO1=`date +"%m"`
DA1=`date +"%d"`
HO1=00
for LN1 in 012
do
LONO=-5.00
LANO=55.00
LOSO=10.00
LASO=60.00
cat > cross.dat << EOF
$LANO $LONO
$LASO $LOSO
EOF

#
MLEVFIL1=FC44_${YR1}${MO1}${DA1}${HO1}00+${LN1}H00M
PLEVFIL1=LARGE_${YR1}${MO1}${DA1}${HO1}00+${LN1}H00M
#
cat > namelist << *EOINFILE
 &namsec
      label       ='Hirlam OPS.'
      lulog       =6,
      logfile     ='cross.log',
      print       =.true.,
      mlevfile_1  ='${MLEVFIL1}'
      plevfile_1  ='${PLEVFIL1}'
      year_1      =$YEAR1,
      month_1     =$MO1,
      day_1       =$DA1,
      hour_1      =$HO1,
      fclen_1     =$LN1,
      rlono       =$LONO,
      rlano       =$LANO,
      rloso       =$LOSO,
      rlaso       =$LASO,
      dxplot      =18.,
      dyplot      =26.,
      pmin        =200,
      pmax        =1050,
      tab         =001,
      par         =11,33,34,39,13
      typ         =109,109,109,100,105,
      diff        =.f.,.f.,.f.,.f.,.f.,
      scale       =1.,1.,1,.1,1.,
      subtract    =0.,0.,0.,0.,0.,
      contours    =.f.,.t.,.f.,.f.,.t.,
      con_int     =5., 2., 2., 2.,5.,
      con_col     =004,001,002,003,007,
      con_thi     =4, 3,  3,  3,2,
      con_lin     =1,  1,  2,  1,2,
      shading     =.f.,.f.,.f.,.f.,.f.,
      sha_num     =7,
      sha_lis     =0.1,0.5,1.0,2.0,3.0,5.0,10.,25.,
      sha_col     =  019,021,023,025,027,029,031,
      sha_pat     =4,4,4,4,4,4,4,4,4,4,
      arrows      =.f.,.f.,.f.,.f.,
      arr_len     = 1.0,
      arr_col      = 002,
      icao        =.f.,
      mountains   =.t.,
      modlev      =.t.,
      mrkfrq      =1,
      hgtmrk      =0.2,
      hgttxt      =0.2,
 &end
*EOINFILE
#
cross_b namelist cross.ps P a4
#
rm namelist
done
 

RED_FILE_PATH=/proj/iflocal/src/metgr/coast/   ; export RED_FILE_PATH
COAST_FILE_PATH=/proj/hirfou/ulf/MetgraF_2/coast/ ; export COAST_FILE_PATH
#
cat > namelists << *EOINFILE
 &namlog
      log_print  =.true.,
      log_unit   = 7,
      log_file   ='metg.log',
 &end
 &namgeo
      scale   = 10.,
      maprep  = 01,
 &end
 &namfil
      gribdir_1  ='$MLEVDIR_1',
      gribfile_1 ='$MLEVFIL1',
      iyr1       = $YEAR1,
      imo1       = $MO1,
      ida1       = $DA1,
      iho1       = $HO1,
      iln1       = $LN1,
 &end
 &namsel
      ipar1   = 33,
      ilev1   = 10,
      ityp1   = 105,
      larrow  = .f.,
 &end
 &namval
      val_file='cross.dat',
      polygon =.t.,
      val_thi =6.0,
      val_hgt =6.0,
      point_hgt=6.0,
      val_col =2,
 &end
*EOINFILE
#
metgraf_b namelists crossline.ps L a4
#
exit
