Context Navigation

grid_noro.F @ 319

Last change on this file since 319 was 228, checked in by lmdzadmin, 23 years ago
detail dans l'initialisation LF
Property svn:eol-style set to `native` Property svn:executable set to ``* Property svn:keywords set to `Author Date Id Revision`
File size: 14.8 KB

Rev	Line
[2]	1
	2	SUBROUTINE grid_noro(imdep, jmdep, xdata, ydata, zdata,
	3	. imar, jmar, x, y,
	4	. zphi,zmea,zstd,zsig,zgam,zthe,
	5	. zpic,zval,mask)
	6	c=======================================================================
	7	c (F. Lott) (voir aussi z.x. Li, A. Harzallah et L. Fairhead)
	8	c
	9	c Compute the Parameters of the SSO scheme as described in
	10	c LOTT & MILLER (1997) and LOTT(1999).
	11	c Target points are on a rectangular grid:
	12	c iim+1 latitudes including North and South Poles;
	13	c jjm+1 longitudes, with periodicity jjm+1=1.
	14	c aux poles. At the poles the fields value is repeated
	15	c jjm+1 time.
	16	c The parameters a,b,c,d represent the limite of the target
	17	c gridpoint region. The means over this region are calculated
	18	c from USN data, ponderated by a weight proportional to the
	19	c surface occupated by the data inside the model gridpoint area.
	20	c In most circumstances, this weight is the ratio between the
	21	c surface of the USN gridpoint area and the surface of the
	22	c model gridpoint area.
	23	c
	24	c (c)
	25	c ----d-----
	26	c \| . . . .\|
	27	c \| \|
	28	c (b)a . * . .b(a)
	29	c \| \|
	30	c \| . . . .\|
	31	c ----c-----
	32	c (d)
	33	C=======================================================================
	34	c INPUT:
	35	c imdep, jmdep: dimensions X and Y input field
	36	c xdata, ydata: coordinates X and Y input field
	37	c zdata: Input field
	38	c In this version it is assumed that the entry data come from
	39	c the USNavy dataset: imdep=iusn=2160, jmdep=jusn=1080.
	40	c OUTPUT:
	41	c imar, jmar: dimensions X and Y Output field
	42	c x, y: ccordinates X and Y Output field.
	43	c zmea: Mean orographie
	44	c zstd: Standard deviation
	45	c zsig: Slope
	46	c zgam: Anisotropy
	47	c zthe: Orientation of the small axis
	48	c zpic: Maximum altitude
	49	c zval: Minimum altitude
	50	C=======================================================================
	51	IMPLICIT INTEGER (I,J)
	52	IMPLICIT REAL(X,Z)
	53
[228]	54	parameter(iusn=2160,jusn=1080,iext=216, epsfra = 1.e-5)
[2]	55	#include "dimensions.h"
	56	REAL xusn(iusn+2*iext),yusn(jusn+2)
	57	REAL zusn(iusn+2*iext,jusn+2)
	58
	59	INTEGER imdep, jmdep
	60	REAL xdata(imdep),ydata(jmdep)
	61	REAL zdata(imdep,jmdep)
	62	c
	63	INTEGER imar, jmar
	64
	65	C INTERMEDIATE FIELDS (CORRELATIONS OF OROGRAPHY GRADIENT)
	66
	67	REAL ztz(iim+1,jjm+1),zxtzx(iim+1,jjm+1)
	68	REAL zytzy(iim+1,jjm+1),zxtzy(iim+1,jjm+1)
	69	REAL weight(iim+1,jjm+1)
	70
	71	C CORRELATIONS OF USN OROGRAPHY GRADIENTS
	72
	73	REAL zxtzxusn(iusn+2iext,jusn+2),zytzyusn(iusn+2iext,jusn+2)
	74	REAL zxtzyusn(iusn+2*iext,jusn+2)
	75	REAL x(imar+1),y(jmar),zphi(imar+1,jmar)
	76	REAL zmea(imar+1,jmar),zstd(imar+1,jmar)
	77	REAL zsig(imar+1,jmar),zgam(imar+1,jmar),zthe(imar+1,jmar)
	78	REAL zpic(imar+1,jmar),zval(imar+1,jmar)
