Index: LMDZ4/branches/V3_test/libf/phylmd/physiq.F
===================================================================
--- LMDZ4/branches/V3_test/libf/phylmd/physiq.F	(revision 719)
+++ LMDZ4/branches/V3_test/libf/phylmd/physiq.F	(revision 720)
@@ -2880,4 +2880,23 @@
 
 c
+cjq - introduce the aerosol direct and first indirect radiative forcings
+cjq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
+      IF (ok_ade.OR.ok_aie) THEN
+         ! Get sulfate aerosol distribution
+         CALL readsulfate(rjourvrai, debut, sulfate)
+         CALL readsulfate_preind(rjourvrai, debut, sulfate_pi)
+
+         ! Calculate aerosol optical properties (Olivier Boucher)
+         CALL aeropt(pplay, paprs, t_seri, sulfate, rhcl,
+     .        tau_ae, piz_ae, cg_ae, aerindex)
+cym
+      ELSE
+        tau_ae(:,:,:)=0.0
+	piz_ae(:,:,:)=0.0
+	cg_ae(:,:,:)=0.0
+cym     
+      ENDIF
+
+c
 cIM calcul nuages par le simulateur ISCCP
 c
@@ -2994,22 +3013,5 @@
      . EXP((Lheat *qsat2m(i))/(RCPD*zt2m(i)))
       ENDDO
-c
-cjq - introduce the aerosol direct and first indirect radiative forcings
-cjq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
-      IF (ok_ade.OR.ok_aie) THEN
-         ! Get sulfate aerosol distribution
-         CALL readsulfate(rjourvrai, debut, sulfate)
-         CALL readsulfate_preind(rjourvrai, debut, sulfate_pi)
-
-         ! Calculate aerosol optical properties (Olivier Boucher)
-         CALL aeropt(pplay, paprs, t_seri, sulfate, rhcl,
-     .        tau_ae, piz_ae, cg_ae, aerindex)
-cym
-      ELSE
-        tau_ae(:,:,:)=0.0
-	piz_ae(:,:,:)=0.0
-	cg_ae(:,:,:)=0.0
-cym     
-      ENDIF
+
 
 #ifdef INCA