source: LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/cm3_to_kg.F

Last change on this file was 2175, checked in by jescribano, 10 years ago

SPLA code included for first time

File size: 563 bytes
Line 
1      SUBROUTINE cm3_to_kg(pplay,t_seri,tr_seri)
2
3      USE dimphy
4      USE infotrac
5      USE indice_sol_mod
6
7      IMPLICIT NONE
8c
9#include "dimensions.h"
10#include "../phylmd/YOMCST.h"
11c #include "../phylmd/dimphy.h"
12c     
13      REAL t_seri(klon,klev), pplay(klon,klev)
14      REAL tr_seri(klon,klev)
15      REAL zrho
16      INTEGER i, k
17c
18      RD = 1000.0 * 1.380658E-23 * 6.0221367E+23 / 28.9644
19      DO k = 1, klev
20      DO i = 1, klon
21        zrho=pplay(i,k)/t_seri(i,k)/RD
22        tr_seri(i,k)=tr_seri(i,k)*1.e6/zrho
23      ENDDO
24      ENDDO
25c
26      END
Note: See TracBrowser for help on using the repository browser.