[99]	79	c$$$ PB integer mask(imar+1,jmar)
[228]	80	real mask(imar+1,jmar), mask_tmp(imar+1,jmar)
[2]	81	real num_tot(2200,1100),num_lan(2200,1100)
	82	c
	83	REAL a(2200),b(2200),c(1100),d(1100)
	84	c
	85	print *,' parametres de l orographie a l echelle sous maille'
	86	xpi=acos(-1.)
	87	rad = 6 371 229.
	88	zdeltay=2.xpi/float(jusn)rad
	89	c
	90	c quelques tests de dimensions:
	91	c
	92	c
	93	if(iim.ne.imar) STOP 'Problem dim. x'
	94	if(jjm.ne.jmar-1) STOP 'Problem dim. y'
	95	IF (imar.GT.2200 .OR. jmar.GT.1100) THEN
	96	PRINT*, 'imar or jmar too big', imar, jmar
	97	CALL ABORT
	98	ENDIF
	99
	100	IF(imdep.ne.iusn.or.jmdep.ne.jusn)then
	101	print *,' imdep or jmdep bad dimensions:',imdep,jmdep
	102	call abort
	103	ENDIF
	104
	105	IF(imar+1.ne.iim+1.or.jmar.ne.jjm+1)THEN
	106	print *,' imar or jmar bad dimensions:',imar,jmar
	107	call abort
	108	ENDIF
	109
	110
	111	c print *,'xdata:',xdata
	112	c print *,'ydata:',ydata
	113	c print *,'x:',x
	114	c print *,'y:',y
	115	c
	116	C EXTENSION OF THE USN DATABASE TO POCEED COMPUTATIONS AT
	117	C BOUNDARIES:
	118	c
	119	DO j=1,jusn
	120	yusn(j+1)=ydata(j)
	121	DO i=1,iusn
	122	zusn(i+iext,j+1)=zdata(i,j)
	123	xusn(i+iext)=xdata(i)
	124	ENDDO
	125	DO i=1,iext
	126	zusn(i,j+1)=zdata(iusn-iext+i,j)
	127	xusn(i)=xdata(iusn-iext+i)-2.*xpi
	128	zusn(iusn+iext+i,j+1)=zdata(i,j)
	129	xusn(iusn+iext+i)=xdata(i)+2.*xpi
	130	ENDDO
	131	ENDDO
	132
	133	yusn(1)=ydata(1)+(ydata(1)-ydata(2))
	134	yusn(jusn+2)=ydata(jusn)+(ydata(jusn)-ydata(jusn-1))
	135	DO i=1,iusn/2+iext
	136	zusn(i,1)=zusn(i+iusn/2,2)
	137	zusn(i+iusn/2+iext,1)=zusn(i,2)
	138	zusn(i,jusn+2)=zusn(i+iusn/2,jusn+1)
	139	zusn(i+iusn/2+iext,jusn+2)=zusn(i,jusn+1)
	140	ENDDO
	141	c
	142	c COMPUTE LIMITS OF MODEL GRIDPOINT AREA
	143	C ( REGULAR GRID)
	144	c
	145	a(1) = x(1) - (x(2)-x(1))/2.0
	146	b(1) = (x(1)+x(2))/2.0
	147	DO i = 2, imar
	148	a(i) = b(i-1)
	149	b(i) = (x(i)+x(i+1))/2.0
	150	ENDDO
	151	a(imar+1) = b(imar)
	152	b(imar+1) = x(imar+1) + (x(imar+1)-x(imar))/2.0
	153
	154	c(1) = y(1) - (y(2)-y(1))/2.0
	155	d(1) = (y(1)+y(2))/2.0
	156	DO j = 2, jmar-1
	157	c(j) = d(j-1)
	158	d(j) = (y(j)+y(j+1))/2.0
	159	ENDDO
	160	c(jmar) = d(jmar-1)
	161	d(jmar) = y(jmar) + (y(jmar)-y(jmar-1))/2.0
	162	c
	163	c initialisations:
	164	c
	165	DO i = 1, imar+1
	166	DO j = 1, jmar
	167	weight(i,j) = 0.0
	168	zxtzx(i,j) = 0.0
	169	zytzy(i,j) = 0.0
	170	zxtzy(i,j) = 0.0
	171	ztz(i,j) = 0.0
	172	zmea(i,j) = 0.0
	173	zpic(i,j) =-1.E+10
	174	zval(i,j) = 1.E+10
	175	ENDDO
	176	ENDDO
	177	c
	178	c COMPUTE SLOPES CORRELATIONS ON USN GRID
	179	c
	180	DO j = 1,jusn+2
	181	DO i = 1, iusn+2*iext
	182	zytzyusn(i,j)=0.0
	183	zxtzxusn(i,j)=0.0
	184	zxtzyusn(i,j)=0.0
	185	ENDDO
	186	ENDDO
	187
	188
	189	DO j = 2,jusn+1
	190	zdeltax=zdeltay*cos(yusn(j))
	191	DO i = 2, iusn+2*iext-1
	192	zytzyusn(i,j)=(zusn(i,j+1)-zusn(i,j-1))2/zdeltay2
	193	zxtzxusn(i,j)=(zusn(i+1,j)-zusn(i-1,j))2/zdeltax2
	194	zxtzyusn(i,j)=(zusn(i,j+1)-zusn(i,j-1))/zdeltay
	195	* *(zusn(i+1,j)-zusn(i-1,j))/zdeltax
	196	ENDDO
	197	ENDDO
	198	c
	199	c SUMMATION OVER GRIDPOINT AREA
	200	c
	201	zleny=xpi/float(jusn)*rad
	202	xincr=xpi/2./float(jusn)
	203	DO ii = 1, imar+1
	204	DO jj = 1, jmar
	205	num_tot(ii,jj)=0.
	206	num_lan(ii,jj)=0.
	207	c PRINT *,' iteration ii jj:',ii,jj
	208	DO j = 2,jusn+1
	209	c DO j = 3,jusn
	210	zlenx=zleny*cos(yusn(j))
	211	zdeltax=zdeltay*cos(yusn(j))
	212	zbordnor=(c(jj)-yusn(j)+xincr)*rad
	213	zbordsud=(yusn(j)-d(jj)+xincr)*rad
	214	weighy=AMAX1(0.,
	215	* amin1(zbordnor,zbordsud,zleny))
	216	IF(weighy.ne.0)THEN
	217	DO i = 2, iusn+2*iext-1
	218	zbordest=(xusn(i)-a(ii)+xincr)radcos(yusn(j))
	219	zbordoue=(b(ii)+xincr-xusn(i))radcos(yusn(j))
	220	weighx=AMAX1(0.,
	221	* amin1(zbordest,zbordoue,zlenx))
	222	IF(weighx.ne.0)THEN
	223	num_tot(ii,jj)=num_tot(ii,jj)+1.0
	224	if(zusn(i,j).ge.1.)num_lan(ii,jj)=num_lan(ii,jj)+1.0
	225	weight(ii,jj)=weight(ii,jj)+weighx*weighy
	226	zxtzx(ii,jj)=zxtzx(ii,jj)+zxtzxusn(i,j)weighxweighy
	227	zytzy(ii,jj)=zytzy(ii,jj)+zytzyusn(i,j)weighxweighy
	228	zxtzy(ii,jj)=zxtzy(ii,jj)+zxtzyusn(i,j)weighxweighy
	229	ztz(ii,jj) =ztz(ii,jj) +zusn(i,j)zusn(i,j)weighx*weighy
	230	c mean
	231	zmea(ii,jj) =zmea(ii,jj)+zusn(i,j)weighxweighy
	232	c peacks
	233	zpic(ii,jj)=amax1(zpic(ii,jj),zusn(i,j))
	234	c valleys
	235	zval(ii,jj)=amin1(zval(ii,jj),zusn(i,j))
	236	ENDIF
	237	ENDDO
	238	ENDIF
	239	ENDDO
	240	ENDDO
	241	ENDDO
	242	c
	243	c COMPUTE PARAMETERS NEEDED BY THE LOTT & MILLER (1997) AND
	244	C LOTT (1999) SSO SCHEME.
	245	c
	246	zllmmea=0.
	247	zllmstd=0.
	248	zllmsig=0.
	249	zllmgam=0.
	250	zllmpic=0.
	251	zllmval=0.
	252	zllmthe=0.
	253	zminthe=0.
	254	c print 100,' '
	255	c100 format(1X,A1,'II JJ',4X,'H',8X,'SD',8X,'SI',3X,'GA',3X,'TH')
	256	DO ii = 1, imar+1
	257	DO jj = 1, jmar
	258	IF (weight(ii,jj) .NE. 0.0) THEN
	259	c Mask
