Changeset 682

Timestamp:

Mar 24, 2006, 4:06:23 PM (19 years ago)

Author:

lmdzadmin

Message:

Ajout de bilans pour INCA - AC
MAF

File:

• LMDZ4/trunk/libf/phylmd/phytrac.F (modified) (1 diff)

LMDZ4/trunk/libf/phylmd/phytrac.F

@@ -552 +552 @@
      $                 jjp1)       !ny
 #ifdef INCAINFO
+#ifdef INCA_AER
+
+c Budget calculation for aerosol species
+CALL tbudget(airephy,pdtphys,nstep,tr_seri,.false.)
+
+c-- summary info----------------------------------------------------------------
+
+if (MOD(nstep,nint(86400./pdtphys)) .eq. 0) then
+print *, "global aerosol optical thickness "
+
+write (form,'(A,I2,A)') "(A,",trnx-trmx+1,"A10)"
+print form,"lamba [nm] ", (solsym(it),it=trmx,trnx)
+
+write (form,'(A,I2,A)') "(I11,",trnx-trmx+1,"F10.4)"
+do i=1,las
+print form,int(lambda(i)),(sum(tausum(:,i,it)*airephy)/sum(airephy),it=trmx,trnx)
+enddo
+
+print *,"global mean angstroem component ", sum(angst*airephy)/sum(airephy)
+endif
+#endif
+#endif
+
+#ifdef INCAINFO
       PRINT *, 'OK.'
       DO it = 1, nbtr