Education Center | Photo Gallery

Data Advanced Data Search More Data Discovery & Access Easy-to-use Data Products Data Collections Data Software & Programming Tools Data Submission Data News Centers & Programs Antarctic Glaciological Data Center GLIMS Glacier Database IARC Frozen Ground Data Center International Polar Year Data & Information Service (IPYDIS) NOAA at NSIDC NASA's Distributed Active Archive Center at NSIDC (NSIDC DAAC) World Data Center for Glaciology, Boulder Science Scientists Research Projects Publications NSIDC Notes Glaciological Data Series Special Reports Posters & Presentations Annual Reports Education Resources News & Events Press Room NSIDC In the News Data News Events About Expertise Data & Services NSIDC Sponsors Contacts Jobs Related Organizations

Coordinates Conversion Program II

The following program converts from latitude and longitude coordinates to polar stereographic SSM/I grid coordinates to pixel sample lines and indexes. Program I converts from pixel sample and line indexes to polar stereographic SSM/I grid coordinates to latitude and longitude coordinates.

c      Convert_11_xy
c
c      This program converts from latitude and longitude
c      to polar stereographic SSM/I grid coordinates (in km)
c      and to Leads ARI grid sample,line indexes.  This program is a
c      version of the program written by C.S. Morris, Jet Propulsion
c      Laboratory,  which appears in the NSIDC SSM/I User's
c      Guide, Section F, "SSMII Polar Grids".
c
c      Program-adapted for Leads ARI grids by F.Fetterer,
c      NOARL Remote Sensing,Stennis Space Center, March 1991

       dimension t (2)
       real*4 x,y,alat,along
       real*8 e,e2,cdr,pi
       character*1 grid

c      Conversion constant, degrees to radians:
       cdr=57.29577951
c      Radius and eccentricity of earth (Hughes ellipsiod, km):
       re=6378.273
       e2=0.006693883
       e=SQRT(e2)
       pi=3.141592654
c      Standard parallel - latitude of no distortion
       slat=70.0
c      Puts 135 deg long. at the top of the grid:
       xlam=-45.

       print*,'Type ''P'' for Pacific grid, ''E'' for European'
       accept*,grid
100    print*,'Type latitude, longitude (9999,9999 to stop)'
       accept*,alat,along
       if (alat.eq.9999.) goto 999
       if(alat.lt.89.995)goto 250
       x=0.
       y=0.
       goto  300
250    rlat=alat
       do 600 II=1,2
        if(II.eq.2)rlat=slat
600    t(II)=TAN((pi/4.)-(rlat/(2.*cdr)))/((1.-e*SIND(rlat))/
     & (1.+e*SIND(rlat)))**(e/2.)

       cm=COSD(slat)/SQRT(1.-e2*(SIND(slat)**2))
       rho-re*cm*t(1)/t(2)
       x=rho*SIND(along-xlam)
       y=rho*COSD(along-xlam)
c      Transform sample, line coordinates to SSMI grid coordinates
300    if(grid.eq.'P'.or.grid.eq.'p')then
        rline=1975-y
        sample=x+2250
       else
        rline=1300-y
        sample=x
       end if
       print*,'    Sample    Line       SSMI x      SSMI y    Lat
     &       Long'
       write (6,310)sample,rline,x,y,alat,along
       print*
       goto  100
310    format(4f1O.2,2f1O.4)
999    continue
       end

The program is taken from:
Fetterer and Hawkins. 1991. An AVHRR data set for the Arctic leads ARI. Naval Oceanographic and Atmospheric Research Laboratory, Technical Note 118. 61 p.