[99]	260	c$$$ if(num_lan(ii,jj)/num_tot(ii,jj).ge.0.5)then
	261	c$$$ mask(ii,jj)=1
	262	c$$$ else
	263	c$$$ mask(ii,jj)=0
	264	c$$$ ENDIF
	265	mask(ii,jj) = num_lan(ii,jj)/num_tot(ii,jj)
[2]	266	c Mean Orography:
	267	zmea (ii,jj)=zmea (ii,jj)/weight(ii,jj)
	268	zxtzx(ii,jj)=zxtzx(ii,jj)/weight(ii,jj)
	269	zytzy(ii,jj)=zytzy(ii,jj)/weight(ii,jj)
	270	zxtzy(ii,jj)=zxtzy(ii,jj)/weight(ii,jj)
	271	ztz(ii,jj) =ztz(ii,jj)/weight(ii,jj)
	272	c Standard deviation:
	273	zstd(ii,jj)=sqrt(AMAX1(0.,ztz(ii,jj)-zmea(ii,jj)**2))
	274	ELSE
	275	PRINT*, 'probleme,ii,jj=', ii,jj
	276	ENDIF
	277	ENDDO
	278	ENDDO
	279
	280	C CORRECT VALUES OF HORIZONTAL SLOPE NEAR THE POLES:
	281
	282	DO ii = 1, imar+1
	283	zxtzx(ii,1)=zxtzx(ii,2)
	284	zxtzx(ii,jmar)=zxtzx(ii,jmar-1)
	285	zxtzy(ii,1)=zxtzy(ii,2)
	286	zxtzy(ii,jmar)=zxtzy(ii,jmar-1)
	287	zytzy(ii,1)=zytzy(ii,2)
	288	zytzy(ii,jmar)=zytzy(ii,jmar-1)
	289	ENDDO
	290
	291	C FILTERS TO SMOOTH OUT FIELDS FOR INPUT INTO SSO SCHEME.
	292
	293	C FIRST FILTER, MOVING AVERAGE OVER 9 POINTS.
	294
	295	CALL MVA9(zmea,iim+1,jjm+1)
	296	CALL MVA9(zstd,iim+1,jjm+1)
	297	CALL MVA9(zpic,iim+1,jjm+1)
	298	CALL MVA9(zval,iim+1,jjm+1)
	299	CALL MVA9(zxtzx,iim+1,jjm+1)
