Changeset 1793 for LMDZ5/trunk/libf/dyn3dpar/comconst.h

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/dyn3dpar/comconst.h (modified) (2 diffs)

LMDZ5/trunk/libf/dyn3dpar/comconst.h

@@ -6 +6 @@
 
       COMMON/comconsti/im,jm,lllm,imp1,jmp1,lllmm1,lllmp1,lcl,          &
-     &                 iflag_top_bound
+     &                 iflag_top_bound,mode_top_bound
       COMMON/comconstr/dtvr,daysec,                                     &
      & pi,dtphys,dtdiss,rad,r,cpp,kappa,cotot,unsim,g,omeg              &
@@ -30 +30 @@
       REAL omeg ! (rad/s) rotation rate of the planet
       REAL dissip_factz,dissip_deltaz,dissip_zref
-      INTEGER iflag_top_bound
-      REAL tau_top_bound
+! top_bound sponge:
+      INTEGER iflag_top_bound ! sponge type
+      INTEGER mode_top_bound  ! sponge mode
+      REAL tau_top_bound ! inverse of sponge characteristic time scale (Hz)
       REAL daylen ! length of solar day, in 'standard' day length
       REAL year_day ! Number of standard days in a year