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

Diff of /trunk/dyn3d/writehist.f

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

revision 261 by guez, Wed Mar 7 13:33:15 2018 UTC revision 262 by guez, Wed Mar 7 13:46:18 2018 UTC
# Line 43  contains Line 43  contains
43    
44      call assert([size(vcov, 1), size(ucov, 1), size(teta, 1), size(phi, 1), &      call assert([size(vcov, 1), size(ucov, 1), size(teta, 1), size(phi, 1), &
45           size(pk, 1), size(ps, 1), size(masse, 1)] == size(q, 1), &           size(pk, 1), size(ps, 1), size(masse, 1)] == size(q, 1), &
46           "writedynav iim")           "writehist iim")
47      call assert([size(vcov, 2) + 1, size(ucov, 2), size(teta, 2), &      call assert([size(vcov, 2) + 1, size(ucov, 2), size(teta, 2), &
48           size(phi, 2), size(pk, 2), size(ps, 2), size(masse, 2)] &           size(phi, 2), size(pk, 2), size(ps, 2), size(masse, 2)] &
49           == size(q, 2), "writedynav jjm")           == size(q, 2), "writehist jjm")
50      call assert([size(vcov, 3), size(ucov, 3), size(teta, 3), size(phi, 3), &      call assert([size(vcov, 3), size(ucov, 3), size(teta, 3), size(phi, 3), &
51           size(pk, 3), size(masse, 3), size(q, 3)] == llm, "writedynav llm")           size(pk, 3), size(masse, 3), size(q, 3)] == llm, "writehist llm")
52    
53      itau_w = itau_dyn + time      itau_w = itau_dyn + time
54      call covnat(llm, ucov, vcov, unat, vnat)      call covnat(llm, ucov, vcov, unat, vnat)

Legend:
Removed from v.261  
changed lines
  Added in v.262

  ViewVC Help
Powered by ViewVC 1.1.21