22 |
use temps, only: day_ref, day_ini, itau_dyn, annee_ref |
use temps, only: day_ref, day_ini, itau_dyn, annee_ref |
23 |
use ener, only: etot0, ang0, ptot0, stot0, ztot0 |
use ener, only: etot0, ang0, ptot0, stot0, ztot0 |
24 |
use advtrac_m, only: tname |
use advtrac_m, only: tname |
25 |
use netcdf95, only: nf95_open, NF90_NOWRITE, nf95_inq_varid, & |
use netcdf95, only: nf95_open, nf95_inq_varid, handle_err, NF95_CLOSE |
26 |
NF90_GET_VAR, handle_err, NF90_NOERR, NF95_CLOSE |
use netcdf, only: NF90_NOWRITE, NF90_GET_VAR, NF90_NOERR |
27 |
use nrutil, only: assert |
use nrutil, only: assert |
28 |
|
|
29 |
! Arguments: |
! Arguments: |
30 |
REAL, intent(out):: vcov(: , :), ucov(:, :), teta(:, :) |
REAL, intent(out):: vcov(: , :), ucov(:, :), teta(:, :) |
31 |
REAL, intent(out):: q(:, :, :), masse(:, :) |
REAL, intent(out):: q(:, :, :), masse(:, :) |
32 |
REAL, intent(out):: ps(:), phis(:) |
REAL, intent(out):: ps(:), phis(:, :) |
33 |
REAL, intent(out):: time |
REAL, intent(out):: time |
34 |
|
|
35 |
! Variables |
! Variables |
44 |
|
|
45 |
call assert(size(vcov, 1) == (iim + 1) * jjm, "dynetat0 vcov 1") |
call assert(size(vcov, 1) == (iim + 1) * jjm, "dynetat0 vcov 1") |
46 |
call assert((/size(ucov, 1), size(teta, 1), size(q, 1), size(masse, 1), & |
call assert((/size(ucov, 1), size(teta, 1), size(q, 1), size(masse, 1), & |
47 |
size(ps), size(phis)/) == (iim + 1) * (jjm + 1), & |
size(ps)/) == (iim + 1) * (jjm + 1), "dynetat0 (iim + 1) * (jjm + 1)") |
48 |
"dynetat0 (iim + 1) * (jjm + 1)") |
call assert(shape(phis) == (/iim + 1, jjm + 1/), "dynetat0 phis") |
49 |
call assert((/size(vcov, 2), size(ucov, 2), size(teta, 2), size(q, 2), & |
call assert((/size(vcov, 2), size(ucov, 2), size(teta, 2), size(q, 2), & |
50 |
size(masse, 2)/) == llm, "dynetat0 llm") |
size(masse, 2)/) == llm, "dynetat0 llm") |
51 |
call assert(size(q, 3) == nqmx, "dynetat0 q 3") |
call assert(size(q, 3) == nqmx, "dynetat0 q 3") |
132 |
call handle_err("dynetat0, aire", ierr, nid) |
call handle_err("dynetat0, aire", ierr, nid) |
133 |
|
|
134 |
call NF95_INQ_VARID (nid, "phisinit", nvarid) |
call NF95_INQ_VARID (nid, "phisinit", nvarid) |
135 |
ierr = NF90_GET_VAR(nid, nvarid, phis, count=(/iim + 1, jjm + 1/)) |
ierr = NF90_GET_VAR(nid, nvarid, phis) |
136 |
call handle_err("dynetat0, phisinit", ierr, nid) |
call handle_err("dynetat0, phisinit", ierr, nid) |
137 |
|
|
138 |
call NF95_INQ_VARID (nid, "temps", nvarid) |
call NF95_INQ_VARID (nid, "temps", nvarid) |