Changeset 184 for LMDZ.3.3/trunk


Ignore:
Timestamp:
Mar 29, 2001, 10:41:34 AM (24 years ago)
Author:
lmdzadmin
Message:

Modifs A.Idelkadi pour le offline
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/trunk/libf/dyn3d/fluxstokenc.F

    r55 r184  
    11      SUBROUTINE fluxstokenc(pbaru,pbarv,masse,teta,phi,phis,
    2      . time_step,itau, fluxid, fluxvid,fluxdid )
     2     . time_step,itau )
    33
    44       USE IOIPSL
     
    3030
    3131      REAL pbarvst(iip1,jjp1,llm),zistdyn
     32        real dtcum
    3233
    33 
    34       INTEGER iadvtr
    35       integer ndex1d(1)
    36       integer ndex2d(ip1jmp1)
    37       integer ndex3dv(ip1jm*llm),ndex3d(ip1jmp1*llm)
     34      INTEGER iadvtr,ndex(1)
    3835      integer nscal
    3936      real tst(1),ist(1),istp(1)
     
    5350     .  time_step,istdyn* time_step,istdyn* time_step,
    5451     . nqmx, fluxid,fluxvid,fluxdid)
     52       
     53        ndex(1) = 0
     54        call histwrite(fluxid, 'phis', 1, phis, iip1*jjp1, ndex)
     55        call histwrite(fluxid, 'aire', 1, aire, iip1*jjp1, ndex)
     56       
     57        ndex(1) = 0
     58        nscal = 1
     59        tst(1) = time_step
     60        call histwrite(fluxdid, 'dtvr', 1, tst, nscal, ndex)
     61        ist(1)=istdyn
     62        call histwrite(fluxdid, 'istdyn', 1, ist, nscal, ndex)
     63        istp(1)= istphy
     64        call histwrite(fluxdid, 'istphy', 1, istp, nscal, ndex)
     65       
    5566        first = .false.
    5667
    5768      endif
    58 
    59 
    60       ndex1d = 0
    61       ndex2d = 0
    62       ndex3dv = 0
    63       ndex3d = 0
    6469
    6570
     
    9398c   Test pour savoir si on advecte a ce pas de temps
    9499      IF ( iadvtr.EQ.istdyn ) THEN
    95 
    96100c    normalisation
    97101      DO l=1,llm
     
    124128
    125129         iadvtr=0
    126 
    127 c     write(*,*)'histwrite phis'
    128       call histwrite(fluxid, 'phis', 1, phis, iip1*jjp1, ndex2d)
    129 c     write(*,*)'histwrite aire'
    130       call histwrite(fluxid, 'aire', 1, aire, iip1*jjp1, ndex2d)
     130        Print*,'ITAU auqel on stoke les fluxmasses',itau
    131131       
    132       nscal = 1
    133       tst(1) = time_step
    134 c     write(*,*)'histwrite dtvr'
    135       call histwrite(fluxdid, 'dtvr', 1, tst, nscal, ndex1d)
    136       ist(1)=istdyn
    137 c     write(*,*)'histwrite istdyn'
    138       call histwrite(fluxdid, 'istdyn', 1, ist, nscal, ndex1d)
    139       istp(1)= istphy
    140 c     write(*,*)'histwrite istphy'
    141       call histwrite(fluxdid, 'istphy', 1, istp, nscal, ndex1d)
    142 
    143 c       write(*,*)'histwrite masse'     
    144         call histwrite(fluxid, 'masse', itau, masse,
    145      .               iip1*jjp1*llm, ndex3d)
     132        call histwrite(fluxid, 'masse', itau, massem,
     133     .               iip1*jjp1*llm, ndex)
    146134       
    147 c       write(*,*)'histwrite pbaru'     
    148135        call histwrite(fluxid, 'pbaru', itau, pbarug,
    149      .               iip1*jjp1*llm, ndex3d)
     136     .               iip1*jjp1*llm, ndex)
    150137       
    151 c       write(*,*)'histwrite pbarv'     
    152         call histwrite(fluxvid, 'pbarv', itau, pbarvst,
    153      .               iim*jjp1*llm, ndex3dv)
     138        call histwrite(fluxvid, 'pbarv', itau, pbarvg,
     139     .               iip1*jjm*llm, ndex)
    154140       
    155 c       write(*,*)'histwrite w'
    156141        call histwrite(fluxid, 'w' ,itau, wg,
    157      .             iip1*jjp1*llm, ndex3d)
     142     .             iip1*jjp1*llm, ndex)
    158143       
    159 c       write(*,*)'histwrite teta'     
    160144        call histwrite(fluxid, 'teta' ,itau, tetac,
    161      .             iip1*jjp1*llm, ndex3d)
     145     .             iip1*jjp1*llm, ndex)
    162146       
    163 c       write(*,*)'histwrite phi'       
    164147        call histwrite(fluxid, 'phi' ,itau, phic,
    165      .             iip1*jjp1*llm, ndex3d)
     148     .             iip1*jjp1*llm, ndex)
    166149       
    167150C
Note: See TracChangeset for help on using the changeset viewer.