+ mkdir temp8358 + cd temp8358 + [ -r fort.1 ] + [ -r fort.2 ] + ln -s /home/npavlis/AC0405/REF02/prdct/und_min1x1_egm2008_Nmax2190_MeanTide_TP_global fort.1 + ln -s /home/npavlis/AC0405/REF02/prdct/dot_min1x1_dot08a_Nmax180_TP_global fort.2 + cat + 1> scr8358.f 0<< program readmin1 implicit real*8(a-h,o-z) parameter(exclud=9999.d0, $ dlatg=01.d0/60.d0,dlong=01.d0/60.d0,nrowsg=10800,ncolsg=21600) real*4 data(ncolsg) dimension sdat(ncolsg,2),stati(22) c----------------------------------------------------------------------- c c Read-in Geoid Undulation 1x1 min. file. c write(6,6001) 6001 format(///15x,'Statistics of Geoid Undulation Values (m)',/) do i = 1, nrowsg read(1) ( data(j),j=1,ncolsg) rlat = 90.d0 - (i-1.d0)*dlatg - dlatg*0.5d0 do j = 1, ncolsg clon = (j-1.d0)*dlong + dlong*0.5d0 c c Print a few values. c if(i.eq. 1.and.j.lt.20) then ! top row write(6,6101) rlat,clon,data(j) 6101 format(5x,2f15.10,f15.5) endif ! top row if(i.eq.5400.and.j.lt.20) then ! mid row write(6,6101) rlat,clon,data(j) endif ! mid row sdat(j,1) = data(j) sdat(j,2) = 1.d0 enddo ! j call stats(90.d0,0.d0,i,dlatg,dlong,nrowsg,ncolsg,sdat,exclud,0, $ stati) enddo ! i c----------------------------------------------------------------------- c c Read-in Dynamic Ocean Topography 1x1 min. file. c write(6,6002) 6002 format(///15x,'Statistics of Dynam. Ocean Topo. Values (m)',/) do i = 1, nrowsg read(2) ( data(j),j=1,ncolsg) rlat = 90.d0 - (i-1.d0)*dlatg - dlatg*0.5d0 do j = 1, ncolsg clon = (j-1.d0)*dlong + dlong*0.5d0 c c Print a few values. c if(i.eq. 1.and.j.lt.20) then ! top row write(6,6101) rlat,clon,data(j) endif ! top row if(i.eq.5400.and.j.lt.20) then ! mid row write(6,6101) rlat,clon,data(j) endif ! mid row sdat(j,1) = data(j) sdat(j,2) = 1.d0 enddo ! j call stats(90.d0,0.d0,i,dlatg,dlong,nrowsg,ncolsg,sdat,exclud,0, $ stati) enddo ! i c----------------------------------------------------------------------- stop end C C C SUBROUTINE STATS(TOPLAT,WSTLON,I,GRDN,GRDE,NROWS,NCOLS,DATA, $ EXCLUD,ISIG,STAT) C----------------------------------------------------------------------- IMPLICIT REAL*8(A-H,O-Z) CHARACTER*20 DLABEL(14) CHARACTER*20 SLABEL( 8) DIMENSION DATA(NCOLS,2) DOUBLE PRECISION DEXCLUD,DG,SD,STAT(22) SAVE DATA PI/3.14159265358979323846D+00/ DATA DLABEL/' Number of Values',' Percentage of Area', $ ' Minimum Value',' Latitude of Minimum', $ 'Longitude of Minimum',' Maximum Value', $ ' Latitude of Maximum','Longitude of Maximum', $ ' Arithmetic Mean',' Area-Weighted Mean', $ ' Arithmetic RMS',' Area-Weighted RMS', $ ' Arithmetic S.Dev.','Area-Weighted S.Dev.'/ DATA SLABEL/' Minimum Sigma',' Latitude of Minimum', $ 'Longitude of Minimum',' Maximum Sigma', $ ' Latitude of Maximum','Longitude of Maximum', $ 'Arithmetic RMS Sigma','Area-wghtd RMS Sigma'/ C----------------------------------------------------------------------- IF(I.EQ.1) THEN DEXCLUD=EXCLUD DTR=PI/180.D0 FOURPI=4.D0*PI DPR=GRDN*DTR DLR=GRDE*DTR CAREA=2.D0*DLR*SIN(DPR/2.D0) DO 10 K=1,22 STAT(K)=0.D0 10 CONTINUE STAT( 3)= DEXCLUD STAT( 6)=-DEXCLUD STAT(15)= DEXCLUD STAT(18)= 0.D0 ENDIF C----------------------------------------------------------------------- DLAT=TOPLAT-(I-1.D0)*GRDN-GRDN/2.D0 COLATC=(90.D0-DLAT)*DTR AREA=CAREA*SIN(COLATC) C----------------------------------------------------------------------- DO 110 J=1,NCOLS DLON=WSTLON+(J-1.D0)*GRDE+GRDE/2.D0 DG=DATA(J,1) SD=DATA(J,2) IF(DG.LT.DEXCLUD) THEN C----------------------------------------------------------------------- STAT( 1)=STAT( 1)+1.D0 STAT( 2)=STAT( 2)+AREA IF(DG.LE.STAT( 3)) THEN STAT( 3)=DG STAT( 4)=DLAT STAT( 5)=DLON ENDIF IF(DG.GE.STAT( 6)) THEN STAT( 6)=DG STAT( 7)=DLAT STAT( 8)=DLON ENDIF STAT( 9)=STAT( 9)+DG STAT(10)=STAT(10)+DG*AREA STAT(11)=STAT(11)+DG**2 STAT(12)=STAT(12)+DG**2*AREA IF(SD.LE.STAT(15)) THEN STAT(15)=SD STAT(16)=DLAT STAT(17)=DLON ENDIF IF(SD.GE.STAT(18)) THEN STAT(18)=SD STAT(19)=DLAT STAT(20)=DLON ENDIF STAT(21)=STAT(21)+SD**2 STAT(22)=STAT(22)+SD**2*AREA C----------------------------------------------------------------------- ENDIF 110 CONTINUE C----------------------------------------------------------------------- IF(I.NE.NROWS) RETURN IF(STAT(1).GT.0.D0) THEN STAT( 9)=STAT( 9)/STAT( 1) STAT(10)=STAT(10)/STAT( 2) STAT(11)=SQRT(STAT(11)/STAT( 1)) STAT(12)=SQRT(STAT(12)/STAT( 2)) STAT(13)=SQRT(STAT(11)**2-STAT( 9)**2) STAT(14)=SQRT(STAT(12)**2-STAT(10)**2) STAT(21)=SQRT(STAT(21)/STAT( 1)) STAT(22)=SQRT(STAT(22)/STAT( 2)) STAT( 2)=STAT( 2)/FOURPI*100.D0 ELSE DO 120 J=3,22 STAT(J)=DEXCLUD 120 CONTINUE ENDIF C======================================================================= NUM=INT(STAT(1)) WRITE(6,6001) DLABEL(1),NUM 6001 FORMAT(/5X,A20,3X,I11) DO 210 K=2,14 WRITE(6,6002) DLABEL(K),STAT(K) 6002 FORMAT(5X,A20,3X,F15.3) 210 CONTINUE WRITE(6,6003) 6003 FORMAT(' ') IF(ISIG.EQ.1) THEN DO 220 K=1,8 WRITE(6,6002) SLABEL(K),STAT(K+14) 220 CONTINUE ENDIF C======================================================================= RETURN END + f90 scr8358.f + a.out Statistics of Geoid Undulation Values (m) 89.9916666667 0.0083333333 15.07745 89.9916666667 0.0250000000 15.07745 89.9916666667 0.0416666667 15.07745 89.9916666667 0.0583333333 15.07745 89.9916666667 0.0750000000 15.07745 89.9916666667 0.0916666667 15.07745 89.9916666667 0.1083333333 15.07745 89.9916666667 0.1250000000 15.07745 89.9916666667 0.1416666667 15.07745 89.9916666667 0.1583333333 15.07745 89.9916666667 0.1750000000 15.07745 89.9916666667 0.1916666667 15.07745 89.9916666667 0.2083333333 15.07745 89.9916666667 0.2250000000 15.07744 89.9916666667 0.2416666667 15.07744 89.9916666667 0.2583333333 15.07744 89.9916666667 0.2750000000 15.07744 89.9916666667 0.2916666667 15.07744 89.9916666667 0.3083333333 15.07744 0.0083333333 0.0083333333 17.75720 0.0083333333 0.0250000000 17.75329 0.0083333333 0.0416666667 17.74838 0.0083333333 0.0583333333 17.74277 0.0083333333 0.0750000000 17.73687 0.0083333333 0.0916666667 17.73118 0.0083333333 0.1083333333 17.72608 0.0083333333 0.1250000000 17.72182 0.0083333333 0.1416666667 17.71840 0.0083333333 0.1583333333 17.71560 0.0083333333 0.1750000000 17.71299 0.0083333333 0.1916666667 17.71008 0.0083333333 0.2083333333 17.70642 0.0083333333 0.2250000000 17.70165 0.0083333333 0.2416666667 17.69561 0.0083333333 0.2583333333 17.68836 0.0083333333 0.2750000000 17.68008 0.0083333333 0.2916666667 17.67109 0.0083333333 0.3083333333 17.66173 Number of Values 233280000 Percentage of Area 100.000 Minimum Value -106.369 Latitude of Minimum 4.692 Longitude of Minimum 78.758 Maximum Value 87.938 Latitude of Maximum -8.392 Longitude of Maximum 147.375 Arithmetic Mean -0.867 Area-Weighted Mean 0.000 Arithmetic RMS 29.225 Area-Weighted RMS 30.576 Arithmetic S.Dev. 29.212 Area-Weighted S.Dev. 30.576 Statistics of Dynam. Ocean Topo. Values (m) 89.9916666667 0.0083333333 -0.39565 89.9916666667 0.0250000000 -0.39565 89.9916666667 0.0416666667 -0.39565 89.9916666667 0.0583333333 -0.39565 89.9916666667 0.0750000000 -0.39565 89.9916666667 0.0916666667 -0.39565 89.9916666667 0.1083333333 -0.39565 89.9916666667 0.1250000000 -0.39565 89.9916666667 0.1416666667 -0.39565 89.9916666667 0.1583333333 -0.39565 89.9916666667 0.1750000000 -0.39565 89.9916666667 0.1916666667 -0.39565 89.9916666667 0.2083333333 -0.39565 89.9916666667 0.2250000000 -0.39565 89.9916666667 0.2416666667 -0.39565 89.9916666667 0.2583333333 -0.39565 89.9916666667 0.2750000000 -0.39565 89.9916666667 0.2916666667 -0.39565 89.9916666667 0.3083333333 -0.39565 0.0083333333 0.0083333333 -0.07798 0.0083333333 0.0250000000 -0.07799 0.0083333333 0.0416666667 -0.07798 0.0083333333 0.0583333333 -0.07798 0.0083333333 0.0750000000 -0.07798 0.0083333333 0.0916666667 -0.07797 0.0083333333 0.1083333333 -0.07796 0.0083333333 0.1250000000 -0.07795 0.0083333333 0.1416666667 -0.07794 0.0083333333 0.1583333333 -0.07792 0.0083333333 0.1750000000 -0.07790 0.0083333333 0.1916666667 -0.07788 0.0083333333 0.2083333333 -0.07786 0.0083333333 0.2250000000 -0.07783 0.0083333333 0.2416666667 -0.07781 0.0083333333 0.2583333333 -0.07778 0.0083333333 0.2750000000 -0.07774 0.0083333333 0.2916666667 -0.07771 0.0083333333 0.3083333333 -0.07767 Number of Values 233280000 Percentage of Area 100.000 Minimum Value -1.955 Latitude of Minimum -60.708 Longitude of Minimum 358.458 Maximum Value 1.129 Latitude of Maximum 31.075 Longitude of Maximum 135.442 Arithmetic Mean -0.196 Area-Weighted Mean -0.026 Arithmetic RMS 0.636 Area-Weighted RMS 0.561 Arithmetic S.Dev. 0.605 Area-Weighted S.Dev. 0.561 + rm scr8358.f a.out fort.1 fort.2 + cd .. + rmdir temp8358