Context Navigation

iniwrite.F @ 1255

Last change on this file since 1255 was 1216, checked in by emillour, 11 years ago

Generic model:
Major cleanup, in order to ease the use of LMDZ.GENERIC with (parallel) dynamics
in LMDZ.COMMON: (NB: this will break LMDZ.UNIVERSAL, which should be thrashed
in the near future)

Updated makegcm_* scripts (and makdim) and added the "-full" (to enforce full recomputation of the model) option
In dyn3d: converted control.h to module control_mod.F90 and converted iniadvtrac.F to module infotrac.F90
Added module mod_const_mpi.F90 in dyn3d (not used in serial mode)
Rearanged input/outputs routines everywhere to handle serial/MPI cases. physdem.F => phyredem.F90 , phyetat0.F => phyetat0.F90 ; all read/write routines for startfi files are gathered in module iostart.F90
added parallelism related routines init_phys_lmdz.F90, comgeomphy.F90, dimphy.F90, iniphysiq.F90, mod_grid_phy_lmdz.F90, mod_phys_lmdz_mpi_data.F90, mod_phys_lmdz_mpi_transfert.F90, mod_phys_lmdz_omp_data.F90, mod_phys_lmdz_omp_transfert.F90, mod_phys_lmdz_para.F90, mod_phys_lmdz_transfert_para.F90 in phymars and mod_const_mpi.F90 in dyn3d (for compliance with parallelism)
added created generic routines 'planetwide_maxval' and 'planetwide_minval', in module "planetwide_mod", that enable obtaining the max and min of a field over the whole planet. This should be further imroved with computation of means (possibly area weighed), etc.

File size: 10.5 KB

