Changeset 1793 for LMDZ5/trunk/libf/dyn3d/leapfrog.F

Timestamp:

Jul 18, 2013, 9:13:18 AM (11 years ago)

Author:

Ehouarn Millour

Message:

Improved sponge layer:

• Sponge tendencies are now computed analytically, instead than using a Forward Euler approximation.
• Sponge tendencies are added within top_bound, and the sponge is applied after physics tendencies have been taken into account.

These changes imply that GCM results (when using sponge layer) will be differentwrt bench test cases using previous revisions.
EM

File:

• LMDZ5/trunk/libf/dyn3d/leapfrog.F (modified) (3 diffs)

LMDZ5/trunk/libf/dyn3d/leapfrog.F

@@ -436 +436 @@
      $               clesphy0, dufi,dvfi,dtetafi,dqfi,dpfi  )
 
-         IF (ok_strato) THEN
-           CALL top_bound( vcov,ucov,teta,masse,dufi,dvfi,dtetafi)
-         ENDIF
-       
 c      ajout des tendances physiques:
 c      ------------------------------
@@ -445 +441 @@
      $                  ucov, vcov, teta , q   ,ps ,
      $                 dufi, dvfi, dtetafi , dqfi ,dpfi  )
+
+         IF (ok_strato) THEN
+           CALL top_bound( vcov,ucov,teta,masse,dtphys)
+         ENDIF
+       
 c
 c  Diagnostique de conservation de l'énergie : difference
@@ -476 +477 @@
         ! Sponge layer (if any)
         IF (ok_strato) THEN
-          dufi(:,:)=0.
-          dvfi(:,:)=0.
-          dtetafi(:,:)=0.
-          dqfi(:,:,:)=0.
-          dpfi(:)=0.
-          CALL top_bound(vcov,ucov,teta,masse,dufi,dvfi,dtetafi)
-          CALL addfi( dtvr, leapf, forward   ,
-     $                  ucov, vcov, teta , q   ,ps ,
-     $                 dufi, dvfi, dtetafi , dqfi ,dpfi  )
+!          dufi(:,:)=0.
+!          dvfi(:,:)=0.
+!          dtetafi(:,:)=0.
+!          dqfi(:,:,:)=0.
+!          dpfi(:)=0.
+!          CALL top_bound(vcov,ucov,teta,masse,dufi,dvfi,dtetafi)
+           CALL top_bound( vcov,ucov,teta,masse,dtvr)
+!          CALL addfi( dtvr, leapf, forward   ,
+!     $                  ucov, vcov, teta , q   ,ps ,
+!     $                 dufi, dvfi, dtetafi , dqfi ,dpfi  )
         ENDIF ! of IF (ok_strato) 
       ENDIF ! of IF (iflag_phys.EQ.2)