SUBROUTINE AEROSOLS
C  NOTE:  only difference between this and the original aerosols.f is the inclusion of the
C         base1 directory name in the read statement.   

#include "com2d.h"
	common/aeros/aerosol_all(2,L$,Z$)
	DIMENSION AEROREAD(2,L$,Z$)
	OPEN (unit=21,NAME='aerosol.dat',
     c		type='OLD', form='FORMATTED')

	do 100 ii=1,2
	do 100 ij=1,L$
		read(21,6663)(aeroread(ii,ij,ik),ik=1,30)
100	continue

	CLOSE(unit=21)

	do 200 ii=1,2
	do 200 ij=1,L$
	  iu=19-ij
	do 300 ik=1,30
	   aerosol_all(ii,iu,ik)=aeroread(ii,ij,ik)
300	continue
200	continue

6663  format(10f6.2)

	IF(Z$ .GT. 30)THEN
		DO 110 II=1,2
		DO 110 IK=31,Z$
		DO 110 IJ=1,L$
			AEROSOL_ALL(II,IJ,IK)=0.0
110		CONTINUE
	ENDIF

	RETURN
	END