Line
1	SUBROUTINE iniwrite(nid,idayref,phis)
2
3	use comsoil_h, only: mlayer, nsoilmx
4	IMPLICIT NONE
5
6	c=======================================================================
7	c
8	c Auteur: L. Fairhead , P. Le Van, Y. Wanherdrick, F. Forget
9	c -------
10	c
11	c Objet:
12	c ------
13	c
14	c 'Initialize' the diagfi.nc file: write down dimensions as well
15	c as time-independent fields (e.g: geopotential, mesh area, ...)
16	c
17	c=======================================================================
18	c-----------------------------------------------------------------------
19	c Declarations:
20	c -------------
21
22	#include "dimensions.h"
23	#include "paramet.h"
24	#include "comcstfi.h"
25	#include "comvert.h"
26	#include "comgeom.h"
27	#include "ener.h"
28	#include "logic.h"
29	#include "netcdf.inc"
30	#include "serre.h"
31
32	c Arguments:
33	c ----------
34
35	integer,intent(in) :: nid ! NetCDF file ID
36	INTEGER*4,intent(in) :: idayref ! date (initial date for this run)
37	real,intent(in) :: phis(ip1jmp1) ! surface geopotential
38
39	c Local:
40	c ------
41	INTEGER length,l
42	parameter (length = 100)
43	REAL tab_cntrl(length) ! run parameters are stored in this array
44	INTEGER ierr
45
46	integer :: nvarid,idim_index,idim_rlonu,idim_rlonv
47	integer :: idim_rlatu,idim_rlatv,idim_llmp1,idim_llm
48	integer :: idim_nsoilmx ! "subsurface_layers" dimension ID #
49	integer, dimension(2) :: id
50	c-----------------------------------------------------------------------
51
52	DO l=1,length
53	tab_cntrl(l)=0.
54	ENDDO
55	tab_cntrl(1) = real(iim)
56	tab_cntrl(2) = real(jjm)
57	tab_cntrl(3) = real(llm)
58	tab_cntrl(4) = real(idayref)
59	tab_cntrl(5) = rad
60	tab_cntrl(6) = omeg
61	tab_cntrl(7) = g
62	tab_cntrl(8) = mugaz
63	tab_cntrl(9) = rcp
64	tab_cntrl(10) = daysec
65	tab_cntrl(11) = dtphys
66	tab_cntrl(12) = etot0
67	tab_cntrl(13) = ptot0
68	tab_cntrl(14) = ztot0
69	tab_cntrl(15) = stot0
70	tab_cntrl(16) = ang0
71	c
72	c .......... P.Le Van ( ajout le 8/04/96 ) .........
73	c ..... parametres pour le zoom ......
74	tab_cntrl(17) = clon
75	tab_cntrl(18) = clat
76	tab_cntrl(19) = grossismx
77	tab_cntrl(20) = grossismy
78	c
79	c ..... ajout le 6/05/97 et le 15/10/97 .......
80	c
81	IF ( fxyhypb ) THEN
82	tab_cntrl(21) = 1.
83	tab_cntrl(22) = dzoomx
84	tab_cntrl(23) = dzoomy
85	ELSE
86	tab_cntrl(21) = 0.
87	tab_cntrl(22) = dzoomx
88	tab_cntrl(23) = dzoomy
89	tab_cntrl(24) = 0.
90	IF( ysinus ) tab_cntrl(24) = 1.
91	ENDIF
92
93	c .........................................................
94
95	! Define dimensions
96
97	ierr = NF_REDEF (nid)
98
99	ierr = NF_DEF_DIM (nid, "index", length, idim_index)
100	ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)
101	ierr = NF_DEF_DIM (nid, "latitude", jjp1, idim_rlatu)
102	ierr = NF_DEF_DIM (nid, "longitude", iip1, idim_rlonv)
103	ierr = NF_DEF_DIM (nid, "rlatv", jjm, idim_rlatv)
104	ierr = NF_DEF_DIM (nid, "interlayer", (llm+1), idim_llmp1)
105	ierr = NF_DEF_DIM (nid, "altitude", llm, idim_llm)
106	ierr = NF_DEF_DIM (nid,"subsurface_layers",nsoilmx,idim_nsoilmx)
107	c
108	ierr = NF_ENDDEF(nid)
109
110	c Contol parameters for this run
111	ierr = NF_REDEF (nid)
112	#ifdef NC_DOUBLE
113	ierr = NF_DEF_VAR (nid, "controle", NF_DOUBLE, 1,
114	. idim_index,nvarid)
115	#else
116	ierr = NF_DEF_VAR (nid, "controle", NF_FLOAT, 1,
117	. idim_index,nvarid)
118	#endif
119	ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 18,
120	. "Control parameters")
121	ierr = NF_ENDDEF(nid)
122	#ifdef NC_DOUBLE
123	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)
124	#else
125	ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl)
126	#endif
127
128	c --------------------------
129	c longitudes and latitudes
130	ierr = NF_REDEF (nid)
131	#ifdef NC_DOUBLE
132	ierr = NF_DEF_VAR (nid, "rlonu", NF_DOUBLE, 1, idim_rlonu,nvarid)
133	#else
134	ierr = NF_DEF_VAR (nid, "rlonu", NF_FLOAT, 1, idim_rlonu,nvarid)
135	#endif
136	ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 21,
137	. "Longitudes at u nodes")
138	ierr = NF_ENDDEF(nid)
139	#ifdef NC_DOUBLE
140	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonu/pi*180)
141	#else
142	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonu/pi*180)
143	#endif
144	c
145	c --------------------------
146	ierr = NF_REDEF (nid)
147	#ifdef NC_DOUBLE
148	ierr =NF_DEF_VAR(nid, "latitude", NF_DOUBLE, 1, idim_rlatu,nvarid)
149	#else
150	ierr =NF_DEF_VAR (nid, "latitude", NF_FLOAT, 1, idim_rlatu,nvarid)
151	#endif
152	ierr =NF_PUT_ATT_TEXT(nid,nvarid,'units',13,"degrees_north")
153	ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 14,
154	. "North latitude")
155	ierr = NF_ENDDEF(nid)
156	#ifdef NC_DOUBLE
157	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatu/pi*180)
158	#else
159	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatu/pi*180)
160	#endif
161	c
162	c --------------------------
163	ierr = NF_REDEF (nid)
164	#ifdef NC_DOUBLE
165	ierr =NF_DEF_VAR(nid,"longitude", NF_DOUBLE, 1, idim_rlonv,nvarid)
166	#else
167	ierr = NF_DEF_VAR(nid,"longitude", NF_FLOAT, 1, idim_rlonv,nvarid)
168	#endif
169	ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 14,
170	. "East longitude")
171	ierr = NF_PUT_ATT_TEXT(nid,nvarid,'units',12,"degrees_east")
172	ierr = NF_ENDDEF(nid)
173	#ifdef NC_DOUBLE
174	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonv/pi*180)
175	#else
176	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonv/pi*180)
177	#endif
178	c
179	c --------------------------
180	ierr = NF_REDEF (nid)
181	#ifdef NC_DOUBLE
182	ierr = NF_DEF_VAR (nid, "altitude", NF_DOUBLE, 1,
183	. idim_llm,nvarid)
184	#else
185	ierr = NF_DEF_VAR (nid, "altitude", NF_FLOAT, 1,
186	. idim_llm,nvarid)
187	#endif
188	ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",10,"pseudo-alt")
189	ierr = NF_PUT_ATT_TEXT (nid,nvarid,'units',2,"km")
190	ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',2,"up")
191
192	ierr = NF_ENDDEF(nid)
193	#ifdef NC_DOUBLE
194	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pseudoalt)
195	#else
196	ierr = NF_PUT_VAR_REAL (nid,nvarid,pseudoalt)
197	#endif
198	c
199	c --------------------------
200	ierr = NF_REDEF (nid)
201	#ifdef NC_DOUBLE
202	ierr = NF_DEF_VAR (nid, "rlatv", NF_DOUBLE, 1, idim_rlatv,nvarid)
203	#else
204	ierr = NF_DEF_VAR (nid, "rlatv", NF_FLOAT, 1, idim_rlatv,nvarid)
205	#endif
206	ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 20,
207	. "Latitudes at v nodes")
208	ierr = NF_ENDDEF(nid)
209	#ifdef NC_DOUBLE
210	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatv/pi*180)
211	#else
212	ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatv/pi*180)
213	#endif
214	c
215	c --------------------------
216	c Vertical levels
217	call def_var(nid,"aps","hybrid pressure at midlayers ","Pa",
218	. 1,idim_llm,nvarid,ierr)
219	#ifdef NC_DOUBLE
220	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aps)
221	#else
222	ierr = NF_PUT_VAR_REAL (nid,nvarid,aps)
223	#endif
224
225	call def_var(nid,"bps","hybrid sigma at midlayers"," ",
226	. 1,idim_llm,nvarid,ierr)
227	#ifdef NC_DOUBLE
228	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bps)
229	#else
230	ierr = NF_PUT_VAR_REAL (nid,nvarid,bps)
231	#endif
232
233	call def_var(nid,"ap","hybrid pressure at interlayers","Pa",
234	. 1,idim_llmp1,nvarid,ierr)
235	#ifdef NC_DOUBLE
236	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ap)
237	#else
238	ierr = NF_PUT_VAR_REAL (nid,nvarid,ap)
239	#endif
240
241	call def_var(nid,"bp","hybrid sigma at interlayers"," ",
242	. 1,idim_llmp1,nvarid,ierr)
243	#ifdef NC_DOUBLE
244	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bp)
245	#else
246	ierr = NF_PUT_VAR_REAL (nid,nvarid,bp)
247	#endif
248
249	!-------------------------------
250	! (soil) depth variable mlayer() (known from comsoil_h)
251	!-------------------------------
252	ierr=NF_REDEF (nid) ! Enter NetCDF (re-)define mode
253	! define variable
254	#ifdef NC_DOUBLE
255	ierr=NF_DEF_VAR(nid,"soildepth",NF_DOUBLE,1,idim_nsoilmx,nvarid)
256	#else
257	ierr=NF_DEF_VAR(nid,"soildepth",NF_FLOAT,1,idim_nsoilmx,nvarid)
258	#endif
259	ierr=NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 20,
260	. "Soil mid-layer depth")
261	ierr=NF_PUT_ATT_TEXT (nid,nvarid,"units",1,"m")
262	ierr=NF_PUT_ATT_TEXT (nid,nvarid,"positive",4,"down")
263	ierr=NF_ENDDEF(nid) ! Leave NetCDF define mode
264	! write variable
265	#ifdef NC_DOUBLE
266	ierr=NF_PUT_VAR_DOUBLE (nid,nvarid,mlayer)
267	#else
268	ierr=NF_PUT_VAR_REAL (nid,nvarid,mlayer)
269	#endif
270
271	c
272	c --------------------------
273	c Mesh area and conversion coefficients cov. <-> contra. <--> natural
274
275	id(1)=idim_rlonu
276	id(2)=idim_rlatu
277	c
278	ierr = NF_REDEF (nid)
279	#ifdef NC_DOUBLE
280	ierr = NF_DEF_VAR (nid, "cu", NF_DOUBLE, 2, id,nvarid)
281	#else
282	ierr = NF_DEF_VAR (nid, "cu", NF_FLOAT, 2, id,nvarid)
283	#endif
284	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 40,
285	. "Conversion coefficients cov <--> natural")
286	ierr = NF_ENDDEF(nid)
287	#ifdef NC_DOUBLE
288	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cu)
289	#else
290	ierr = NF_PUT_VAR_REAL (nid,nvarid,cu)
291	#endif
292	c
293	id(1)=idim_rlonv
294	id(2)=idim_rlatv
295	c
296	c --------------------------
297	ierr = NF_REDEF (nid)
298	#ifdef NC_DOUBLE
299	ierr = NF_DEF_VAR (nid, "cv", NF_DOUBLE, 2, id,nvarid)
300	#else
301	ierr = NF_DEF_VAR (nid, "cv", NF_FLOAT, 2, id,nvarid)
302	#endif
303	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 40,
304	. "Conversion coefficients cov <--> natural")
305	ierr = NF_ENDDEF(nid)
306	#ifdef NC_DOUBLE
307	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cv)
308	#else
309	ierr = NF_PUT_VAR_REAL (nid,nvarid,cv)
310	#endif
311	c
312	id(1)=idim_rlonv
313	id(2)=idim_rlatu
314	c
315	c --------------------------
316	ierr = NF_REDEF (nid)
317	#ifdef NC_DOUBLE
318	ierr = NF_DEF_VAR (nid, "aire", NF_DOUBLE, 2, id,nvarid)
319	#else
320	ierr = NF_DEF_VAR (nid, "aire", NF_FLOAT, 2, id,nvarid)
321	#endif
322	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9,
323	. "Mesh area")
324	ierr = NF_ENDDEF(nid)
325	#ifdef NC_DOUBLE
326	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aire)
327	#else
328	ierr = NF_PUT_VAR_REAL (nid,nvarid,aire)
329	#endif
330	c
331	c Surface geopotential
332	id(1)=idim_rlonv
333	id(2)=idim_rlatu
334	c
335	ierr = NF_REDEF (nid)
336	#ifdef NC_DOUBLE
337	ierr = NF_DEF_VAR (nid, "phisinit", NF_DOUBLE, 2, id,nvarid)
338	#else
339	ierr = NF_DEF_VAR (nid, "phisinit", NF_FLOAT, 2, id,nvarid)
340	#endif
341	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 27,
342	. "Geopotential at the surface")
343	ierr = NF_ENDDEF(nid)
344	#ifdef NC_DOUBLE
345	ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,phis)
346	#else
347	ierr = NF_PUT_VAR_REAL (nid,nvarid,phis)
348	#endif
349	c
350
351	write(,)'iniwrite: iim,jjm,llm,idayref',iim,jjm,llm,idayref
352	write(,)'iniwrite: rad,omeg,g,mugaz,rcp',
353	s rad,omeg,g,mugaz,rcp
354	write(,)'iniwrite: daysec,dtphys',daysec,dtphys
355
356	END

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format