Changeset 1216 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90

Timestamp:

Jul 28, 2009, 10:05:59 AM (15 years ago)

Author:

jghattas

Message:

Bug-fix dans le calcul de aerosol charge(load) dans la colonne pour
l'interpolation des aerosols.

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90 (modified) (8 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90

@@ -1 +1 @@
 ! $Id$
 !
-SUBROUTINE readaerosol_interp(id_aero, r_day, first, pplay, paprs, mass_out, pi_mass_out)
+SUBROUTINE readaerosol_interp(id_aero, r_day, first, pplay, paprs, t_seri, mass_out, pi_mass_out)
 !
 ! This routine will return the mass concentration at actual day(mass_out) and 
@@ -36 +36 @@
   REAL, DIMENSION(klon,klev), INTENT(IN) :: pplay  ! pression at model mid-layers
   REAL, DIMENSION(klon,klev+1),INTENT(IN):: paprs  ! pression between model layers
+  REAL, DIMENSION(klon,klev), INTENT(IN) :: t_seri ! air temperature
 !      
 ! Output:      
@@ -51 +52 @@
 !$OMP THREADPRIVATE(klev_src)
 
+  REAL                              :: zrho      ! Air density [kg/m3]
+  REAL                              :: volm      ! Volyme de melange [kg/kg]
   REAL, DIMENSION(klon)             :: psurf_day, pi_psurf_day
   REAL, DIMENSION(klon)             :: load_src, pi_load_src  ! Mass load at source grid
@@ -149 +152 @@
     
      IF (debug) THEN
+        CALL writefield_phy('var_year_jan',var_year(:,:,1,id_aero),klev_src)
+        CALL writefield_phy('var_year_dec',var_year(:,:,12,id_aero),klev_src)
         CALL writefield_phy('psurf_src',psurf_year(:,:,id_aero),1)
         CALL writefield_phy('pi_psurf_src',pi_psurf_year(:,:,id_aero),1)
@@ -300 +305 @@
         DO k= 1, klev
            DO i = 1, klon
-              load_tgt(i) = load_tgt(i) + 1/RG * var_day(i,k,id_aero)*delp(i,k)
+              zrho = pplay(i,k)/t_seri(i,k)/RD       ! [kg/m3]
+              volm = var_day(i,k,id_aero)*1.E-9/zrho ! [kg/kg]
+              load_tgt(i) = load_tgt(i) + 1/RG * volm *delp(i,k)
            END DO
         END DO
@@ -315 +322 @@
            DO k= 1, klev
               DO i = 1, klon
-                 load_tgt_test(i) = load_tgt_test(i) + 1/RG * var_day(i,k,id_aero)*delp(i,k)
+                 zrho = pplay(i,k)/t_seri(i,k)/RD       ! [kg/m3]
+                 volm = var_day(i,k,id_aero)*1.E-9/zrho ! [kg/kg]
+                 load_tgt_test(i) = load_tgt_test(i) + 1/RG * volm*delp(i,k)
               END DO
            END DO
@@ -353 +362 @@
         DO k = 1, klev
            DO i = 1, klon
-              load_tgt(i) = load_tgt(i) + 1/RG * pi_var_day(i,k,id_aero)*delp(i,k)
+              zrho = pplay(i,k)/t_seri(i,k)/RD          ! [kg/m3]
+              volm = pi_var_day(i,k,id_aero)*1.E-9/zrho ! [kg/kg]
+              load_tgt(i) = load_tgt(i) + 1/RG * volm * delp(i,k)
            END DO
         END DO
@@ -367 +378 @@
            DO k = 1, klev
               DO i = 1, klon
-                 load_tgt_test(i) = load_tgt_test(i) + 1/RG * pi_var_day(i,k,id_aero)*delp(i,k)
+                 zrho = pplay(i,k)/t_seri(i,k)/RD          ! [kg/m3]
+                 volm = pi_var_day(i,k,id_aero)*1.E-9/zrho ! [kg/kg]
+                 load_tgt_test(i) = load_tgt_test(i) + 1/RG * volm * delp(i,k)
               END DO
            END DO