Context Navigation

massdair_loc.F @ 3322

Last change on this file since 3322 was 1632, checked in by Laurent Fairhead, 12 years ago

Import initial du répertoire dyn3dmem

Attention! ceci n'est qu'une version préliminaire du code "basse mémoire":
le code contenu dans ce répertoire est basé sur la r1320 et a donc besoin
d'être mis à jour par rapport à la dynamique parallèle d'aujourd'hui.
Ce code est toutefois mis à disposition pour circonvenir à des problèmes
de mémoire que certaines configurations du modèle pourraient rencontrer.
Dans l'état, il compile et tourne sur vargas et au CCRT

Initial import of dyn3dmem

Warning! this is just a preliminary version of the memory light code:
it is based on r1320 of the code and thus needs to be updated before
it can replace the present dyn3dpar code. It is nevertheless put at your
disposal to circumvent some memory problems some LMDZ configurations may
encounter. In its present state, it will compile and run on vargas and CCRT

File size: 3.2 KB

Line
1	SUBROUTINE massdair_loc( p, masse )
2	USE parallel
3	c
4	c *********************************************************************
5	c .... Calcule la masse d'air dans chaque maille ....
6	c *********************************************************************
7	c
8	c Auteurs : P. Le Van , Fr. Hourdin .
9	c ..........
10	c
11	c .. p est un argum. d'entree pour le s-pg ...
12	c .. masse est un argum.de sortie pour le s-pg ...
13	c
14	c .... p est defini aux interfaces des llm couches .....
15	c
16	IMPLICIT NONE
17	c
18	#include "dimensions.h"
19	#include "paramet.h"
20	#include "comconst.h"
21	#include "comgeom.h"
22	c
23	c ..... arguments ....
24	c
25	REAL p(ijb_u:ije_u,llmp1), masse(ijb_u:ije_u,llm)
26
27	c .... Variables locales .....
28
29	INTEGER l,ij
30	INTEGER ijb,ije
31	REAL massemoyn, massemoys
32
33	REAL SSUM
34	EXTERNAL SSUM
35	c
36	c
37	c Methode pour calculer massebx et masseby .
38	c ----------------------------------------
39	c
40	c A chaque point scalaire P (i,j) est affecte 4 coefficients d'aires
41	c alpha1(i,j) calcule au point ( i+1/4,j-1/4 )
42	c alpha2(i,j) calcule au point ( i+1/4,j+1/4 )
43	c alpha3(i,j) calcule au point ( i-1/4,j+1/4 )
44	c alpha4(i,j) calcule au point ( i-1/4,j-1/4 )
45	c
46	c Avec alpha1(i,j) = aire(i+1/4,j-1/4)/ aire(i,j)
47	c
48	c N.B . Pour plus de details, voir s-pg ... iniconst ...
49	c
50	c
51	c
52	c alpha4 . . alpha1 . alpha4
53	c (i,j) (i,j) (i+1,j)
54	c
55	c P . U . . P
56	c (i,j) (i,j) (i+1,j)
57	c
58	c alpha3 . . alpha2 .alpha3
59	c (i,j) (i,j) (i+1,j)
60	c
61	c V . Z . . V
62	c (i,j)
63	c
64	c alpha4 . . alpha1 .alpha4
65	c (i,j+1) (i,j+1) (i+1,j+1)
66	c
67	c P . U . . P
68	c (i,j+1) (i+1,j+1)
69	c
70	c
71	c
72	c On a :
73	c
74	c massebx(i,j) = masse(i ,j) * ( alpha1(i ,j) + alpha2(i,j)) +
75	c masse(i+1,j) * ( alpha3(i+1,j) + alpha4(i+1,j) )
76	c localise au point ... U (i,j) ...
77	c
78	c masseby(i,j) = masse(i,j ) * ( alpha2(i,j ) + alpha3(i,j ) +
79	c masse(i,j+1) * ( alpha1(i,j+1) + alpha4(i,j+1)
80	c localise au point ... V (i,j) ...
81	c
82	c
83	c=======================================================================
84
85
86
87
88	ijb=ij_begin-iip1
89	ije=ij_end+2*iip1
90
91	if (pole_nord) ijb=ij_begin
92	if (pole_sud) ije=ij_end
93
94	c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
95	DO 100 l = 1 , llm
96	c
97	DO ij = ijb, ije
98	masse(ij,l) = airesurg(ij) * ( p(ij,l) - p(ij,l+1) )
99	ENDDO
100	c
101	DO ij = ijb, ije,iip1
102	masse(ij+ iim,l) = masse(ij,l)
103	ENDDO
104	c
105	c DO ij = 1, iim
106	c masse( ij ,l) = masse( ij ,l) * aire( ij )
107	c masse(ij+ip1jm,l) = masse(ij+ip1jm,l) * aire(ij+ip1jm)
108	c ENDDO
109	c massemoyn = SSUM(iim,masse( 1 ,l),1)/ apoln
110	c massemoys = SSUM(iim,masse(ip1jm+1,l),1)/ apols
111	c DO ij = 1, iip1
112	c masse( ij ,l ) = massemoyn
113	c masse(ij+ip1jm,l ) = massemoys
114	c ENDDO
115
116	100 CONTINUE
117	c$OMP END DO NOWAIT
118	c
119	RETURN
120	END

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

Download in other formats:

Original Format