Home

Context Navigation

source: palm/trunk/SOURCE/swap_timelevel.f90 @ 103

Last change on this file since 103 was 102, checked in by raasch, 17 years ago
preliminary version for coupled runs
Property svn:keywords set to `Id`
File size: 6.2 KB

Line
1	SUBROUTINE swap_timelevel
2
3	!------------------------------------------------------------------------------!
4	! Actual revisions:
5	! -----------------
6	! swaping of uswst, vswst included
7	!
8	! Former revisions:
9	! -----------------
10	! $Id: swap_timelevel.f90 102 2007-07-27 09:09:17Z raasch $
11	!
12	! 95 2007-06-02 16:48:38Z raasch
13	! Swaping of salinity
14	!
15	! 75 2007-03-22 09:54:05Z raasch
16	! moisture renamed humidity
17	!
18	! 19 2007-02-23 04:53:48Z raasch
19	! Swaping of top fluxes
20	!
21	! RCS Log replace by Id keyword, revision history cleaned up
22	!
23	! Revision 1.8 2004/01/28 15:28:18 raasch
24	! Swaping for Runge-Kutta schemes implemented
25	!
26	! Revision 1.1 2000/01/10 10:08:58 10:08:58 raasch (Siegfried Raasch)
27	! Initial revision
28	!
29	!
30	! Description:
31	! ------------
32	! Swap of timelevels of variables after each timestep
33	!------------------------------------------------------------------------------!
34
35	USE arrays_3d
36	USE cpulog
37	USE interfaces
38	USE control_parameters
39
40	IMPLICIT NONE
41
42
43	CALL cpu_log( log_point(28), 'swap_timelevel', 'start' )
44
45	!
46	!-- Incrementing timestep counter
47	timestep_count = timestep_count + 1
48
49	!
50	!-- Swap of 3-level variables
51	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
52
53	SELECT CASE ( MOD( timestep_count, 3 ) )
54
55	CASE ( 0 )
56
57	u_m => u_1; u => u_2; u_p => u_3
58	v_m => v_1; v => v_2; v_p => v_3
59	w_m => w_1; w => w_2; w_p => w_3
60	pt_m => pt_1; pt => pt_2; pt_p => pt_3
61	IF ( .NOT. constant_diffusion ) THEN
62	e_m => e_1; e => e_2; e_p => e_3
63	ENDIF
64	IF ( humidity .OR. passive_scalar ) THEN
65	q_m => q_1; q => q_2; q_p => q_3
66	ENDIF
67
68	CASE ( 1 )
69
70	u_m => u_2; u => u_3; u_p => u_1
71	v_m => v_2; v => v_3; v_p => v_1
72	w_m => w_2; w => w_3; w_p => w_1
73	pt_m => pt_2; pt => pt_3; pt_p => pt_1
74	IF ( .NOT. constant_diffusion ) THEN
75	e_m => e_2; e => e_3; e_p => e_1
76	ENDIF
77	IF ( humidity .OR. passive_scalar ) THEN
78	q_m => q_2; q => q_3; q_p => q_1
79	ENDIF
80
81	CASE ( 2 )
82
83	u_m => u_3; u => u_1; u_p => u_2
84	v_m => v_3; v => v_1; v_p => v_2
85	w_m => w_3; w => w_1; w_p => w_2
86	pt_m => pt_3; pt => pt_1; pt_p => pt_2
87	IF ( .NOT. constant_diffusion ) THEN
88	e_m => e_3; e => e_1; e_p => e_2
89	ENDIF
90	IF ( humidity .OR. passive_scalar ) THEN
91	q_m => q_3; q => q_1; q_p => q_2
92	ENDIF
93
94	END SELECT
95
96	ENDIF
97
98	!
99	!-- Swap of 2-level variables
100	SELECT CASE ( MOD( timestep_count, 2 ) )
101
102	CASE ( 0 )
103
104	IF ( timestep_scheme(1:5) == 'runge' ) THEN
105
106	u => u_1; u_p => u_2
107	v => v_1; v_p => v_2
108	w => w_1; w_p => w_2
109	pt => pt_1; pt_p => pt_2
110	IF ( .NOT. constant_diffusion ) THEN
111	e => e_1; e_p => e_2
112	ENDIF
113	IF ( ocean ) THEN
114	sa => sa_1; sa_p => sa_2
115	ENDIF
116	IF ( humidity .OR. passive_scalar ) THEN
117	q => q_1; q_p => q_2
118	ENDIF
119
120	ELSE
121	!
122	!-- Old timelevels are needed for explicit diffusion within leapfrog
123	IF ( .NOT. constant_diffusion ) THEN
124	kh_m => kh_1; kh => kh_2
125	km_m => km_1; km => km_2
126	IF ( use_surface_fluxes ) THEN
127	usws_m => usws_1; usws => usws_2
128	vsws_m => vsws_1; vsws => vsws_2
129	shf_m => shf_1; shf => shf_2
130	IF ( humidity .OR. passive_scalar ) THEN
131	qsws_m => qsws_1; qsws => qsws_2
132	ENDIF
133	ENDIF
134	IF ( prandtl_layer ) THEN
135	rif_m => rif_1; rif => rif_2
136	ENDIF
137	IF ( use_top_fluxes ) THEN
138	uswst_m => uswst_1; uswst => uswst_2
139	vswst_m => vswst_1; vswst => vswst_2
140	tswst_m => tswst_1; tswst => tswst_2
141	IF ( humidity .OR. passive_scalar ) THEN
142	qswst_m => qswst_1; qswst => qswst_2
143	ENDIF
144	ENDIF
145	ENDIF
146
147	IF ( humidity ) THEN
148	vpt_m => vpt_1; vpt => vpt_2
149	ENDIF
150
151	ENDIF
152
153	CASE ( 1 )
154
155	IF ( timestep_scheme(1:5) == 'runge' ) THEN
156
157	u => u_2; u_p => u_1
158	v => v_2; v_p => v_1
159	w => w_2; w_p => w_1
160	pt => pt_2; pt_p => pt_1
161	IF ( .NOT. constant_diffusion ) THEN
162	e => e_2; e_p => e_1
163	ENDIF
164	IF ( ocean ) THEN
165	sa => sa_2; sa_p => sa_1
166	ENDIF
167	IF ( humidity .OR. passive_scalar ) THEN
168	q => q_2; q_p => q_1
169	ENDIF
170
171	ELSE
172
173	IF ( .NOT. constant_diffusion ) THEN
174	kh_m => kh_2; kh => kh_1
175	km_m => km_2; km => km_1
176	IF ( use_surface_fluxes ) THEN
177	usws_m => usws_2; usws => usws_1
178	vsws_m => vsws_2; vsws => vsws_1
179	shf_m => shf_2; shf => shf_1
180	IF ( humidity .OR. passive_scalar ) THEN
181	qsws_m => qsws_2; qsws => qsws_1
182	ENDIF
183	ENDIF
184	IF ( prandtl_layer ) THEN
185	rif_m => rif_2; rif => rif_1
186	ENDIF
187	IF ( use_top_fluxes ) THEN
188	uswst_m => uswst_2; uswst => uswst_1
189	vswst_m => vswst_2; vswst => vswst_1
190	tswst_m => tswst_2; tswst => tswst_1
191	IF ( humidity .OR. passive_scalar ) THEN
192	qswst_m => qswst_2; qswst => qswst_1
193	ENDIF
194	ENDIF
195	ENDIF
196
197	IF ( humidity ) THEN
198	vpt_m => vpt_2; vpt => vpt_1
199	ENDIF
200
201	ENDIF
202
203	END SELECT
204
205	CALL cpu_log( log_point(28), 'swap_timelevel', 'stop' )
206
207	END SUBROUTINE swap_timelevel
208
209

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |