parameter (nHtrn=39,nvar=14,nfix=5,ncfc=5)
      character*100 fn
      character*1 ca,cl
      character*39 cc
      character*20 cin1,cin2,cin3,cin4
      real df
      real r(nHtrn),r_(ncfc)
      real xCont(7)
      real cfcm(ncfc)
      real, allocatable :: plev(:),tlev(:),w(:,:),alt(:),w_(:,:)
      integer, allocatable :: molID(:),indx(:)
      integer molvarID(nvar),molFixID(nfix),molCFCID(ncfc)
      data molvarID/1,2,3,4,5,6,11,12,13,15,16,17,18,19/
      data molFixID/7,8,9,10,14/
      data molCFCID/20,21,22,23,24/
      data drym,cfcm/28.964,153.820,88.000,137.37, 120.910, 86.470/
      narg=iargc()
      call getarg(1,cin3)
      read(cin3,*)v1
      call getarg(2,cin4)
      read(cin4,*)v2
      call getarg(3,cin1)
      read(cin1,*)imol
      ca='A'
      cl='L'
      cc='CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC'
      xCont=(/0.0,0.0,0.0,0.0,0.0,0.0,0.0/)
      fn='tmpJG.dat'
      open(11,file=fn,status='old')
      read(11,*)nlev
      allocate (plev(nlev),tlev(nlev),alt(nlev))
      read(11,*)plev(1:nlev)
      read(11,*)alt(1:nlev)
      read(11,*)tlev(1:nlev)
      close(11)
      fn='JGuserDflt.dat'
      open(11,file=fn,status='old')
      read(11,*)nl,nmol
      allocate (molid(nmol))
      read(11,*)molid(1:nmol)
      allocate(w(nlev,nmol),w_(nlev,6))
      do n=1,nmol
         read(11,*)w(1:nlev,n)
      enddo
      if (imol == 7) w(1:nlev,1)=0.
      fn='TAPE5'
      open(21,file=fn,status='unknown')
      
      alt(:)=0.0 !fictitious set , corresponds to IMMAX < 0
      nlevm=-nlev
      v1=v1-1.0
      v2=v2+1.0
      df=0.0002
      if (v2 <= 1561.1) then
         df=0.0001
      endif
      if (v2 <= 781.1) then
         df=0.00005
      endif
      if (v2 <= 391.1) then
         df=0.00002
      endif
      
      write(21,'(a1)')'$'
      write(21,'(14(a4,i1),2(1x,i4))')' HI=',1,' F4=',1,' CN=',6,     
     &      'AE=',0,' EM=',1,' SC=',0,' FI=',0,' PL=',0,' TS=',0,
     &      ' AM=',1,' MG=',0,' LA=',0,' OD=',1,' XS=',1,0,0
      if (imol == 1) xCont(1:4)=1.0
      if (imol == 2) xCont(3:4)=1.0
      write(21,'(7(f10.3))')xCont(1:7)
      dvset=0.0
      write(21,'(3f10.3,f10.6,4f10.3,4x,i1,5x,f10.5)')v1,v2,0.0,dvset,
     &      0.0,0.0,0.0,0.0,0,df
      write(21,'(7f10.3,4x,a1)')0.0,1.0,0.0,0.0
      write(21,'(7i5,i2,1x,i2,3f10.3,10x,f10.3)')0,2,nlevm,1,1,nHtrn,
     &      1,0,0,0.0,120.0,0.5*(v1+v2),45.0
      write(21,'(5f10.3,i5,5x,f10.3)')plev(nlev),plev(1),0.0,0.0,0.0,
     &      0,0.0
      write(21,'(8F10.3)')plev(nlev:1:-1) 
      write(21,'(i5)')nlevm
      do n=nlev,1,-1
         write(21,'(f10.3,f10.4,f10.3,5x,2a1,1x,a1,1x,a39)')alt(n),
     &         plev(n),tlev(n),ca,ca,cl,cc
         w_(n,1:6)=0.0
         if (imol == 2) w_(n,2:6)=w(n,2:6)
         if (imol == 1) w_(n,1:6)=w(n,1:6)
         write(21,'(8e15.7)')w_(n,1:6)*1.e3,0.0,0.0
         write(21,'(8e15.7)')0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
         write(21,'(8e15.7)')0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
         write(21,'(8e15.7)')0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
         write(21,'(8e15.7)')0.0,0.0,0.0,0.0,0.0,0.0,0.0
      enddo
      write(21,'(3(i5))')ncfc,0,0
      write(21,'(5a10)')'CCL4','CF4','F11','F12','CHCLF2'
      write(21,'(2i5)')nlev,1
      r_(:)=0.0
      do n=nlev,1,-1
!         r_(1:ncfc)=w(n,nvar+nfix+1:nvar+nfix+ncfc)*1.e6*drym/
!     &        cfcm(1:ncfc)
         write(21,'(f10.3,5x,5a1)')plev(n),ca,ca,ca,ca,ca
         write(21,'(5(E10.3))')r_(1:ncfc)
      enddo
      write(21,'(a1)')'%'
      end