Changeset 2699 for LMDZ5/trunk

Timestamp:

Nov 7, 2016, 10:23:46 PM (8 years ago)

Author:

oboucher

Message:

Fixing a small bug on the vertical discretization of the source term and removing a useless USE module

File:

• LMDZ5/trunk/libf/phylmd/StratAer/traccoag_mod.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/StratAer/traccoag_mod.F90

@@ -24 +24 @@
     USE phys_cal_mod
     USE sulfate_aer_mod
-    USE nucleate_mod
     USE phys_local_var_mod, ONLY: stratomask
     USE sulfate_aer_mod, ONLY : stracomp, denh2sa 
@@ -88 +87 @@
     REAL,DIMENSION(klev)                   :: f_lay_emiss         ! fraction of emission for every vertical layer
     REAL                                   :: f_lay_sum           ! sum of layer emission fractions
-    REAL                                   :: alt_lower           ! integral boundary for Gaussian emission profile in altitude
-    REAL                                   :: alt_upper           ! integral boundary for Gaussian emission profile in altitude
+    REAL                                   :: alt                 ! altitude for integral calculation
     INTEGER,PARAMETER                      :: n_int_alt=10        ! number of subintervals for integration over Gaussian emission profile
     REAL,DIMENSION(nbtr_bin)               :: r_bin               ! particle radius in size bin [m]
@@ -184 +182 @@
             DO k=1, klev
               f_lay_emiss(k)=0.0
-              alt_lower=altLMDz(k)
               DO i_int=1, n_int_alt
-                alt_upper=alt_lower+(altLMDz(k+1)-altLMDz(k))/float(n_int_alt)
+                alt=altLMDz(k)+float(i_int)*(altLMDz(k+1)-altLMDz(k))/float(n_int_alt)
                 f_lay_emiss(k)=f_lay_emiss(k)+1./(sqrt(2.*RPI)*sigma_alt_vol)* &
-                & (exp(-0.5*((altemiss_vol-alt_lower)/sigma_alt_vol)**2.)+     &
-                &  exp(-0.5*((alt_upper-altemiss_vol)/sigma_alt_vol)**2.))*    &
-                & (alt_upper-alt_lower)/2.
+                &              exp(-0.5*((alt-altemiss_vol)/sigma_alt_vol)**2.)*   & 
+                &              (altLMDz(k+1)-altLMDz(k))/float(n_int_alt) 
               ENDDO
               f_lay_sum=f_lay_sum+f_lay_emiss(k)
@@ -227 +223 @@
           DO k=1, klev
             f_lay_emiss(k)=0.0
-            alt_lower=altLMDz(k)
             DO i_int=1, n_int_alt
-            alt_upper=alt_lower+(altLMDz(k+1)-altLMDz(k))/float(n_int_alt)
-            f_lay_emiss(k)=f_lay_emiss(k)+1./(sqrt(2.*RPI)*sigma_alt_sai)* &
-              & (exp(-0.5*((altemiss_sai-alt_lower)/sigma_alt_sai)**2.)+   &
-              &  exp(-0.5*((alt_upper-altemiss_sai)/sigma_alt_sai)**2.))*  &
-              & (alt_upper-alt_lower)/2.
+              alt=altLMDz(k)+float(i_int)*(altLMDz(k+1)-altLMDz(k))/float(n_int_alt)
+              f_lay_emiss(k)=f_lay_emiss(k)+1./(sqrt(2.*RPI)*sigma_alt_sai)* &
+              &              exp(-0.5*((alt-altemiss_sai)/sigma_alt_sai)**2.)*   & 
+              &              (altLMDz(k+1)-altLMDz(k))/float(n_int_alt) 
             ENDDO
             f_lay_sum=f_lay_sum+f_lay_emiss(k)