/[lmdze]/trunk/dyn3d/dynetat0.f
ViewVC logotype

Diff of /trunk/dyn3d/dynetat0.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 10 by guez, Fri Apr 18 14:45:53 2008 UTC
# Line 22  contains Line 22  contains
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(:) ! in Pa
33        REAL, intent(out):: phis(:, :)
34      REAL, intent(out):: time      REAL, intent(out):: time
35    
36      !   Variables      !   Variables
# Line 44  contains Line 45  contains
45    
46      call assert(size(vcov, 1) == (iim + 1) * jjm, "dynetat0 vcov 1")      call assert(size(vcov, 1) == (iim + 1) * jjm, "dynetat0 vcov 1")
47      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), &
48           size(ps), size(phis)/) == (iim + 1) * (jjm + 1), &           size(ps)/) == (iim + 1) * (jjm + 1), "dynetat0 (iim + 1) * (jjm + 1)")
49           "dynetat0 (iim + 1) * (jjm + 1)")      call assert(shape(phis) == (/iim + 1, jjm + 1/), "dynetat0 phis")
50      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), &
51           size(masse, 2)/) == llm, "dynetat0 llm")           size(masse, 2)/) == llm, "dynetat0 llm")
52      call assert(size(q, 3) == nqmx, "dynetat0 q 3")      call assert(size(q, 3) == nqmx, "dynetat0 q 3")
# Line 132  contains Line 133  contains
133      call handle_err("dynetat0, aire", ierr, nid)      call handle_err("dynetat0, aire", ierr, nid)
134    
135      call NF95_INQ_VARID (nid, "phisinit", nvarid)      call NF95_INQ_VARID (nid, "phisinit", nvarid)
136      ierr = NF90_GET_VAR(nid, nvarid, phis, count=(/iim + 1, jjm + 1/))      ierr = NF90_GET_VAR(nid, nvarid, phis)
137      call handle_err("dynetat0, phisinit", ierr, nid)      call handle_err("dynetat0, phisinit", ierr, nid)
138    
139      call NF95_INQ_VARID (nid, "temps", nvarid)      call NF95_INQ_VARID (nid, "temps", nvarid)
Legend:
Removed from v.3  
changed lines
  Added in v.10
  ViewVC Help
Powered by ViewVC 1.1.21