source: LMDZ5/trunk/libf/phylmd/Dust/checkmass.F90 @ 5098

Last change on this file since 5098 was 2630, checked in by fhourdin, 8 years ago

Importation du modèle d'aérosols de Boucher, Escribano et al.

File size: 1.2 KB
Line 
1SUBROUTINE checkmass(zq,RNAVOG,mass,zdz,pplay,t_seri,iscm3,comment)
2  USE dimphy
3USE geometry_mod , ONLY:cell_area
4  IMPLICIT NONE
5
6#include "YOMCST.h"
7
8! Entrees
9  REAL,DIMENSION(klon,klev), INTENT(IN)   :: zq
10  REAL,DIMENSION(klon,klev), INTENT(IN)   :: zdz
11  REAL,INTENT(IN)                         :: mass
12  REAL,INTENT(IN)                         :: RNAVOG
13  CHARACTER(LEN=*),INTENT(IN)             :: comment
14REAL,DIMENSION(klon,klev) ,INTENT(IN)  :: pplay
15REAL,DIMENSION(klon,klev) ,INTENT(IN) :: t_seri
16LOGICAL :: iscm3
17
18
19! Local 
20!  INTEGER,DIMENSION(klon)     :: jadrs
21  INTEGER                     :: i, k
22  REAL                              :: totmass
23  REAL,DIMENSION(klon)              ::  burden
24  REAL,DIMENSION(klon,klev)   :: zqv
25
26
27DO i=1, klon
28burden(i)=0.0
29  DO k=1, klev
30zqv(i,k)=zq(i,k)
31
32ENDDO
33
34ENDDO
35
36totmass=0.0
37
38
39IF (iscm3) THEN
40
41ELSE
42
43CALL kg_to_cm3(pplay,t_seri,zqv)
44
45ENDIF
46
47  DO k=1, klev
48      DO i=1, klon
49        burden(i)=burden(i)+(zqv(i,k)*1.e6*zdz(i,k)* &
50                 mass*1.e3/RNAVOG)     !--mg S/m2
51    ENDDO
52  ENDDO
53    DO i=1, klon
54        totmass=totmass+burden(i)*cell_area(i)
55ENDDO
56WRITE(*,9999) comment," totmass= ",totmass
579999  format(a20,a10,e15.8)
58
59END SUBROUTINE checkmass
Note: See TracBrowser for help on using the repository browser.