LMDZ
ini_bilKP_ins.h
Go to the documentation of this file.
1 c
2 c $Id: ini_bilKP_ins.h 1907 2013-11-26 13:10:46Z lguez $
3 c
4  IF (ok_journe) THEN
5 c
6  zsto = dtime
7  zout = dtime
8  typeval=tinst
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 c
25 cIM 280405 BEG
26 c
27 cIM cf. AM 081204 BEG region
28  imin_ins=1
29  imax_ins=iim
30  jmin_ins=1
31  jmax_ins=jjmp1
32 cym do i=1,iim-1
33 cym if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i
34 cym if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1
35 cym enddo
36 cym do j=1,jjmp1
37 cym if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j
38 cym if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j
39 cym enddo
40 c
41  print*,'On stoke le fichier bilKP instantanne sur ',
42  s imin_ins,imax_ins,jmin_ins,jmax_ins
43  print*,'On stoke le fichier bilKP instantanne sur ',
44  s zx_lon(imin_ins,1),zx_lon(imax_ins,1),
45  s zx_lat(1,jmin_ins),zx_lat(1,jmax_ins)
46 cIM cf. AM 081204 END region
47 c
48 cIM 280405 END
49 c
50 cym IF(1.EQ.0) THEN
51 cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1,
52 cym . zx_lat(1,:),
53 cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
54 cym . nhori, nid_bilKPins)
55  ENDIF
56 c
57 cIM 280405 BEG
58 c
59 cIM cf. AM 081204 BEG region
60 cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1),
61 cym . jjmp1,zx_lat(1,:),
62 cym . imin_ins,imax_ins-imin_ins+1,
63 cym . jmin_ins,jmax_ins-jmin_ins+1,
64 cym . itau_phy, zjulian, dtime,
65 cym . nhori, nid_bilKPins)
66  CALL histbeg_phy("histbilKP_ins", itau_phy, zjulian, dtime,
67  . nhori, nid_bilKPins)
68 cIM 081204 END
69 c
70 cIM 280405 END
71 c
72  write(*,*)'Journee ', itau_phy, zjulian
73  CALL histvert(nid_bilKPins, "presnivs",
74  . "Vertical levels","mb",
75  . klev, presnivs/100., nvert)
76 c
77 c Champs 3D:
78 c
79  CALL histdef(nid_bilKPins,"ue",
80  . "Zonal energy transport","-",
81  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
82  . typeval, zsto,zout)
83 c
84  CALL histdef(nid_bilKPins,"ve",
85  . "Merid energy transport","-",
86  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
87  . typeval, zsto,zout)
88 c
89  CALL histdef(nid_bilKPins,"uq",
90  . "Zonal humidity transport","-",
91  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
92  . typeval, zsto,zout)
93 c
94  CALL histdef(nid_bilKPins,"vq",
95  . "Merid humidity transport","-",
96  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
97  . typeval, zsto,zout)
98 c
99 c Champs 3D:
100 c
101  CALL histdef(nid_bilKPins, "temp",
102  . "Air temperature", "K",
103  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
104  . typeval, zsto,zout)
105 c
106  CALL histdef(nid_bilKPins,"ovap",
107  . "Specific humidity","Kg/Kg",
108  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
109  . typeval, zsto,zout)
110 c
111  CALL histdef(nid_bilKPins,"geop",
112  . "Geopotential height", "m",
113  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
114  . typeval, zsto,zout)
115 c
116  CALL histdef(nid_bilKPins,"vitu",
117  . "Zonal wind", "m/s",
118  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
119  . typeval, zsto,zout)
120 c
121  CALL histdef(nid_bilKPins,"vitv",
122  . "Meridional wind", "m/s",
123  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
124  . typeval, zsto,zout)
125 c
126  CALL histdef(nid_bilKPins, "vitw",
127  . "Vertical wind", "m/s",
128  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
129  . typeval, zsto,zout)
130 c
131  CALL histdef(nid_bilKPins, "pres",
132  . "Inter-Layer Air pressure",
133  . "Pa",
134  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
135  . typeval, zsto,zout)
136 c
137  CALL histdef(nid_bilKPins, "play",
138  . "Mean-Layer Air pressure",
139  . "Pa",
140  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
141  . typeval, zsto,zout)
142 c
143  CALL histdef(nid_bilKPins, "oliq",
144  . "Liquid water content",
145  . "kg/kg",
146  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
147  . typeval, zsto,zout)
148 c
149  CALL histdef(nid_bilKPins, "dtdyn",
150  . "Dynamics dT", "K/s",
151  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
152  . typeval, zsto,zout)
153 c
154  CALL histdef(nid_bilKPins, "dqdyn",
155  . "Dynamics dQ", "Kg/Kg/s",
156  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
157  . typeval, zsto,zout)
158 c
159  CALL histdef(nid_bilKPins, "dtcon",
160  . "Convection dT", "K/s",
161  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
162  . typeval, zsto,zout)
163 c
164  CALL histdef(nid_bilKPins, "ducon",
165  . "Convection du", "m/s2",
166  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
167  . typeval, zsto,zout)
168 c
169  CALL histdef(nid_bilKPins, "dvcon",
170  . "Convection dv", "m/s2",
171  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
172  . typeval, zsto,zout)
173 c
174  CALL histdef(nid_bilKPins,"dqcon",
175  . "Convection dQ","Kg/Kg/s",
176  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
177  . typeval, zsto,zout)
178 c
179  CALL histdef(nid_bilKPins, "dtlsc",
180  . "Condensation dT", "K/s",
181  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
182  . typeval, zsto,zout)
183 c
184  CALL histdef(nid_bilKPins,"dqlsc",
185  . "Condensation dQ","Kg/Kg/s",
186  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
187  . typeval, zsto,zout)
188 c
189  CALL histdef(nid_bilKPins,"dtvdf",
190  . "Boundary-layer dT","K/s",
191  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
192  . typeval, zsto,zout)
193 c
194  CALL histdef(nid_bilKPins, "dqvdf",
195  . "Boundary-layer dQ",
196  . "Kg/Kg/s",
197  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
198  . typeval, zsto,zout)
199 c
200  CALL histdef(nid_bilKPins,"dtajs",
201  . "Ajustement sec dT","K/s",
202  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
203  . typeval, zsto,zout)
204 c
205  CALL histdef(nid_bilKPins,"dqajs",
206  . "Ajustement sec dQ",
207  . "Kg/Kg/s",
208  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
209  . typeval, zsto,zout)
210 c
211  CALL histdef(nid_bilKPins,"dteva",
212  . "Reevaporation dT","K/s",
213  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
214  . typeval, zsto,zout)
215 c
216  CALL histdef(nid_bilKPins,"dqeva",
217  . "Reevaporation dQ",
218  . "Kg/Kg/s",
219  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
220  . typeval, zsto,zout)
221 
222 c
223  CALL histdef(nid_bilKPins, "dtswr",
224  . "SW radiation dT", "K/s",
225  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
226  . typeval, zsto,zout)
227 c
228  CALL histdef(nid_bilKPins, "dtsw0",
229  . "SW radiation dT", "K/s",
230  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
231  . typeval, zsto,zout)
232 c
233  CALL histdef(nid_bilKPins, "dtlwr",
234  . "LW radiation dT", "K/s",
235  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
236  . typeval, zsto,zout)
237 c
238  CALL histdef(nid_bilKPins, "dtlw0",
239  . "LW radiation dT", "K/s",
240  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
241  . typeval, zsto,zout)
242 c
243  CALL histdef(nid_bilKPins,"duvdf",
244  . "Boundary-layer dU","m/s2",
245  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
246  . typeval, zsto,zout)
247 c
248  CALL histdef(nid_bilKPins,"dvvdf",
249  . "Boundary-layer dV","m/s2",
250  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
251  . typeval, zsto,zout)
252 c
253  IF (ok_orodr) THEN
254  IF (ok_orolf) THEN
255  CALL histdef(nid_bilKPins, "duoli",
256  . "Orography dU", "m/s2",
257  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
258  . typeval, zsto,zout)
259 c
260  CALL histdef(nid_bilKPins, "dvoli",
261  . "Orography dV", "m/s2",
262  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
263  . typeval, zsto,zout)
264 c
265  ENDIF
266  ENDIF
267 C
268  CALL histdef(nid_bilKPins, "duphy",
269  . "Physiq dU","-",
270  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
271  . typeval, zsto,zout)
272 C
273  CALL histdef(nid_bilKPins, "dvphy",
274  . "Physiq dV","-",
275  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
276  . typeval, zsto,zout)
277 C
278  CALL histdef(nid_bilKPins, "dtphy",
279  . "Physiq dT","-",
280  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
281  . typeval, zsto,zout)
282 C
283  CALL histdef(nid_bilKPins, "dqphy",
284  . "Physiq dQ","-",
285  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
286  . typeval, zsto,zout)
287 C
288  CALL histdef(nid_bilKPins, "dqlphy",
289  . "Physiq dQl","-",
290  . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
291  . typeval, zsto,zout)
292 cIM 280405 BEG
293 c
294 c Champs 2D:
295 c
296 c u850, v850
297 c DO k=1, nlevSTD
298  DO k=1, 12
299 c
300  IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
301  IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
302 c
303  IF(bb2.EQ."850") THEN
304 c
305  CALL histdef(nid_bilKPins, "u"//bb2,
306  . "Zonal wind "//bb2//"mb","m/s",
307  . iim,jjphy_nb,nhori, 1,1,1, -99, 32,
308  . typeval, zsto,zout)
309 c
310  CALL histdef(nid_bilKPins, "v"//bb2,
311  . "Meridional wind "//bb2//"mb","m/s",
312  . iim,jjphy_nb,nhori, 1,1,1, -99, 32,
313  . typeval, zsto,zout)
314 c
315  ENDIF !(bb2.EQ."850")
316 c
317  ENDDO !k=1, 12
318 c
319 cIM 280405 END
320 c
321  CALL histend(nid_bilKPins)
322 c
323  ndex2d = 0
324  ndex3d = 0
325 c
326  ENDIF ! fin de test sur ok_journe
idayref
c c idayref
Definition: ini_bilKP_ins.h:11
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
zjulian
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 zjulian
Definition: ini_bilKP_ins.h:53
dtime
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 dtime
Definition: ini_bilKP_ins.h:53
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
znivsig
c c zjulian c cym CALL iim cym klev znivsig(ll)
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
imin_ins
c c zjulian c cym CALL iim cym klev iim cym jjmp1 cym On stoke le fichier bilKP instantanne s imin_ins
Definition: ini_bilKP_ins.h:41
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
jjmp1
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 jjmp1
Definition: ini_bilKP_ins.h:51
fichier
!$Header!integer nvarmx s s s fichier
Definition: gradsdef.h:20
k
c c $Id c c calculs statistiques distribution nuage ftion du regime dynamique c c Ce calcul doit etre fait a partir de valeurs mensuelles CALL nbregdyn DO k
Definition: calcul_REGDYN.h:12
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
nhori
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 nhori
Definition: ini_bilKP_ins.h:53
nlevSTD
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL nlevSTD
Definition: calcul_STDlev.h:26
imax_ins
c c zjulian c cym CALL iim cym klev iim cym jjmp1 cym On stoke le fichier bilKP instantanne s imax_ins
Definition: ini_bilKP_ins.h:41
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)
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
jmin_ins
c c zjulian c cym CALL iim cym klev iim cym jjmp1 cym On stoke le fichier bilKP instantanne s jmin_ins
Definition: ini_bilKP_ins.h:41
j
!$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 j
Definition: 1Dconv.h:27
cym
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
Definition: ini_bilKP_ins.h:53
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
u
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm u(l)
gr_fi_ecrit
c c zjulian c cym CALL gr_fi_ecrit(1, klon, iim, jjmp1, rlon, zx_lon) cym DO i
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
zx_lon
c c zjulian c cym CALL iim cym zx_lon(i, 1)
zx_lat
c c zjulian c cym CALL iim cym klev * zx_lat
Definition: ini_bilKP_ins.h:23
$Id
c c $Id
Definition: ini_bilKP_ins.h:11
iim
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 iim
Definition: ini_bilKP_ins.h:51
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
if
c c zjulian c cym CALL iim cym klev iim cym if(zx_lon(i, 1).lt.lonmin_ins) imin_ins
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"
itau_phy
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 itau_phy
Definition: ini_bilKP_ins.h:53
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