SUBROUTINE PS2LNPS(PS, ALNPS, NLON, NLAT)
      DIMENSION ALNPS(NLON, NLAT), PS(NLON, NLAT)
      ALNPS = LOG(PS)    ! Compute natural logarithm of surface pressure
      RETURN
      END