Changeset 541 for LMDZ4/trunk/libf/dyn3d


Ignore:
Timestamp:
Jun 22, 2004, 1:45:36 PM (20 years ago)
Author:
lmdzadmin
Message:

Convergence avec la version d'Olivia Coindreau incluant:

  • le offline
  • les thermiques
  • mellor & yamada dans la couche limite

LF

Location:
LMDZ4/trunk/libf/dyn3d
Files:
1 added
8 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk/libf/dyn3d/bilan_dyn.F

    r524 r541  
    22! $Header$
    33!
    4 c
    5 c $Header$
    6 c
    74      SUBROUTINE bilan_dyn (ntrac,dt_app,dt_cum,
    85     s  ps,masse,pk,flux_u,flux_v,teta,phi,ucov,vcov,trac)

  • LMDZ4/trunk/libf/dyn3d/conf_gcm.F

    r527 r541  
    574574      ENDIF
    575575c
     576!Config  Key  = offline
     577!Config  Desc = Nouvelle eau liquide
     578!Config  Def  = n
     579!Config  Help = Permet de mettre en route la
     580!Config         nouvelle parametrisation de l'eau liquide !
     581       offline = .FALSE.
     582       CALL getin('offline',offline)
     583
    576584
    577585      write(lunout,*)' #########################################'
    … …  
    615623      write(lunout,*)' tauxx = ', tauxx
    616624      write(lunout,*)' tauyy = ', tauyy
     625      write(lunout,*)' offline = ', offline
    617626
    618627      RETURN
    … …  
    720729       CALL getin('ysinus',ysinus)
    721730c
     731!Config  Key  = offline
     732!Config  Desc = Nouvelle eau liquide
     733!Config  Def  = n
     734!Config  Help = Permet de mettre en route la
     735!Config         nouvelle parametrisation de l'eau liquide !
     736       offline = .FALSE.
     737       CALL getin('offline',offline)
     738      write(lunout,*)' offline = ', offline
     739
    722740
    723741      write(lunout,*)' #########################################'
    … …  
    761779      write(lunout,*)' taux = ', taux
    762780      write(lunout,*)' tauy = ', tauy
     781      write(lunout,*)' offline = ', offline
    763782c
    764783      RETURN

  • LMDZ4/trunk/libf/dyn3d/control.h

    r524 r541  
    66
    77      COMMON/control/nday,day_step,
    8      $              iperiod,iapp_tracvl,iconser,iecri,idissip,iphysiq ,
    9      $              periodav,ecritphy,iecrimoy,dayref,anneeref,
    10      $              raz_date
     8     .              iperiod,iapp_tracvl,iconser,iecri,idissip,iphysiq ,
     9     .              periodav,ecritphy,iecrimoy,dayref,anneeref,
     10     .              raz_date,offline
    1111
    1212      INTEGER   nday,day_step,iperiod,iapp_tracvl,iconser,iecri,
    13      $          idissip,iphysiq,iecrimoy,dayref,anneeref, raz_date
     13     .          idissip,iphysiq,iecrimoy,dayref,anneeref, raz_date
    1414      REAL periodav, ecritphy
     15      logical offline
    1516
    1617c-----------------------------------------------------------------------

  • LMDZ4/trunk/libf/dyn3d/dynetat0.F

    r524 r541  
    106106      IF ( tab_cntrl(24).EQ.1. )  THEN
    107107        fxyhypb  = . TRUE .
    108         dzoomx   = tab_cntrl(25)
    109         dzoomy   = tab_cntrl(26)
    110         taux     = tab_cntrl(28)
    111         tauy     = tab_cntrl(29)
     108c        dzoomx   = tab_cntrl(25)
     109c        dzoomy   = tab_cntrl(26)
     110c        taux     = tab_cntrl(28)
     111c        tauy     = tab_cntrl(29)
    112112      ELSE
    113113        fxyhypb = . FALSE .

  • LMDZ4/trunk/libf/dyn3d/dynredem.F

    r524 r541  
    394394#include "comgeom.h"
    395395#include "advtrac.h"
     396#include "temps.h"
    396397
    397398      INTEGER nq, l
    … …  
    406407      INTEGER ierr
    407408      INTEGER iq
     409      INTEGER length
     410      PARAMETER (length = 100)
     411      REAL tab_cntrl(length) ! tableau des parametres du run
    408412      character*20 modname
    409413      character*80 abort_message
    … …  
    435439#endif
    436440      PRINT*, "Enregistrement pour ", nb, time
     441
     442c
     443c  Re-ecriture du tableau de controle, itaufin n'est plus defini quand
     444c  on passe dans dynredem0
     445      ierr = NF_INQ_VARID (nid, "controle", nvarid)
     446      IF (ierr .NE. NF_NOERR) THEN
     447         abort_message="dynredem1: Le champ <controle> est absent"
     448         ierr = 1
     449         CALL abort_gcm(modname,abort_message,ierr)
     450      ENDIF
     451#ifdef NC_DOUBLE
     452      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tab_cntrl)
     453#else
     454      ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)
     455#endif
     456       tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
     457#ifdef NC_DOUBLE
     458      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)
     459#else
     460      ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl)
     461#endif
    437462
    438463c  Ecriture des champs

  • LMDZ4/trunk/libf/dyn3d/gcm.F

    r524 r541  
    5656#include "com_io_dyn.h"
    5757#include "iniprint.h"
    58 
    59 c#include "tracstoke.h"
     58#include "tracstoke.h"
    6059
    6160
    … …  
    118117c-jld
    119118
    120 
    121       LOGICAL offline  ! Controle du stockage ds "fluxmass"
    122       PARAMETER (offline=.false.)
    123119
    124120      character*80 dynhist_file, dynhistave_file
    … …  
    368364#endif
    369365
     366c  Choix des frequences de stokage pour le offline
     367c      istdyn=day_step/4     ! stockage toutes les 6h=1jour/4
     368c      istdyn=day_step/12     ! stockage toutes les 2h=1jour/12
     369      istdyn=day_step/4     ! stockage toutes les 6h=1jour/12
     370      istphy=istdyn/iphysiq     
     371
     372
    370373c
    371374c-----------------------------------------------------------------------

  • LMDZ4/trunk/libf/dyn3d/leapfrog.F

    r524 r541  
    149149c-jld
    150150
    151       LOGICAL offline  ! Controle du stockage ds "fluxmass"
    152       PARAMETER (offline=.false.)
    153 
    154151      character*80 dynhist_file, dynhistave_file
    155152      character*20 modname
    … …  
    199196
    200197#ifdef CPP_IOIPSL
    201       if (ok_guide) then
     198      if (ok_guide.and.(itaufin-itau-1)*dtvr.gt.21600) then
    202199        call guide(itau,ucov,vcov,teta,q,masse,ps)
     200      else
     201        print*,'attention on ne guide pas les 6 dernieres heures'
    203202      endif
    204203#endif
    … …  
    289288
    290289#ifdef CPP_IOIPSL
    291 c           CALL fluxstokenc(pbaru,pbarv,masse,teta,phi,phis,
    292 c    .   time_step, itau)
     290           CALL fluxstokenc(pbaru,pbarv,masse,teta,phi,phis,
     291     .   dtvr, itau)
    293292#endif
    294293
    … …  
    486485
    487486            IF( itau. EQ. itaufinp1 ) then 
    488 c      write(79,*) 'ucov',ucov
    489 c      write(80,*) 'vcov',vcov
    490 c      write(81,*) 'teta',teta
    491 c      write(82,*) 'ps',ps
    492 c      write(83,*) 'q',q
    493 c      WRITE(85,*) 'q1 = ',q(:,:,1)
    494 c      WRITE(86,*) 'q3 = ',q(:,:,3)
     487       write(79,*) 'ucov',ucov
     488       write(80,*) 'vcov',vcov
     489       write(81,*) 'teta',teta
     490       write(82,*) 'ps',ps
     491       write(83,*) 'q',q
     492       WRITE(85,*) 'q1 = ',q(:,:,1)
     493       WRITE(86,*) 'q3 = ',q(:,:,3)
    495494
    496495              abort_message = 'Simulation finished'

  • LMDZ4/trunk/libf/dyn3d/wrgrads.F

    r524 r541  
    2929      writectl=.false.
    3030
    31 c     print*,if,iid(if),jid(if),ifd(if),jfd(if)
     31      print*,if,iid(if),jid(if),ifd(if),jfd(if)
    3232      iii=iid(if)
    3333      iji=jid(if)
    … …  
    3838      lm=lmd(if)
    3939
    40 c      print*,'im,jm,lm,name,firsttime(if)'
    41 c      print*,im,jm,lm,name,firsttime(if)
     40      print*,'im,jm,lm,name,firsttime(if)'
     41      print*,im,jm,lm,name,firsttime(if)
    4242
    4343      if(firsttime(if)) then
    … …  
    8181      endif
    8282
    83 c     print*,'ivar(if),nvar(if),var(ivar(if),if),writectl'
    84 c     print*,ivar(if),nvar(if),var(ivar(if),if),writectl
     83      print*,'ivar(if),nvar(if),var(ivar(if),if),writectl'
     84      print*,ivar(if),nvar(if),var(ivar(if),if),writectl
    8585      do l=1,nl
    8686         irec(if)=irec(if)+1
Note: See TracChangeset for help on using the changeset viewer.