utils.F90

Go to the documentation of this file.

 module utils

 implicit none

 public

 character(len=512) :: topo_file = "orog"
 character(len=128) :: topo_field = "orog_filt"
 character(len=128) :: mask_field = "slmsk"
 character(len=512) :: grid_file = "atmos_mosaic.nc"

 logical :: zero_ocean = .true.
 logical :: nested = .false.
 logical :: regional = .false.

 integer :: grid_type = 0
 
 real    :: stretch_fac = 1.0
 real    :: res = 48.

 contains

 subroutine read_namelist

 implicit none

 integer :: stdunit = 6, unit=7, io_status
 logical :: opened

 namelist /filter_topo_nml/ topo_file, topo_field, mask_field, grid_file, zero_ocean, &
        stretch_fac, res, nested, grid_type, regional

 do
   inquire( unit=unit, opened=opened )
   if( .NOT.opened )exit
   unit = unit + 1
   if( unit.EQ.100 )call handle_err(-1, 'Unable to locate unit number.' )
 end do

 open( unit=unit, file='input.nml', iostat=io_status )
 read( unit,filter_topo_nml, iostat=io_status )
 close(unit)

 if (io_status > 0) call handle_err(-1, 'Error reading input.nml')

 write (stdunit, nml=filter_topo_nml)

 end subroutine read_namelist

 subroutine handle_err(status, string)

 implicit none

#include <netcdf.inc>

 integer,          intent(in) :: status
 character(len=*), intent(in) :: string
 character(len=256) :: errmsg

 if (status .ne. nf_noerr) then
   errmsg = nf_strerror(status)
   errmsg = trim(errmsg) // " " // trim(string)
   print *, "FATAL ERROR:"
   print *, trim(errmsg)
   error stop 'Stopped'
 endif

 end subroutine handle_err

 end module utils