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

Timestamp:

Jul 24, 2013, 1:36:44 PM (11 years ago)

Author:

emillour

Message:

Common dynamics: Improved sponge layer scheme (top_bound):

• Sponge tendencies are now computed analytically, instead of using a Forward Euler approximation.
• Sponge tendencies are now added within top_bound.

EM

File:

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

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

@@ -103 +103 @@
 
 c   tendances de la couche superieure */s
-      REAL dvtop(ip1jm,llm),dutop(ip1jmp1,llm)
-      REAL dtetatop(ip1jmp1,llm)
-      REAL dqtop(ip1jmp1,llm,nqtot),dptop(ip1jmp1)
+c      REAL dvtop(ip1jm,llm)
+      REAL dutop(ip1jmp1,llm)
+c      REAL dtetatop(ip1jmp1,llm)
+c      REAL dqtop(ip1jmp1,llm,nqtot),dptop(ip1jmp1)
 
 c   TITAN : tendances due au forces de marees */s
@@ -216 +217 @@
         dtetadis(:,:)=0.
         dutop(:,:)   =0.
-        dvtop(:,:)   =0.
-        dtetatop(:,:)=0.
-        dqtop(:,:,:) =0.
-        dptop(:)     =0.
+c        dvtop(:,:)   =0.
+c        dtetatop(:,:)=0.
+c        dqtop(:,:,:) =0.
+c        dptop(:)     =0.
         dufi(:,:)   =0.
         dvfi(:,:)   =0.
@@ -504 +505 @@
 c      -------------------
          IF (ok_strato) THEN
-           CALL top_bound( vcov,ucov,teta,phi,masse,
-     $                     dutop,dvtop,dtetatop)
-c dqtop=0, dptop=0
-           CALL addfi( dtphys, leapf, forward   ,
-     $                  ucov, vcov, teta , q   ,ps ,
-     $                 dutop, dvtop, dtetatop , dqtop ,dptop  )
+           CALL top_bound(vcov,ucov,teta,masse,dtphys,dutop)
+           dutop(:,:)=dutop(:,:)/dtphys   ! convert to a tendency in (m/s)/s
          ENDIF
 
@@ -543 +540 @@
         ! Sponge layer (if any)
         IF (ok_strato) THEN
-          CALL top_bound(vcov,ucov,teta,phi,
-     $                   masse,dutop,dvtop,dtetatop)
-c dqtop=0, dptop=0
-          CALL addfi( dtvr, leapf, forward   ,
-     $                  ucov, vcov, teta , q   ,ps ,
-     $                 dutop, dvtop, dtetatop , dqtop ,dptop  )
+           CALL top_bound(vcov,ucov,teta,masse,dtvr,dutop)
+           dutop(:,:)=dutop(:,:)/dtvr   ! convert to a tendency in (m/s)/s
         ENDIF ! of IF (ok_strato) 
       ENDIF ! of IF (iflag_phys.EQ.2)