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

Timestamp:

Aug 22, 2013, 4:02:07 PM (11 years ago)

Author:

emillour

Message:

Common dynamics: (and collateral adaptations in Venus physics)
Improved cpdet routines in and additional sponge mode:

• Additionnal sponge mode (trigered with "callsponge" flag), in line with the one used in the Generic and Martian GCM. This sponge is called whenever there is a dissipation step.
• Improvement of the cpdet routines : created routines tpot2t_glo_p and t2tpot_glo_p which handle fields on the whole dynamics (scaler) grid, which are more efficient than calling tpot2t_p or t2tpot_p with slabs of data (generated use of intermediate copies of these chunks of data at every call)
• Turned cpdet.F into a module cpdet_mod.F90 (and correspondingly adapted all routines in the Venus physics).

EM

File:

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

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

@@ -16 +16 @@
       USE write_field
       USE control_mod
+      use cpdet_mod, only: cpdet,tpot2t,t2tpot
+      use sponge_mod, only: callsponge,mode_sponge,sponge
       IMPLICIT NONE
 
@@ -197 +199 @@
       ! for CP(T)
       real :: dtec
-      real,external :: cpdet
       real :: ztetaec(ip1jmp1,llm)
 
 c dummy: sinon cette routine n'est jamais compilee...
       if(1.eq.0) then
+#ifdef CPP_PHYS
         CALL init_phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
+#endif
       endif
 
@@ -562 +565 @@
       IF(apdiss) THEN
 
+        ! sponge layer
+        if (callsponge) then
+          CALL sponge(ucov,vcov,teta,ps,dtdiss,mode_sponge)
+        endif
 
 c   calcul de l'energie cinetique avant dissipation