Page History

            izlat = int( zlat * 100.0_JPRB)
            izlon = int( zlon * 100.0_JPRB)
            i_refcoord = sign(1,izlat) * (abs(izlat) * 100000 + izlon)
            OBS_RAWTABLE(jpraw_refcoord, jtotal) =&
              &sign(1,i_refcoord) * (abs(i_refcoord) +&
              &1 * 0.01_JPRB)

!**   Used to compute averages (always in the same order)

call KEYSORT(irc, POBS, KLEN,key=jpraw_refcoord,transposed=.TRUE.)

!**   (Lat, Lon) of observations

jlatN = NDGL + 1
LoopAllObs: do j=1,KLEN
  iflag = 0
  ZOBSLAT = POBS(kplat,j)
  ZOBSLON = mod(POBS(kplon,j) + 360.0_JPRB, 360.0_JPRB)

!**   Calculate latitudinal NW-pivot values for each observation (i.e. point#1)


  WHILE_LOOP: do while (jlatN > 0.and. ZOBSLAT > ALAT(jlatN))
    jlatN = jlatN - 1
  enddo WHILE_LOOP