Context Navigation

tracinca_mod.F90 @ 4582

Last change on this file since 4582 was 4500, checked in by Laurent Fairhead, 21 months ago
Some modifications needed for DYNAMICO interface
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory
File size: 7.7 KB

Line
1	!$Id $
2	!
3	MODULE tracinca_mod
4	!
5	! This module prepares and calls the INCA main subroutines.
6	!
7	IMPLICIT NONE
8
9	CHARACTER(len=4),SAVE :: config_inca
10	!$OMP THREADPRIVATE(config_inca)
11	! config_inca='none' => without INCA
12	! config_inca='chem' => INCA with chemistry
13	! config_inca='aero' => INCA with aerosols
14	! config_inca='aeNP' => INCA with aerosols NP (?)
15	CONTAINS
16
17	SUBROUTINE tracinca_init(aerosol,lessivage)
18	! This subroutine initialize some control varaibles.
19
20	USE infotrac_phy, ONLY: nbtr, type_trac
21	IMPLICIT NONE
22
23	! Output variables
24	LOGICAL,DIMENSION(nbtr), INTENT(OUT) :: aerosol
25	LOGICAL,INTENT(OUT) :: lessivage
26
27	! Initialization
28	lessivage =.FALSE.
29	aerosol(:) = .FALSE.
30
31	!--- COHERENCE TEST BETWEEN "type_trac" AND "config_inca"
32	IF(ANY(type_trac == ['inca','inco']) .AND. ALL(config_inca /= ['aero','aeNP','chem'])) CALL abort_physic('tracinca_init', &
33	'INCA enabled, but unknown config_inca = "'//TRIM(config_inca)//'". Please modify "run.def"', 1)
34
35	!--- PROBLEM IF "config_inca" DIFFERS FROM "none" AND INCA HAS NOT BEEN ACTIVATED
36	IF(ALL(type_trac /= ['inca','inco']) .AND. config_inca /= 'none') CALL abort_physic('tracinca_init', &
37	'INCA disabled, but config_inca = "'//TRIM(config_inca)//'" should be "none". Please modify "run.def"', 1)
38
39	END SUBROUTINE tracinca_init
40
41	SUBROUTINE tracinca( &
42	nstep, julien, gmtime, lafin, &
43	pdtphys, t_seri, paprs, pplay, &
44	pmfu, upwd, ftsol, pctsrf, pphis, &
45	pphi, albsol, sh, ch, rh, &
46	cldfra, rneb, diafra, cldliq, &
47	itop_con, ibas_con, pmflxr, pmflxs, &
48	prfl, psfl, aerosol_couple, flxmass_w, &
49	tau_aero, piz_aero, cg_aero, ccm, &
50	rfname, &
51	tr_seri, source)
52
53	!========================================================
54	! -- CHIMIE INCA --
55	!========================================================
56
57	USE dimphy
58	USE infotrac_phy, ONLY: nbtr
59	USE vampir
60	USE indice_sol_mod
61	USE geometry_mod, ONLY: cell_area
62	USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat
63	USE aero_mod, ONLY : naero_grp
64	IMPLICIT NONE
65
66	!==========================================================================
67	! -- DESCRIPTION DES ARGUMENTS --
68	!==========================================================================
69
70
71	! EN ENTREE ...
72	!
73	!Configuration grille,temps:
74	INTEGER,INTENT(IN) :: nstep ! Appel physique
75	INTEGER,INTENT(IN) :: julien ! Jour julien
76	REAL,INTENT(IN) :: gmtime
77	REAL,INTENT(IN) :: pdtphys ! Pas d'integration pour la physique (seconde)
78	LOGICAL,INTENT(IN) :: lafin ! le flag de la fin de la physique
79
80
81	!Physique:
82	!--------
83	REAL,DIMENSION(klon,klev),INTENT(IN) :: t_seri ! Temperature
84	REAL,DIMENSION(klon,klev),INTENT(IN) :: sh ! humidite specifique
85	REAL,DIMENSION(klon,klev),INTENT(IN) :: ch ! eau liquide
86	REAL,DIMENSION(klon,klev),INTENT(IN) :: rh ! humidite relative
87	REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs ! pression pour chaque inter-couche (en Pa)
88	REAL,DIMENSION(klon,klev),INTENT(IN) :: pplay ! pression pour le mileu de chaque couche (en Pa)
89	REAL,DIMENSION(klon,klev),INTENT(IN) :: pphi ! geopotentiel
90	REAL,DIMENSION(klon),INTENT(IN) :: pphis
91	REAL,DIMENSION(klon,klev),INTENT(IN) :: cldliq ! eau condensee pour le radiatif
92	REAL,DIMENSION(klon,klev),INTENT(IN) :: cldfra ! fraction nuageuse (tous les nuages)
93	REAL,DIMENSION(klon,klev),INTENT(IN) :: diafra ! fraction nuageuse (convection ou stratus artificiels)
94	REAL,DIMENSION(klon,klev),INTENT(IN) :: rneb ! fraction nuageuse (grande echelle)
95	INTEGER,DIMENSION(klon),INTENT(IN) :: itop_con
96	INTEGER,DIMENSION(klon),INTENT(IN) :: ibas_con
97	REAL,DIMENSION(klon),INTENT(IN) :: albsol ! albedo surface
98	!
99	!Convection:
100	!----------
101	REAL,DIMENSION(klon,klev),INTENT(IN) :: pmfu ! flux de masse dans le panache montant - Tiedtke
102	REAL,DIMENSION(klon,klev),INTENT(IN) :: upwd ! flux de masse dans le panache montant - Emanuel
103
104	!...Tiedke
105	REAL,DIMENSION(klon,klev+1),INTENT(IN) :: pmflxr, pmflxs ! Flux precipitant de pluie, neige aux interfaces [convection]
106	REAL,DIMENSION(klon,klev+1),INTENT(IN) :: prfl, psfl ! Flux precipitant de pluie, neige aux interfaces [large-scale]
107
108	LOGICAL,INTENT(IN) :: aerosol_couple
109	REAL,DIMENSION(klon,klev),INTENT(IN) :: flxmass_w
110	REAL,DIMENSION(klon,klev,naero_grp,2),INTENT(IN) :: tau_aero
111	REAL,DIMENSION(klon,klev,naero_grp,2),INTENT(IN) :: piz_aero
112	REAL,DIMENSION(klon,klev,naero_grp,2),INTENT(IN) :: cg_aero
113	CHARACTER(len=4),DIMENSION(naero_grp),INTENT(IN) :: rfname
114	REAL,DIMENSION(klon,klev,2),INTENT(IN) :: ccm
115
116	! Arguments necessaires pour les sources et puits de traceur:
117	REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: ftsol ! Temperature du sol (surf)(Kelvin)
118	REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf ! Pourcentage de sol f(nature du sol)
119
120
121	! InOutput argument
122	REAL,DIMENSION(klon,klev,nbtr),INTENT(INOUT) :: tr_seri ! Concentration Traceur [U/KgA]
123
124	! Output arguments
125	REAL,DIMENSION(klon,nbtr), INTENT(OUT) :: source ! a voir lorsque le flux de surface est prescrit
126
127	!=======================================================================================
128	! -- VARIABLES LOCALES TRACEURS --
129	!=======================================================================================
130
131	INTEGER :: k
132	REAL,DIMENSION(klon,klev) :: pdel
133	REAL,DIMENSION(klon,klev) :: zpmfu ! flux de masse dans le panache montant
134	REAL :: calday
135	INTEGER :: ncsec
136
137	CALL VTe(VTphysiq)
138	CALL VTb(VTinca)
139
140	calday = REAL(julien) + gmtime
141	ncsec = NINT (86400.*gmtime)
142
143	DO k = 1, klev
144	pdel(:,k) = paprs(:,k) - paprs (:,k+1)
145	END DO
146
147	#ifdef INCA
148	IF (config_inca == 'aero' .OR. config_inca == 'chem') THEN
149	zpmfu(:,:)=pmfu(:,:)
150	ELSE IF (config_inca == 'aeNP') THEN
151	zpmfu(:,:)=upwd(:,:)
152	ENDIF
153
154	CALL aerosolmain( &
155	aerosol_couple,tr_seri,pdtphys, &
156	pplay,pdel,prfl,pmflxr,psfl, &
157	pmflxs,zpmfu,itop_con,ibas_con, &
158	pphi,cell_area,nstep,rneb,t_seri, &
159	rh,tau_aero,piz_aero,cg_aero, &
160	rfname,ccm,lafin)
161	#endif
162
163
164	#ifdef INCA
165	CALL chemmain (tr_seri, & !mmr
166	nstep, & !nstep
167	calday, & !calday
168	julien, & !ncdate
169	ncsec, & !ncsec
170	1, & !lat
171	pdtphys, & !delt
172	paprs(1,1), & !ps
173	pplay, & !pmid
174	pdel, & !pdel
175	cell_area, &
176	pctsrf(1,1),& !oro
177	ftsol, & !tsurf
178	albsol, & !albs
179	pphi, & !zma
180	pphis, & !phis
181	cldfra, & !cldfr
182	rneb, & !cldfr_st
183	diafra, & !cldfr_cv
184	itop_con, & !cldtop
185	ibas_con, & !cldbot
186	cldliq, & !cwat
187	prfl, & !flxrst
188	pmflxr, & !flxrcv
189	psfl, & !flxsst
190	pmflxs, & !flxscv
191	zpmfu, & !flxupd !--now depends on whether AP or NP
192	flxmass_w, & !flxmass_w
193	t_seri, & !tfld
194	sh, & !sh
195	ch, & !ql
196	rh, & !rh
197	nbp_lon, & !nx
198	nbp_lat, & !ny
199	source )
200	#endif
201
202	CALL VTe(VTinca)
203	CALL VTb(VTphysiq)
204
205
206	END SUBROUTINE tracinca
207
208
209	END MODULE tracinca_mod

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

Download in other formats:

Original Format