#!/bin/ksh
DIR=/usr_home/rozante/orog/layers/20km_50L_op
set -x
# 
# This script takes an existing deta distribution w/50 mb top and
# redistributes the layers as the top is raised to 25 mb
#
#   FIND THE ETA VALUE ON TOP OF A LOWEST ETA LAYER
#   OVER THE OCEAN THAT IS 20 METERS THICK IN THE
#   STANDARD ATMOSPHERE

cd $DIR
F_FF12=deta.38
F_FF14=hold38
export F_FF12 F_FF14

#f90 -o lyrs.new.x  lyrs.new.f
#lyrs.new.x > holdl

F_FF12=deta.38
F_FF23=eta.co2
F_FF24=deta.50.bin
export F_FF12 F_FF23 F_FF24

#f90  -o insert.lyrs.x insert.lyrs.f
#insert.lyrs.x > hold50

# 
# Get files needed by gesprep, etc with new deta 
# distribution
#

#f90 -float2  -o deltaeta.x  deltaeta.f

F_FF11=deta.50.bin
F_FF51=si2data.l50
F_FF52=deta.50
export F_FF11 F_FF51 F_FF52
#deltaeta.x > deltaeta.out

f90 -float2  -o gfdlco2.x gfdlco2.f

E2D=/usr_home/rozante/orog/layers
F_FF23=eta.co2                     ### -Fcos -Nibm -Cascii  fort.23
F_FF61=${E2D}/tr49t85.ascii        ### -Fcos -Cebcdic -Nibm fort.61
F_FF62=${E2D}/tr49t67.ascii        ### -Fcos -Cebcdic -Nibm fort.62
F_FF63=${E2D}/tr67t85.ascii        ### -Fcos -Cebcdic -Nibm fort.63
F_FF66=outco2 
export F_FF23 F_FF61 F_FF62 F_FF63 F_FF66 
gfdlco2.x < modtop.co2 > gfdlco2.out

mv outco2 radco2.50_25mb
\rm ru*e.e*  *e.o* *.x