スプライン補間 (Spline Interpolation)

[Tue Mar 15 11:29:37 JST 2011]

[am@aofd30 processor=x86_64]

[~/11.Work10/15.HeatBudjet/51.jcope2/32.NCEP_FLUX/23.Spline_Test]

$ ls

griddeg_cut.asc input/ output/ para_org.txt ps/ sp*

griddeg_cut_org.asc obj/ para.txt plot/ run.sh* src/

$ cd src

$ ls

0.Readme.txt makefile splnd2.f srcbackup/

griddeg_cut.asc sp.f90 splnd2.f90

$ make

if [ ! -d ../obj ]; then \

mkdir -p ../obj ; \

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/sp.o sp.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/splnd2.o splnd2.f

ifort -o ../sp ../obj/sp.o ../obj/splnd2.o -module ../obj

real 0m1.049s

user 0m0.919s

sys 0m0.110s

$ run.sh

Program sp starts.

Input: input/lhtfl.sfc.gauss.2000.raw.asc

Output: output/lhtfl.sfc.gauss.2000.spl.asc

sig= 1.00000

np= 192

......

nx= 237

ny= 307

Done program sp.

$ cd plot

$ ls ../output

air2m.gauss.2000.spl.asc lhtfl.sfc.gauss.2000.spl.asc

$ plflux.sh lhtfl.sfc.gauss.2000.spl.asc

Bash script ./plflux.sh starts.

Input : ../output/lhtfl.sfc.gauss.2000.spl.asc

Output : ../ps/lhtfl.sfc.gauss.2000.spl.ps

115.042 134.709 18.7917 44.2917 1.70472 255.712

115.042

134.709

18.7917

44.2917

1.70472

255.712

pstext: Record 6 is incomplete (skipped)

Done ./plflux.sh

プログラム一式: spline.tar.gz

以前のバージョン。改行コードがwindows用になっている。