Changeset 371 for codes/icosagcm

Timestamp:

04/08/16 09:18:36 (8 years ago)

Author:

dubos

Message:

Bugfix : NH slow tendencies for W,Phi were ignored

File:

• codes/icosagcm/trunk/src/hevi_scheme.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/hevi_scheme.f90

@@ -19 +19 @@
     ! See Weller et al., 2013 - ARK2 scheme Fig. 2
     REAL(rstd) :: dt
-    REAL(rstd), PARAMETER :: alpha=(3.+SQRT(8.))/6., delta=.5/SQRT(2.), gamma=1.-2.*delta
+    REAL(rstd), PARAMETER :: delta=.5/SQRT(2.), gamma=1.-2.*delta
+!    REAL(rstd), PARAMETER :: alpha=(3.+SQRT(8.))/6. ! original value in GKC2013
+    REAL(rstd), PARAMETER :: alpha=0.7
     REAL(rstd), PARAMETER, DIMENSION(3) :: wj = (/delta,delta,gamma/)
     CALL set_coefs_hevi(dt, &
@@ -84 +86 @@
           CALL update(cjl(l,j), f_u, f_du_fast(:,l))
           IF(.NOT. hydrostatic) THEN
- !            CALL update(bjl(l,j), f_W, f_dW_slow(:,l)) ! slow tendencies of w, Phi not implemented yet
+             CALL update(bjl(l,j), f_W, f_dW_slow(:,l))
              CALL update(cjl(l,j), f_W, f_dW_fast(:,l))
- !            CALL update(bjl(l,j), f_geopot, f_dPhi_slow(:,l))
+             CALL update(bjl(l,j), f_geopot, f_dPhi_slow(:,l))
              CALL update(cjl(l,j), f_geopot, f_dPhi_fast(:,l))
           END IF