/[lmdze]/trunk/libf/dyn3d/dynetat0.f90
ViewVC logotype

Diff of /trunk/libf/dyn3d/dynetat0.f90

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

revision 6 by guez, Wed Feb 27 13:16:39 2008 UTC revision 7 by guez, Mon Mar 31 12:24:17 2008 UTC
# Line 29  contains Line 29  contains
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
# Line 44  contains Line 44  contains
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")
# Line 132  contains Line 132  contains
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)
Legend:
Removed from v.6  
changed lines
  Added in v.7
  ViewVC Help
Powered by ViewVC 1.1.21