14 |
|
|
15 |
! This is the main procedure for the "physics" part of the program. |
! This is the main procedure for the "physics" part of the program. |
16 |
|
|
17 |
|
! Libraries: |
18 |
|
use netcdf95, only: NF95_CLOSE |
19 |
|
use nr_util, only: assert |
20 |
|
|
21 |
use aaam_bud_m, only: aaam_bud |
use aaam_bud_m, only: aaam_bud |
22 |
USE abort_gcm_m, ONLY: abort_gcm |
USE abort_gcm_m, ONLY: abort_gcm |
23 |
use ajsec_m, only: ajsec |
use ajsec_m, only: ajsec |
47 |
USE indicesol, ONLY: clnsurf, epsfra, nbsrf |
USE indicesol, ONLY: clnsurf, epsfra, nbsrf |
48 |
USE ini_histins_m, ONLY: ini_histins, nid_ins |
USE ini_histins_m, ONLY: ini_histins, nid_ins |
49 |
use lift_noro_m, only: lift_noro |
use lift_noro_m, only: lift_noro |
|
use netcdf95, only: NF95_CLOSE |
|
50 |
use newmicro_m, only: newmicro |
use newmicro_m, only: newmicro |
|
use nr_util, only: assert |
|
51 |
USE orbite_m, ONLY: orbite |
USE orbite_m, ONLY: orbite |
52 |
USE ozonecm_m, ONLY: ozonecm |
USE ozonecm_m, ONLY: ozonecm |
53 |
USE phyetat0_m, ONLY: phyetat0 |
USE phyetat0_m, ONLY: phyetat0 |
82 |
! géopotentiel de chaque couche (référence sol) |
! géopotentiel de chaque couche (référence sol) |
83 |
|
|
84 |
REAL, intent(in):: pphis(:) ! (klon) géopotentiel du sol |
REAL, intent(in):: pphis(:) ! (klon) géopotentiel du sol |
85 |
|
REAL, intent(in):: u(:, :) ! (klon, llm) zonal wind, in m / s |
86 |
REAL, intent(in):: u(:, :) ! (klon, llm) |
REAL, intent(in):: v(:, :) ! (klon, llm) meridional wind, in m / s |
|
! vitesse dans la direction X (de O a E) en m / s |
|
|
|
|
|
REAL, intent(in):: v(:, :) ! (klon, llm) vitesse Y (de S a N) en m / s |
|
87 |
REAL, intent(in):: t(:, :) ! (klon, llm) temperature (K) |
REAL, intent(in):: t(:, :) ! (klon, llm) temperature (K) |
88 |
|
|
89 |
REAL, intent(in):: qx(:, :, :) ! (klon, llm, nqmx) |
REAL, intent(in):: qx(:, :, :) ! (klon, llm, nqmx) |
356 |
REAL t_seri(klon, llm) |
REAL t_seri(klon, llm) |
357 |
real q_seri(klon, llm) ! mass fraction of water vapor |
real q_seri(klon, llm) ! mass fraction of water vapor |
358 |
REAL ql_seri(klon, llm) |
REAL ql_seri(klon, llm) |
359 |
REAL u_seri(klon, llm), v_seri(klon, llm) |
REAL u_seri(klon, llm), v_seri(klon, llm) ! wind, in m s-1 |
360 |
REAL tr_seri(klon, llm, nqmx - 2) |
REAL tr_seri(klon, llm, nqmx - 2) |
361 |
|
|
362 |
REAL zx_rh(klon, llm) |
REAL zx_rh(klon, llm) |
984 |
CALL histwrite_phy("msnow", sum(fsnow * pctsrf, dim = 2)) |
CALL histwrite_phy("msnow", sum(fsnow * pctsrf, dim = 2)) |
985 |
call histwrite_phy("qsurf", sum(fqsurf * pctsrf, dim = 2)) |
call histwrite_phy("qsurf", sum(fqsurf * pctsrf, dim = 2)) |
986 |
call histwrite_phy("flat", zxfluxlat) |
call histwrite_phy("flat", zxfluxlat) |
987 |
|
call histwrite_phy("rld", lwdn) |
988 |
|
call histwrite_phy("rldcs", lwdn0) |
989 |
|
|
990 |
DO nsrf = 1, nbsrf |
DO nsrf = 1, nbsrf |
991 |
CALL histwrite_phy("fract_"//clnsurf(nsrf), pctsrf(:, nsrf)) |
CALL histwrite_phy("fract_"//clnsurf(nsrf), pctsrf(:, nsrf)) |