Changeset 270 for trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F

Timestamp:

Aug 18, 2011, 12:09:27 PM (13 years ago)

Author:

emillour

Message:

Ehouarn: Mise a jour des dynamiques (seq et ) pour suivre
les changements et developpements de LMDZ5 terrestre
(mise a niveau avec LMDZ5 trunk, rev 1560). Ce qui ne devrais pas changer grand chose au fonctionnement de base du GCM).
Modification importante: correction du bug sur le cumul des flux de masse pour le transport des traceurs (mais dans les faits semble avoir peu d'impact).
(voir "commit_importants.log" pour les details des ajouts et modifications).

File:

• trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F (modified) (5 diffs)

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

@@ -152 +152 @@
       logical ok_sync
       parameter (ok_sync = .true.) 
+      logical physic
 
       data callinigrads/.true./
@@ -222 +223 @@
 
       itau = 0
+      physic=.true.
+      if (iflag_phys==0.or.iflag_phys==2) physic=.false.
+
 c      iday = day_ini+itau/day_step
 c      time = REAL(itau-(iday-day_ini)*day_step)/day_step+time_0
@@ -322 +326 @@
       ! Purely Matsuno time stepping
          IF( MOD(itau,iconser) .EQ.0.AND.  forward    ) conser = .TRUE.
-         IF( MOD(itau,idissip ).EQ.0.AND..NOT.forward ) apdiss = .TRUE.
+         IF( MOD(itau,dissip_period ).EQ.0.AND..NOT.forward ) 
+     s        apdiss = .TRUE.
          IF( MOD(itau,iphysiq ).EQ.0.AND..NOT.forward 
-     s          .and. iflag_phys.EQ.1                 ) apphys = .TRUE.
+     s          .and. physic                        ) apphys = .TRUE.
       ELSE
       ! Leapfrog/Matsuno time stepping 
          IF( MOD(itau   ,iconser) .EQ. 0              ) conser = .TRUE.
-         IF( MOD(itau+1,idissip)  .EQ. 0              ) apdiss = .TRUE.
-         IF( MOD(itau+1,iphysiq).EQ.0.AND.iflag_phys.EQ.1) apphys=.TRUE.
+         IF( MOD(itau+1,dissip_period).EQ.0 .AND. .NOT. forward )
+     s        apdiss = .TRUE.
+         IF( MOD(itau+1,iphysiq).EQ.0.AND.physic       ) apphys=.TRUE.
       END IF
 
@@ -357 +363 @@
 c   -------------------------------------------------------------
 
-      IF( forward. OR . leapf )  THEN
-
+!      IF( forward. OR . leapf )  THEN
+      IF((.not.forward).OR. leapf )  THEN
+        ! Ehouarn: gather mass fluxes during backward Matsuno or LF step
          CALL caladvtrac(q,pbaru,pbarv,
      *        p, masse, dq,  teta,
@@ -497 +504 @@
           ENDDO
         ENDDO ! of DO l=1,llm 
-          call friction(ucov,vcov,dtvr)
     
-       if (planet_type.eq."giant") then
-          ! Intrinsic heat flux
-          ! Aymeric -- for giant planets
-          if (ihf .gt. 1.e-6) then
-          !print *, '**** INTRINSIC HEAT FLUX ****', ihf 
-            DO ij=1,ip1jmp1
-              teta(ij,1) = teta(ij,1)
-     &        + dtvr * aire(ij) * ihf / cpp / masse(ij,1)
-            ENDDO
-          !print *, '**** d teta ' 
-          !print *, dtvr * aire(ijb:ije) * ihf / cpp / masse(ijb:ije,1)
-          endif
-       endif
+        if (planet_type.eq."giant") then
+          ! add an intrinsic heat flux at the base of the atmosphere
+          teta(:,1)=teta(:,1)+dtvr*aire(:)*ihf/cpp/masse(:,1)
+        endif
+
+        call friction(ucov,vcov,dtvr)