sfctest.f
76 lines
| 2.6 KiB
| text/x-fortran
|
FortranFixedLexer
| r0 | PROGRAM SFCCONV | |||
| c | ||||
| c $Revision: 1.1 $ | ||||
| c | ||||
| csfcconv - program to test the Geographic to/from Geomagnetic Coordinate | ||||
| c subprogram to be run at AFSFC | ||||
| c | ||||
| c | ||||
| c Comments: This program and the SFC subroutine have been developed on a | ||||
| c DEC VAX 7000 using VMS FORTRAN and trying to maintain the | ||||
| c SFC style of coding. The calling sequence is identical to the | ||||
| c previous version of the subprogram. | ||||
| c | ||||
| c This version has been ported to the unix environment. The original | ||||
| c program used a loop and ran through a set of altitudes from | ||||
| c 0 km to 2000 km. It opened the file SFCOUT.DAT for each height | ||||
| c and made use of the fact the VMS keeps multiple versions of | ||||
| c files to allow you to have a separate file for each altitude. | ||||
| c | ||||
| c For the UNIX implementation the automatic loop through the | ||||
| c altitudes has been replaced by a request for the user to | ||||
| c input an arbitrary altitude to be used. | ||||
| c | ||||
| c $Log: not supported by cvs2svn $ | ||||
| c Revision 1.1 94/10/14 10:59:07 10:59:07 baker (Kile Baker S1G) | ||||
| c Initial revision | ||||
| c | ||||
| c | ||||
| C BEGIN PROGRAM | ||||
| C RUN FOR ALTITUDES FROM 0 KM TO 2000 KM | ||||
| C DO 200 JH = 0, 20 | ||||
| C .. Local Scalars .. | ||||
| DOUBLE PRECISION CGLAT,CGLON,GLAT,GLON,GMLAT_,GMLON_,R_HEIGHT_IN | ||||
| INTEGER I_ERROR,LAT,LON | ||||
| C .. | ||||
| C .. External Subroutines .. | ||||
| EXTERNAL CONVERT_GEO_COORD | ||||
| C .. | ||||
| C .. Intrinsic Functions .. | ||||
| INTRINSIC DBLE | ||||
| C .. | ||||
| R_HEIGHT_IN = 0.0D0 | ||||
| C type *,'enter the height to be used' | ||||
| C accept *,r_height_in | ||||
| C | ||||
| C | ||||
| OPEN (1,FILE='SFCOUT.DAT',FORM='FORMATTED',STATUS='UNKNOWN') | ||||
| CRUN FOR ALL GEOGRAPHICS LATITUDES (STEPPING BY 2 DEGREES) AND LONGITUDES | ||||
| C (STEPPING BY TEN DEGREES) | ||||
| DO 20 LAT = 67,69,2 | ||||
| GLAT = DBLE(LAT) | ||||
| DO 10 LON = 309,319,10 | ||||
| GLON = DBLE(LON) | ||||
| C Convert from Geograhic (Geocentric) to Geomagnetic | ||||
| CALL CONVERT_GEO_COORD(GLAT,GLON,R_HEIGHT_IN,CGLAT,CGLON,1, | ||||
| * I_ERROR) | ||||
| IF (I_ERROR.NE.0) THEN | ||||
| PRINT *,'err: glat, glong, i_err',GLAT,GLON,I_ERROR | ||||
| END IF | ||||
| C Convert back to Geographics as a sanity check | ||||
| CALL CONVERT_GEO_COORD(CGLAT,CGLON,R_HEIGHT_IN,GMLAT_, | ||||
| * GMLON_,2,I_ERROR) | ||||
| IF (I_ERROR.NE.0) THEN | ||||
| PRINT *,'err: mlat, mlong, i_err',CGLAT,CGLON,I_ERROR | ||||
| END IF | ||||
| C | ||||
| WRITE (1,FMT=9000) R_HEIGHT_IN,GLAT,GLON,CGLAT,CGLON, | ||||
| * GMLAT_,GMLON_ | ||||
| 10 CONTINUE | ||||
| C | ||||
| 20 CONTINUE | ||||
| CLOSE (1) | ||||
| 30 CONTINUE | ||||
| 9000 FORMAT (3f8.2,4F13.5) | ||||
| STOP | ||||
| END | ||||
