cloud0.f

Go to the documentation of this file.

      SUBROUTINE cloud0
*
*
*       Cloud initialization.
*       ---------------------
*
      include 'common6.h'
      common/clouds/  xcl(3,mcl),xdotcl(3,mcl),bodycl(mcl),rcl2(mcl),
     &                clm(mcl),clmdot(mcl),cldot,vcl,sigma,rb2,pcl2,
     &                tcl,stepcl,ncl,newcl
*
*
*       Initialize cloud variables.
      ncl = 0
      tcl = 0.0d0
      newcl = 0
      pcl2 = 0.0
*
*       Read the cloud parameters (pc, km/sec, Msun & pc).
      READ (5,*)  ncl, rb2, vcl, sigma, (clm(j),j=1,ncl),
     &            (rcl2(j),j=1,ncl)
      WRITE (6,1)  ncl, rb2, vcl, sigma
    1 FORMAT (/,12x,'CLOUDS:    NCL =',i4,'  RB =',f5.1,
     &                      '   MEAN VELOC =',f5.1,'  DISP =',f5.1)
      WRITE (6,2)  (clm(j),j=1,ncl)
    2 FORMAT (/,12x,'CLOUD MASSES:   ',1p,10e9.1)
      WRITE (6,3)  (rcl2(j),j=1,ncl)
    3 FORMAT (/,12x,'PLUMMER RADII:  ',1p,10e9.1)
      rbar1 = rbar
      IF (rbar.EQ.0.0) rbar1 = 1.0
*       Set cloud parameters in scaled units.
      rb2 = rb2/rbar1
      a1 = 0.047*sqrt(zmass*zmbar/rbar1)
*       Rms velocity of cluster members in km/sec.
      a2 = a1/sqrt(0.5d0*zmass)
*       Velocity unit.
      vcl = vcl/a2
*       Cloud velocity in scaled units.
      sigma = sigma/a2
*       Specify conservative cloud integration step using crossing time.
      stepcl = 0.002*tcr*rb2/vcl
*
*       Adopt a quantized value.
      CALL stepk(stepcl,dtn)
      stepcl = dtn
*
*       Scale radii & masses to model units.
      DO 10 j = 1,ncl
          rcl2(j) = rcl2(j)/rbar1
          clm(j) = clm(j)/zmbar
   10 CONTINUE
*
      WRITE (6,15)  rb2, vcl, sigma, stepcl
   15 FORMAT (/,12x,'SCALING:     RB =',f6.1,'  VCL =',f5.1,
     &                        '  SIGMA =',f5.1,'  STEP =',1p,e10.2)
      WRITE (6,20)  (clm(j),j=1,ncl)
   20 FORMAT (/,12x,'SCALED MASSES:  ',10f7.2)
      WRITE (6,25)  (rcl2(j),j=1,ncl)
   25 FORMAT (/,12x,'SCALED RADII:   ',10f7.2)
      cldot = 0.1*rb2/vcl
*       Time scale for 'sun-rise' is 0.05 of the cloud crossing time.
      cldot = 1.0/cldot
*
*       Define the square of cloud half-mass radii & growth times.
      DO 30 j = 1,ncl
          rcl2(j) = rcl2(j)**2
          clmdot(j) = clm(j)*cldot
   30 CONTINUE
*
*       Set square boundary radius & impact parameter.
      rb2 = rb2**2
      pcl2 = rb2
*       Define density centre for routine CLOUD.
      DO 40 k = 1,3
          rdens(k) = 0.0
   40 CONTINUE
*
*       Initialize new clouds on the boundary.
      DO 50 icl = 1,ncl
          CALL cloud(icl)
   50 CONTINUE
*
      RETURN
*
      END