clint.f

Go to the documentation of this file.

      SUBROUTINE clint
*
*
*       Cloud integration.
*       ------------------
*
      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
      REAL*8  fc(3)
*
*
*       Check the time for advancing clouds.
      IF (tcl + stepcl.GT.time + toff) go to 30 
*
*       Set time-step and update the cloud reference time.
      dt = time + toff - tcl
      dt1 = 0.5d0*dt
      tcl = time + toff
*
*       Integrate cloud orbits in rotating coordinates with tidal effects.
      DO 20 icl = 1,ncl
          fc(1) = xcl(1,icl)*tidal(1) + xdotcl(2,icl)*tidal(4)
          fc(2) = -xdotcl(1,icl)*tidal(4)
          fc(3) = xcl(3,icl)*tidal(3)
          xcl2 = 0.0
*
          DO 10 k = 1,3
              xcl(k,icl) = (fc(k)*dt1 + xdotcl(k,icl))*dt + xcl(k,icl)
              xdotcl(k,icl) = fc(k)*dt + xdotcl(k,icl)
              xcl2 = xcl2 + (xcl(k,icl) - rdens(k))**2
   10     CONTINUE
*
*       Determine the minimum impact parameter (squared).
          pcl2 = min(xcl2,pcl2)
*
*       Check for modification of the cloud mass (increase or reduction).
          IF (xcl2.LT.rb2) THEN
              IF (bodycl(icl).LT.clm(icl)) THEN
                  bodycl(icl) = bodycl(icl) + clmdot(icl)*dt
              END IF
          ELSE
*       Reduce cloud mass gradually by 'sun-set' procedure.
              bodycl(icl) = bodycl(icl) - clmdot(icl)*dt
          END IF
*
*       Initialize a new cloud when current mass becomes negative.
          IF (bodycl(icl).LE.0.0) CALL cloud(icl)
   20 CONTINUE
*
   30 RETURN
*
      END