source: LMDZ5/branches/LMDZ6_rc0/libf/phylmd/tracinca_mod.F90

Last change on this file was 2575, checked in by acozic, 8 years ago

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