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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 139 - (show annotations)
Tue May 26 17:46:03 2015 UTC (9 years ago) by guez
File size: 5823 byte(s) 
dynetat0 read rlonu, rlatu, rlonv, rlatv, cu_2d, cv_2d, aire_2d from
"start.nc" and then these variables were overwritten by
inigeom. Corrected this. Now, inigeom does not compute rlonu, rlatu,
rlonv and rlatv. Moreover, cu_2d, cv_2d, aire_2d are not written to
"restart.nc". Since xprimu, xprimv, xprimm025, xprimp025, rlatu1,
rlatu2, yprimu1, yprimu2 are computed at the same time as rlonu,
rlatu, rlonv, rlatv, and since it would not be convenient to separate
those computations, we decide to write xprimu, xprimv, xprimm025,
xprimp025, rlatu1, rlatu2, yprimu1, yprimu2 into "restart.nc", read
them from "start.nc" and not compute them in inigeom. So, in summary,
"start.nc" contains all the coordinates and their derivatives, and
inigeom only computes the 2D-variables.

Technical details:

Moved variables rlatu, rlonv, rlonu, rlatv, xprimu, xprimv from module
comgeom to module dynetat0_m. Upgraded local variables rlatu1,
yprimu1, rlatu2, yprimu2, xprimm025, xprimp025 of procedure inigeom to
variables of module dynetat0_m.

Removed unused local variable yprimu of procedure inigeom and
corresponding argument yyprimu of fyhyp.

Moved variables clat, clon, grossismx, grossismy, dzoomx, dzoomy,
taux, tauy from module serre to module dynetat0_m (since they are read
from "start.nc"). The default values are now defined in read_serre
instead of in the declarations. Changed name of module serre to
read_serre_m, no more module variable here.

The calls to fxhyp and fyhyp are moved from inigeom to etat0.

Side effects in programs other than gcm: etat0 and read_serre write
variables of module dynetat0; the programs test_fxyp and
test_inter_barxy need more source files.

Removed unused arguments len and nd of cv3_tracer. Removed unused
argument PPSOL of LWU.

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