/[lmdze]/trunk/phylmd/initphysto.f
ViewVC logotype

Annotation of /trunk/phylmd/initphysto.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 12 - (hide annotations)
Mon Jul 21 16:05:07 2008 UTC (15 years, 10 months ago) by guez
Original Path: trunk/libf/phylmd/initphysto.f
File size: 7892 byte(s) 
-- Minor modification of input/output:

Created procedure "read_logic". Variables of module "logic" are read
by "read_logic" instead of "conf_gcm". Variable "offline" of module
"conf_gcm" is read from namelist instead of "*.def".

Deleted arguments "dtime", "co2_ppm_etat0", "solaire_etat0",
"tabcntr0" and local variables "radpas", "tab_cntrl" of
"phyetat0". "phyetat0" does not read "controle" in "startphy.nc" any
longer. "phyetat0" now reads global attribute "itau_phy" from
"startphy.nc". "phyredem" does not create variable "controle" in
"startphy.nc" any longer. "phyredem" now writes global attribute
"itau_phy" of "startphy.nc". Deleted argument "tabcntr0" of
"printflag". Removed diagnostic messages written by "printflag" for
comparison of the variable "controle" of "startphy.nc" and the
variables read from "*.def" or namelist input.

-- Removing unwanted functionality:

Removed variable "lunout" from module "iniprint", replaced everywhere
by standard output.

Removed case "ocean == 'couple'" in "clmain", "interfsurf_hq" and
"physiq". Removed procedure "interfoce_cpl".

-- Should not change anything at run time:

Automated creation of graphs in documentation. More documentation on
input files.

Converted Fortran files to free format: "phyredem.f90", "printflag.f90".

Split module "clesphy" into "clesphys" and "clesphys2".

Removed variables "conser", "leapf", "forward", "apphys", "apdiss" and
"statcl" from module "logic". Added arguments "conser" to "advect",
"leapf" to "integrd". Added local variables "forward", "leapf",
"apphys", "conser", "apdiss" in "leapfrog".

Added intent attributes.

Deleted arguments "dtime" of "phyredem", "pdtime" of "flxdtdq", "sh"
of "phytrac", "dt" of "yamada".

