LMDZ
write_bilKP_ins.h
Go to the documentation of this file.
1  c
2 c $Id$
3 c
4  IF (ok_journe) THEN
5 c
6  ndex2d = 0
7  ndex3d = 0
8 c
9  itau_w = itau_phy + itap + start_time * day_step / iphysiq
10 c
11 c Champs 3D:
12 c
13 cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, ue_lay,zx_tmp_3d)
14  CALL histwrite_phy(nid_bilKPins,"ue",itau_w,ue_lay)
15 c
16 cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, ve_lay,zx_tmp_3d)
17  CALL histwrite_phy(nid_bilKPins,"ve",itau_w,ve_lay)
18 c
19 cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, uq_lay,zx_tmp_3d)
20  CALL histwrite_phy(nid_bilKPins,"uq",itau_w,uq_lay)
21 c
22 cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, vq_lay,zx_tmp_3d)
23  CALL histwrite_phy(nid_bilKPins,"vq",itau_w,vq_lay)
24 c
25 c Champs 3D:
26 C
27 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
28  CALL histwrite_phy(nid_bilKPins,"temp",itau_w,t_seri)
29 c
30 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
31  CALL histwrite_phy(nid_bilKPins,"ovap",itau_w,qx(:,:,ivap))
32 c
33 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
34  CALL histwrite_phy(nid_bilKPins,"geop",itau_w,zphi)
35 c
36 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
37  CALL histwrite_phy(nid_bilKPins,"vitu",itau_w,u_seri)
38 c
39 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
40  CALL histwrite_phy(nid_bilKPins,"vitv",itau_w,v_seri)
41 c
42 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
43  CALL histwrite_phy(nid_bilKPins,"vitw",itau_w,omega)
44 c
45 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
46  CALL histwrite_phy(nid_bilKPins,"pres",itau_w,pplay)
47 c
48 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, paprs, zx_tmp_3d)
49  CALL histwrite_phy(nid_bilKPins,"play",itau_w,paprs)
50 c
51 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)
52  CALL histwrite_phy(nid_bilKPins,"oliq",itau_w,cldliq)
53 c
54 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
55  CALL histwrite_phy(nid_bilKPins,"dtdyn",itau_w,d_t_dyn)
56 c
57 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
58  CALL histwrite_phy(nid_bilKPins,"dqdyn",itau_w,d_q_dyn)
59 c
60 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d)
61  CALL histwrite_phy(nid_bilKPins,"dtcon",itau_w,d_t_con)
62 c
63 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_con, zx_tmp_3d)
64  CALL histwrite_phy(nid_bilKPins,"ducon",itau_w,d_u_con)
65 c
66 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_con, zx_tmp_3d)
67  CALL histwrite_phy(nid_bilKPins,"dvcon",itau_w,d_v_con)
68 c
69 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
70  CALL histwrite_phy(nid_bilKPins,"dqcon",itau_w,d_q_con)
71 c
72 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d)
73  CALL histwrite_phy(nid_bilKPins,"dtlsc",itau_w,d_t_lsc)
74 c
75 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
76  CALL histwrite_phy(nid_bilKPins,"dqlsc",itau_w,d_q_lsc)
77 c
78 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d)
79  CALL histwrite_phy(nid_bilKPins,"dtvdf",itau_w,d_t_vdf)
80 c
81 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
82  CALL histwrite_phy(nid_bilKPins,"dqvdf",itau_w,d_q_vdf)
83 c
84 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d)
85  CALL histwrite_phy(nid_bilKPins,"dtajs",itau_w,d_t_ajs)
86 c
87 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
88  CALL histwrite_phy(nid_bilKPins,"dqajs",itau_w,d_q_ajs)
89 c
90 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d)
91  CALL histwrite_phy(nid_bilKPins,"dteva",itau_w,d_t_eva)
92 c
93 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
94  CALL histwrite_phy(nid_bilKPins,"dqeva",itau_w,d_q_eva)
95 c
96 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d)
97  CALL histwrite_phy(nid_bilKPins,"dtswr",itau_w,heat)
98 c
99 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d)
100  CALL histwrite_phy(nid_bilKPins,"dtsw0",itau_w,heat0)
101 c
102 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d)
103  CALL histwrite_phy(nid_bilKPins,"dtlwr",itau_w,cool)
104 c
105 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d)
106  CALL histwrite_phy(nid_bilKPins,"dtlw0",itau_w,cool0)
107 c
108 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
109  CALL histwrite_phy(nid_bilKPins,"duvdf",itau_w,d_u_vdf)
110 c
111 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
112  CALL histwrite_phy(nid_bilKPins,"dvvdf",itau_w,d_v_vdf)
113 c
114  IF (ok_orodr) THEN
115  IF (ok_orolf) THEN
116 c
117  DO k = 1, klev
118  DO i = 1, klon
119  d_u_oli(i,k) = d_u_oro(i,k) + d_u_lif(i,k)
120  d_v_oli(i,k) = d_v_oro(i,k) + d_v_lif(i,k)
121  ENDDO
122  ENDDO
123 c
124 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oli, zx_tmp_3d)
125  CALL histwrite_phy(nid_bilKPins,"duoli",itau_w,d_u_oli)
126 c
127 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oli, zx_tmp_3d)
128  CALL histwrite_phy(nid_bilKPins,"dvoli",itau_w,d_v_oli)
129 c
130  ENDIF
131  ENDIF
132 C
133 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u, zx_tmp_3d)
134  CALL histwrite_phy(nid_bilKPins,"duphy",itau_w,d_u)
135 c
136 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v, zx_tmp_3d)
137  CALL histwrite_phy(nid_bilKPins,"dvphy",itau_w,d_v)
138 c
139 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
140  CALL histwrite_phy(nid_bilKPins,"dtphy",itau_w,d_t)
141 c
142 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,1),
143 cym .zx_tmp_3d)
144  CALL histwrite_phy(nid_bilKPins,"dqphy",itau_w,d_qx(:,:,1))
145 c
146 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,2),
147 cym .zx_tmp_3d)
148  CALL histwrite_phy(nid_bilKPins,"dqlphy",itau_w,d_qx(:,:,2))
149 c
150 cIM 280405 BEG
151 c
152 c Champs 2D:
153 c
154 c Ecriture de champs dynamiques sur des niveaux de pression
155 c DO k=1, nlevSTD
156  DO k=1, 12
157 c
158  IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
159  IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
160 c
161  IF(bb2.EQ."850") THEN
162 c
163 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,1),zx_tmp_2d)
164  CALL histwrite_phy(nid_bilKPins,"u"//bb2,itau_w,usumSTD(:,k,1))
165 c
166 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,1),zx_tmp_2d)
167  CALL histwrite_phy(nid_bilKPins,"v"//bb2,itau_w,vsumSTD(:,k,1))
168 c
169  ENDIF !(bb2.EQ."850")
170 c
171  ENDDO !k=1, 12
172 c
173 cIM 280405 END
174 C
175  if (ok_sync) then
176  call histsync(nid_bilKPins)
177  endif
178  ENDIF
179 
ok_orolf
!$Header!c include clesph0 h c COMMON clesph0 ok_orolf
Definition: clesph0.h:6
if
!$Id Turb_fcg_gcssold if(prt_level.ge.1) then print *
phys_local_var_mod::d_t_con
real, dimension(:,:), allocatable, save d_t_con
Definition: phys_local_var_mod.F90:29
klon
INTERFACE SUBROUTINE RRTM_ECRT_140GP klon
Definition: rrtm_ecrt_140gp.intfb.h:3
$Id
c c $Id
Definition: ini_bilKP_ave.h:11
phys_local_var_mod::d_u_oro
real, dimension(:,:), allocatable, save d_u_oro
Definition: phys_local_var_mod.F90:76
phys_state_var_mod::heat0
real, dimension(:,:), allocatable, save heat0
Definition: phys_state_var_mod.F90:302
ok_orodr
!$Header!c include clesph0 h c COMMON clesph0 ok_orodr
Definition: clesph0.h:6
phys_local_var_mod::d_q_lsc
real, dimension(:,:), allocatable, save d_q_lsc
Definition: phys_local_var_mod.F90:35
phys_state_var_mod::heat
real, dimension(:,:), allocatable, save heat
Definition: phys_state_var_mod.F90:300
phys_local_var_mod::d_v_oli
real, dimension(:,:), allocatable, save d_v_oli
Definition: phys_local_var_mod.F90:62
itap
!$Id klon initialisation mois suivants day_rain itap
Definition: calcul_divers.h:18
phys_local_var_mod::d_t_ajs
real, dimension(:,:), allocatable, save d_t_ajs
Definition: phys_local_var_mod.F90:41
phys_local_var_mod::zphi
real, dimension(:,:), allocatable, save zphi
Definition: phys_local_var_mod.F90:349
phys_local_var_mod::d_t_dyn
real, dimension(:,:), allocatable, save d_t_dyn
Definition: phys_local_var_mod.F90:21
phys_state_var_mod::cool
real, dimension(:,:), allocatable, save cool
Definition: phys_state_var_mod.F90:304
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
phys_local_var_mod::ve
real, dimension(:), allocatable, save ve
Definition: phys_local_var_mod.F90:273
control_mod::day_step
integer, save day_step
Definition: control_mod.F90:15
control_mod::iphysiq
integer, save iphysiq
Definition: control_mod.F90:24
itau_phy
!$Id itau_phy
Definition: temps.h:15
IF
c c $Id c IF(ok_journe) THEN c ndex2d=0ndex3d=0c itau_w
k
c c $Id c klev DO klon nlevSTD DO k
Definition: write_bilKP_ins.h:156
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
phys_local_var_mod::d_v_con
real, dimension(:,:), allocatable, save d_v_con
Definition: phys_local_var_mod.F90:31
phys_local_var_mod::d_q_ajs
real, dimension(:,:), allocatable, save d_q_ajs
Definition: phys_local_var_mod.F90:41
phys_state_var_mod::cool0
real, dimension(:,:), allocatable, save cool0
Definition: phys_state_var_mod.F90:306
pression
subroutine pression(ngrid, ap, bp, ps, p)
Definition: pression.F90:2
phys_local_var_mod::d_u_lif
real, dimension(:,:), allocatable, save d_u_lif
Definition: phys_local_var_mod.F90:80
d_u_oli
c c $Id c klev DO klon d_u_oli(i, k)
phys_local_var_mod::d_q_dyn
real, dimension(:,:), allocatable, save d_q_dyn
Definition: phys_local_var_mod.F90:21
cym
c c zjulian c cym CALL iim cym klev cym cym
Definition: ini_bilKP_ave.h:26
zx_tmp_2d
!FH On elimine toutes les clefs physiques dans la dynamique zx_tmp_2d
Definition: write_paramLMDZ_dyn.h:9
phys_local_var_mod::d_u_con
real, dimension(:,:), allocatable, save d_u_con
Definition: phys_local_var_mod.F90:31
phys_local_var_mod::uq
real, dimension(:), allocatable, save uq
Definition: phys_local_var_mod.F90:273
pplay
!$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 pplay
Definition: calcul_STDlev.h:26
paprs
!$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 ulevSTD CALL &zphi philevSTD CALL &zx_rh rhlevSTD!DO klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon du jour ou toutes les read_climoz CALL &zphi geo500!IM on interpole a chaque pas de temps le paprs
Definition: calcul_STDlev.h:144
phys_local_var_mod::d_v_lif
real, dimension(:,:), allocatable, save d_v_lif
Definition: phys_local_var_mod.F90:80
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)
jjmp1
c c zjulian c cym CALL iim cym klev jjmp1
Definition: ini_bilKP_ave.h:24
phys_local_var_mod::d_q_con
real, dimension(:,:), allocatable, save d_q_con
Definition: phys_local_var_mod.F90:29
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
phys_local_var_mod::d_v_oro
real, dimension(:,:), allocatable, save d_v_oro
Definition: phys_local_var_mod.F90:76
phys_local_var_mod::ue
real, dimension(:), allocatable, save ue
Definition: phys_local_var_mod.F90:273
phys_local_var_mod::d_q_vdf
real, dimension(:,:), allocatable, save d_q_vdf
Definition: phys_local_var_mod.F90:64
phys_local_var_mod::d_t_eva
real, dimension(:,:), allocatable, save d_t_eva
nrlmd
Definition: phys_local_var_mod.F90:51
u_seri
!$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 & u_seri
Definition: calcul_STDlev.h:26
iim
c c zjulian c cym CALL iim cym klev iim
Definition: ini_bilKP_ave.h:24
phys_local_var_mod::d_u_vdf
real, dimension(:,:), allocatable, save d_u_vdf
Definition: phys_local_var_mod.F90:66
gr_fi_ecrit
subroutine gr_fi_ecrit(nfield, nlon, iim, jjmp1, fi, ecrit)
Definition: gr_fi_ecrit.F90:5
phys_local_var_mod::t_seri
real, dimension(:,:), allocatable, save t_seri
Definition: phys_local_var_mod.F90:12
phys_local_var_mod::d_v_vdf
real, dimension(:,:), allocatable, save d_v_vdf
Definition: phys_local_var_mod.F90:66
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
phys_local_var_mod::d_t_vdf
real, dimension(:,:), allocatable, save d_t_vdf
Definition: phys_local_var_mod.F90:64
de
do llm!au dessus de
Definition: 1D_nudge_sandu_astex.h:1
histwrite_phy
!$Id ***************************************!ECRITURE DU pphis CALL histwrite_phy(nid_tra,.FALSE.,"aire", itau_w, airephy) CALL histwrite_phy(nid_tra
phys_local_var_mod::d_q_eva
real, dimension(:,:), allocatable, save d_q_eva
Definition: phys_local_var_mod.F90:51
itau_w
!FH On elimine toutes les clefs physiques dans la dynamique itau_w
Definition: write_paramLMDZ_dyn.h:9
phys_local_var_mod::omega
real, dimension(:,:), allocatable, save omega
Definition: phys_local_var_mod.F90:342
start_time
!$Id start_time
Definition: temps.h:15
i
c c $Id c klev DO i
Definition: write_bilKP_ins.h:118
phys_local_var_mod::v_seri
real, dimension(:,:), allocatable, save v_seri
Definition: phys_local_var_mod.F90:16
phys_local_var_mod::d_t_lsc
real, dimension(:,:), allocatable, save d_t_lsc
Definition: phys_local_var_mod.F90:35