plume_model.f90

Go to the documentation of this file.

!********************************************************************************
!*********************************************************************



PROGRAM plume_model
  
  USE inpout, ONLY: initialize , read_inp , check_hysplit
  
  USE inpout, ONLY: open_file_units , close_file_units

  USE inversion, ONLY: invert_height
  
  USE rise, ONLY: plumerise 

  USE sw_umbrella, ONLY : solve_umbrella
  
  USE solver_module, ONLY: allocate_matrix

  USE variables, ONLY: hysplit_flag , inversion_flag , umbrella_flag

  USE variables, ONLY : wp

  IMPLICIT NONE
 
  REAL(wp) :: t1 , t2

  WRITE(*,*)
  WRITE(*,*)
  WRITE(*,*) '-------------------- PlumeMoM V.2 ---------------------'
  WRITE(*,*)
  WRITE(*,*) 'Created by M.de'' Michieli Vitturi(1) and F.Pardini (1)'
  WRITE(*,*) ''
  WRITE(*,*) '(1) Istituto Nazionale di Geofisica e Vulcanologia'
  WRITE(*,*) '    Sezione Pisa, Pisa, Italy'
  WRITE(*,*)
  
  CALL cpu_time(t1)

  !
  !***  Initialize the input variables
  !
  CALL initialize
  !
  !***  Read from file the input parameters
  !
  CALL read_inp
  !
  !***  Open the units for output
  !
  CALL open_file_units
  !
  !***  Allocate varaibles for the colum model
  !
  CALL allocate_matrix
  
  IF ( inversion_flag ) THEN
     
     !***  Solve the plume model
     CALL invert_height

     IF ( umbrella_flag ) THEN

        CALL solve_umbrella
        
     END IF

  ELSE

     !***  Solve the plume model
     CALL plumerise

     IF ( umbrella_flag ) THEN

        CALL solve_umbrella
        
     END IF
     
  END IF
    
  CALL close_file_units

  IF ( hysplit_flag ) CALL check_hysplit
  
  CALL cpu_time(t2)

  WRITE(*,*) 'Time taken by the code was',t2-t1,'seconds'

END PROGRAM plume_model