Changeset 720 for LMDZ4/branches/V3_test/libf/phylmd

Timestamp:

Oct 6, 2006, 12:38:30 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Corrections bugs divers YM
LF

File:

• LMDZ4/branches/V3_test/libf/phylmd/physiq.F (modified) (2 diffs)

LMDZ4/branches/V3_test/libf/phylmd/physiq.F

@@ -2880 +2880 @@
 
 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 +3013 @@
      . 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