Changeset 495 for trunk/LMDZ.COMMON/libf


Ignore:
Timestamp:
Jan 10, 2012, 5:59:10 PM (13 years ago)
Author:
slebonnois
Message:

Mise a jour physique Titan, ajout des forces de marees (dans la dynamique, sous flag titan). SL.

Location:
trunk/LMDZ.COMMON/libf
Files:
2 added
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F

    r492 r495  
    338338       CALL getin('mode_top_bound',mode_top_bound)
    339339       CALL getin('tau_top_bound',tau_top_bound)
     340
     341! FOR TITAN: tidal forces
     342       tidal=.TRUE.
     343       CALL getin('tidal',tidal)
    340344
    341345!Config  Key  = coefdis

  • trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F

    r492 r495  
    106106      REAL dtetatop(ip1jmp1,llm)
    107107      REAL dqtop(ip1jmp1,llm,nqtot),dptop(ip1jmp1)
     108
     109c   TITAN : tendances due au forces de marees */s
     110      REAL dvtidal(ip1jm,llm),dutidal(ip1jmp1,llm)
    108111
    109112c   tendances physiques */s
    … …  
    396399     $              finvmaold                                    )
    397400
     401       IF ((planet_type.eq.titan).and.(tidal)) then
     402c-----------------------------------------------------------------------
     403c   Marées gravitationnelles causées par Saturne
     404c   B. Charnay (28/10/2010)
     405c   ----------------------------------------------------------
     406            CALL tidal_forces(rdaym_ini, dutidal, dvtidal)
     407            ucov=ucov+dutidal*dt
     408            vcov=vcov+dvtidal*dt
     409       ENDIF
    398410
    399411c .P.Le Van (26/04/94  ajout de  finvpold dans l'appel d'integrd)

  • trunk/LMDZ.COMMON/libf/dyn3d/logic.h

    r127 r495  
    1010      COMMON/logicl/ purmats,forward,leapf,apphys,                      &
    1111     &  statcl,conser,apdiss,apdelq,saison,ecripar,fxyhypb,ysinus       &
    12      &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
     12     &  ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile               &
    1313     &  ,ok_limit,ok_etat0,grilles_gcm_netcdf,hybrid
    1414
    … …  
    1717      LOGICAL purmats,forward,leapf,apphys,statcl,conser,               &
    1818     & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus                      &
    19      &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
     19     &  ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile               &
    2020     &  ,ok_limit,ok_etat0,grilles_gcm_netcdf
    2121      logical hybrid ! vertical coordinate is hybrid if true (sigma otherwise)

  • trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F

    r492 r495  
    352352       CALL getin('tau_top_bound',tau_top_bound)
    353353
    354 !
     354! FOR TITAN: tidal forces
     355       tidal=.TRUE.
     356       CALL getin('tidal',tidal)
     357
    355358!Config  Key  = coefdis
    356359!Config  Desc = coefficient pour gamdissip

  • trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F

    r492 r495  
    117117      REAL,SAVE :: dptop(ip1jmp1)
    118118      REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: dqtop
     119
     120c   TITAN : tendances due au forces de marees */s
     121      REAL,SAVE :: dvtidal(ip1jm,llm),dutidal(ip1jmp1,llm)
    119122
    120123c   variables pour le fichier histoire
    … …  
    696699     $              finvmaold                                    )
    697700
     701       IF ((planet_type.eq.titan).and.(tidal)) then
     702c-----------------------------------------------------------------------
     703c   Marées gravitationnelles causées par Saturne
     704c   B. Charnay (28/10/2010)
     705c   ----------------------------------------------------------
     706            CALL tidal_forces(rdaym_ini, dutidal, dvtidal)
     707            ucov=ucov+dutidal*dt
     708            vcov=vcov+dvtidal*dt
     709       ENDIF
     710
    698711!       CALL FTRACE_REGION_END("integrd")
    699712c$OMP BARRIER

  • trunk/LMDZ.COMMON/libf/dyn3dpar/logic.h

    r127 r495  
    1010      COMMON/logicl/ purmats,forward,leapf,apphys,                      &
    1111     &  statcl,conser,apdiss,apdelq,saison,ecripar,fxyhypb,ysinus       &
    12      &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
     12     &  ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile               &
    1313     &  ,ok_limit,ok_etat0,grilles_gcm_netcdf,hybrid
    1414
    … …  
    1717      LOGICAL purmats,forward,leapf,apphys,statcl,conser,               &
    1818     & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus                      &
    19      &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
     19     &  ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile               &
    2020     &  ,ok_limit,ok_etat0,grilles_gcm_netcdf
    2121      logical hybrid ! vertical coordinate is hybrid if true (sigma otherwise)
Note: See TracChangeset for help on using the changeset viewer.