Changeset 1183 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/nuage.F

Timestamp:

Jun 16, 2009, 5:38:46 PM (15 years ago)

Author:

jghattas

Message:

• Ajoute test de coherence entre les longitudes et latitudes des fichiers de forcage des aerosols et le model.

• Corection pour la variable mass_solu_aero : changement du nom, commentaire et contenu(il y avait l'aerosol dust en trop)

• Changement du nom de module car dans INCA le meme nom existe deja : aerosol_mod => aero_mod

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/nuage.F (modified) (5 diffs)

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

@@ -5 +5 @@
      .                  pch, pcl, pcm, pct, pctlwp,
      e                  ok_aie,
-     e                  mass_ins_aero, mass_ins_aero_pi, 
+     e                  mass_solu_aero, mass_solu_aero_pi, 
      e                  bl95_b0, bl95_b1,
      s                  cldtaupi, re, fl)
@@ -19 +19 @@
 c pclc----input-R-couverture nuageuse pour le rayonnement (0 a 1)
 c ok_aie--input-L-apply aerosol indirect effect or not
-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 mass_solu_aero-----input-R-total mass concentration for all soluble aerosols[ug/m^3]
+c mass_solu_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 (    -"-      )
@@ -73 +73 @@
       LOGICAL ok_aie            ! Apply AIE or not?
       
-      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 mass_solu_aero(klon, klev)    ! total mass concentration for all soluble aerosols[ug m-3]
+      REAL mass_solu_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]
@@ -107 +107 @@
             !             
             cdnc(i,k) = 10.**(bl95_b0+bl95_b1*
-     .           log(MAX(mass_ins_aero(i,k),1.e-4))/log(10.))*1.e6 !-m-3
+     .           log(MAX(mass_solu_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]
@@ -113 +113 @@
             cdnc(i,k)=MIN(1000.e6,MAX(20.e6,cdnc(i,k)))
             cdnc_pi(i,k) = 10.**(bl95_b0+bl95_b1*
-     .           log(MAX(mass_ins_aero_pi(i,k),1.e-4))/log(10.))*1.e6 !-m-3
+     .           log(MAX(mass_solu_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)))
             !