Changeset 339

Timestamp:

07/15/15 16:06:42 (9 years ago)

Author:

dubos

Message:

Kinnmark (1984) RK2.5 scheme

File:

• codes/icosagcm/trunk/src/timeloop_gcm.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/timeloop_gcm.f90

@@ -6 +6 @@
   PUBLIC :: init_timeloop, timeloop
 
-  INTEGER, PARAMETER :: euler=1, rk4=2, mlf=3
+  INTEGER, PARAMETER :: euler=1, rk4=2, mlf=3, rk25=4
   INTEGER, PARAMETER :: itau_sync=10
+  REAL(rstd), DIMENSION(4), PARAMETER :: coef_rk4 = (/ .25, 1./3., .5, 1. /)
+  REAL(rstd), DIMENSION(5), PARAMETER :: coef_rk25 = (/ .25, 1./6., 3./8., .5, 1. /)
 
   TYPE(t_message),SAVE :: req_ps0, req_mass0, req_theta_rhodz0, req_u0, req_q0
@@ -95 +97 @@
         scheme=rk4
         nb_stage=4
+      CASE ('RK2.5')
+        scheme=rk25
+        nb_stage=5
       CASE ('leapfrog_matsuno')
         scheme=mlf
@@ -248 +253 @@
           CALL euler_scheme(.TRUE.)
        CASE (rk4)
-          CALL rk_scheme(stage)
+          CALL rk_scheme(stage, coef_rk4)
+       CASE (rk25)
+          CALL rk_scheme(stage, coef_rk25)
        CASE (mlf)
           CALL  leapfrog_matsuno_scheme(stage)
@@ -388 +395 @@
     END SUBROUTINE Euler_scheme
 
-    SUBROUTINE RK_scheme(stage)
+    SUBROUTINE RK_scheme(stage,coef)
       IMPLICIT NONE
       INTEGER :: ind, stage
-      REAL(rstd), DIMENSION(4), PARAMETER :: coef = (/ .25, 1./3., .5, 1. /)
+      REAL(rstd), INTENT(IN) :: coef(:)
       REAL(rstd) :: tau
       INTEGER :: i,j,ij,l