[99]	300	CALL MVA9(zxtzy,iim+1,jjm+1)
	301	c$$$ CALL MVA9(zytzy,iim+1,jjm+1)
[228]	302	C$$$ Masque prenant en compte maximum de terre
	303	mask_tmp= 0.0
	304	WHERE(mask .GE. epsfra) mask_tmp = 1.
[2]	305
	306	DO ii = 1, imar
	307	DO jj = 1, jmar
	308	IF (weight(ii,jj) .NE. 0.0) THEN
	309	c Coefficients K, L et M:
	310	xk=(zxtzx(ii,jj)+zytzy(ii,jj))/2.
	311	xl=(zxtzx(ii,jj)-zytzy(ii,jj))/2.
	312	xm=zxtzy(ii,jj)
	313	xp=xk-sqrt(xl2+xm2)
	314	xq=xk+sqrt(xl2+xm2)
	315	xw=1.e-8
	316	if(xp.le.xw) xp=0.
	317	if(xq.le.xw) xq=xw
	318	if(abs(xm).le.xw) xm=xw*sign(1.,xm)
[228]	319	c slope:
[99]	320	c$$$ zsig(ii,jj)=sqrt(xq)*mask(ii,jj)
	321	c$$$c isotropy:
	322	c$$$ zgam(ii,jj)=xp/xq*mask(ii,jj)
	323	c$$$c angle theta:
	324	c$$$ zthe(ii,jj)=57.29577951atan2(xm,xl)/2.mask(ii,jj)
	325	c$$$ zphi(ii,jj)=zmea(ii,jj)*mask(ii,jj)
	326	c$$$ zmea(ii,jj)=zmea(ii,jj)*mask(ii,jj)
	327	c$$$ zpic(ii,jj)=zpic(ii,jj)*mask(ii,jj)
	328	c$$$ zval(ii,jj)=zval(ii,jj)*mask(ii,jj)
	329	c$$$ zstd(ii,jj)=zstd(ii,jj)*mask(ii,jj)
	330	C$$* PB modif pour maque de terre fractionnaire
