Changeset 1179 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/newmicro.F

Timestamp:

Jun 11, 2009, 4:18:47 PM (15 years ago)

Author:

jghattas

Message:

• Ajout de l'interpolation vertical pour les nouveaux fichiers de forcage des aerosols. Utilisant les anciennes fichiers de SO4 pas d'interpolation possible. Convergence numerique avec la version precedente en utilisant les anciens fichiers des SO4. aerosol_optic.F90 change du nom pour readaerosol_optic.F90 (lecture d'aerosol + optic) Les fichiers de forcage aerosol doit maintenant avoir le suffix .nc.

• Correction des bugs pour inca et certain diagnostiques optionelles de radlwsw.
• Ajout de test pour le choix advection schema.

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/newmicro.F (modified) (6 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/newmicro.F

@@ -1 @@
-!
-! $Header$
-!
+! $Id$
+!     
       SUBROUTINE newmicro (paprs, pplay,ok_newmicro,
      .                  t, pqlwp, pclc, pcltau, pclemi,
@@ -7 +6 @@
      s                  xflwp, xfiwp, xflwc, xfiwc,
      e                  ok_aie, 
-     e                  sulfate, sulfate_pi, 
+     e                  mass_ins_aero, mass_ins_aero_pi, 
      e                  bl95_b0, bl95_b1,
      s                  cldtaupi, re, fl)
@@ -22 +21 @@
 c 
 c ok_aie--input-L-apply aerosol indirect effect or not
-c sulfate-input-R-sulfate aerosol mass concentration [um/m^3]
-c sulfate_pi-input-R-dito, pre-industrial value
+c mass_ins_aero-----input-R-total mass concentration for all indissoluble aerosols[ug/m^3]
+c mass_ins_aero_pi--input-R-dito, pre-industrial value
 c bl95_b0-input-R-a parameter, may be varied for tests (s-sea, l-land)
 c bl95_b1-input-R-a parameter, may be varied for tests (    -"-      )
@@ -94 +93 @@
       LOGICAL ok_a1lwpdep       ! a1 LWP dependent?
       
-      REAL sulfate(klon, klev)  ! sulfate aerosol mass concentration [ug m-3]
+      REAL mass_ins_aero(klon, klev)    ! total mass concentration for all indissoluble aerosols [ug m-3]
+      REAL mass_ins_aero_pi(klon, klev) ! - " - (pre-industrial value)
       REAL cdnc(klon, klev)     ! cloud droplet number concentration [m-3]
       REAL re(klon, klev)       ! cloud droplet effective radius [um]
-      REAL sulfate_pi(klon, klev)  ! sulfate aerosol mass concentration [ug m-3] (pre-industrial value)
       REAL cdnc_pi(klon, klev)     ! cloud droplet number concentration [m-3] (pi value)
       REAL re_pi(klon, klev)       ! cloud droplet effective radius [um] (pi value)
@@ -157 +156 @@
                                 !             
                   cdnc(i,k) = 10.**(bl95_b0+bl95_b1*
-     &                 log(MAX(sulfate(i,k),1.e-4))/log(10.))*1.e6 !-m-3
+     &                 log(MAX(mass_ins_aero(i,k),1.e-4))/log(10.))*1.e6 !-m-3
                                 ! Cloud droplet number concentration (CDNC) is restricted
                                 ! to be within [20, 1000 cm^3]
@@ -165 +164 @@
                                 !
                   cdnc_pi(i,k) = 10.**(bl95_b0+bl95_b1*
-     &                 log(MAX(sulfate_pi(i,k),1.e-4))/log(10.))*1.e6 !-m-3
+     &                 log(MAX(mass_ins_aero_pi(i,k),1.e-4))/log(10.))*1.e6 !-m-3
                   cdnc_pi(i,k)=MIN(1000.e6,MAX(20.e6,cdnc_pi(i,k)))
                ENDDO