Context Navigation

massbar_loc.F @ 1745

Last change on this file since 1745 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.1 KB

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

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

Download in other formats:

Original Format