source: LMDZ5/trunk/libf/phylmd/tracinca_mod.F90 @ 2364

Last change on this file since 2364 was 2363, checked in by acozic, 9 years ago

make some modification to allows compilation with inca model

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