LMDZ
write_paramLMDZ_dyn.h
Go to the documentation of this file.
1 ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
2 ! Attention : il n'y a aucune raison pour ecrire ces constantes
3 ! comme des champs 2D. A corriger un jour ...
4 
5 c
6  ndex2d = 0
7  itau_w=itau_dyn+itau
8 c
9  zx_tmp_2d(1:iip1,1:jjp1)=REAL(prt_level)
10  CALL histwrite(nid_ctesGCM, "prt_level", itau_w,
11  . zx_tmp_2d,iip1*jjp1,ndex2d)
12 c
13  zx_tmp_2d(1:iip1,1:jjp1)=REAL(dayref)
14  CALL histwrite(nid_ctesGCM, "dayref", itau_w,
15  . zx_tmp_2d,iip1*jjp1,ndex2d)
16 c
17  zx_tmp_2d(1:iip1,1:jjp1)=REAL(anneeref)
18  CALL histwrite(nid_ctesGCM, "anneeref", itau_w,
19  . zx_tmp_2d,iip1*jjp1,ndex2d)
20 c
21  zx_tmp_2d(1:iip1,1:jjp1)=REAL(raz_date)
22  CALL histwrite(nid_ctesGCM, "raz_date", itau_w,
23  . zx_tmp_2d,iip1*jjp1,ndex2d)
24 c
25  zx_tmp_2d(1:iip1,1:jjp1)=REAL(nday)
26  CALL histwrite(nid_ctesGCM, "nday", itau_w,
27  . zx_tmp_2d,iip1*jjp1,ndex2d)
28 c
29  zx_tmp_2d(1:iip1,1:jjp1)=REAL(day_step)
30  CALL histwrite(nid_ctesGCM, "day_step", itau_w,
31  . zx_tmp_2d,iip1*jjp1,ndex2d)
32 c
33  zx_tmp_2d(1:iip1,1:jjp1)=REAL(iperiod)
34  CALL histwrite(nid_ctesGCM, "iperiod", itau_w,
35  . zx_tmp_2d,iip1*jjp1,ndex2d)
36 c
37  zx_tmp_2d(1:iip1,1:jjp1)=REAL(iapp_tracvl)
38  CALL histwrite(nid_ctesGCM, "iapp_tracvl", itau_w,
39  . zx_tmp_2d,iip1*jjp1,ndex2d)
40 c
41  zx_tmp_2d(1:iip1,1:jjp1)=REAL(iconser)
42  CALL histwrite(nid_ctesGCM, "iconser", itau_w,
43  . zx_tmp_2d,iip1*jjp1,ndex2d)
44 c
45  zx_tmp_2d(1:iip1,1:jjp1)=REAL(iecri)
46  CALL histwrite(nid_ctesGCM, "iecri", itau_w,
47  . zx_tmp_2d,iip1*jjp1,ndex2d)
48 c
49  zx_tmp_2d(1:iip1,1:jjp1)=periodav
50  CALL histwrite(nid_ctesGCM, "periodav", itau_w,
51  . zx_tmp_2d,iip1*jjp1,ndex2d)
52 c
53  zx_tmp_2d(1:iip1,1:jjp1)=REAL(dissip_period)
54  CALL histwrite(nid_ctesGCM, "dissip_period", itau_w,
55  . zx_tmp_2d,iip1*jjp1,ndex2d)
56 c
57  IF(lstardis) THEN
58  zx_tmp_2d(1:iip1,1:jjp1)=1.
59  ELSE
60  zx_tmp_2d(1:iip1,1:jjp1)=0.
61  ENDIF
62  CALL histwrite(nid_ctesGCM, "lstardis", itau_w,
63  . zx_tmp_2d,iip1*jjp1,ndex2d)
64 c
65  zx_tmp_2d(1:iip1,1:jjp1)=REAL(nitergdiv)
66  CALL histwrite(nid_ctesGCM, "nitergdiv", itau_w,
67  . zx_tmp_2d,iip1*jjp1,ndex2d)
68 c
69  zx_tmp_2d(1:iip1,1:jjp1)=REAL(nitergrot)
70  CALL histwrite(nid_ctesGCM, "nitergrot", itau_w,
71  . zx_tmp_2d,iip1*jjp1,ndex2d)
72 c
73  zx_tmp_2d(1:iip1,1:jjp1)=REAL(niterh)
74  CALL histwrite(nid_ctesGCM, "niterh", itau_w,
75  . zx_tmp_2d,iip1*jjp1,ndex2d)
76 c
77  zx_tmp_2d(1:iip1,1:jjp1)=tetagdiv
78  CALL histwrite(nid_ctesGCM, "tetagdiv", itau_w,
79  . zx_tmp_2d,iip1*jjp1,ndex2d)
80 c
81  zx_tmp_2d(1:iip1,1:jjp1)=tetagrot
82  CALL histwrite(nid_ctesGCM, "tetagrot", itau_w,
83  . zx_tmp_2d,iip1*jjp1,ndex2d)
84 c
85  zx_tmp_2d(1:iip1,1:jjp1)=tetatemp
86  CALL histwrite(nid_ctesGCM, "tetatemp", itau_w,
87  . zx_tmp_2d,iip1*jjp1,ndex2d)
88 c
89  zx_tmp_2d(1:iip1,1:jjp1)=coefdis
90  CALL histwrite(nid_ctesGCM, "coefdis", itau_w,
91  . zx_tmp_2d,iip1*jjp1,ndex2d)
92 c
93  IF(purmats) THEN
94  zx_tmp_2d(1:iip1,1:jjp1)=1.
95  ELSE
96  zx_tmp_2d(1:iip1,1:jjp1)=0.
97  ENDIF
98  CALL histwrite(nid_ctesGCM, "purmats", itau_w,
99  . zx_tmp_2d,iip1*jjp1,ndex2d)
100 c
101  IF(ok_guide) THEN
102  zx_tmp_2d(1:iip1,1:jjp1)=1.
103  ELSE
104  zx_tmp_2d(1:iip1,1:jjp1)=0.
105  ENDIF
106  CALL histwrite(nid_ctesGCM, "ok_guide", itau_w,
107  . zx_tmp_2d,iip1*jjp1,ndex2d)
108 c
109  if (calend == 'earth_360d') then
110  zx_tmp_2d(1:iip1,1:jjp1)=1.
111  else if (calend == 'earth_365d') then
112  zx_tmp_2d(1:iip1,1:jjp1)=2.
113  else if (calend == 'earth_366d') then
114  zx_tmp_2d(1:iip1,1:jjp1)=3.
115  endif
116 
117  CALL histwrite(nid_ctesGCM, "true_calendar", itau_w,
118  . zx_tmp_2d,iip1*jjp1,ndex2d)
119 c
120  zx_tmp_2d(1:iip1,1:jjp1)=REAL(iflag_phys)
121  CALL histwrite(nid_ctesGCM, "iflag_phys", itau_w,
122  . zx_tmp_2d,iip1*jjp1,ndex2d)
123 c
124  zx_tmp_2d(1:iip1,1:jjp1)=REAL(iphysiq)
125  CALL histwrite(nid_ctesGCM, "iphysiq", itau_w,
126  . zx_tmp_2d,iip1*jjp1,ndex2d)
127 c
128 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
129 ! FH 2008/05/02
130 ! La variable cycle_diurne n'est pas vue par la dynamique
131 ! IF(cycle_diurne) THEN
132 ! zx_tmp_2d(1:iip1,1:jjp1)=1.
133 ! ELSE
134 ! zx_tmp_2d(1:iip1,1:jjp1)=0.
135 ! ENDIF
136 ! CALL histwrite(nid_ctesGCM, "cycle_diurne", itau_w,
137 ! . zx_tmp_2d,iip1*jjp1,ndex2d)
138 !
139 ! IF(soil_model) THEN
140 ! zx_tmp_2d(1:iip1,1:jjp1)=1.
141 ! ELSE
142 ! zx_tmp_2d(1:iip1,1:jjp1)=0.
143 ! ENDIF
144 ! CALL histwrite(nid_ctesGCM, "soil_model", itau_w,
145 ! . zx_tmp_2d,iip1*jjp1,ndex2d)
146 !
147 ! IF(new_oliq) THEN
148 ! zx_tmp_2d(1:iip1,1:jjp1)=1.
149 ! ELSE
150 ! zx_tmp_2d(1:iip1,1:jjp1)=0.
151 ! ENDIF
152 ! CALL histwrite(nid_ctesGCM, "new_oliq", itau_w,
153 ! . zx_tmp_2d,iip1*jjp1,ndex2d)
154 !
155 ! IF(ok_orodr) THEN
156 ! zx_tmp_2d(1:iip1,1:jjp1)=1.
157 ! ELSE
158 ! zx_tmp_2d(1:iip1,1:jjp1)=0.
159 ! ENDIF
160 ! CALL histwrite(nid_ctesGCM, "ok_orodr", itau_w,
161 ! . zx_tmp_2d,iip1*jjp1,ndex2d)
162 !
163 ! IF(ok_orolf) THEN
164 ! zx_tmp_2d(1:iip1,1:jjp1)=1.
165 ! ELSE
166 ! zx_tmp_2d(1:iip1,1:jjp1)=0.
167 ! ENDIF
168 ! CALL histwrite(nid_ctesGCM, "ok_orolf", itau_w,
169 ! . zx_tmp_2d,iip1*jjp1,ndex2d)
170 !
171 ! IF(ok_limitvrai) THEN
172 ! zx_tmp_2d(1:iip1,1:jjp1)=1.
173 ! ELSE
174 ! zx_tmp_2d(1:iip1,1:jjp1)=0.
175 ! ENDIF
176 ! CALL histwrite(nid_ctesGCM, "ok_limitvrai", itau_w,
177 ! . zx_tmp_2d,iip1*jjp1,ndex2d)
178 !
179 ! zx_tmp_2d(1:iip1,1:jjp1)=nbapp_rad
180 ! CALL histwrite(nid_ctesGCM, "nbapp_rad", itau_w,
181 ! . zx_tmp_2d,iip1*jjp1,ndex2d)
182 !
183 ! zx_tmp_2d(1:iip1,1:jjp1)=iflag_con
184 ! CALL histwrite(nid_ctesGCM, "iflag_con", itau_w,
185 ! . zx_tmp_2d,iip1*jjp1,ndex2d)
186 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
187 c
188  zx_tmp_2d(1:iip1,1:jjp1)=clon
189  CALL histwrite(nid_ctesGCM, "clon", itau_w,
190  . zx_tmp_2d,iip1*jjp1,ndex2d)
191 c
192  zx_tmp_2d(1:iip1,1:jjp1)=clat
193  CALL histwrite(nid_ctesGCM, "clat", itau_w,
194  . zx_tmp_2d,iip1*jjp1,ndex2d)
195 c
196  zx_tmp_2d(1:iip1,1:jjp1)=grossismx
197  CALL histwrite(nid_ctesGCM, "grossismx", itau_w,
198  . zx_tmp_2d,iip1*jjp1,ndex2d)
199 c
200  zx_tmp_2d(1:iip1,1:jjp1)=grossismy
201  CALL histwrite(nid_ctesGCM, "grossismy", itau_w,
202  . zx_tmp_2d,iip1*jjp1,ndex2d)
203 c
204  IF(fxyhypb) THEN
205  zx_tmp_2d(1:iip1,1:jjp1)=1.
206  ELSE
207  zx_tmp_2d(1:iip1,1:jjp1)=0.
208  ENDIF
209  CALL histwrite(nid_ctesGCM, "fxyhypb", itau_w,
210  . zx_tmp_2d,iip1*jjp1,ndex2d)
211 c
212  zx_tmp_2d(1:iip1,1:jjp1)=dzoomx
213  CALL histwrite(nid_ctesGCM, "dzoomx", itau_w,
214  . zx_tmp_2d,iip1*jjp1,ndex2d)
215 c
216  zx_tmp_2d(1:iip1,1:jjp1)=dzoomy
217  CALL histwrite(nid_ctesGCM, "dzoomy", itau_w,
218  . zx_tmp_2d,iip1*jjp1,ndex2d)
219 c
220  zx_tmp_2d(1:iip1,1:jjp1)=taux
221  CALL histwrite(nid_ctesGCM, "taux", itau_w,
222  . zx_tmp_2d,iip1*jjp1,ndex2d)
223 c
224  zx_tmp_2d(1:iip1,1:jjp1)=tauy
225  CALL histwrite(nid_ctesGCM, "tauy", itau_w,
226  . zx_tmp_2d,iip1*jjp1,ndex2d)
227 c
228  IF(ysinus) THEN
229  zx_tmp_2d(1:iip1,1:jjp1)=1.
230  ELSE
231  zx_tmp_2d(1:iip1,1:jjp1)=0.
232  ENDIF
233  CALL histwrite(nid_ctesGCM, "ysinus", itau_w,
234  . zx_tmp_2d,iip1*jjp1,ndex2d)
235 c
236  zx_tmp_2d(1:iip1,1:jjp1)=ip_ebil_dyn
237  CALL histwrite(nid_ctesGCM, "ip_ebil_dyn", itau_w,
238  . zx_tmp_2d,iip1*jjp1,ndex2d)
239 c
240 c=================================================================
241 c
242  if (ok_sync) then
243  call histsync(nid_ctesGCM)
244  endif
245 c
246 c=================================================================
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 *
dzoomy
!$Header!c!c!c include serre h!c REAL dzoomy
Definition: serre.h:8
grossismx
!$Header!c!c!c include serre h!c REAL && grossismx
Definition: serre.h:8
soil_model
!$Header!c include clesph0 h c COMMON clesph0 soil_model
Definition: clesph0.h:6
ok_orodr
!$Header!c include clesph0 h c COMMON clesph0 ok_orodr
Definition: clesph0.h:6
iflag_con
!$Id ok_orolf LOGICAL ok_limitvrai LOGICAL ok_all_xml INTEGER iflag_con
Definition: clesphys.h:12
cycle_diurne
!$Header!c include clesph0 h c COMMON clesph0 cycle_diurne
Definition: clesph0.h:6
ok_limitvrai
!$Header!c include clesph0 h c COMMON clesph0 ok_limitvrai
Definition: clesph0.h:6
clon
!$Header!c!c!c include serre h!c REAL clon
Definition: serre.h:8
new_oliq
!$Header!c include clesph0 h c COMMON clesph0 new_oliq
Definition: clesph0.h:6
nbapp_rad
!$Header!c include clesph0 h c COMMON clesph0 nbapp_rad
Definition: clesph0.h:6
IF
!$Id klon IF(pctsrf(i, is_ter).GT.0.) THEN paire_ter(i)
zx_tmp_2d
!FH On elimine toutes les clefs physiques dans la dynamique zx_tmp_2d
Definition: write_paramLMDZ_dyn.h:9
fxyhypb
!$Id fxyhypb
Definition: logic.h:10
dzoomx
!$Header!c!c!c include serre h!c REAL dzoomx
Definition: serre.h:8
grossismy
!$Header!c!c!c include serre h!c REAL grossismy
Definition: serre.h:8
taux
!$Header!c!c!c include serre h!c REAL taux
Definition: serre.h:8
control_mod::ip_ebil_dyn
integer, save ip_ebil_dyn
Definition: control_mod.F90:29
clat
!$Header!c!c!c include serre h!c REAL clat
Definition: serre.h:8
FH
!$Id!Parameters for nlm real spfac!IM cf FH
Definition: cv30param.h:5
itau_w
!FH On elimine toutes les clefs physiques dans la dynamique itau_w
Definition: write_paramLMDZ_dyn.h:9
jjp1
!FH On elimine toutes les clefs physiques dans la dynamique iip1 * jjp1
Definition: write_paramLMDZ_dyn.h:9