Changes between Version 8 and Version 9 of 2021WP/KNL-02_Jerome_RK3_stage1_tsplit

Timestamp:

2021-03-04T10:31:34+01:00 (3 years ago)

Author:

jchanut

Comment:

--

2021WP/KNL-02_Jerome_RK3_stage1_tsplit

@@ -97 +97 @@
       Nnn_e  =  Naa_e
    END DO
+   !
+   ! ----------------------Leapfrog only--------------------------------!
+   IF( .NOT.( kt == nit000 .AND. l_1st_euler ) ) THEN
+      DO_2D( 1, 1, 1, 1 )
+         zun_save = un_adv(ji,jj)
+         zvn_save = vn_adv(ji,jj)
+
+         un_adv(ji,jj) = r1_2 * ( ub2_b(ji,jj) + zun_save - rn_atfp * un_bf(ji,jj) )
+         vn_adv(ji,jj) = r1_2 * ( vb2_b(ji,jj) + zvn_save - rn_atfp * vn_bf(ji,jj) )
+!                          ! Update corrective fluxes for next time step
+         un_bf(ji,jj)  = rn_atfp * un_bf(ji,jj) + ( zun_save - ub2_b(ji,jj) )
+         vn_bf(ji,jj)  = rn_atfp * vn_bf(ji,jj) + ( zvn_save - vb2_b(ji,jj) )
+!                          ! Save integrated transport for next computation
+         ub2_b(ji,jj) = zun_save
+         vb2_b(ji,jj) = zvn_save
+      END_2D
+    END IF
 }}}