Changeset 1793 for LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.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/dyn3dmem/leapfrog_loc.F (modified) (1 diff)

LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F

@@ -1125 +1125 @@
         ! Sponge layer (if any)
         IF (ok_strato) THEN
-          ! set dufi,dvfi,... to zero
-          ijb=ij_begin
-          ije=ij_end
-!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-          do l=1,llm
-            dufi(ijb:ije,l)=0
-            dtetafi(ijb:ije,l)=0
-            dqfi(ijb:ije,l,1:nqtot)=0
-          enddo
-!$OMP END DO
-!$OMP MASTER
-          dpfi(ijb:ije)=0
-!$OMP END MASTER
-          ijb=ij_begin
-          ije=ij_end
-          if (pole_sud) ije=ije-iip1
-!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-          do l=1,llm
-            dvfi(ijb:ije,l)=0
-          enddo
-!$OMP END DO
-
-          CALL top_bound_loc(vcov,ucov,teta,masse,dufi,dvfi,dtetafi)
-          CALL addfi_loc( dtvr, leapf, forward   ,
-     $                  ucov, vcov, teta , q   ,ps ,
-     $                 dufi, dvfi, dtetafi , dqfi ,dpfi  )
+          CALL top_bound_loc(vcov,ucov,teta,masse,dtvr)
 !$OMP BARRIER
         ENDIF ! of IF (ok_strato)