subroutine anthrP

      implicit none


      character(120) :: fname, new_file_name

      integer iou
      logical exists

      include "stdunits.h"
      include "size.h"
      include "isleperim.h"
      include "anthrP.h"

      real c1, c0
      real tmpijk(101)
      c1 = 1.
      c0 = 0.


#  if defined O_anthrP
      fname = new_file_name ("FutureofP.nc")
      inquire (file=trim(fname), exist=exists)
      if (.not. exists) then
         print*, "Error => ", trim(fname), " does not exist."
         stop 'topog in topog.f'
      endif
      call openfile (fname, iou)
      if (.not. exists) then
         print*, "==> Warning: FutureofP data does not exist."
      else
         call getvara ('anthrP', iou, 101, 1, 101
     &,             tmpijk, c1, c0)
      endif
      anthrP_data = tmpijk

         call getvara ('time', iou, 101, 1, 101
     &,             tmpijk, c1, c0)
      anthrP_time = tmpijk      
!      print*, anthrP_time

#  endif
      return
      end