c*********************************************************************** subroutine mztvc_626fh(ig) c jul 2011 malv+fgg c*********************************************************************** implicit none !!!!!!!!!!!!!!!!!!!!!!! ! common variables & constants include 'nltedefs.h' include 'nlte_atm.h' include 'nlte_data.h' include 'tcr_15um.h' include 'nlte_matrix.h' include 'nlte_curtis.h' !!!!!!!!!!!!!!!!!!!!!!! ! arguments integer ig ! ADDED FOR TRACEBACK !!!!!!!!!!!!!!!!!!!!!!! ! local variables real*4 cdummy(nl,nl), csngl(nl,nl) real*8 cax1(nl,nl), cax2(nl,nl), cax3(nl,nl) real*8 v1(nl), v2(nl), v3(nl), cm_factor, vc_factor integer itauout,icfout,itableout, interpol,ismooth, isngldble integer i,j,ik,ist,isot,ib,itt !character bandcode*2 character isotcode*2 !character codmatrx_hot*5 !!!!!!!!!!!!!!!!!!!!!!! ! external functions external bandid character*2 bandid !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! subroutines called: ! mz4sub, dmzout, readc_mz4, readcupdw, mztf !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! formatos 132 format(i2) ************************************************************************ ************************************************************************ isngldble = 1 ! =1 --> dble precission fileroot = 'cfl' ist = 1 isot = 26 write (isotcode,132) isot call zerov( vc121, nl ) do 11, ik=1,3 ib=ik+1 call mztvc (ig,v1, ib, 1, irw_mztf, imu, 0,0,0 ) do i=1,nl if(ik.eq.1)then vc_factor = dble(667.75/618.03) elseif(ik.eq.2)then vc_factor = 1.d0 elseif(ik.eq.3)then vc_factor = dble(667.75/720.806) end if vc121(i) = vc121(i) + v1(i) * vc_factor end do 11 continue return end