Deleted local variables "dtime", "co2_ppm_etat0", "solaire_etat0",
"length", "tabcntr0" in "physiq". Replaced all references to "dtime"
by references to "pdtphys".
1 guez 3 !
2     ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/initphysto.F,v 1.2 2004/06/22 11:45:32 lmdzadmin Exp $
3     !
4     C
5     C
6     subroutine initphysto
7     . (infile,
8     . rlon, rlat, tstep,t_ops,t_wrt,nq,fileid)
9    
10     USE IOIPSL
11    
12     use dimens_m
13     use paramet_m
14     use comconst
15     use indicesol
16     use dimphy
17     use logic
18     use comgeom
19     use serre
20     use temps
21     use ener
22     implicit none
23    
24     C
25     C Routine d'initialisation des ecritures des fichiers histoires LMDZ
26     C au format IOIPSL
27     C
28     C Appels succesifs des routines: histbeg
29     C histhori
30     C histver
31     C histdef
32     C histend
33     C
34     C Entree:
35     C
36     C infile: nom du fichier histoire a creer
37     C day0,anne0: date de reference
38     C tstep: duree du pas de temps en seconde
39     C t_ops: frequence de l'operation pour IOIPSL
40     C t_wrt: frequence d'ecriture sur le fichier
41     C nq: nombre de traceurs
42     C
43     C Sortie:
44     C fileid: ID du fichier netcdf cree
45     C filevid:ID du fichier netcdf pour la grille v
46     C
47     C L. Fairhead, LMD, 03/99
48     C
49     C =====================================================================
50     C
51     C Declarations
52    
53     C Arguments
54     character*(*) infile
55     integer*4 nhoriid, i
56 guez 12 real, intent(in):: tstep
57     real t_ops, t_wrt
58 guez 3 integer fileid, filevid
59     integer nq,l
60     real nivsigs(llm)
61    
62     C Variables locales
63     C
64     integer tau0
65     real zjulian
66     character*3 str
67     character*10 ctrac
68     integer iq
69     integer uhoriid, vhoriid, thoriid, zvertiid
70     integer ii,jj
71     integer zan, idayref
72     logical ok_sync
73     REAL zx_lon(iim,jjm+1), zx_lat(iim,jjm+1)
74     C
75     REAL, intent(in):: rlon(klon), rlat(klon)
76    
77     C Initialisations
78     C
79     pi = 4. * atan (1.)
80     str='q '
81     ctrac = 'traceur '
82     ok_sync= .true.
83     C
84     C Appel a histbeg: creation du fichier netcdf et initialisations diverses
85     C
86    
87     zan = annee_ref
88     idayref = day_ref
89     CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
90     tau0 = 0
91    
92     CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon)
93     DO i = 1, iim
94     zx_lon(i,1) = rlon(i+1)
95     zx_lon(i,jjm+1) = rlon(i+1)
96     ENDDO
97     CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat)
98    
99    
100     call histbeg_totreg(infile, iim, zx_lon(:,1), jjm+1, zx_lat(1,:),
101     . 1, iim, 1, jjm+1,
102     . tau0, zjulian, tstep, nhoriid, fileid)
103    
104     C Appel a histvert pour la grille verticale
105     C
106     DO l=1,llm
107     nivsigs(l)=float(l)
108     ENDDO
109    
110     write(*,*) 'avant histvert ds initphysto'
111    
112     call histvert(fileid, 'sig_s', 'Niveaux sigma',
113     . 'sigma_level',
114     . llm, nivsigs, zvertiid)
115     C
116     C Appels a histdef pour la definition des variables a sauvegarder
117     C
118     write(*,*) 'apres histvert ds initphysto'
119    
120     CALL histdef(fileid, "phis", "Surface geop. height", "-",
121     . iim,jjm+1,nhoriid, 1,1,1, -99, 32,
122     . "once", t_ops, t_wrt)
123     c
124     write(*,*) 'apres phis ds initphysto'
125    
126     CALL histdef(fileid, "aire", "Grid area", "-",
127     . iim,jjm+1,nhoriid, 1,1,1, -99, 32,
128     . "once", t_ops, t_wrt)
129     write(*,*) 'apres aire ds initphysto'
130    
131     CALL histdef(fileid, "dtime", "tps phys ", "s",
132     . 1,1,nhoriid, 1,1,1, -99, 32,
133     . "once", t_ops, t_wrt)
134    
135     CALL histdef(fileid, "istphy", "tps stock", "s",
136     . 1,1,nhoriid, 1,1,1, -99, 32,
137     . "once", t_ops, t_wrt)
138    
139     C T
140     C
141     call histdef(fileid, 't', 'Temperature', 'K',
142     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
143     . 32, 'inst(X)', t_ops, t_wrt)
144     write(*,*) 'apres t ds initphysto'
145     C mfu
146     C
147     call histdef(fileid, 'mfu', 'flx m. pan. mt', 'kg m/s',
148     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
149     . 32, 'inst(X)', t_ops, t_wrt)
150     write(*,*) 'apres mfu ds initphysto'
151     C
152     C mfd
153     C
154     call histdef(fileid, 'mfd', 'flx m. pan. des', 'kg m/s',
155     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
156     . 32, 'inst(X)', t_ops, t_wrt)
157    
158     C
159     C en_u
160     C
161     call histdef(fileid, 'en_u', 'flx ent pan mt', 'kg m/s',
162     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
163     . 32, 'inst(X)', t_ops, t_wrt)
164     write(*,*) 'apres en_u ds initphysto'
165     C
166     C de_u
167     C
168     call histdef(fileid, 'de_u', 'flx det pan mt', 'kg m/s',
169     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
170     . 32, 'inst(X)', t_ops, t_wrt)
171    
172     C
173     C en_d
174     C
175     call histdef(fileid, 'en_d', 'flx ent pan dt', 'kg m/s',
176     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
177     . 32, 'inst(X)', t_ops, t_wrt)
178     C
179    
180     C
181     C de_d
182     C
183     call histdef(fileid, 'de_d', 'flx det pan dt', 'kg m/s',
184     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
185     . 32, 'inst(X)', t_ops, t_wrt)
186    
187     c coefh frac_impa,frac_nucl
188    
189     call histdef(fileid, "coefh", " ", " ",
190     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
191     . 32, "inst(X)", t_ops, t_wrt)
192    
193     c abderrahmane le 16 09 02
194     call histdef(fileid, "fm_th", " ", " ",
195     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
196     . 32, "inst(X)", t_ops, t_wrt)
197    
198     call histdef(fileid, "en_th", " ", " ",
199     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
200     . 32, "inst(X)", t_ops, t_wrt)
201     c fin aj
202    
203     write(*,*) 'apres coefh ds initphysto'
204    
205     call histdef(fileid, 'frac_impa', ' ', ' ',
206     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
207     . 32, 'inst(X)', t_ops, t_wrt)
208    
209     call histdef(fileid, 'frac_nucl', ' ', ' ',
210     . iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
211     . 32, 'inst(X)', t_ops, t_wrt)
212    
213     c
214     c pyu1
215     c
216     CALL histdef(fileid, "pyu1", " ", " ",
217     . iim,jjm+1,nhoriid, 1,1,1, -99, 32,
218     . "inst(X)", t_ops, t_wrt)
219    
220     c
221     c pyv1
222     c
223     CALL histdef(fileid, "pyv1", " ", " ",
224     . iim,jjm+1,nhoriid, 1,1,1, -99, 32,
225     . "inst(X)", t_ops, t_wrt)
226    
227     write(*,*) 'apres pyv1 ds initphysto'
228     c
229     c ftsol1
230     c
231     call histdef(fileid, "ftsol1", " ", " ",
232     . iim, jjm+1, nhoriid, 1, 1,1, -99,32,
233     . "inst(X)", t_ops, t_wrt)
234    
235     c
236     c ftsol2
237     c
238     call histdef(fileid, "ftsol2", " ", " ",
239     . iim, jjm+1, nhoriid, 1, 1,1, -99,32,
240     . "inst(X)", t_ops, t_wrt)
241    
242     c
243     c ftsol3
244     c
245     call histdef(fileid, "ftsol3", " ", " ",
246     . iim, jjm+1, nhoriid, 1, 1,1, -99,
247     . 32, "inst(X)", t_ops, t_wrt)
248    
249     c
250     c ftsol4
251     c
252     call histdef(fileid, "ftsol4", " ", " ",
253     . iim, jjm+1, nhoriid, 1, 1,1, -99,
254     . 32, "inst(X)", t_ops, t_wrt)
255    
256     c
257     c rain
258     c
259     call histdef(fileid, "rain", " ", " ",
260     . iim, jjm+1, nhoriid, 1, 1,1, -99,
261     . 32, "inst(X)", t_ops, t_wrt)
262    
263     c
264     c psrf1
265     c
266     call histdef(fileid, "psrf1", " ", " ",
267     . iim, jjm+1, nhoriid, 1, 1, 1, -99,
268     . 32, "inst(X)", t_ops, t_wrt)
269    
270     c
271     c psrf2
272     c
273     call histdef(fileid, "psrf2", " ", " ",
274     . iim, jjm+1, nhoriid, 1, 1, 1, -99,
275     . 32, "inst(X)", t_ops, t_wrt)
276    
277     c
278     c psrf3
279     c
280     call histdef(fileid, "psrf3", " ", " ",
281     . iim, jjm+1, nhoriid, 1, 1, 1, -99,
282     . 32, "inst(X)", t_ops, t_wrt)
283    
284     c
285     c psrf4
286     c
287     call histdef(fileid, "psrf4", " ", " ",
288     . iim, jjm+1, nhoriid, 1, 1, 1, -99,
289     . 32, "inst(X)", t_ops, t_wrt)
290    
291     write(*,*) 'avant histend ds initphysto'
292    
293     call histend(fileid)
294     c if (ok_sync) call histsync(fileid)
295     if (ok_sync) call histsync
296    
297    
298    
299     return
300     end
  ViewVC Help
Powered by ViewVC 1.1.21