[2]	331	c slope:
[228]	332	zsig(ii,jj)=sqrt(xq)*mask_tmp(ii,jj)
[2]	333	c isotropy:
[228]	334	zgam(ii,jj)=xp/xq*mask_tmp(ii,jj)
[2]	335	c angle theta:
[228]	336	zthe(ii,jj)=57.29577951atan2(xm,xl)/2.mask_tmp(ii,jj)
	337	zphi(ii,jj)=zmea(ii,jj)*mask_tmp(ii,jj)
	338	zmea(ii,jj)=zmea(ii,jj)*mask_tmp(ii,jj)
	339	zpic(ii,jj)=zpic(ii,jj)*mask_tmp(ii,jj)
	340	zval(ii,jj)=zval(ii,jj)*mask_tmp(ii,jj)
	341	zstd(ii,jj)=zstd(ii,jj)*mask_tmp(ii,jj)
[2]	342	c print 101,ii,jj,
	343	c * zmea(ii,jj),zstd(ii,jj),zsig(ii,jj),zgam(ii,jj),
	344	c * zthe(ii,jj)
	345	c101 format(1x,2(1x,i2),2(1x,f7.1),1x,f7.4,2x,f4.2,1x,f5.1)
	346	ELSE
	347	c PRINT*, 'probleme,ii,jj=', ii,jj
	348	ENDIF
	349	zllmmea=AMAX1(zmea(ii,jj),zllmmea)
	350	zllmstd=AMAX1(zstd(ii,jj),zllmstd)
	351	zllmsig=AMAX1(zsig(ii,jj),zllmsig)
	352	zllmgam=AMAX1(zgam(ii,jj),zllmgam)
	353	zllmthe=AMAX1(zthe(ii,jj),zllmthe)
	354	zminthe=amin1(zthe(ii,jj),zminthe)
	355	zllmpic=AMAX1(zpic(ii,jj),zllmpic)
	356	zllmval=AMAX1(zval(ii,jj),zllmval)
	357	ENDDO
	358	ENDDO
	359
	360	print *,' MEAN ORO:',zllmmea
	361	print *,' ST. DEV.:',zllmstd
	362	print *,' PENTE:',zllmsig
	363	print *,' ANISOTROP:',zllmgam
	364	print *,' ANGLE:',zminthe,zllmthe
	365	print *,' pic:',zllmpic
	366	print *,' val:',zllmval
	367
	368	C
	369	c gamma and theta a 1. and 0. at poles
	370	c
	371	DO jj=1,jmar
	372	zmea(imar+1,jj)=zmea(1,jj)
