source: LMDZ6/trunk/libf/phylmd/Dust/cm3_to_kg.f90

Last change on this file was 5274, checked in by abarral, 9 hours ago

Replace yomcst.h by existing module

File size: 1.3 KB
RevLine 
[5246]1SUBROUTINE cm3_to_kg(pplay,t_seri,tr_seri)
[2630]2
[5246]3  USE dimphy
4  USE infotrac
5  USE indice_sol_mod
[2630]6
[5271]7  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
[5274]8USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
9          , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
10          , R_ecc, R_peri, R_incl                                      &
11          , RA, RG, R1SA                                         &
12          , RSIGMA                                                     &
13          , R, RMD, RMV, RD, RV, RCPD                    &
14          , RMO3, RMCO2, RMC, RMCH4, RMN2O, RMCFC11, RMCFC12        &
15          , RCPV, RCVD, RCVV, RKAPPA, RETV, eps_w                    &
16          , RCW, RCS                                                 &
17          , RLVTT, RLSTT, RLMLT, RTT, RATM                           &
18          , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS            &
19          , RALPD, RBETD, RGAMD
[5271]20IMPLICIT NONE
[5246]21  !
[5271]22
[5274]23
[5246]24  !
25  REAL :: t_seri(klon,klev), pplay(klon,klev)
26  REAL :: tr_seri(klon,klev)
27  REAL :: zrho
28  INTEGER :: i, k
29  !
30  !JE20150707      RD = 1000.0 * 1.380658E-23 * 6.0221367E+23 / 28.9644
31  DO k = 1, klev
32  DO i = 1, klon
33    zrho=pplay(i,k)/t_seri(i,k)/RD
34    tr_seri(i,k)=tr_seri(i,k)*1.e6/zrho
35  ENDDO
36  ENDDO
37  !
38END SUBROUTINE cm3_to_kg
Note: See TracBrowser for help on using the repository browser.