Coastal & Marine Geology InfoBank

Home FACS Activities Atlas Geology School Related Sites More

USGS InfoBank program -- magredgeomag

Skip navigational links
Loading
InfoBank Programs: by Name   by Topic  
Expanded Description
Topic Description
Name
magredgeomag
Synopsis
/infobank/programs/mag/geomag/magredgeomag.for
Description
  Program to calculate IGRF regional magnetic field and residual.

       
PROMPTS
  First  input file is magnetics  (.200 or .240).

  Second input file is navigation (.050 or .070).
  Enter name of input file #1:
  Enter name of input file #2:
  Enter output file name:
  Enter last 2 digits of year:
    
INCLUDE_FILES
  drsvt

    
SUBROUTINES
  vescinit            (shared PMG external)         

  getsymbol           (shared PMG external)         
  percentage          (shared PMG external)         
  openem              (shared PMG external)         
  converttime         (shared PMG external)         
  field               (internal)                    
  forerror            (shared PMG external)         
    
FUNCTIONS
  abs                 (FORTRAN intrinsic function)  

    
FILES_USED
  unit=*(write)

  unit=(p)iinunit(read)
  unit=(p)ifilmag(read,bcksp,close)
  unit=(p)ifilnav(read,close)
  unit=(p)ifilout(write,close)
    
AUTHOR
  Clint Steele               12/26/83          .for

    based on program MRMFR.
    used subroutine named GMRF75 to calculate regional field
  modified by Carolyn Degnan, 09/14/84
    to accept ida,ihr,imin,isec,itens
  modified by Carolyn Degnan, 10/20/84
    to check for values with sensitivity to tenths of a gamma
  modified by Carolyn Degnan, 02/28/85
    to only accept new times of ida,ihr,imin,isec,itens
    to output new standard formats without navigation
    ran it thru MAT
  modified by Carolyn Degnan, 04/15/85
    to use GEOMAG subroutines and IGRF80 data
  modified by Carolyn Degnan, 11/06/85
    to use IGRFDEMO subroutines and
    IGRF45, IGRF50, IGRF55, IGRF60, DGRF65, DGRF70, DGRF75,
    DGRF80, IGRF85, & IGRF85S data
  converted to subroutine, CDegnan, 2/17/87
  converted for and ran MAT, CDegnan, 3/4/87
  added calls to PERCENTAGE, CDegnan, 6/10/87
  modified to read new gradiometer format (...,t20,f7.1,...), CDegnan, 9/15/87
  changed reference from DRS$ to PMG$, CDegnan, 2/19/92
  modified by Carolyn Degnan, 07/30/92
    to use subroutine FIELD to calculate regional value.
    FIELD is from NGDC and has all IGRF models.

    
Reads
      read ( iinunit, cfmtiyr ) iyr
      read ( ifilmag, cfmtmag2, iostat = ierr, err = 100,
     *    end = 290 ) ida, ihr, imin, isec, iten, rmstobs, rslvobs
              read ( ifilmag, cfmtmag1, iostat = ierr, err = 290,
     *            end = 800 ) ida, ihr, imin, isec, iten, observed
                  read ( ifilmag, cfmtmag2, iostat = ierr, err = 290,
     *                end = 800 ) ida, ihr, imin, isec, iten,
     *                rmstobs, rslvobs
              read ( ifilnav, cfmtnav, iostat = ierr, err = 291,
     *            end = 700 ) ida, ihr, imin, isec, iten,
     *            rlatnew, rlonnew
Writes
                      write ( ifilout, cfmtout2, iostat = ierr,
     *                    err = 92 ) ida, ihr, imin, isec, iten,
                      write ( ifilout, cfmtout1, iostat = ierr,
     *                    err = 92 ) ida, ihr, imin, isec, iten,
     *                    observed, regional, residual
Opens
Calls
      call getcentury (iyr, Century)
              call converttime ('COMPRESS', ida, ihr, imin, isec, 
     *                          iten, itmmag, ErrorFlag)
                  call converttime ('COMPRESS', ida, ihr, imin, isec, 
     *                              iten, itmmag, ErrorFlag)
          call percentage ( ifilmag )
              call converttime ('COMPRESS', ida, ihr, imin, isec, 
     *                          iten, itmnavnew, ErrorFlag)
              call field (flat, flon, elev, date,
     *                    d, i, h, x, y, z, f,
     *                    ddot, idot, hdot, xdot, ydot, zdot, fdot,
     *                    jerr)
                      call converttime ('EXPAND', ida, ihr, imin, 
     *                          isec, iten, itmmag, ErrorFlag)
                      call converttime ('EXPAND', ida, ihr, imin, 
     *                             isec, iten, itmmag, ErrorFlag)
      call forerror ( ierr )
      call forerror ( ierr )

Skip footer navigational links

Coastal and Marine Science Centers:  Pacific   St. Petersburg   Woods Hole  
InfoBank   Coastal and Marine Geology Program   Geologic Information   Ask-A-Geologist   USGS Disclaimer  


Accessibility FOIA Privacy Policies and Notices

Take Pride in America logo USA.gov logo U.S. Department of the Interior | U.S. Geological Survey
URL: http://walrus.wr.usgs.gov/infobank/programs/mag/geomag/magredgeomag.doc.html
Page Contact Information: InfoBank staff
Page Last Modified: Mon Sep 16 03:37:07 PDT 2013  (chd)