[7]	373	zphi(imar+1,jj)=zphi(1,jj)
[2]	374	zpic(imar+1,jj)=zpic(1,jj)
	375	zval(imar+1,jj)=zval(1,jj)
	376	zstd(imar+1,jj)=zstd(1,jj)
	377	zsig(imar+1,jj)=zsig(1,jj)
	378	zgam(imar+1,jj)=zgam(1,jj)
	379	zthe(imar+1,jj)=zthe(1,jj)
	380	ENDDO
	381
	382
	383	zmeanor=0.0
	384	zmeasud=0.0
	385	zstdnor=0.0
	386	zstdsud=0.0
	387	zsignor=0.0
	388	zsigsud=0.0
	389	zweinor=0.0
	390	zweisud=0.0
	391	zpicnor=0.0
	392	zpicsud=0.0
	393	zvalnor=0.0
	394	zvalsud=0.0
	395
	396	DO ii=1,imar
	397	zweinor=zweinor+ weight(ii, 1)
	398	zweisud=zweisud+ weight(ii,jmar)
	399	zmeanor=zmeanor+zmea(ii, 1)*weight(ii, 1)
	400	zmeasud=zmeasud+zmea(ii,jmar)*weight(ii,jmar)
	401	zstdnor=zstdnor+zstd(ii, 1)*weight(ii, 1)
	402	zstdsud=zstdsud+zstd(ii,jmar)*weight(ii,jmar)
	403	zsignor=zsignor+zsig(ii, 1)*weight(ii, 1)
	404	zsigsud=zsigsud+zsig(ii,jmar)*weight(ii,jmar)
	405	zpicnor=zpicnor+zpic(ii, 1)*weight(ii, 1)
	406	zpicsud=zpicsud+zpic(ii,jmar)*weight(ii,jmar)
	407	zvalnor=zvalnor+zval(ii, 1)*weight(ii, 1)
	408	zvalsud=zvalsud+zval(ii,jmar)*weight(ii,jmar)
	409	ENDDO
	410
	411	DO ii=1,imar+1
	412	zmea(ii, 1)=zmeanor/zweinor
	413	zmea(ii,jmar)=zmeasud/zweisud
	414	zphi(ii, 1)=zmeanor/zweinor
	415	zphi(ii,jmar)=zmeasud/zweisud
	416	zpic(ii, 1)=zpicnor/zweinor
	417	zpic(ii,jmar)=zpicsud/zweisud
	418	zval(ii, 1)=zvalnor/zweinor
	419	zval(ii,jmar)=zvalsud/zweisud
	420	zstd(ii, 1)=zstdnor/zweinor
	421	zstd(ii,jmar)=zstdsud/zweisud
	422	zsig(ii, 1)=zsignor/zweinor
	423	zsig(ii,jmar)=zsigsud/zweisud
	424	zgam(ii, 1)=1.
	425	zgam(ii,jmar)=1.
	426	zthe(ii, 1)=0.
	427	zthe(ii,jmar)=0.
	428	ENDDO
	429
	430	RETURN
	431	END
	432
	433	SUBROUTINE MVA9(X,IMAR,JMAR)
	434
	435	C MAKE A MOVING AVERAGE OVER 9 GRIDPOINTS OF THE X FIELDS
	436
	437	PARAMETER (ISMo=300,JSMo=200)
	438	REAL X(IMAR,JMAR),XF(ISMo,JSMo)
[99]	439	real WEIGHTpb(-1:1,-1:1)
[2]	440
	441	if(imar.gt.ismo) stop'surdimensionner ismo dans mva9 (grid_noro)'
	442	if(jmar.gt.jsmo) stop'surdimensionner jsmo dans mva9 (grid_noro)'
	443
	444	SUM=0.
	445	DO IS=-1,1
[99]	446	DO JS=-1,1
	447	WEIGHTpb(IS,JS)=1./FLOAT((1+IS*2)(1+JS**2))
	448	SUM=SUM+WEIGHTpb(IS,JS)
	449	ENDDO
[2]	450	ENDDO
[99]	451
	452	WRITE(,) 'MVA9 ', IMAR, JMAR
	453	WRITE(,) 'MVA9 ', WEIGHTpb
	454	WRITE(,) 'MVA9 SUM ', SUM
[2]	455	DO IS=-1,1
[99]	456	DO JS=-1,1
	457	WEIGHTpb(IS,JS)=WEIGHTpb(IS,JS)/SUM
	458	ENDDO
[2]	459	ENDDO
	460
	461	DO J=2,JMAR-1
[99]	462	DO I=2,IMAR-1
	463	XF(I,J)=0.
	464	DO IS=-1,1
	465	DO JS=-1,1
	466	XF(I,J)=XF(I,J)+X(I+IS,J+JS)*WEIGHTpb(IS,JS)
	467	ENDDO
	468	ENDDO
	469	ENDDO
[2]	470	ENDDO
	471
	472	DO J=2,JMAR-1
[99]	473	XF(1,J)=0.
	474	IS=IMAR-1
	475	DO JS=-1,1
	476	XF(1,J)=XF(1,J)+X(IS,J+JS)*WEIGHTpb(-1,JS)
	477	ENDDO
	478	DO IS=0,1
	479	DO JS=-1,1
	480	XF(1,J)=XF(1,J)+X(1+IS,J+JS)*WEIGHTpb(IS,JS)
	481	ENDDO
	482	ENDDO
	483	XF(IMAR,J)=XF(1,J)
[2]	484	ENDDO
	485
	486	DO I=1,IMAR
[99]	487	XF(I,1)=XF(I,2)
	488	XF(I,JMAR)=XF(I,JMAR-1)
[2]	489	ENDDO
	490
	491	DO I=1,IMAR
[99]	492	DO J=1,JMAR
	493	X(I,J)=XF(I,J)
	494	ENDDO
[2]	495	ENDDO
	496
	497	RETURN
	498	END
	499
	500
	501

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format