#!/bin/csh
#
set echo
#
# --- extract 2-d fields from a single HYCOM archive file.
# --- configured for 28 layers.
#
# --- output can be formatted, unformatted (BINARY), .[ab] (HYCOM).
# --- or use ${E}archv2ncdf2d for netCDF output.
#
# --- output is HYCOM .a files, converted to "raw" .A files.
# --- this is an example, customize it for your datafile needs.
#
#
touch regional.depth.a regional.depth.b
if (-z regional.depth.a) then
  /bin/rm regional.depth.a
  /bin/ln -s ../topo/depth_GLBt0.72_15.a regional.depth.a
endif
if (-z regional.depth.b) then
  /bin/rm regional.depth.b
  /bin/ln -s ../topo/depth_GLBt0.72_15.b regional.depth.b
endif
#
touch regional.grid.a regional.grid.b
if (-z regional.grid.a) then
  /bin/rm regional.grid.a
  /bin/ln -s ../topo/regional.grid.a regional.grid.a
endif
if (-z regional.grid.b) then
  /bin/rm regional.grid.b
  /bin/ln -s ../topo/regional.grid.b regional.grid.b
endif
#
# --- D,y,d select the archive files.
setenv CDF_TITLE        "HYCOM GLBt0.72"
setenv CDF_INST         "Naval Research Laboratory"
#
setenv D /nexsan/people/abozec/HYCOM_CESM/GLBt0.72/BHY/expt_04.1/
setenv O /nexsan/people/abozec/HYCOM_CESM/GLBt0.72/BHY/expt_04.1/
setenv E 041

setenv Y 0009
setenv m 1
setenv d 350
setenv h 19
#
#    foreach d ( 016 046 077 107 138 168 199 229 260 290 321 351 )
    setenv d1 `echo $m | awk '{printf("%02d", $1)}'`
   setenv CDF033  ${O}/${E}_archm.${Y}_${d}_${h}_3dw.nc
   /bin/rm $CDF033

  /Net/yucatan/abozec/PACIFIC/HYCOM/hycom/HYCOM-tools/archive/src/archv2ncdf2d  <<E-o-D
${D}/archm.${Y}_${d}_${h}.a
netCDF
000	'iexpt ' = experiment number x10 (000=from archive file)
  3	'yrflag' = days in year flag (0=360J16,1=366J16,2=366J01,3-actual)
500 	'idm   ' = longitudinal array size
382	'jdm   ' = latitudinal  array size
 32	'kdm   ' = number of layers
 34.	'thbase' = reference density (sigma units)
  0	'smooth' = smooth fields before plotting (0=F,1=T)
  0	'mthin ' = mask thin layers from plots   (0=F,1=T)
  1	'iorign' = i-origin of plotted subregion
  1	'jorign' = j-origin of plotted subregion
  0	'idmp  ' = i-extent of plotted subregion (<=idm; 0 implies idm)
  0	'jdmp  ' = j-extent of plotted subregion (<=jdm; 0 implies jdm)
  0	'botio ' = bathymetry       I/O unit (0 no I/O)
 0	'flxio ' = surf. heat flux  I/O unit (0 no I/O)
 0	'empio ' = surf. evap-pcip  I/O unit (0 no I/O)
  0	'ttrio ' = surf. temp trend I/O unit (0 no I/O)
  0	'strio ' = surf. saln trend I/O unit (0 no I/O)
  0	'icvio ' = ice coverage     I/O unit (0 no I/O)
  0	'ithio ' = ice thickness    I/O unit (0 no I/O)
  0	'ictio ' = ice temperature  I/O unit (0 no I/O)
  0	'sshio ' = sea surf. height I/O unit (0 no I/O)
  0     'bkeio ' = baro. k.e.       I/O unit (0 no I/O)
  0	'bsfio ' = baro. strmfn.    I/O unit (0 no I/O)
  0	'uvmio ' = mix. lay. u-vel. I/O unit (0 no I/O)
 0	'vvmio ' = mix. lay. v-vel. I/O unit (0 no I/O)
  0	'spmio ' = mix. lay. speed  I/O unit (0 no I/O)
 0	'bltio ' = bnd. lay. thick. I/O unit (0 no I/O)
 0	'mltio ' = mix. lay. thick. I/O unit (0 no I/O)
 0	'sstio ' = mix. lay. temp.  I/O unit (0 no I/O)
 0	'sssio ' = mix. lay. saln.  I/O unit (0 no I/O)
 0	'ssdio ' = mix. lay. dens.  I/O unit (0 no I/O)
  0     'mkeio ' = mix. lay. k.e.   I/O unit (0 no I/O)
-1	'kf    ' = first output layer (=0 end output; <0 label with layer #)
  32	'kl    ' = last  output layer
  0	'uvlio ' = layer k   u-vel. I/O unit (0 no I/O)
  0	'vvlio ' = layer k   v-vel. I/O unit (0 no I/O)
  0	'splio ' = layer k   speed. I/O unit (0 no I/O)
  33	'wvlio ' = layer k   w-vel. I/O unit (0 no I/O)
  33	'infio ' = layer k   i.dep. I/O unit (0 no I/O)
  33	'thkio ' = layer k   thick. I/O unit (0 no I/O)
  0	'temio ' = layer k   temp   I/O unit (0 no I/O)
  0	'salio ' = layer k   saln.  I/O unit (0 no I/O)
  0	'tthio ' = layer k   dens,  I/O unit (0 no I/O)
  0     'keio  ' = kinetic energy   I/O unit (0 no I/O)
  0	'sfnio ' = layer k  strmfn. I/O unit (0 no I/O)
  0	'kf    ' = first output layer (=0 end output; <0 label with layer #)
E-o-D
#    set m=`expr $m + 1`
#end