|
Last change
on this file since 5159 was
5159,
checked in by abarral, 7 weeks ago
|
|
Put dimensions.h and paramet.h into modules
|
|
File size:
538 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 | USE lmdz_yomcst |
|---|
| 7 | |
|---|
| 8 | USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm |
|---|
| 9 | IMPLICIT NONE |
|---|
| 10 | |
|---|
| 11 | |
|---|
| 12 | |
|---|
| 13 | REAL :: t_seri(klon, klev), pplay(klon, klev) |
|---|
| 14 | REAL :: tr_seri(klon, klev) |
|---|
| 15 | REAL :: zrho |
|---|
| 16 | INTEGER :: i, k |
|---|
| 17 | |
|---|
| 18 | !JE20150707 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 |
|---|
| 25 | |
|---|
| 26 | END SUBROUTINE cm3_to_kg |
|---|
Note: See
TracBrowser
for help on using the repository browser.
