Ignore:
Timestamp:
Mar 22, 2016, 1:19:13 PM (9 years ago)
Author:
Ehouarn Millour
Message:

Reinstate writehist, writedyn and bilan_dyn in dyn3dmem so that ouputs in the dynamics (dyn_hist* and dynzon files) may be generated when in parallel.
EM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/dyn3dmem/initdynav_loc.F

    r1907 r2475  
    154154
    155155      ddid=(/ 1,2 /)
    156       dsg=(/ iip1,jjp1 /)
     156      dsg=(/ iip1,jjm /)
    157157      dsl=(/ iip1,jjn /)
    158158      dpf=(/ 1,jjb /)
     
    171171     
    172172! Grille U
     173
     174      do jj = 1, jjp1
     175        do ii = 1, iip1
     176          rlong(ii,jj) = rlonu(ii) * 180. / pi
     177          rlat(ii,jj) = rlatu(jj) * 180. / pi
     178        enddo
     179      enddo
    173180
    174181      jjb=jj_begin
     
    209216C  Vents U
    210217C
     218      jjn=jj_nb
    211219      call histdef(histuaveid, 'u', 'vent u moyen ',
    212      .             'm/s', iip1, jjp1, uhoriid, llm, 1, llm, zvertiidu,
     220     .             'm/s', iip1, jjn, uhoriid, llm, 1, llm, zvertiidu,
    213221     .             32, 'ave(X)', t_ops, t_wrt)
    214222
     
    216224C  Vents V
    217225C
     226      if (pole_sud) jjn=jj_nb-1
    218227      call histdef(histvaveid, 'v', 'vent v moyen',
    219      .             'm/s', iip1, jjm, vhoriid, llm, 1, llm, zvertiidv,
     228     .             'm/s', iip1, jjn, vhoriid, llm, 1, llm, zvertiidv,
    220229     .             32, 'ave(X)', t_ops, t_wrt)
    221230
     
    223232C  Temperature
    224233C
     234      jjn=jj_nb
    225235      call histdef(histaveid, 'temp', 'temperature moyenne', 'K',
    226      .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
     236     .             iip1, jjn, thoriid, llm, 1, llm, zvertiid,
    227237     .             32, 'ave(X)', t_ops, t_wrt)
    228238C
     
    230240C
    231241      call histdef(histaveid, 'theta', 'temperature potentielle', 'K',
    232      .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
     242     .             iip1, jjn, thoriid, llm, 1, llm, zvertiid,
    233243     .             32, 'ave(X)', t_ops, t_wrt)
    234244
     
    238248C
    239249      call histdef(histaveid, 'phi', 'geopotentiel moyen', '-',
    240      .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
     250     .             iip1, jjn, thoriid, llm, 1, llm, zvertiid,
    241251     .             32, 'ave(X)', t_ops, t_wrt)
    242252C
     
    251261C  Masse
    252262C
    253       call histdef(histaveid, 'masse', 'masse', 'kg',
    254      .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
     263      call histdef(histaveid, 'masse', 'masse moyenne', 'kg',
     264     .             iip1, jjn, thoriid, llm, 1, llm, zvertiid,
    255265     .             32, 'ave(X)', t_ops, t_wrt)
    256266C
     
    258268C
    259269      call histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa',
    260      .             iip1, jjp1, thoriid, 1, 1, 1, -99,
    261      .             32, 'ave(X)', t_ops, t_wrt)
    262 C
    263 Pression au sol
     270     .             iip1, jjn, thoriid, 1, 1, 1, -99,
     271     .             32, 'ave(X)', t_ops, t_wrt)
     272C
     273Geopotentiel au sol
    264274C
    265275!      call histdef(histaveid, 'phis', 'geopotentiel au sol', '-',
     
    273283      call histend(histvaveid)
    274284#else
    275       write(lunout,*)'initdynav_p: Needs IOIPSL to function'
     285      write(lunout,*)'initdynav_loc: Needs IOIPSL to function'
    276286#endif
    277287! #endif of #ifdef CPP_IOIPSL
    278       return
    279288      end
Note: See TracChangeset for help on using the changeset viewer.