Ignore:
Timestamp:
Mar 21, 2002, 4:09:32 PM (23 years ago)
Author:
lmdzadmin
Message:

2 changements pour les fichiers histoire:

  • utilisation de l'entree "rectilineaire" de IOIPSL pour ne plus avoir

a

lancer ncregular a chaque fois

  • le calendrier des fichiers histoire est maintenant base sur la date d'initialisation de la simulation plutot que sur la date de depart du

job

en cours

LF

Location:
LMDZ.3.3/branches/rel-LF/libf/bibio
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/branches/rel-LF/libf/bibio/initdynav.F

    r350 r353  
     1c
     2c $Header$
     3c
    14      subroutine initdynav(infile,day0,anne0,tstep,t_ops,t_wrt
    25     .                     ,nq,fileid)

  • LMDZ.3.3/branches/rel-LF/libf/bibio/initfluxsto.F

    r177 r353  
    6666      integer uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid
    6767      integer ii,jj
    68       integer zan
     68      integer zan, idayref
    6969      logical ok_sync
    7070C
    … …  
    7979C         
    8080
    81       zan = anne_ini
    82       CALL ymds2ju(zan, 1, 1, 0.0, zjulian)
    83       zjulian = zjulian + day_ini
    84       tau0 = 0
     81      zan = annee_ref
     82      idayref = day_ref
     83      CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
     84      tau0 = itau_dyn
    8585       
    8686        do jj = 1, jjp1
    … …  
    9191      enddo
    9292 
    93       call histbeg(infile, iip1, rlong, jjp1, rlat,
     93      call histbeg(infile, iip1, rlong(:,1), jjp1, rlat(1,:),
    9494     .             1, iip1, 1, jjp1,
    9595     .             tau0, zjulian, tstep, uhoriid, fileid)
    … …  
    107107      enddo
    108108
    109       call histbeg('fluxstokev.nc', iip1, rlong, jjm, rlat,
     109      call histbeg('fluxstokev.nc', iip1, rlong(:,1), jjm, rlat(1,:),
    110110     .             1, iip1, 1, jjm,
    111111     .             tau0, zjulian, tstep, vhoriid, filevid)

  • LMDZ.3.3/branches/rel-LF/libf/bibio/inithist.F

    r177 r353  
    6464      integer uhoriid, vhoriid, thoriid, zvertiid
    6565      integer ii,jj
    66       integer zan
     66      integer zan, dayref
    6767C
    6868C  Initialisations
    … …  
    7676
    7777      zan = anne0
    78       CALL ymds2ju(zan, 1, 1, 0.0, zjulian)
    79       zjulian = zjulian + day0
    80       tau0 = 0
     78      dayref = day0
     79      CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
     80      tau0 = itau_dyn
    8181     
    8282      do jj = 1, jjp1
    … …  
    8787      enddo
    8888       
    89       call histbeg(infile, iip1, rlong, jjp1, rlat,
     89      call histbeg(infile, iip1, rlong(:,1), jjp1, rlat(1,:),
    9090     .             1, iip1, 1, jjp1,
    9191     .             tau0, zjulian, tstep, uhoriid, fileid)
    … …  
    102102      enddo
    103103
    104       call histbeg('dyn_histv.nc', iip1, rlong, jjm, rlat,
     104      call histbeg('dyn_histv.nc', iip1, rlong(:,1), jjm, rlat(1,:),
    105105     .             1, iip1, 1, jjm,
    106106     .             tau0, zjulian, tstep, vhoriid, filevid)

  • LMDZ.3.3/branches/rel-LF/libf/bibio/writedynav.F

    r295 r353  
    6464      REAL vnat(ip1jm,llm),unat(ip1jmp1,llm)
    6565      logical ok_sync
     66      integer itau_w
    6667C
    6768C  Initialisations
    … …  
    7677      vnat = 999.999
    7778      unat = 999.999
     79      itau_w = itau_dyn + time
    7880
    7981C Passage aux composantes naturelles du vent
    … …  
    8688C
    8789      call gr_u_scal(llm, unat, us)
    88       call histwrite(histid, 'u', time, us,
     90      call histwrite(histid, 'u', itau_w, us,
    8991     .               iip1*jjp1*llm, ndex3d)
    9092C
    … …  
    9294C
    9395      call gr_v_scal(llm, vnat, vs)
    94       call histwrite(histid, 'v', time, vs,
     96      call histwrite(histid, 'v', itau_w, vs,
    9597     .               iip1*jjp1*llm, ndex3d)
    9698C
    9799C  Temperature potentielle moyennee
    98100C
    99       call histwrite(histid, 'theta', time, teta,
     101      call histwrite(histid, 'theta', itau_w, teta,
    100102     .                iip1*jjp1*llm, ndex3d)
    101103C
    … …  
    105107        tm(ii) = teta(ii) * ppk(ii)/cpp
    106108      enddo
    107       call histwrite(histid, 'temp', time, tm,
     109      call histwrite(histid, 'temp', itau_w, tm,
    108110     .                iip1*jjp1*llm, ndex3d)
    109111C
    110112C  Geopotentiel
    111113C
    112       call histwrite(histid, 'phi', time, phi,
     114      call histwrite(histid, 'phi', itau_w, phi,
    113115     .                iip1*jjp1*llm, ndex3d)
    114116C
    … …  
    122124            WRITE(str(2:3),'(i2.2)') iq
    123125          ENDIF
    124           call histwrite(histid, str, time, q(:,:,iq),
     126          call histwrite(histid, str, itau_w, q(:,:,iq),
    125127     .                   iip1*jjp1*llm, ndex3d)
    126128        enddo
    … …  
    129131C  Masse
    130132C
    131        call histwrite(histid, 'masse', time, masse, iip1*jjp1, ndex2d)
     133       call histwrite(histid, 'masse', itau_w, masse, iip1*jjp1, ndex2d)
    132134C
    133135C  Pression au sol
    134136C
    135        call histwrite(histid, 'ps', time, ps, iip1*jjp1, ndex2d)
     137       call histwrite(histid, 'ps', itau_w, ps, iip1*jjp1, ndex2d)
    136138C
    137139C  Geopotentiel au sol
    138140C
    139        call histwrite(histid, 'phis', time, phis, iip1*jjp1, ndex2d)
     141       call histwrite(histid, 'phis', itau_w, phis, iip1*jjp1, ndex2d)
    140142C
    141143C  Fin

  • LMDZ.3.3/branches/rel-LF/libf/bibio/writehist.F

    r30 r353  
    6363      character*3 str
    6464      logical ok_sync
     65      integer itau_w
    6566C
    6667C  Initialisations
    … …  
    7172      ndex2d = 0
    7273      ok_sync =.TRUE.
     74      itau_w = itau_dyn + time
    7375C
    7476C  Appels a histwrite pour l'ecriture des variables a sauvegarder
    … …  
    7678C  Vents U
    7779C
    78       call histwrite(histid, 'ucov', time, ucov,
     80      call histwrite(histid, 'ucov', itau_w, ucov,
    7981     .               iip1*jjp1*llm, ndexu)
    8082
    … …  
    8284C  Vents V
    8385C
    84       call histwrite(histvid, 'vcov', time, vcov,
     86      call histwrite(histvid, 'vcov', itau_w, vcov,
    8587     .               iip1*jjm*llm, ndexv)
    8688
    … …  
    8890C  Temperature potentielle
    8991C
    90       call histwrite(histid, 'teta', time, teta,
     92      call histwrite(histid, 'teta', itau_w, teta,
    9193     .                iip1*jjp1*llm, ndexu)
    9294C
    9395C  Geopotentiel
    9496C
    95       call histwrite(histid, 'phi', time, phi,
     97      call histwrite(histid, 'phi', itau_w, phi,
    9698     .                iip1*jjp1*llm, ndexu)
    9799C
    … …  
    105107            WRITE(str(2:3),'(i2.2)') iq
    106108          ENDIF
    107           call histwrite(histid, str, time, q(:,:,iq),
     109          call histwrite(histid, str, itau_w, q(:,:,iq),
    108110     .                   iip1*jjp1*llm, ndexu)
    109111        enddo
    … …  
    112114C  Masse
    113115C
    114       call histwrite(histid, 'masse', time, masse, iip1*jjp1, ndex2d)
     116      call histwrite(histid, 'masse', itau_w, masse, iip1*jjp1, ndex2d)
    115117C
    116118C  Pression au sol
    117119C
    118       call histwrite(histid, 'ps', time, ps, iip1*jjp1, ndex2d)
     120      call histwrite(histid, 'ps', itau_w, ps, iip1*jjp1, ndex2d)
    119121C
    120122C  Geopotentiel au sol
    121123C
    122       call histwrite(histid, 'phis', time, phis, iip1*jjp1, ndex2d)
     124      call histwrite(histid, 'phis', itau_w, phis, iip1*jjp1, ndex2d)
    123125C
    124126C  Fin

  • LMDZ.3.3/branches/rel-LF/libf/bibio/writephys.F90

    r284 r353  
    9393  real, dimension(llm) :: zlev
    9494  real                 :: date0, dtime
     95  integer              :: itau_deb
    9596!
    9697! Variables locales
    … …  
    121122! Initialisation du fichier
    122123!
    123   call histbeg(nom_fichier, iim, temp_lon, jjm, temp_lat, &
    124      &                 1, iim, 1, jjm, 0, date0, dtime, &
     124  call histbeg(nom_fichier, iim, temp_lon(:,1), jjm, temp_lat(1,:), &
     125     &                 1, iim, 1, jjm, itau_deb, date0, dtime, &
    125126     &                 nhori, file_id)
    126127  call histvert(file_id, "presnivs", "Vertical levels", "mb", &
Note: See TracChangeset for help on using the changeset viewer.