*ID VOLVDC
*/
*/ Optional changes for double call to go with volvar_right added by SFBT 17/11/99.
*/ note if double call wanted then this mod must be added AFTER RFORC405
*DECLARE RAD_CTL1
*B ARE2F404.87 
     &      VOLCMASS_DIAG(P_FIELDDA),      ! Volcanic fields set to 
!     background levels by setting variable to negative value
!! Q for John Edwards -- is this necessary (i.e. could I just force 2nd call
!! to use background aerolol???
*B  ADB1F400.91
	volcmass_diag=-10 !! F90 notn set all values -ve.
*/ change to 2nd call to SW.
*/ running out of continuation lines so need to add on to a line.
*D  RFSWV.1
     &  DISS_SULPHATE(FIRST_POINT_SULPC, 1),VOLCMASS_DIAG(FIRST_POINT),     
*/ change to 2nd call to LW to add volcanic aerosol to the call.
*I RFLWV.1
     &   volcmass_diag(FP_LOCAL(I)),
*DECLARE FILL3A
*I ADB2F404.466
!! if the volcanic aerosol is < 0 then use the background values
!! done like this so as to keep bit reproducability.
!! NB 5 is hard wired as I=5 means volcanoes in this case.
	DO L=1, N_PROFILE
          LG=I_GATHER(L)
          IF ( VOLCMASS(LG) .LT. 0. ) THEN
             AEROSOL_MIX_RATIO_CLIM(L,NLEVS+1-TRINDX(LG),5) =
     &           STRAT_MASS(5) * G/
     &           ((AB(TRINDX(LG))+BB(TRINDX(LG))*PSTAR(LG))       
     &              -(AB(NLEVS+1)+BB(NLEVS+1)*PSTAR(LG)))      
	  ENDIF
        ENDDO
