LMDZ
ini_bilKP_ave.h
Go to the documentation of this file.
1 c
2 c $Id: ini_bilKP_ave.h 1907 2013-11-26 13:10:46Z lguez $
3 c
4  IF (ok_journe) THEN
5 c
6  zsto = dtime
7  zout = ecrit_day
8  typeval=tave
9 c
10  idayref = day_ref
11  CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12 c
13 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
14 cym DO i = 1, iim
15 cym zx_lon(i,1) = rlon(i+1)
16 cym zx_lon(i,jjmp1) = rlon(i+1)
17 cym ENDDO
18  DO ll=1,klev
19  znivsig(ll)=REAL(ll)
20  ENDDO
21 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
22 cym write(*,*)'zx_lon = ',zx_lon(:,1)
23 cym write(*,*)'zx_lat = ',zx_lat(1,:)
24 cym CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1,
25 cym . zx_lat(1,:),
26 cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
27 cym . nhori, nid_bilKPave)
28  CALL histbeg_phy("histbilKP_ave", itau_phy, zjulian, dtime,
29  . nhori, nid_bilKPave)
30 
31  write(*,*)'Journee ', itau_phy, zjulian
32  CALL histvert(nid_bilKPave, "presnivs",
33  . "Vertical levels","mb",
34  . klev, presnivs/100., nvert)
35 c
36 c
37 c Champs 3D:
38 c
39  CALL histdef(nid_bilKPave,"ue",
40  . "Zonal energy transport","-",
41  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
42  . typeval, zsto,zout)
43 c
44  CALL histdef(nid_bilKPave,"ve",
45  . "Merid energy transport","-",
46  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
47  . typeval, zsto,zout)
48 c
49  CALL histdef(nid_bilKPave,"uq",
50  . "Zonal humidity transport","-",
51  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
52  . typeval, zsto,zout)
53 c
54  CALL histdef(nid_bilKPave,"vq",
55  . "Merid humidity transport","-",
56  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
57  . typeval, zsto,zout)
58 c
59 c Champs 3D:
60 c
61  CALL histdef(nid_bilKPave,"temp",
62  . "Air temperature","K",
63  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
64  . typeval, zsto,zout)
65 c
66  CALL histdef(nid_bilKPave,"ovap",
67  . "Specific humidity","Kg/Kg",
68  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
69  . typeval, zsto,zout)
70 c
71  CALL histdef(nid_bilKPave,"geop",
72  . "Geopotential height","m",
73  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
74  . typeval, zsto,zout)
75 c
76  CALL histdef(nid_bilKPave,"vitu",
77  . "Zonal wind","m/s",
78  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
79  . typeval, zsto,zout)
80 c
81  CALL histdef(nid_bilKPave,"vitv",
82  . "Meridional wind","m/s",
83  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
84  . typeval, zsto,zout)
85 c
86  CALL histdef(nid_bilKPave, "vitw",
87  . "Vertical wind", "m/s",
88  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
89  . typeval, zsto,zout)
90 c
91  CALL histdef(nid_bilKPave, "pres",
92  . "Inter-Layer Air pressure",
93  . "Pa",
94  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
95  . typeval, zsto,zout)
96 c
97  CALL histdef(nid_bilKPave, "play",
98  . "Mean-Layer Air pressure",
99  . "Pa",
100  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
101  . typeval, zsto,zout)
102 c
103  CALL histdef(nid_bilKPave, "oliq",
104  . "Liquid water content",
105  . "kg/kg",
106  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
107  . typeval, zsto,zout)
108 c
109  CALL histdef(nid_bilKPave, "dtdyn",
110  . "Dynamics dT", "K/s",
111  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
112  . typeval, zsto,zout)
113 c
114  CALL histdef(nid_bilKPave, "dqdyn",
115  . "Dynamics dQ", "Kg/Kg/s",
116  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
117  . typeval, zsto,zout)
118 c
119  CALL histdef(nid_bilKPave, "dtcon",
120  . "Convection dT", "K/s",
121  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
122  . typeval, zsto,zout)
123 c
124  CALL histdef(nid_bilKPave, "ducon",
125  . "Convection du", "m/s2",
126  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
127  . typeval, zsto,zout)
128 c
129  CALL histdef(nid_bilKPave, "dvcon",
130  . "Convection dv", "m/s2",
131  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
132  . typeval, zsto,zout)
133 c
134  CALL histdef(nid_bilKPave,"dqcon",
135  . "Convection dQ","Kg/Kg/s",
136  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
137  . typeval, zsto,zout)
138 c
139  CALL histdef(nid_bilKPave, "dtlsc",
140  . "Condensation dT", "K/s",
141  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
142  . typeval, zsto,zout)
143 c
144  CALL histdef(nid_bilKPave,"dqlsc",
145  . "Condensation dQ","Kg/Kg/s",
146  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
147  . typeval, zsto,zout)
148 c
149  CALL histdef(nid_bilKPave,"dtvdf",
150  . "Boundary-layer dT","K/s",
151  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
152  . typeval, zsto,zout)
153 c
154  CALL histdef(nid_bilKPave, "dqvdf",
155  . "Boundary-layer dQ",
156  . "Kg/Kg/s",
157  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
158  . typeval, zsto,zout)
159 c
160  CALL histdef(nid_bilKPave,"dtajs",
161  . "Ajustement sec dT","K/s",
162  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
163  . typeval, zsto,zout)
164 c
165  CALL histdef(nid_bilKPave, "dqajs",
166  . "Ajustement sec dQ",
167  . "Kg/Kg/s",
168  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
169  . typeval, zsto,zout)
170 c
171  CALL histdef(nid_bilKPave, "dteva",
172  . "Reevaporation dT", "K/s",
173  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
174  . typeval, zsto,zout)
175 c
176  CALL histdef(nid_bilKPave,"dqeva",
177  . "Reevaporation dQ",
178  . "Kg/Kg/s",
179  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
180  . typeval, zsto,zout)
181 
182 c
183  CALL histdef(nid_bilKPave, "dtswr",
184  . "SW radiation dT", "K/s",
185  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
186  . typeval, zsto,zout)
187 c
188  CALL histdef(nid_bilKPave, "dtsw0",
189  . "SW radiation dT", "K/s",
190  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
191  . typeval, zsto,zout)
192 c
193  CALL histdef(nid_bilKPave, "dtlwr",
194  . "LW radiation dT", "K/s",
195  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
196  . typeval, zsto,zout)
197 c
198  CALL histdef(nid_bilKPave, "dtlw0",
199  . "LW radiation dT", "K/s",
200  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
201  . typeval, zsto,zout)
202 c
203  CALL histdef(nid_bilKPave,"duvdf",
204  . "Boundary-layer dU","m/s2",
205  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
206  . typeval, zsto,zout)
207 c
208  CALL histdef(nid_bilKPave,"dvvdf",
209  . "Boundary-layer dV","m/s2",
210  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
211  . typeval, zsto,zout)
212 c
213  IF (ok_orodr) THEN
214  IF (ok_orolf) THEN
215  CALL histdef(nid_bilKPave, "duoli",
216  . "Orography dU", "m/s2",
217  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
218  . typeval, zsto,zout)
219 c
220  CALL histdef(nid_bilKPave, "dvoli",
221  . "Orography dV", "m/s2",
222  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
223  . typeval, zsto,zout)
224 c
225  ENDIF
226  ENDIF
227 C
228  CALL histdef(nid_bilKPave, "duphy",
229  . "Physiq dU","-",
230  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
231  . typeval, zsto,zout)
232 C
233  CALL histdef(nid_bilKPave, "dvphy",
234  . "Physiq dV","-",
235  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
236  . typeval, zsto,zout)
237 C
238  CALL histdef(nid_bilKPave, "dtphy",
239  . "Physiq dT","-",
240  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
241  . typeval, zsto,zout)
242 C
243  CALL histdef(nid_bilKPave, "dqphy",
244  . "Physiq dQ","-",
245  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
246  . typeval, zsto,zout)
247 C
248  CALL histdef(nid_bilKPave, "dqlphy",
249  . "Physiq dQl","-",
250  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
251  . typeval, zsto,zout)
252 C
253 C
254  CALL histend(nid_bilKPave)
255 c
256  ndex2d = 0
257  ndex3d = 0
258 c
259  ENDIF ! fin de test sur ok_journe
ok_orolf
!$Header!c include clesph0 h c COMMON clesph0 ok_orolf
Definition: clesph0.h:6
zsto
!$Id La frequence de sortie REGDYN est ecrit_mth c iwmax il faut mettre la meme valeur pour c zout et zsto c dtime est passe par ailleurs a histbeg c c zsto
Definition: ini_histREGDYN.h:33
klon
INTERFACE SUBROUTINE RRTM_ECRT_140GP klon
Definition: rrtm_ecrt_140gp.intfb.h:3
$Id
c c $Id
Definition: ini_bilKP_ave.h:11
ymds2ju
!$Id zjulian!correction pour l heure initiale!jyg!jyg CALL ymds2ju(annee_ref, 1, day_ref, hour, zjulian)!jyg CALL histbeg_phy("histrac"
ok_orodr
!$Header!c include clesph0 h c COMMON clesph0 ok_orodr
Definition: clesph0.h:6
zjulian
c c zjulian c cym CALL iim cym klev cym zjulian
Definition: ini_bilKP_ave.h:26
histbeg_phy
c c zjulian c cym CALL iim cym klev cym cym nid_bilKPave CALL histbeg_phy("histbilKP_ave", itau_phy, zjulian, dtime,.nhori, nid_bilKPave) write(*
cosp_output_mod::nvert
integer, dimension(3), save nvert
Definition: cosp_output_mod.F90:17
REAL
!$Id ok_orolf LOGICAL ok_limitvrai LOGICAL ok_all_xml INTEGER iflag_ener_conserv REAL solaire REAL(kind=8) RCO2
phys_local_var_mod::s2
real, dimension(:), allocatable, save s2
Definition: phys_local_var_mod.F90:334
nhori
c c zjulian c cym CALL iim cym klev cym cym nhori
Definition: ini_bilKP_ave.h:26
Journee
c c zjulian c cym CALL iim cym klev cym cym nid_bilKPave CALL * Journee
Definition: ini_bilKP_ave.h:31
dimphy::klev
integer, save klev
Definition: dimphy.F90:7
phys_local_var_mod::vq
real, dimension(:), allocatable, save vq
Definition: phys_local_var_mod.F90:273
itau_phy
c c zjulian c cym CALL iim cym klev cym itau_phy
Definition: ini_bilKP_ave.h:26
phys_local_var_mod::ve
real, dimension(:), allocatable, save ve
Definition: phys_local_var_mod.F90:273
i
!$Id klon initialisation mois suivants day_rain itap ENDIF!Calcul fin de nday_rain calcul nday_rain itap DO i
Definition: calcul_divers.h:24
pressure
!$Id vert_prof_dissip LOGICAL lstardis INTEGER niterh integer vert_prof_dissip!vertical profile of horizontal dissipation!Allowed function of pressure
Definition: comdissnew.h:13
IF
!$Id klon IF(pctsrf(i, is_ter).GT.0.) THEN paire_ter(i)
cym
c c zjulian c cym CALL iim cym klev cym cym
Definition: ini_bilKP_ave.h:26
znivsig
c c zjulian c cym CALL iim cym klev znivsig(ll)
histvert
!$Id zjulian!correction pour l heure initiale!jyg!jyg CALL nid_tra CALL histvert(nid_tra,"presnivs","Vertical levels","Pa", klev, presnivs, nvert,"down") zsto
phys_local_var_mod::uq
real, dimension(:), allocatable, save uq
Definition: phys_local_var_mod.F90:273
day_ref
!$Id day_ref
Definition: temps.h:15
height
!$Id nt_fire height
Definition: 1D_read_forc_cases.h:21
zout
!$Id La frequence de sortie REGDYN est ecrit_mth c iwmax il faut mettre la meme valeur pour c zout et zsto c dtime est passe par ailleurs a histbeg c c zout
Definition: ini_histREGDYN.h:31
idayref
c c idayref
Definition: ini_bilKP_ave.h:11
jjmp1
c c zjulian c cym CALL iim cym klev jjmp1
Definition: ini_bilKP_ave.h:24
rlat
c c $Id c nbregdyn DO klon c rlat(i) c ENDIF!lon c ENDIF!lat ENDIF!pctsrf ENDDO!klon ENDDO!nbregdyn cIM 190504 ENDIF!ok_regdyn cIM somme de toutes les nhistoW BEG IF(debut) THEN DO nreg
sur
c c zjulian c cym CALL iim cym klev iim cym jjmp1 cym On stoke le fichier bilKP instantanne sur
Definition: ini_bilKP_ins.h:41
ll
!$Id Turb_fcg!implicit none!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc!cette routine permet d obtenir hq et ainsi de!pouvoir calculer la convergence et le cisaillement dans la physiq!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc INTEGER klev REAL ll
Definition: 1Dconv.h:27
phys_local_var_mod::ue
real, dimension(:), allocatable, save ue
Definition: phys_local_var_mod.F90:273
call_calfis_mod::dv
real, dimension(:,:), pointer, save dv
Definition: call_calfis_mod.F90:19
call_calfis_mod::du
real, dimension(:,:), pointer, save du
Definition: call_calfis_mod.F90:18
gr_fi_ecrit
c c zjulian c cym CALL gr_fi_ecrit(1, klon, iim, jjmp1, rlon, zx_lon) cym DO i
iim
c c zjulian c cym CALL iim cym klev iim
Definition: ini_bilKP_ave.h:24
zx_lat
c c zjulian c cym CALL iim cym klev * zx_lat
Definition: ini_bilKP_ave.h:23
play
!$Id Turb_fcg!implicit none!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc!cette routine permet d obtenir hq et ainsi de!pouvoir calculer la convergence et le cisaillement dans la physiq!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc INTEGER klev REAL in CHARACTER file_fordat COMMON com1_phys_gcss play
Definition: 1Dconv.h:27
de
do llm!au dessus de
Definition: 1D_nudge_sandu_astex.h:1
histbeg
c c zjulian c cym CALL iim cym klev iim cym jjmp1 cym On stoke le fichier bilKP instantanne s jmax_ins print On stoke le fichier bilKP instantanne s s cym cym nid_bilKPins ENDIF c cIM BEG c cIM cf AM BEG region cym CALL histbeg("histbilKP_ins", iim, zx_lon(:, 1), cym.jjmp1, zx_lat(1,:), cym.imin_ins, imax_ins-imin_ins+1, cym.jmin_ins, jmax_ins-jmin_ins+1, cym.itau_phy, zjulian, dtime, cym.nhori, nid_bilKPins) CALL histbeg_phy("histbilKP_ins"
zx_lon
c c zjulian c cym CALL iim cym zx_lon(i, 1)
dtime
c c zjulian c cym CALL iim cym klev cym dtime
Definition: ini_bilKP_ave.h:26
rlon
c c $Id c nbregdyn DO klon c rlon(i)
vertical_layers_mod::presnivs
real, dimension(:), allocatable, save presnivs
Definition: vertical_layers_mod.F90:11
annee_ref
!$Id annee_ref
Definition: temps.h:15