Index: /trunk/LMDZ.MARS/util/aeroptical.F90
===================================================================
--- /trunk/LMDZ.MARS/util/aeroptical.F90	(revision 2317)
+++ /trunk/LMDZ.MARS/util/aeroptical.F90	(revision 2318)
@@ -1,23 +1,22 @@
 program aeroptical
 ! program for computation of aerosol opacities
-! author : Antoine Bierjon, April 2020
+! author : Antoine Bierjon, April-May 2020
 !
 !===============================================================================
 !     PREFACE
 !===============================================================================
-! This program takes in input a GCM output file (diagfi,stats,concat) that
+! This program takes as inputs a GCM output file (diagfi,stats,concat) that
 ! contains:
 !      - the Mass Mixing Ratios of dust (dustq) and water ice (h2o_ice)
 !      - their effective radius (reffdust, reffice(*))
 !      - the atmospheric density (rho)
-! as well as the 2 boundaries of a wavelength interval and the directory
-! containing the ASCII files of the optical properties of the aerosols.
+! as well as a wavelength to calibrate the opacities, and the directory
+! containing the ASCII files of the aerosols' optical properties.
 !
 ! It outputs a netcdf file containing the opacities [1/km] of the dust and
-! water ice aerosols as well as the combined opacity of the 2 aerosols,
-! averaged within the prescribed wavelength interval.
+! water ice aerosols calibrated at the prescribed wavelength.
 !
 !
-! (*) due to a high uncertainty of reffice in the gcm, the value is asked
+! (*) due to a high uncertainty of reffice in the GCM, the value is asked
 ! directly to the user (if the user returns 0, then the program reads the GCM
 ! file to get reffice)
@@ -26,5 +25,5 @@
 ! 1) if one wanted to add another aerosol to compute, one should look for
 ! the comments + NEW AER? that are disseminated all along the code to indicate
-! the parts of the code that must be modified.
+! the parts of the code that should be modified.
 ! 2) another enhancement of this program could be the possibility to read its
 ! own product files and recalibrate the opacities at another wavelength
@@ -46,5 +45,5 @@
 integer :: ierr                               ! [netcdf] subroutine returned error code
 character(len=256) :: error_text              ! text to display in case of error
-integer :: tmpvarid                           ! temporary store a variable ID number
+integer :: tmpvarid                           ! temporarily stores a variable ID number
 integer :: lonlen,latlen,altlen,timelen       ! nb of grid points along longitude,latitude,altitude,time
 integer :: GCM_layers                         ! number of GCM layers
@@ -57,5 +56,5 @@
 real,dimension(:,:,:,:),allocatable :: rho    ! atmospheric density [kg.m-3]
 integer :: naerkind                           ! nb of aerosols to consider
-integer :: iaer                               ! aerosol kind index (1=dust, 2=water ice)
+integer :: iaer                               ! aerosol kind index (1=dust, 2=water ice) ! + NEW AER?
 integer :: ilon,ilat,ialt,it                  ! Loop indices for coordinates
 
@@ -69,10 +68,10 @@
                                                  ! altdimout: stores altitude dimension ID number
                                                  ! timedimout: stores time dimension ID number
-integer :: tmpvaridout                           ! temporary stores a variable ID number
+integer :: tmpvaridout                           ! temporarily stores a variable ID number
 
 real :: wvl_val                                  ! reference wavelength of the output opacities (given by the user) [m]
 integer :: varshape(4)                           ! stores a variable shape (of dimensions' IDs)
-character(len=16) :: tmpvarname                  ! temporary stores a variable name
-real,dimension(:,:,:,:,:),allocatable :: opa_aer ! Opacity of the aerosols [1/km]
+character(len=16) :: tmpvarname                  ! temporarily stores a variable name
+real,dimension(:,:,:,:,:),allocatable :: opa_aer ! Aerosols opacities [1/km]
 real :: missval                                  ! Value to put in outfile when the reff is out of bounds
 PARAMETER(missval=1E+20)
@@ -100,5 +99,5 @@
 ! Opacity computation
 integer :: iwvl1,iwvl2,isize1,isize2     ! Wavelength and Particle size indices for the interpolation
-real :: coeff                            ! Interpolation coeficient
+real :: coeff                            ! Interpolation coefficient
 real,dimension(2) :: reffQext            ! Qext after reff interpolation
 real :: Qext_val                         ! Qext after both interpolations