fabm_tavg_def_4d Subroutine

      subroutine fabm_tavg_def_4d(iou,it)
         !! define FABM variables in NetCDF context
      integer, intent(in) :: iou
         !! NetCDF file id
      integer, intent(in) :: it(:)
         !! array of dimension ids

         integer :: j,k,n

         real, parameter :: c0 = 0.
         real, parameter :: c1 = 1.
         real, parameter :: c100 = 100.
         real, parameter :: c500 = 500.

         associate(VAR => model%interior_state_variables)
         do n=1,npelagic
!            call defvar ('O_'//trim(mapt(2+n)), iou, 4, it,
            print*, 'O_'//trim(VAR(n)%name)
            call defvar ('O_'//trim(VAR(n)%name), iou, 4, it,
     &                 VAR(n)%minimum, VAR(n)%maximum, ' ', 'F',
     &                 trim(VAR(n)%long_name), ' ', trim(VAR(n)%units))
         end do
         end associate

      end subroutine fabm_tavg_def_4d