c*********************************************************************** subroutine mzescape_fh(ig) c jul 2011 malv+fgg c*********************************************************************** implicit none include 'nltedefs.h' include 'nlte_matrix.h' include 'nlte_atm.h' include 'tcr_15um.h' c local variables integer i, ib, ik, istyle integer ig ! ADDED FOR TRACEBACK real*8 tau_factor real*8 aux(nl), aux2(nl), aux3(nl) c*********************************************************************** call zero4v( aux, taustar12,tauinf121,tauii121, nl) do ik=1,3 ib=ik+1 call mzescape ( ig,aux,aux2,aux3, ib, 1,irw_mztf,imu ) tau_factor = 1.d0 if (ik.eq.1) tau_factor = dble(667.75/618.03) if (ik.eq.3) tau_factor = dble(667.75/720.806) do i=1,nl taustar12(i) = taustar12(i) + aux(i) * tau_factor tauinf121(i) = tauinf121(i) + aux2(i) * tau_factor tauii121(i) = tauii121(i) + aux3(i) * tau_factor enddo enddo istyle = 2 call mzescape_normaliz ( taustar12, istyle ) c end return end