var/obsproc/MAP_plot/Map.KMA.csh

   1 #! /bin/csh
   2 ##################################################################
   3 #  This deck is used to produce a map of mesoscale observations
   4 #  sites, and produce the data files for 4DVAR.
   5 #
   6 #                                    Yong-Run Guo
   7 #                                     12/06/2001
   8 ###################################################################
   9 #set echo
  10  unalias mv rm
  11
  12
  13  set TIME_WINDOW_MIN = '2001071323'
  14  set TIME_ANALYSIS   = '2001071400'
  15  set TIME_WINDOW_MAX = '2001071401'
  16
  17  set OBSDATA  = ../obs_gts.3dvar.2001071400.syn
  18
  19  echo " "
  20
  21  if (-e $OBSDATA) then
  22
  23  if ( -e gmeta )         rm gmeta
  24  if ( -e namelist.file ) rm namelist.file
  25
  26 cat >! namelist.file << EOF
  27
  28  &MAPBG
  29  PHIC  =  38.0 ,
  30  XLONC = 126.0,
  31  IEXP  =  0,          ; Domain expanded (1) or not (0) ?
  32  AEXP  =   120.,      ; Distance of expansion in km on each side of domain
  33  IPROJ = 'LAMCON',
  34  ;IPROJ = 'POLSTR',   ; MAP PROJECTION
  35  ;IPROJ = 'MERCAT',   ; MAP PROJECTION
  36  TRUELAT1=30.,
  37  TRUELAT2=60.,
  38  &
  39
  40  &DOMAINS
  41  MAXNES =   3,
  42  NESTIX =  171,  160,  141,  130,  194,
  43  NESTJX =  191,  178,  141,  146,  146,
  44  DIS    =  30.,  10.,   5.,  15.,   5.,
  45  NUMC   =    1,    1,   2,     3,    4,
  46  NESTI  =    1,   56,  34,    32,   30,
  47  NESTJ  =    1,   64,  67,    29,   50,
  48  &
  49
  50  &TIME_WINDOW
  51   TIME_WINDOW_MIN = '$TIME_WINDOW_MIN', ; beginning of time window
  52   TIME_WINDOW_MAX = '$TIME_WINDOW_MAX', ; end of time window
  53  &
  54
  55  &SKEW_PLOTS
  56   skewt_plot = .FALSE.,
  57   is_plot = .FALSE.,  ; .TRUE. to plot Tb19v innovation with different symbol.
  58  &
  59 EOF
  60
  61  sed -f no_comment.sed  namelist.file >! tmp.file
  62  mv tmp.file  namelist.file
  63
  64 # make clean
  65   make
  66
  67  ln -s -f $OBSDATA        fort.99
  68  ln -s -f namelist.file   fort.15
  69
  70  echo " "
  71  echo "Map.exe >&! Map.out.${TIME_ANALYSIS}"
  72        Map.exe >&! Map.out.${TIME_ANALYSIS}
  73  echo " "
  74
  75  if (-e gmeta) \
  76  echo "Generate plots in metafile gmeta.${TIME_ANALYSIS}"
  77  mv gmeta gmeta.${TIME_ANALYSIS}
  78
  79 else
  80  echo "Cannot find input file $OBSDATA"
  81 endif
  82  echo " "