#
# Runscript for pom08:
#
# Note: "params" replaces original values set in pom08.f
#
# ------------------ PARAMETER FILES FOR EACH CASE -----------
#
# -- SEAMOUNT TEST CASE:
#
#echo '      parameter(im=65, jm=49, kb=21)' >  grid
#echo '      iproblem=  1                  ' >  params
#echo '      nwad=      0                  ' >> params
#echo '      nsmolar=   0                  ' >> params
#echo '      hc=        0.05               ' >> params
#echo '      wadsmoth=  0.0                ' >> params
#echo '      dte=       6.                 ' >> params
#echo '      isplit=    30                 ' >> params
#echo '      hhi0=      20.0               ' >> params
#echo '      days=      50.0               ' >> params
#echo '      prtd1=      1.0               ' >> params
#echo '      horcon=    0.1                ' >> params
#echo '      aam_init=  50.e0              ' >> params
#echo '      kl1=       2                  ' >> params
#echo '      kl2=       kb-1               ' >> params
#
#
# -- BOX TEST CASE:
#
#echo '      parameter(im=65, jm=49, kb=21)' >  grid
#echo '      iproblem=  2                  ' >  params
#echo '      nwad=      0                  ' >> params
#echo '      nsmolar=   0                  ' >> params
#echo '      hc=        0.05               ' >> params
#echo '      wadsmoth=  0.0                ' >> params
#echo '      dte=       6.                 ' >> params
#echo '      isplit=    30                 ' >> params
#echo '      hhi0=      20.0               ' >> params
#echo '      days=      50.0               ' >> params
#echo '      prtd1=      1.0               ' >> params
#echo '      horcon=    0.1                ' >> params
#echo '      aam_init=  50.e0              ' >> params
#echo '      kl1=       2                  ' >> params
#echo '      kl2=       kb-1               ' >> params
#
# -- COAST TEST CASE
#	(from pom2k - not tested, but should be ok)
#
#echo '      parameter(im=41, jm=61, kb=16)' >  grid
#echo '      iproblem=  3                  ' >  params
#echo '      nwad=      0                  ' >> params
#echo '      days=      0.50               ' >> params
#echo '      prtd1=     0.25               ' >> params
#echo '      dte=      12.                 ' >> params
#cp IC.dat fort.40            # Initial Condition file
#
# -- WADSEAMOUNT TEST CASE (initial cond. calc. in pom08.f)
#
 echo '      parameter(im=65, jm=49, kb=21)' >  grid
 echo '      iproblem=  41                 ' >  params
 echo '      nwad=      1                  ' >> params
 echo '      nsmolar=   0                  ' >> params
 echo '      hc=        0.05               ' >> params
 echo '      wadsmoth=  0.05               ' >> params
 echo '      dte=       12.                ' >> params
 echo '      isplit=    10                 ' >> params
 echo '      mode=       3                 ' >> params
 echo '      hhi0=      20.0               ' >> params
 echo '      days=      50.0               ' >> params
 echo '      prtd1=     0.125              ' >> params
 echo '      horcon=    0.1                ' >> params
 echo '      aam_init=  50.e0              ' >> params
 echo '      kl1=       2                  ' >> params
 echo '      kl2=       kb-1               ' >> params
#
# ------------------ COMPILE & RUN ----------------------------
#
 ifort pom08.f -o a.out /usr/local/lib/libnetcdf.a
#
# ------------------ OUTPUT FILES -----------------------------
#
a.out > pom08.out        # printout file
#
#       pom08.nc         # netCDF file
#
#Run-timing info: ~30min per 50d run on 1 Intel processor
#1825.976u 53.476s 32:34.44 96.1%        0+0k 0+0io 78pf+0w
#