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

Annotation of /trunk/Sources/phylmd/initphysto.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 186 - (hide annotations)
Mon Mar 21 15:36:26 2016 UTC (8 years, 2 months ago) by guez
File size: 5992 byte(s) 
Removed variables nlm and nlp of module cv30_param_m. We do not
believe much in the benefit of these intermediary variables so we go
for clarity.

Removed variable noff of module cv30_param_m. Never used anywhere
else. Just set the value of nl explicitly in cv30_param.

Removed argument nd of cv30_param. Only called with nd = klev.

Replaced calls to zilch by array assignments. There was a strange
double call to zilch with the same arguments in cv30_mixing.

Removed procedure cv_flag. Just set the value of variable cvflag_grav
of module cvflag at declaration.
1 guez 155 module initphysto_m
2 guez 3
3 guez 155 IMPLICIT NONE
4 guez 3
5 guez 155 contains
6 guez 3
7 guez 155 SUBROUTINE initphysto(infile,rlon,rlat,tstep,t_ops,t_wrt,fileid)
8 guez 3
9 guez 155 ! From phylmd/initphysto.F,v 1.2 2004/06/22 11:45:32
10     ! Routine d'initialisation des ecritures des fichiers histoires LMDZ
11     ! au format IOIPSL
12 guez 3
13 guez 155 ! Appels succesifs des routines: histbeg
14     ! histhori
15     ! histver
16     ! histdef
17     ! histend
18 guez 3
19 guez 155 ! Entree:
20 guez 3
21 guez 155 ! infile: nom du fichier histoire a creer
22     ! day0,anne0: date de reference
23     ! tstep: duree du pas de temps en seconde
24     ! t_ops: frequence de l'operation pour IOIPSL
25     ! t_wrt: frequence d'ecriture sur le fichier
26 guez 3
27 guez 155 ! Sortie:
28     ! fileid: ID du fichier netcdf cree
29     ! filevid:ID du fichier netcdf pour la grille v
30 guez 3
31 guez 155 ! L. Fairhead, 03/99
32 guez 3
33 guez 155 use dynetat0_m, only: day_ref, annee_ref
34 guez 186 use gr_fi_ecrit_m, only: gr_fi_ecrit
35 guez 155 USE histbeg_totreg_m, ONLY : histbeg_totreg
36     USE histdef_m, ONLY : histdef
37     USE histend_m, ONLY : histend
38     use histsync_m, only: histsync
39     USE histvert_m, ONLY : histvert
40     USE dimens_m
41     USE paramet_m
42     USE comconst
43     USE indicesol
44     USE dimphy
45     use conf_gcm_m
46     USE comgeom
47     USE ymds2ju_m
48 guez 3
49 guez 155 ! Arguments
50     CHARACTER(len=*) infile
51     INTEGER nhoriid, i
52     REAL, INTENT (IN) :: tstep
53     REAL t_ops, t_wrt
54     INTEGER fileid
55     INTEGER l
56     REAL nivsigs(llm)
57 guez 3
58 guez 155 ! Variables locales
59 guez 3
60 guez 155 INTEGER tau0
61     REAL zjulian
62     INTEGER zvertiid
63     LOGICAL ok_sync
64     REAL zx_lon(iim,jjm+1), zx_lat(iim,jjm+1)
65 guez 3
66 guez 155 REAL, INTENT (IN) :: rlon(klon), rlat(klon)
67 guez 3
68 guez 155 !-----------------------------------------------------
69 guez 3
70 guez 155 ! Initialisations
71     ok_sync = .TRUE.
72 guez 3
73 guez 155 ! Appel a histbeg: creation du fichier netcdf et initialisations
74     ! diverses
75 guez 3
76 guez 155 CALL ymds2ju(annee_ref,1,day_ref,0.0,zjulian)
77     tau0 = 0
78 guez 3
79 guez 155 CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon)
80     DO i = 1, iim
81     zx_lon(i,1) = rlon(i+1)
82     zx_lon(i,jjm+1) = rlon(i+1)
83     END DO
84     CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat)
85 guez 3
86    
87 guez 155 CALL histbeg_totreg(infile,zx_lon(:,1),zx_lat(1,:),1,iim,1,jjm+1,tau0, &
88     zjulian,tstep,nhoriid,fileid)
89 guez 3
90 guez 155 ! Appel a histvert pour la grille verticale
91 guez 3
92 guez 155 DO l = 1, llm
93     nivsigs(l) = float(l)
94     END DO
95 guez 3
96 guez 155 CALL histvert(fileid,'sig_s','Niveaux sigma','sigma_level',nivsigs, &
97     zvertiid)
98 guez 3
99 guez 155 ! Appels a histdef pour la definition des variables a sauvegarder
100 guez 3
101 guez 155 CALL histdef(fileid,'phis','Surface geop. height','-',iim,jjm+1,nhoriid, &
102     1,1,1,-99,'once',t_ops,t_wrt)
103 guez 3
104 guez 155 CALL histdef(fileid,'aire','Grid area','-',iim,jjm+1,nhoriid,1,1,1,-99, &
105     'once',t_ops,t_wrt)
106 guez 3
107 guez 155 CALL histdef(fileid,'dtime','tps phys ','s',1,1,nhoriid,1,1,1,-99, &
108     'once',t_ops,t_wrt)
109 guez 3
110 guez 155 CALL histdef(fileid,'istphy','tps stock','s',1,1,nhoriid,1,1,1,-99, &
111     'once',t_ops,t_wrt)
112 guez 3
113 guez 155 ! T
114 guez 3
115 guez 155 CALL histdef(fileid,'t','Temperature','K',iim,jjm+1,nhoriid,llm,1,llm, &
116     zvertiid,'inst(X)',t_ops,t_wrt)
117 guez 31
118 guez 155 CALL histdef(fileid,'mfu','flx m. pan. mt','kg m/s',iim,jjm+1,nhoriid, &
119     llm,1,llm,zvertiid,'inst(X)',t_ops,t_wrt)
120 guez 31
121 guez 155 CALL histdef(fileid,'mfd','flx m. pan. des','kg m/s',iim,jjm+1,nhoriid, &
122     llm,1,llm,zvertiid,'inst(X)',t_ops,t_wrt)
123 guez 31
124    
125 guez 155 ! en_u
126 guez 31
127 guez 155 CALL histdef(fileid,'en_u','flx ent pan mt','kg m/s',iim,jjm+1,nhoriid, &
128     llm,1,llm,zvertiid,'inst(X)',t_ops,t_wrt)
129 guez 31
130 guez 155 CALL histdef(fileid,'de_u','flx det pan mt','kg m/s',iim,jjm+1,nhoriid, &
131     llm,1,llm,zvertiid,'inst(X)',t_ops,t_wrt)
132 guez 31
133    
134 guez 155 ! en_d
135 guez 31
136 guez 155 CALL histdef(fileid,'en_d','flx ent pan dt','kg m/s',iim,jjm+1,nhoriid, &
137     llm,1,llm,zvertiid,'inst(X)',t_ops,t_wrt)
138 guez 31
139    
140    
141 guez 155 ! de_d
142 guez 31
143 guez 155 CALL histdef(fileid,'de_d','flx det pan dt','kg m/s',iim,jjm+1,nhoriid, &
144     llm,1,llm,zvertiid,'inst(X)',t_ops,t_wrt)
145 guez 31
146 guez 155 ! coefh frac_impa,frac_nucl
147 guez 31
148 guez 155 CALL histdef(fileid,'coefh',' ',' ',iim,jjm+1,nhoriid,llm,1,llm, &
149     zvertiid,'inst(X)',t_ops,t_wrt)
150 guez 31
151 guez 155 ! abderrahmane le 16 09 02
152     CALL histdef(fileid,'fm_th',' ',' ',iim,jjm+1,nhoriid,llm,1,llm, &
153     zvertiid,'inst(X)',t_ops,t_wrt)
154 guez 31
155 guez 155 CALL histdef(fileid,'en_th',' ',' ',iim,jjm+1,nhoriid,llm,1,llm, &
156     zvertiid,'inst(X)',t_ops,t_wrt)
157     ! fin aj
158 guez 31
159 guez 155 CALL histdef(fileid,'frac_impa',' ',' ',iim,jjm+1,nhoriid,llm,1,llm, &
160     zvertiid,'inst(X)',t_ops,t_wrt)
161 guez 31
162 guez 155 CALL histdef(fileid,'frac_nucl',' ',' ',iim,jjm+1,nhoriid,llm,1,llm, &
163     zvertiid,'inst(X)',t_ops,t_wrt)
164 guez 31
165    
166 guez 155 ! pyu1
167 guez 31
168 guez 155 CALL histdef(fileid,'pyu1',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
169     'inst(X)',t_ops,t_wrt)
170 guez 31
171    
172 guez 155 ! pyv1
173 guez 31
174 guez 155 CALL histdef(fileid,'pyv1',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
175     'inst(X)',t_ops,t_wrt)
176 guez 31
177 guez 155 CALL histdef(fileid,'ftsol1',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
178     'inst(X)',t_ops,t_wrt)
179 guez 31
180    
181 guez 155 ! ftsol2
182 guez 31
183 guez 155 CALL histdef(fileid,'ftsol2',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
184     'inst(X)',t_ops,t_wrt)
185 guez 31
186    
187 guez 155 ! ftsol3
188 guez 31
189 guez 155 CALL histdef(fileid,'ftsol3',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
190     'inst(X)',t_ops,t_wrt)
191 guez 31
192    
193 guez 155 ! ftsol4
194 guez 31
195 guez 155 CALL histdef(fileid,'ftsol4',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
196     'inst(X)',t_ops,t_wrt)
197 guez 31
198    
199 guez 155 ! rain
200 guez 31
201 guez 155 CALL histdef(fileid,'rain',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
202     'inst(X)',t_ops,t_wrt)
203 guez 31
204    
205 guez 155 ! psrf1
206 guez 31
207 guez 155 CALL histdef(fileid,'psrf1',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
208     'inst(X)',t_ops,t_wrt)
209 guez 31
210    
211 guez 155 ! psrf2
212 guez 31
213 guez 155 CALL histdef(fileid,'psrf2',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
214     'inst(X)',t_ops,t_wrt)
215 guez 31
216    
217 guez 155 ! psrf3
218 guez 31
219 guez 155 CALL histdef(fileid,'psrf3',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
220     'inst(X)',t_ops,t_wrt)
221 guez 31
222    
223 guez 155 ! psrf4
224 guez 31
225 guez 155 CALL histdef(fileid,'psrf4',' ',' ',iim,jjm+1,nhoriid,1,1,1,-99, &
226     'inst(X)',t_ops,t_wrt)
227    
228     CALL histend(fileid)
229     IF (ok_sync) CALL histsync
230    
231     END SUBROUTINE initphysto
232    
233     end module initphysto_m
  ViewVC Help
Powered by ViewVC 1.1.21