source: trunk/LMDZ.GENERIC/libf/phystd/phys_state_var_mod.F90 @ 2613

Last change on this file since 2613 was 2550, checked in by emillour, 3 years ago

Generic GCM:
Some OpenMP fixes in routines initracer.F, nonoro_gwd_ran_mod.F90,
phys_state_var_mod.F90 and sugas_corrk.F90
EM
File size: 13.7 KB
RevLine 
[1835]1!
2! $Id: phys_state_var_mod.F90 1670 2012-10-17 08:42:04Z idelkadi $
3!
4      MODULE phys_state_var_mod
5! Variables sauvegardees pour le startphy.nc
6!======================================================================
7!
8!
9!======================================================================
10! Declaration des variables
11      USE dimphy, only : klon,klev
12      USE comsoil_h, only : nsoilmx
13      use comsaison_h, only: mu0, fract
14      use radcommon_h, only: glat
15      use slab_ice_h, only : noceanmx
16      USE radinc_h, only : L_NSPECTI, L_NSPECTV,naerkind
17      use surfdat_h, only: phisfi, albedodat,  &
18                        zmea, zstd, zsig, zgam, zthe
[1883]19      use turb_mod, only: q2,sensibFlux,wstar,ustar,tstar,hfmax_th,zmax_th
[2299]20      use nonoro_gwd_ran_mod, only: ini_nonoro_gwd_ran, end_nonoro_gwd_ran
[1883]21
[1835]22      real,allocatable,dimension(:,:),save :: ztprevious ! Previous loop Atmospheric Temperature (K)
23! Useful for Dynamical Heating calculation.
24      real,allocatable,dimension(:,:),save :: zuprevious
[1883]25!$OMP THREADPRIVATE(ztprevious,zuprevious)
26
[1835]27      real, dimension(:),allocatable,save ::  tsurf                ! Surface temperature (K).
28      real, dimension(:,:),allocatable,save ::  tsoil              ! Sub-surface temperatures (K).
29      real, dimension(:,:),allocatable,save :: albedo              ! Surface Spectral albedo. By MT2015.
30      real, dimension(:),allocatable,save :: albedo_equivalent     ! Spectral Mean albedo.
[2485]31      real, dimension(:),allocatable,save :: albedo_equivalent1     ! Spectral Mean albedo clear sky.
[1835]32      real, dimension(:),allocatable,save :: albedo_snow_SPECTV    ! Snow Spectral albedo.
33      real, dimension(:),allocatable,save :: albedo_co2_ice_SPECTV ! CO2 ice Spectral albedo.
[2485]34!$OMP THREADPRIVATE(tsurf,tsoil,albedo,albedo_equivalent,albedo_equivalent1,albedo_snow_SPECTV,albedo_co2_ice_SPECTV)
[1835]35
36      real,dimension(:),allocatable,save :: albedo_bareground ! Bare Ground Albedo. By MT 2015.
37      real,dimension(:),allocatable,save :: rnat              ! Defines the type of the grid (ocean,continent,...). By BC.
38!$OMP THREADPRIVATE(albedo_bareground,rnat)
39
40      real,dimension(:),allocatable,save :: emis        ! Thermal IR surface emissivity.
41      real,dimension(:,:),allocatable,save :: dtrad     ! Net atmospheric radiative heating rate (K.s-1).
42      real,dimension(:),allocatable,save :: fluxrad_sky ! Radiative flux from sky absorbed by surface (W.m-2).
43      real,dimension(:),allocatable,save :: fluxrad     ! Net radiative surface flux (W.m-2).
44      real,dimension(:),allocatable,save :: capcal      ! Surface heat capacity (J m-2 K-1).
45      real,dimension(:),allocatable,save :: fluxgrd     ! Surface conduction flux (W.m-2).
46      real,dimension(:,:),allocatable,save :: qsurf     ! Tracer on surface (e.g. kg.m-2).
[1883]47!$OMP THREADPRIVATE(emis,dtrad,fluxrad_sky,fluxrad,capcal,fluxgrd,qsurf)
[1835]48
49      ! FOR DIAGNOSTIC :
50
51      real,dimension(:),allocatable,save :: fluxsurf_lw     ! Incident Long Wave (IR) surface flux (W.m-2).
[2485]52      real,dimension(:),allocatable,save :: fluxsurf_lw1     ! Incident Long Wave (IR) surface flux (W.m-2) clear sky.
[1835]53      real,dimension(:),allocatable,save :: fluxsurf_sw     ! Incident Short Wave (stellar) surface flux (W.m-2).
[2485]54      real,dimension(:),allocatable,save :: fluxsurf_sw1     ! Incident Short Wave (stellar) surface flux (W.m-2) clear sky.
[1835]55      real,dimension(:),allocatable,save :: fluxsurfabs_sw  ! Absorbed Short Wave (stellar) flux by the surface (W.m-2).
[2485]56!$OMP THREADPRIVATE(fluxsurf_lw,fluxsurf_lw1,fluxsurf_sw,fluxsurf_sw1,fluxsurfabs_sw)
[1883]57
[1835]58      real,dimension(:),allocatable,save :: fluxtop_lw      ! Outgoing LW (IR) flux to space (W.m-2).
[2485]59      real,dimension(:),allocatable,save :: fluxtop_lw1      ! Outgoing LW (IR) flux to space (W.m-2) clear sky.
[1835]60      real,dimension(:),allocatable,save :: fluxabs_sw      ! Absorbed SW (stellar) flux (W.m-2).
[2485]61      real,dimension(:),allocatable,save :: fluxabs_sw1      ! Absorbed SW (stellar) flux (W.m-2) clear sky.
[1835]62      real,dimension(:),allocatable,save :: fluxtop_dn      ! Incoming SW (stellar) radiation at the top of the atmosphere (W.m-2).
63      real,dimension(:),allocatable,save :: fluxdyn         ! Horizontal heat transport by dynamics (W.m-2).
[2485]64!$OMP THREADPRIVATE(fluxtop_lw,fluxtop_lw1,fluxabs_sw,fluxabs_sw1,fluxtop_dn,fluxdyn)
[1883]65
[1835]66      real,dimension(:,:),allocatable,save :: OLR_nu        ! Outgoing LW radiation in each band (Normalized to the band width (W/m2/cm-1)).
67      real,dimension(:,:),allocatable,save :: OSR_nu        ! Outgoing SW radiation in each band (Normalized to the band width (W/m2/cm-1)).
[2537]68      real,dimension(:,:),allocatable,save :: GSR_nu        ! Surface SW (or 'Ground' SW) radiation in each band (Normalized to the band width (W/m2/cm-1)).
[1835]69      real,dimension(:,:),allocatable,save :: zdtlw         ! LW heating tendencies (K/s).
70      real,dimension(:,:),allocatable,save :: zdtsw         ! SW heating tendencies (K/s).
71      !real,dimension(:),allocatable,save :: sensibFlux      ! Turbulent flux given by the atmosphere to the surface (W.m-2).
[2537]72!$OMP THREADPRIVATE(OLR_nu,OSR_nu,GSR_nu,zdtlw,zdtsw)
[1835]73
[2133]74      real,dimension(:,:,:),allocatable,save :: int_dtauv   ! VI optical thickness of layers within narrowbands for diags ().
75      real,dimension(:,:,:),allocatable,save :: int_dtaui   ! IR optical thickness of layers within narrowbands for diags ().
76!$OMP THREADPRIVATE(int_dtaui,int_dtauv)
77
[1835]78      real,allocatable,dimension(:),save :: tau_col ! Total Aerosol Optical Depth.
79      real,allocatable,save :: hice(:) ! Oceanic Ice height. by BC
80!$OMP THREADPRIVATE(tau_col,hice)
81
82      real,allocatable,dimension(:,:),save :: cloudfrac  ! Fraction of clouds (%).
83      real,allocatable,dimension(:),save :: totcloudfrac ! Column fraction of clouds (%).
84!$OMP THREADPRIVATE(cloudfrac,totcloudfrac)
85
86      real,allocatable,dimension(:,:),save :: qsurf_hist
87      real,allocatable,dimension(:,:,:),save :: nueffrad ! Aerosol effective radius variance. By RW
[1883]88!$OMP THREADPRIVATE(qsurf_hist,nueffrad)
[1835]89
90      real,allocatable,dimension(:),save :: ice_initial
91      real,allocatable,dimension(:),save :: ice_min
[1883]92!$OMP THREADPRIVATE(ice_initial,ice_min)
[1835]93
[2485]94      real,dimension(:),allocatable,save ::  pctsrf_sic
95      real,dimension(:,:),allocatable,save :: tslab        ! Slab_ocean temperature (K)
96      real,dimension(:),allocatable,save ::  tsea_ice      ! Sea ice temperature (K)
97      real,dimension(:),allocatable,save :: sea_ice
98      real,dimension(:),allocatable,save :: zmasq
99      integer,dimension(:),allocatable,save ::knindex
[1835]100      real,allocatable,dimension(:,:,:),save :: reffrad
[2485]101      real,dimension(:), allocatable,save :: fluxgrdocean  ! Surface conduction flux (W.m-2) with slab_ocean.
102!$OMP THREADPRIVATE(pctsrf_sic,tslab,tsea_ice,sea_ice,zmasq,knindex,reffrad,fluxgrdocean)
103
104      real,dimension(:,:),allocatable,save :: dEzdiff   ! Turbulent diffusion heating (W.m-2)
105      real,dimension(:),allocatable,save :: dEdiff      ! Integrated turbulent diffusion heating (W.m-2)
106      real,dimension(:),allocatable,save :: dEdiffs
107      real,dimension(:,:),allocatable,save :: dEzRadsw  ! Radiative heating (W.m-2)
108      real,dimension(:,:),allocatable,save :: dEzRadlw  ! Radiative heating (W.m-2)
[2550]109!$OMP THREADPRIVATE(dEzdiff,dEdiff,dEdiffs,dEzRadsw,dEzRadlw)
[2485]110
111      real,dimension(:),allocatable,save :: madjdE      ! Heat from moistadj (W.m-2)
112      real,dimension(:,:),allocatable,save :: madjdEz   ! Heat from moistadj (W.m-2)
113      real,dimension(:),allocatable,save :: lscaledE    ! Heat from largescale (W.m-2)
114      real,dimension(:,:),allocatable,save :: lscaledEz ! Heat from largescale (W.m-2)
115!$OMP THREADPRIVATE(madjdE,madjdEz,lscaledE,lscaledEz)
[1842]116     
[2485]117      real,dimension(:),allocatable,save :: H2Omaxcol ! Maximum possible H2O column amount (at 100% saturation) (kg.m-2).
118!$OMP THREADPRIVATE(H2Omaxcol)
119     
[1835]120CONTAINS
121
122!======================================================================
[1842]123SUBROUTINE phys_state_var_init(nqtot)
[1835]124
125IMPLICIT NONE
126
[1883]127        integer,intent(in) :: nqtot
[1842]128
[1835]129!  Parametres de l'Orographie a l'Echelle Sous-Maille (OESM):
130!
131!zmea(:)   ! orographie moyenne
132!zstd(:)   ! deviation standard de l'OESM
133!zsig(:)   ! pente de l'OESM
134!zgam(:)   ! anisotropie de l'OESM
135!zthe(:)   ! orientation de l'OESM
136!zpic(:)   ! Maximum de l'OESM
137!zval(:)   ! Minimum de l'OESM
138!rugoro(:) ! longueur de rugosite de l'OESM
[2485]139
[1835]140        ALLOCATE(phisfi(klon))
141        ALLOCATE(tsurf(klon))
142        ALLOCATE(tsoil(klon,nsoilmx))
143        ALLOCATE(albedo(klon,L_NSPECTV))
144        ALLOCATE(albedo_equivalent(klon))
[2485]145        ALLOCATE(albedo_equivalent1(klon))
[1835]146        ALLOCATE(albedo_snow_SPECTV(L_NSPECTV))
147        ALLOCATE(albedo_co2_ice_SPECTV(L_NSPECTV))
148        ALLOCATE(albedo_bareground(klon))
149        ALLOCATE(rnat(klon))
150        ALLOCATE(emis(klon))
151        ALLOCATE(dtrad(klon,klev))
152        ALLOCATE(fluxrad_sky(klon))
153        ALLOCATE(fluxrad(klon))
154        ALLOCATE(capcal(klon))
155        ALLOCATE(fluxgrd(klon))
[2485]156        ALLOCATE(fluxgrdocean(klon))
[1835]157        ALLOCATE(qsurf(klon,nqtot))
158        ALLOCATE(q2(klon,klev+1))
159        ALLOCATE(ztprevious(klon,klev))
160        ALLOCATE(zuprevious(klon,klev))
161        ALLOCATE(cloudfrac(klon,klev))
162        ALLOCATE(totcloudfrac(klon))
163        ALLOCATE(hice(klon))
164        ALLOCATE(qsurf_hist(klon,nqtot))
165        ALLOCATE(reffrad(klon,klev,naerkind))
166        ALLOCATE(nueffrad(klon,klev,naerkind))
167        ALLOCATE(ice_initial(klon))
168        ALLOCATE(ice_min(klon))
169        ALLOCATE(fluxsurf_lw(klon))
[2485]170        ALLOCATE(fluxsurf_lw1(klon))
[1835]171        ALLOCATE(fluxsurf_sw(klon))
[2485]172        ALLOCATE(fluxsurf_sw1(klon))
[1835]173        ALLOCATE(fluxsurfabs_sw(klon))
174        ALLOCATE(fluxtop_lw(klon))
[2485]175        ALLOCATE(fluxtop_lw1(klon))
[1835]176        ALLOCATE(fluxabs_sw(klon))
[2485]177        ALLOCATE(fluxabs_sw1(klon))
[1835]178        ALLOCATE(fluxtop_dn(klon))
179        ALLOCATE(fluxdyn(klon))
180        ALLOCATE(OLR_nu(klon,L_NSPECTI))
181        ALLOCATE(OSR_nu(klon,L_NSPECTV))
[2537]182        ALLOCATE(GSR_nu(klon,L_NSPECTV))
[2133]183        ALLOCATE(int_dtaui(klon,klev,L_NSPECTI))
184        ALLOCATE(int_dtauv(klon,klev,L_NSPECTV))
[1835]185        ALLOCATE(sensibFlux(klon))
186        ALLOCATE(zdtlw(klon,klev))
187        ALLOCATE(zdtsw(klon,klev))
188        ALLOCATE(tau_col(klon))
189        ALLOCATE(pctsrf_sic(klon))
190        ALLOCATE(tslab(klon,noceanmx))
191        ALLOCATE(tsea_ice(klon))
192        ALLOCATE(sea_ice(klon))
193        ALLOCATE(zmasq(klon))
194        ALLOCATE(knindex(klon))
[2485]195        ALLOCATE(dEzdiff(klon,klev))
196        ALLOCATE(dEdiff(klon))
197        ALLOCATE(dEdiffs(klon))
198        ALLOCATE(dEzRadsw(klon,klev))
199        ALLOCATE(dEzRadlw(klon,klev))
200        ALLOCATE(madjdE(klon))
201        ALLOCATE(madjdEz(klon,klev))
202        ALLOCATE(lscaledE(klon))
203        ALLOCATE(lscaledEz(klon,klev))
204        ALLOCATE(H2Omaxcol(klon))
[1835]205        ! This is defined in comsaison_h
206        ALLOCATE(mu0(klon))
207        ALLOCATE(fract(klon))
208         ! This is defined in radcommon_h
209        ALLOCATE(glat(klon))
210        ALLOCATE(albedodat(klon))
211        ALLOCATE(zmea(klon))
212        ALLOCATE(zstd(klon))
213        ALLOCATE(zsig(klon))
214        ALLOCATE(zgam(klon))
215        ALLOCATE(zthe(klon))
216        !allocate(l0(klon))
217        allocate(wstar(klon))
218        allocate(ustar(klon))
219        allocate(tstar(klon))
220        allocate(hfmax_th(klon))
221        allocate(zmax_th(klon))
[2299]222        ! allocate arrays in "nonoro_gwd_ran_mod"
223        call ini_nonoro_gwd_ran(klon,klev)
[1835]224END SUBROUTINE phys_state_var_init
225
226!======================================================================
227SUBROUTINE phys_state_var_end
228
229IMPLICIT NONE
230
231        DEALLOCATE(tsurf)
232        DEALLOCATE(tsoil)
233        DEALLOCATE(albedo)
234        DEALLOCATE(albedo_equivalent)
[2485]235        DEALLOCATE(albedo_equivalent1)
[1835]236        DEALLOCATE(albedo_snow_SPECTV)
237        DEALLOCATE(albedo_co2_ice_SPECTV)
238        DEALLOCATE(albedo_bareground)
239        DEALLOCATE(rnat)
240        DEALLOCATE(emis)
241        DEALLOCATE(dtrad)
242        DEALLOCATE(fluxrad_sky)
243        DEALLOCATE(fluxrad)
244        DEALLOCATE(capcal)
245
246        DEALLOCATE(fluxgrd)
[2485]247        DEALLOCATE(fluxgrdocean)
[1835]248        DEALLOCATE(qsurf)
249        DEALLOCATE(q2)
250        DEALLOCATE(ztprevious)
251        DEALLOCATE(zuprevious)
252        DEALLOCATE(cloudfrac)
253        DEALLOCATE(totcloudfrac)
254        DEALLOCATE(hice)
255        DEALLOCATE(qsurf_hist)
256        DEALLOCATE(reffrad)
257        DEALLOCATE(nueffrad)
258        DEALLOCATE(ice_initial)
259        DEALLOCATE(ice_min)
260        DEALLOCATE(fluxsurf_lw)
[2485]261        DEALLOCATE(fluxsurf_lw1)
[1835]262        DEALLOCATE(fluxsurf_sw)
[2485]263        DEALLOCATE(fluxsurf_sw1)
[1835]264        DEALLOCATE(fluxsurfabs_sw)
265        DEALLOCATE(fluxtop_lw)
[2485]266        DEALLOCATE(fluxtop_lw1)
[1835]267        DEALLOCATE(fluxabs_sw)
[2485]268        DEALLOCATE(fluxabs_sw1)
[1835]269        DEALLOCATE(fluxtop_dn)
270        DEALLOCATE(fluxdyn)
271        DEALLOCATE(OLR_nu)
272        DEALLOCATE(OSR_nu)
[2537]273        DEALLOCATE(GSR_nu)
[2133]274        DEALLOCATE(int_dtaui)
275        DEALLOCATE(int_dtauv)
[1835]276        DEALLOCATE(sensibFlux)
277        DEALLOCATE(zdtlw)
278        DEALLOCATE(zdtsw)
279        DEALLOCATE(tau_col)
280        DEALLOCATE(pctsrf_sic)
281        DEALLOCATE(tslab)
282        DEALLOCATE(tsea_ice)
283        DEALLOCATE(sea_ice)
284        DEALLOCATE(zmasq)
285        DEALLOCATE(knindex)
[2485]286        DEALLOCATE(dEzdiff)
287        DEALLOCATE(dEdiff)
288        DEALLOCATE(dEdiffs)
289        DEALLOCATE(dEzRadsw)
290        DEALLOCATE(dEzRadlw)
291        DEALLOCATE(madjdE)
292        DEALLOCATE(madjdEz)
293        DEALLOCATE(lscaledE)
294        DEALLOCATE(lscaledEz)
295        DEALLOCATE(H2Omaxcol)       
[1835]296        DEALLOCATE(mu0)
297        DEALLOCATE(fract)
298        DEALLOCATE(glat)
299        DEALLOCATE(phisfi)
300        DEALLOCATE(albedodat)
301        DEALLOCATE(zmea)
302        DEALLOCATE(zstd)
303        DEALLOCATE(zsig)
304        DEALLOCATE(zgam)
305        DEALLOCATE(zthe)
306        !deallocate(l0)
307        deallocate(wstar)
308        deallocate(ustar)
309        deallocate(tstar)
310        deallocate(hfmax_th)
311        deallocate(zmax_th)
[2299]312        ! deallocate arrays in "nonoro_gwd_ran_mod"
313        call end_nonoro_gwd_ran
[1835]314
315END SUBROUTINE phys_state_var_end
316
317      END MODULE phys_state_var_mod
Note: See TracBrowser for help on using the repository browser.