subroutine check_tvib_qual(zt,tt,tvt,ml,mgb, nl, vibsts, & idgb, ngb ) implicit none external mp_status character*100 mp_message logical mp_error common/mp_status/ mp_error, mp_message cccccccccccccccccccccccccccccccccccccccccccccccccccccccc integer*4 mgb, ml, nl,i,j,k, idgb, ngb real*4 zt, tt, tvt character*8 vibsts dimension zt(ml), tt(ml), tvt(ml,mgb,2), & vibsts(2,mgb), idgb(ngb) cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c if(mp_error) then ! tvib calculation barfed print*, 'error returned from tvib calculation' print*, mp_message c we will print out the temps and vibrational temps write(73,*) '# kinetic' do i =1, nl write(73,*) tt(i),zt(i) enddo write(73,*) '&' do k=1,2 do j=1,ngb write(73,*) '# ',idgb(j), vibsts(k,j) do i =1, nl write(73,*) tvt(i,j,k ),zt(i) enddo write(73,*) '&' enddo enddo stop endif return end