Changeset 495 for trunk/LMDZ.COMMON/libf/dyn3d

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/dyn3d

Files:

• conf_gcm.F (modified) (1 diff)
• leapfrog.F (modified) (2 diffs)
• logic.h (modified) (2 diffs)
• tidal_forces.F (added)

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

@@ -338 +338 @@
        CALL getin('mode_top_bound',mode_top_bound)
        CALL getin('tau_top_bound',tau_top_bound)
+
+! FOR TITAN: tidal forces
+       tidal=.TRUE.
+       CALL getin('tidal',tidal)
 
 !Config  Key  = coefdis

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

@@ -106 +106 @@
       REAL dtetatop(ip1jmp1,llm)
       REAL dqtop(ip1jmp1,llm,nqtot),dptop(ip1jmp1)
+
+c   TITAN : tendances due au forces de marees */s
+      REAL dvtidal(ip1jm,llm),dutidal(ip1jmp1,llm)
 
 c   tendances physiques */s
@@ -396 +399 @@
      $              finvmaold                                    )
 
+       IF ((planet_type.eq.titan).and.(tidal)) then
+c-----------------------------------------------------------------------
+c   Marées gravitationnelles causées par Saturne
+c   B. Charnay (28/10/2010)
+c   ----------------------------------------------------------
+            CALL tidal_forces(rdaym_ini, dutidal, dvtidal)
+            ucov=ucov+dutidal*dt
+            vcov=vcov+dvtidal*dt
+       ENDIF
 
 c .P.Le Van (26/04/94  ajout de  finvpold dans l'appel d'integrd)

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

@@ -10 +10 @@
       COMMON/logicl/ purmats,forward,leapf,apphys,                      &
      &  statcl,conser,apdiss,apdelq,saison,ecripar,fxyhypb,ysinus       &
-     &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
+     &  ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile               &
      &  ,ok_limit,ok_etat0,grilles_gcm_netcdf,hybrid
 
@@ -17 +17 @@
       LOGICAL purmats,forward,leapf,apphys,statcl,conser,               &
      & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus                      &
-     &  ,read_start,ok_guide,ok_strato,ok_gradsfile                     &
+     &  ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile               &
      &  ,ok_limit,ok_etat0,grilles_gcm_netcdf
       logical hybrid ! vertical coordinate is hybrid if true (sigma otherwise)