Commit ac5bf535 authored by agebhard's avatar agebhard

Changes for R 2.14.x

parent 1bb6e43a
Package: akima
Version: 0.5-5
Date: 2010-01-14
Version: 0.5-6
Date: 2012-01-06
Title: Interpolation of irregularly spaced data
Author: Fortran code by H. Akima
R port by Albrecht Gebhardt <albrecht.gebhardt@uni-klu.ac.at>
aspline function by Thomas Petzoldt <petzoldt@rcs.urz.tu-dresden.de>
aspline function by Thomas Petzoldt <thomas.petzoldt@tu-dresden.de>
enhancements and corrections by Martin Maechler <maechler@stat.math.ethz.ch>
Maintainer: Albrecht Gebhardt <albrecht.gebhardt@uni-klu.ac.at>
Description: Linear or cubic spline interpolation for irregular gridded data
......
useDynLib(akima)
export(aspline,interp,interpp,interp.old,interp.new)
.First.lib <- function(lib, pkg) {
if(getRversion()<="0.62")
stop("This version for R 0.62 or later")
library.dynam("akima", pkg, lib)
}
akima Waveform Distortion Data for Bivariate
Interpolation
"akima" <-
list(x = c(11.16, 24.2, 12.85, 19.85, 10.35, 24.65, 19.72, 15.91,
0, 20.87, 6.71, 3.45, 19.99, 14.26, 10.28, 4.51, 17.43, 22.8,
0, 7.58, 16.7, 6.08, 1.99, 25, 14.9, 3.22, 0, 9.66, 2.56, 5.22,
11.77, 17.25, 15.1, 25, 12.13, 25, 22.33, 11.52, 14.59, 15.2,
7.54, 5.23, 17.32, 2.14, 0.51, 22.69, 25, 5.47, 21.67, 3.31), y = c(1.24,
16.23, 3.06, 10.72, 4.11, 2.4, 1.39, 7.74, 20, 20, 6.26, 12.78,
4.62, 17.87, 15.16, 20, 3.46, 12.39, 4.48, 1.98, 19.65, 4.58,
5.6, 11.87, 3.12, 16.78, 0, 20, 3.02, 14.66, 10.47, 19.57, 17.19,
3.87, 10.79, 0, 6.21, 8.53, 8.71, 0, 10.69, 10.72, 13.78, 15.03,
8.37, 19.63, 20, 17.13, 14.36, 0.13), z = c(22.15, 2.83, 22.11,
7.97, 22.33, 10.25, 16.83, 15.3, 34.6, 7.54, 30.97, 41.24, 14.72,
10.74, 21.59, 15.61, 18.6, 5.47, 61.77, 29.87, 6.31, 35.74, 51.81,
4.4, 21.7, 39.93, 58.2, 4.73, 50.55, 40.36, 13.62, 6.43, 12.57,
8.74, 13.71, 12, 10.25, 15.74, 14.81, 21.6, 19.31, 26.5, 12.11,
53.1, 49.43, 3.25, 0.6, 28.63, 5.52, 44.08))
File added
......@@ -164,9 +164,10 @@ image (akima.smooth, main = "interp(<akima data>, *) on finer grid")
contour(akima.smooth, add = TRUE, col = "thistle")
points(akima, pch = 3, cex = 2, col = "blue")
# use triangulation package to show underlying triangulation:
\dontrun{
if(library(tripack, logical.return=TRUE))
plot(tri.mesh(akima), add=TRUE, lty="dashed")
}
# use only 15 points (interpolation only within convex hull!)
akima.part <- with(akima, interp(x[1:15], y[1:15], z[1:15]))
image(akima.part)
......
......@@ -144,16 +144,16 @@
IER = 0
RETURN
* Error exit
20 WRITE (*,FMT=9000) MD,NDP
20 CONTINUE WRITE (*,FMT=9000) MD,NDP
IER = 1
RETURN
30 WRITE (*,FMT=9010) MD,NDP,NDPPV
30 CONTINUE WRITE (*,FMT=9010) MD,NDP,NDPPV
IER = 2
RETURN
40 WRITE (*,FMT=9020) MD,NDP,NIP
40 CONTINUE WRITE (*,FMT=9020) MD,NDP,NIP
IER = 3
RETURN
50 WRITE (*,FMT=9030)
50 CONTINUE WRITE (*,FMT=9030)
IER = 9
RETURN
* Format statement for error message
......@@ -324,19 +324,19 @@
IER = 0
RETURN
* Error exit
40 WRITE (*,FMT=9000) MD,NDP
40 CONTINUE WRITE (*,FMT=9000) MD,NDP
IER = 1
RETURN
50 WRITE (*,FMT=9010) MD,NDP,NDPPV
50 CONTINUE WRITE (*,FMT=9010) MD,NDP,NDPPV
IER = 2
RETURN
60 WRITE (*,FMT=9020) MD,NDP,NXI,NYI
60 CONTINUE WRITE (*,FMT=9020) MD,NDP,NXI,NYI
IER = 3
RETURN
70 WRITE (*,FMT=9030) MD,NDP,NXI,NYI
70 CONTINUE WRITE (*,FMT=9030) MD,NDP,NXI,NYI
IER = 4
RETURN
80 WRITE (*,FMT=9040)
80 CONTINUE WRITE (*,FMT=9040)
IER = 9
RETURN
* Format statement for error message
......@@ -445,22 +445,22 @@
* Error exit
10 IF (IERTM.EQ.-1) THEN
IERT = 1
WRITE (*,FMT=9000) NDP
CONTINUE WRITE (*,FMT=9000) NDP
ELSE IF (IERTM.EQ.-2) THEN
IERT = 2
WRITE (*,FMT=9010)
CONTINUE WRITE (*,FMT=9010)
ELSE
IERT = 3
IP1 = IERTM
WRITE (*,FMT=9020) NDP,IP1,XD(IP1),YD(IP1)
CONTINUE WRITE (*,FMT=9020) NDP,IP1,XD(IP1),YD(IP1)
END IF
RETURN
20 IF (IERTL.EQ.1) THEN
IERT = 4
WRITE (*,FMT=9030) NDP
CONTINUE WRITE (*,FMT=9030) NDP
ELSE IF (IERTL.EQ.2) THEN
IERT = 5
WRITE (*,FMT=9040)
CONTINUE WRITE (*,FMT=9040)
END IF
RETURN
* Format statements
......
......@@ -137,7 +137,7 @@ C INTERPOLATES THE ZI VALUES. (FOR MD=1,2,3)
81 CONTINUE
RETURN
C ERROR EXIT
90 WRITE (LUN,2090) MD0,NCP0,NDP0,NIP0
90 CONTINUE WRITE (LUN,2090) MD0,NCP0,NDP0,NIP0
RETURN
C FORMAT STATEMENT FOR ERROR MESSAGE
2090 FORMAT(1X/41H *** IMPROPER INPUT PARAMETER VALUE(S)./
......
......@@ -101,10 +101,10 @@ C - REPLACES THE LOCAL ARRAY FOR THE OUTPUT ARRAY.
59 CONTINUE
RETURN
C ERROR EXIT
90 WRITE (LUN,2090)
90 CONTINUE WRITE (LUN,2090)
GO TO 92
91 WRITE (LUN,2091)
92 WRITE (LUN,2092) NDP0,NCP0
91 CONTINUE WRITE (LUN,2091)
92 CONTINUE WRITE (LUN,2092) NDP0,NCP0
IPC(1)=0
RETURN
C FORMAT STATEMENTS FOR ERROR MESSAGES
......
......@@ -176,7 +176,7 @@ C INTERPOLATES THE ZI VALUES. (FOR MD=1,2,3)
89 CONTINUE
RETURN
C ERROR EXIT
90 WRITE (LUN,2090) MD0,NCP0,NDP0,NXI0,NYI0
90 CONTINUE WRITE (LUN,2090) MD0,NCP0,NDP0,NXI0,NYI0
RETURN
C FORMAT STATEMENT FOR ERROR MESSAGE
2090 FORMAT(1X/41H *** IMPROPER INPUT PARAMETER VALUE(S)./
......
......@@ -355,12 +355,12 @@ C ARE LISTED COUNTER-CLOCKWISE.
NL=NL0
RETURN
C ERROR EXIT
90 WRITE (LUN,2090) NDP0
90 CONTINUE WRITE (LUN,2090) NDP0
GO TO 93
91 WRITE (LUN,2091) NDP0,IP1,IP2,X1,Y1
91 CONTINUE WRITE (LUN,2091) NDP0,IP1,IP2,X1,Y1
GO TO 93
92 WRITE (LUN,2092) NDP0
93 WRITE (LUN,2093)
92 CONTINUE WRITE (LUN,2092) NDP0
93 CONTINUE WRITE (LUN,2093)
NT=0
RETURN
C FORMAT STATEMENTS
......
......@@ -1609,7 +1609,7 @@ C
C Error flag returned by OPTIM.
C
25 IER = 5
WRITE (*,100) NIT, IERR
CONTINUE WRITE (*,100) NIT, IERR
RETURN
100 FORMAT (//5X,'*** Error in OPTIM: NIT = ',I4,
. ', IER = ',I1,' ***'/)
......@@ -2105,7 +2105,7 @@ C Invalid triangulation data structure or collinear nodes
C on convex hull boundary.
C
33 IER = 3
WRITE (*,130) IN1, IN2
CONTINUE WRITE (*,130) IN1, IN2
130 FORMAT (//5X,'*** Error in EDGE: Invalid triangula',
. 'tion or null triangles on boundary'/
. 9X,'IN1 =',I4,', IN2=',I4/)
......@@ -2114,7 +2114,7 @@ C
C Error flag returned by OPTIM.
C
34 IER = 4
WRITE (*,140) NIT, IERR
CONTINUE WRITE (*,140) NIT, IERR
140 FORMAT (//5X,'*** Error in OPTIM: NIT = ',I4,
. ', IER = ',I1,' ***'/)
RETURN
......@@ -4423,7 +4423,7 @@ C
C
C Print a heading and test for invalid input.
C
WRITE (LUN,100)
CONTINUE WRITE (LUN,100)
NL = 1
IF (N .LT. 3 .OR. N .GT. NMAX .OR.
. (NROW .NE. 6 .AND. NROW .NE. 9) .OR.
......@@ -4431,21 +4431,21 @@ C
C
C Print an error message and bypass the loops.
C
WRITE (LUN,110) N, NROW, NT
CONTINUE WRITE (LUN,110) N, NROW, NT
GO TO 3
ENDIF
IF (PRNTX) THEN
C
C Print X and Y.
C
WRITE (LUN,101)
CONTINUE WRITE (LUN,101)
NL = 6
DO 1 I = 1,N
IF (NL .GE. NLMAX) THEN
WRITE (LUN,106)
CONTINUE WRITE (LUN,106)
NL = 0
ENDIF
WRITE (LUN,102) I, X(I), Y(I)
CONTINUE WRITE (LUN,102) I, X(I), Y(I)
NL = NL + 1
1 CONTINUE
ENDIF
......@@ -4453,21 +4453,21 @@ C
C Print the triangulation LTRI.
C
IF (NL .GT. NLMAX/2) THEN
WRITE (LUN,106)
CONTINUE WRITE (LUN,106)
NL = 0
ENDIF
IF (NROW .EQ. 6) THEN
WRITE (LUN,103)
CONTINUE WRITE (LUN,103)
ELSE
WRITE (LUN,104)
CONTINUE WRITE (LUN,104)
ENDIF
NL = NL + 5
DO 2 K = 1,NT
IF (NL .GE. NLMAX) THEN
WRITE (LUN,106)
CONTINUE WRITE (LUN,106)
NL = 0
ENDIF
WRITE (LUN,105) K, (LTRI(I,K), I = 1,NROW)
CONTINUE WRITE (LUN,105) K, (LTRI(I,K), I = 1,NROW)
NL = NL + 1
2 CONTINUE
C
......@@ -4476,13 +4476,13 @@ C triangles).
C
NB = 2*N - NT - 2
NA = NT + N - 1
IF (NL .GT. NLMAX-6) WRITE (LUN,106)
WRITE (LUN,107) NB, NA, NT
IF (NL .GT. NLMAX-6) CONTINUE WRITE (LUN,106)
CONTINUE WRITE (LUN,107) NB, NA, NT
C
C Print NCC and LCT.
C
3 WRITE (LUN,108) NCC
IF (NCC .GT. 0) WRITE (LUN,109) (LCT(I), I = 1,NCC)
3 CONTINUE WRITE (LUN,108) NCC
IF (NCC .GT. 0) CONTINUE WRITE (LUN,109) (LCT(I), I = 1,NCC)
RETURN
C
C Print formats:
......@@ -5561,7 +5561,7 @@ C
C
C Output header comments.
C
WRITE (LUN,100,ERR=13) IPX1, IPY1, IPX2, IPY2
CONTINUE WRITE (LUN,100,ERR=13) IPX1, IPY1, IPX2, IPY2
100 FORMAT ('%!PS-Adobe-3.0 EPSF-3.0'/
. '%%BoundingBox:',4I4/
. '%%Title: Triangulation'/
......@@ -5583,13 +5583,13 @@ C Set the line thickness to 2 points, and draw the
C viewport boundary.
C
T = 2.0
WRITE (LUN,110,ERR=13) T
WRITE (LUN,120,ERR=13) IPX1, IPY1
WRITE (LUN,130,ERR=13) IPX1, IPY2
WRITE (LUN,130,ERR=13) IPX2, IPY2
WRITE (LUN,130,ERR=13) IPX2, IPY1
WRITE (LUN,140,ERR=13)
WRITE (LUN,150,ERR=13)
CONTINUE WRITE (LUN,110,ERR=13) T
CONTINUE WRITE (LUN,120,ERR=13) IPX1, IPY1
CONTINUE WRITE (LUN,130,ERR=13) IPX1, IPY2
CONTINUE WRITE (LUN,130,ERR=13) IPX2, IPY2
CONTINUE WRITE (LUN,130,ERR=13) IPX2, IPY1
CONTINUE WRITE (LUN,140,ERR=13)
CONTINUE WRITE (LUN,150,ERR=13)
110 FORMAT (F12.6,' setlinewidth')
120 FORMAT (2I4,' moveto')
130 FORMAT (2I4,' lineto')
......@@ -5602,7 +5602,7 @@ C
SFY = DBLE(IH)/DY
TX = IPX1 - SFX*WX1
TY = IPY1 - SFY*WY1
WRITE (LUN,160,ERR=13) TX, TY, SFX, SFY
CONTINUE WRITE (LUN,160,ERR=13) TX, TY, SFX, SFY
160 FORMAT (2F12.6,' translate'/
. 2F12.6,' scale')
C
......@@ -5611,17 +5611,17 @@ C changed to reflect the new scaling which is applied to
C all subsequent output. Set it to 1.0 point.
C
T = 2.0/(SFX+SFY)
WRITE (LUN,110,ERR=13) T
CONTINUE WRITE (LUN,110,ERR=13) T
C
C Save the current graphics state, and set the clip path to
C the boundary of the window.
C
WRITE (LUN,170,ERR=13)
WRITE (LUN,180,ERR=13) WX1, WY1
WRITE (LUN,190,ERR=13) WX2, WY1
WRITE (LUN,190,ERR=13) WX2, WY2
WRITE (LUN,190,ERR=13) WX1, WY2
WRITE (LUN,200,ERR=13)
CONTINUE WRITE (LUN,170,ERR=13)
CONTINUE WRITE (LUN,180,ERR=13) WX1, WY1
CONTINUE WRITE (LUN,190,ERR=13) WX2, WY1
CONTINUE WRITE (LUN,190,ERR=13) WX2, WY2
CONTINUE WRITE (LUN,190,ERR=13) WX1, WY2
CONTINUE WRITE (LUN,200,ERR=13)
170 FORMAT ('gsave')
180 FORMAT (2F12.6,' moveto')
190 FORMAT (2F12.6,' lineto')
......@@ -5645,7 +5645,7 @@ C
C
C Add the edge to the path.
C
WRITE (LUN,210,ERR=13) X0, Y0, X(N1), Y(N1)
CONTINUE WRITE (LUN,210,ERR=13) X0, Y0, X(N1), Y(N1)
210 FORMAT (2F12.6,' moveto',2F12.6,' lineto')
ENDIF
IF (LP .NE. LPL) GO TO 2
......@@ -5702,7 +5702,7 @@ C (PASS1=FALSE and CNSTR=TRUE); i.e., CNSTR and PASS1
C have opposite values.
C
IF (CNSTR .NEQV. PASS1)
. WRITE (LUN,210,ERR=13) X0, Y0, X(N1), Y(N1)
. CONTINUE WRITE (LUN,210,ERR=13) X0, Y0, X(N1), Y(N1)
ENDIF
IF (N1 .EQ. N0BAK) CNSTR = .FALSE.
C
......@@ -5720,9 +5720,9 @@ C are applied to everything, the dash length must be
C specified in world coordinates.
C
PASS1 = .FALSE.
WRITE (LUN,150,ERR=13)
CONTINUE WRITE (LUN,150,ERR=13)
T = DASHL*2.0/(SFX+SFY)
WRITE (LUN,220,ERR=13) T
CONTINUE WRITE (LUN,220,ERR=13) T
220 FORMAT ('[',F12.6,'] 0 setdash')
GO TO 4
ENDIF
......@@ -5730,8 +5730,8 @@ C
C Paint the path and restore the saved graphics state (with
C no clip path).
C
WRITE (LUN,150,ERR=13)
WRITE (LUN,230,ERR=13)
CONTINUE WRITE (LUN,150,ERR=13)
CONTINUE WRITE (LUN,230,ERR=13)
230 FORMAT ('grestore')
IF (NUMBR) THEN
C
......@@ -5740,7 +5740,7 @@ C Convert FSIZN from points to world coordinates, and
C output the commands to select a font and scale it.
C
T = FSIZN*2.0/(SFX+SFY)
WRITE (LUN,240,ERR=13) T
CONTINUE WRITE (LUN,240,ERR=13) T
240 FORMAT ('/Helvetica findfont'/
. F12.6,' scalefont setfont')
C
......@@ -5756,8 +5756,8 @@ C Move to (X0,Y0), and draw the label N0. The first char-
C acter will have its lower left corner about one
C character width to the right of the nodal position.
C
WRITE (LUN,180,ERR=13) X0, Y0
WRITE (LUN,250,ERR=13) N0
CONTINUE WRITE (LUN,180,ERR=13) X0, Y0
CONTINUE WRITE (LUN,250,ERR=13) N0
250 FORMAT ('(',I3,') show')
9 CONTINUE
ENDIF
......@@ -5766,15 +5766,15 @@ C Convert FSIZT from points to world coordinates, and output
C the commands to select a font and scale it.
C
T = FSIZT*2.0/(SFX+SFY)
WRITE (LUN,240,ERR=13) T
CONTINUE WRITE (LUN,240,ERR=13) T
C
C Display TITLE centered above the plot:
C
Y0 = WY2 + 3.0*T
WRITE (LUN,260,ERR=13) TITLE, (WX1+WX2)/2.0, Y0
CONTINUE WRITE (LUN,260,ERR=13) TITLE, (WX1+WX2)/2.0, Y0
260 FORMAT (A80/' stringwidth pop 2 div neg ',F12.6,
. ' add ',F12.6,' moveto')
WRITE (LUN,270,ERR=13) TITLE
CONTINUE WRITE (LUN,270,ERR=13) TITLE
270 FORMAT (A80/' show')
IF (ANNOT) THEN
C
......@@ -5782,10 +5782,10 @@ C Display the window extrema below the plot.
C
X0 = WX1
Y0 = WY1 - 100.0/(SFX+SFY)
WRITE (LUN,180,ERR=13) X0, Y0
WRITE (LUN,280,ERR=13) WX1, WX2
CONTINUE WRITE (LUN,180,ERR=13) X0, Y0
CONTINUE WRITE (LUN,280,ERR=13) WX1, WX2
Y0 = Y0 - 2.0*T
WRITE (LUN,290,ERR=13) X0, Y0, WY1, WY2
CONTINUE WRITE (LUN,290,ERR=13) X0, Y0, WY1, WY2
280 FORMAT ('(Window: WX1 = ',E9.3,', WX2 = ',E9.3,
. ') show')
290 FORMAT ('(Window: ) stringwidth pop ',F12.6,' add',
......@@ -5796,7 +5796,7 @@ C
C Paint the path and output the showpage command and
C end-of-file indicator.
C
WRITE (LUN,300,ERR=13)
CONTINUE WRITE (LUN,300,ERR=13)
300 FORMAT ('stroke'/
. 'showpage'/
. '%%EOF')
......@@ -5805,7 +5805,7 @@ C HP's interpreters require a one-byte End-of-PostScript-Job
C indicator (to eliminate a timeout error message):
C ASCII 4.
C
WRITE (LUN,310,ERR=13) CHAR(4)
CONTINUE WRITE (LUN,310,ERR=13) CHAR(4)
310 FORMAT (A1)
C
C No error encountered.
......@@ -5895,12 +5895,12 @@ C
C
C Print a heading and test the range of N.
C
WRITE (LUN,100) NN
CONTINUE WRITE (LUN,100) NN
IF (NN .LT. 3 .OR. NN .GT. NMAX) THEN
C
C N is outside its valid range.
C
WRITE (LUN,110)
CONTINUE WRITE (LUN,110)
GO TO 5
ENDIF
C
......@@ -5914,7 +5914,7 @@ C
C Print LIST only. K is the number of neighbors of NODE
C which are stored in NABOR.
C
WRITE (LUN,101)
CONTINUE WRITE (LUN,101)
DO 2 NODE = 1,NN
LPL = LEND(NODE)
LP = LPL
......@@ -5942,17 +5942,17 @@ C
INC = (K-1)/14 + 2
NL = NL + INC
IF (NL .GT. NLMAX) THEN
WRITE (LUN,106)
CONTINUE WRITE (LUN,106)
NL = INC
ENDIF
WRITE (LUN,103) NODE, (NABOR(I), I = 1,K)
IF (K .NE. 14) WRITE (LUN,105)
CONTINUE WRITE (LUN,103) NODE, (NABOR(I), I = 1,K)
IF (K .NE. 14) CONTINUE WRITE (LUN,105)
2 CONTINUE
ELSE
C
C Print X, Y, and LIST.
C
WRITE (LUN,102)
CONTINUE WRITE (LUN,102)
DO 4 NODE = 1,NN
LPL = LEND(NODE)
LP = LPL
......@@ -5977,12 +5977,12 @@ C
INC = (K-1)/8 + 2
NL = NL + INC
IF (NL .GT. NLMAX) THEN
WRITE (LUN,106)
CONTINUE WRITE (LUN,106)
NL = INC
ENDIF
WRITE (LUN,104) NODE, X(NODE), Y(NODE),
. (NABOR(I), I = 1,K)
IF (K .NE. 8) WRITE (LUN,105)
CONTINUE WRITE (LUN,104) NODE, X(NODE), Y(NODE),
. (NABOR(I), I = 1,K)
IF (K .NE. 8) CONTINUE WRITE (LUN,105)
4 CONTINUE
ENDIF
C
......@@ -5991,13 +5991,13 @@ C triangles).
C
NT = 2*NN - NB - 2
NA = NT + NN - 1
IF (NL .GT. NLMAX-6) WRITE (LUN,106)
WRITE (LUN,107) NB, NA, NT
IF (NL .GT. NLMAX-6) CONTINUE WRITE (LUN,106)
CONTINUE WRITE (LUN,107) NB, NA, NT
C
C Print NCC and LCC.
C
5 WRITE (LUN,108) NCC
IF (NCC .GT. 0) WRITE (LUN,109) (LCC(I), I = 1,NCC)
5 CONTINUE WRITE (LUN,108) NCC
IF (NCC .GT. 0) CONTINUE WRITE (LUN,109) (LCC(I), I = 1,NCC)
RETURN
C
C Print formats:
......
......@@ -79,31 +79,31 @@ C CALCULATION ID000670
16 CONTINUE ID000810
17 CONTINUE ID000820
C PRINTING OF INPUT DATA ID000830
20 WRITE (LUN,2020) NDP ID000840
20 CONTINUE WRITE (LUN,2020) NDP ID000840
DO 23 IDP=1,NDP ID000850
IF(MOD(IDP,5).EQ.1) WRITE (LUN,2021) ID000860
WRITE (LUN,2022) IDP,XD(IDP),YD(IDP),ZD(IDP) ID000870
IF(MOD(IDP,5).EQ.1) CONTINUE WRITE (LUN,2021) ID000860
CONTINUE WRITE (LUN,2022) IDP,XD(IDP),YD(IDP),ZD(IDP) ID000870
23 CONTINUE ID000880
C PRINTING OF OUTPUT RESULTS ID000890
30 WRITE (LUN,2030) ID000900
WRITE (LUN,2031) YI ID000910
30 CONTINUE WRITE (LUN,2030) ID000900
CONTINUE WRITE (LUN,2031) YI ID000910
DO 33 IXI=1,NXI ID000920
WRITE (LUN,2032) XI(IXI),(ZI1(IXI,IYI),IYI=1,NYI) ID000930
CONTINUE WRITE (LUN,2032) XI(IXI),(ZI1(IXI,IYI),IYI=1,NYI) ID000930
33 CONTINUE ID000940
40 WRITE (LUN,2040) ID000950
WRITE (LUN,2031) YI ID000960
40 CONTINUE WRITE (LUN,2040) ID000950
CONTINUE WRITE (LUN,2031) YI ID000960
DO 43 IXI=1,NXI ID000970
WRITE (LUN,2032) XI(IXI),(DZI1(IXI,IYI),IYI=1,NYI) ID000980
CONTINUE WRITE (LUN,2032) XI(IXI),(DZI1(IXI,IYI),IYI=1,NYI) ID000980
43 CONTINUE ID000990
50 WRITE (LUN,2050) ID001000
WRITE (LUN,2031) YI ID001010
50 CONTINUE WRITE (LUN,2050) ID001000
CONTINUE WRITE (LUN,2031) YI ID001010
DO 53 IXI=1,NXI ID001020
WRITE (LUN,2032) XI(IXI),(ZI2(IXI,IYI),IYI=1,NYI) ID001030
CONTINUE WRITE (LUN,2032) XI(IXI),(ZI2(IXI,IYI),IYI=1,NYI) ID001030
53 CONTINUE ID001040
60 WRITE (LUN,2060) ID001050
WRITE (LUN,2031) YI ID001060
60 CONTINUE WRITE (LUN,2060) ID001050
CONTINUE WRITE (LUN,2031) YI ID001060
DO 63 IXI=1,NXI ID001070
WRITE (LUN,2032) XI(IXI),(DZI2(IXI,IYI),IYI=1,NYI) ID001080
CONTINUE WRITE (LUN,2032) XI(IXI),(DZI2(IXI,IYI),IYI=1,NYI) ID001080
63 CONTINUE ID001090
STOP ID001100
C FORMAT STATEMENTS ID001110
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment