Changeset 108 for trunk/libf/dyn3d/leapfrog.F

Timestamp:

Apr 12, 2011, 11:16:02 AM (14 years ago)

Author:

slebonnois

Message:

Sebastien Lebonnois: sponge layer et dissip horizontale.

File:

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

trunk/libf/dyn3d/leapfrog.F

@@ -111 +111 @@
       REAL dvfi(ip1jm,llm),dufi(ip1jmp1,llm)
       REAL dtetafi(ip1jmp1,llm),dqfi(ip1jmp1,llm,nqtot),dpfi(ip1jmp1)
-
-      real :: duspg(ip1jmp1,llm) ! for bilan_dyn
 
 c   variables pour le fichier histoire
@@ -463 +461 @@
 c      -------------------
          IF (ok_strato) THEN
-           CALL top_bound( vcov,ucov,teta,masse,dufi,dvfi,dtetafi)
+           CALL top_bound( vcov,ucov,teta,masse,dutop,dvtop,dtetatop)
+c dqtop=0, dptop=0
+           CALL addfi( dtphys, leapf, forward   ,
+     $                  ucov, vcov, teta , q   ,ps ,
+     $                 dutop, dvtop, dtetatop , dqtop ,dptop  )
          ENDIF
-c dqtop=0, dptop=0
-       
-c      ajout des tendances physiques:
-c      ------------------------------
-          CALL addfi( dtphys, leapf, forward   ,
-     $                  ucov, vcov, teta , q   ,ps ,
-     $                 dufi, dvfi, dtetafi , dqfi ,dpfi  )
-c
+
 c  Diagnostique de conservation de l'énergie : difference
          IF (ip_ebil_dyn.ge.1 ) THEN 
@@ -512 +507 @@
         ! Sponge layer (if any)
         IF (ok_strato) THEN
-          dutop(:,:)=0.
-          dvtop(:,:)=0.
-          dtetatop(:,:)=0.
-          dqtop(:,:,:)=0.
-          dptop(:)=0.
           CALL top_bound(vcov,ucov,teta,masse,dutop,dvtop,dtetatop)
+c dqtop=0, dptop=0
           CALL addfi( dtvr, leapf, forward   ,
      $                  ucov, vcov, teta , q   ,ps ,
@@ -669 +660 @@
                  CALL bilan_dyn(dtvr*iperiod,dtvr*day_step*periodav,
      &                 ps,masse,pk,pbaru,pbarv,teta,phi,ucov,vcov,
-     &                 du,dudis,duspg,dufi)
+     &                 du,dudis,dutop,dufi)
 #endif
                END IF
@@ -801 +792 @@
                  CALL bilan_dyn(dtvr*iperiod,dtvr*day_step*periodav,
      &                 ps,masse,pk,pbaru,pbarv,teta,phi,ucov,vcov,
-     &                 du,dudis,duspg,dufi)
+     &                 du,dudis,dutop,dufi)
 #endif
                ENDIF