Context Navigation

tracinca_mod.F90 @ 5185

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

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

Download in other formats:

Original Format