/[lmdze]/trunk/dyn3d/dynredem0.f
ViewVC logotype

Annotation of /trunk/dyn3d/dynredem0.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 10 - (hide annotations)
Fri Apr 18 14:45:53 2008 UTC (16 years, 1 month ago) by guez
Original Path: trunk/libf/dyn3d/dynredem0.f90
File size: 9137 byte(s) 
Added NetCDF directory "/home/guez/include" in "g95.mk" and
"nag_tools.mk".

Added some "intent" attributes in "PVtheta", "advtrac", "caladvtrac",
"calfis", "diagedyn", "dissip", "vlspltqs", "aeropt", "ajsec",
"calltherm", "clmain", "cltrac", "cltracrn", "concvl", "conema3",
"conflx", "fisrtilp", "newmicro", "nuage", "diagcld1", "diagcld2",
"drag_noro", "lift_noro", "SUGWD", "physiq", "phytrac", "radlwsw", "thermcell".

Removed the case "ierr == 0" in "abort_gcm"; moved call to "histclo"
and messages for end of run from "abort_gcm" to "gcm"; replaced call
to "abort_gcm" in "leapfrog" by exit from outer loop.

In "calfis": removed argument "pp" and variable "unskap"; changed
"pksurcp" from scalar to rank 2; use "pressure_var"; rewrote
computation of "zplev", "zplay", "ztfi", "pcvgt" using "dyn_phy";
added computation of "pls".

Removed unused variable in "dynredem0".

In "exner_hyb": changed "dellta" from scalar to rank 1; replaced call
to "ssum" by call to "sum"; removed variables "xpn" and "xps";
replaced some loops by array expressions.

In "leapfrog": use "pressure_var"; deleted variables "p", "longcles".

Converted common blocks "YOECUMF", "YOEGWD" to modules.

Removed argument "pplay" in "cvltr", "diagetpq", "nflxtr".

Created module "raddimlw" from include file "raddimlw.h".

Corrected call to "new_unit" in "test_disvert".
1 guez 7 module dynredem0_m
2 guez 3
3 guez 7 IMPLICIT NONE
4 guez 3
5 guez 7 contains
6 guez 3
7 guez 7 SUBROUTINE dynredem0(fichnom, iday_end, phis)
8 guez 3
9 guez 7 ! From dyn3d/dynredem.F, v 1.2 2004/06/22 11:45:30
10 guez 3
11 guez 7 ! Ecriture du fichier de redémarrage au format NetCDF (initialisation)
12 guez 3
13 guez 7 USE IOIPSL, only: ymds2ju, ju2ymds
14     use dimens_m, only: iim, jjm, llm, nqmx
15 guez 10 use paramet_m, only: iip1, jjp1, llmp1
16 guez 7 use comconst, only: rad, cpp, daysec, dtvr, kappa, g, omeg
17     use comvert, only: pa, bp, ap, nivsigs, preff, presnivs, nivsig
18     use logic
19     use comgeom
20     use serre
21     use temps, only: annee_ref, day_ref, itaufin, itau_dyn
22     use ener
23     use advtrac_m, only: tname, ttext
24     use netcdf95, only: nf95_create, nf95_put_att, nf95_def_dim, &
25     nf95_def_var, NF95_ENDDEF, NF95_PUT_VAR
26     use netcdf, only: NF90_CLOBBER, NF90_GLOBAL, NF90_UNLIMITED, nf90_float
27 guez 5
28 guez 7 CHARACTER(len=*), intent(in):: fichnom
29     INTEGER, intent(in):: iday_end
30     REAL, intent(in):: phis(:, :)
31 guez 5
32 guez 7 ! Local:
33 guez 3
34 guez 7 include "netcdf.inc"
35 guez 5
36 guez 7 INTEGER iq, l
37     INTEGER length
38     PARAMETER (length = 100)
39     REAL tab_cntrl(length) ! tableau des parametres du run
40     INTEGER ierr
41 guez 3
42 guez 7 ! Variables locales pour NetCDF:
43 guez 3
44 guez 7 INTEGER dims2(2), dims3(3), dims4(4)
45     INTEGER idim_index
46     INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
47     INTEGER idim_s, idim_sig
48     INTEGER idim_tim
49     INTEGER nid, nvarid
50 guez 3
51 guez 7 REAL zjulian, hours
52     INTEGER yyears0, jjour0, mmois0
53     character(len=30) unites
54 guez 3
55 guez 7 !-----------------------------------------------------------------------
56 guez 3
57 guez 7 print *, "Call sequence information: dynredem0"
58 guez 3
59 guez 7 call ymds2ju(annee_ref, 1, iday_end, 0.0, zjulian)
60     call ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)
61 guez 5
62 guez 7 DO l=1, length
63     tab_cntrl(l) = 0.
64     ENDDO
65     tab_cntrl(1) = REAL(iim)
66     tab_cntrl(2) = REAL(jjm)
67     tab_cntrl(3) = REAL(llm)
68     tab_cntrl(4) = REAL(day_ref)
69     tab_cntrl(5) = REAL(annee_ref)
70     tab_cntrl(6) = rad
71     tab_cntrl(7) = omeg
72     tab_cntrl(8) = g
73     tab_cntrl(9) = cpp
74     tab_cntrl(10) = kappa
75     tab_cntrl(11) = daysec
76     tab_cntrl(12) = dtvr
77     tab_cntrl(13) = etot0
78     tab_cntrl(14) = ptot0
79     tab_cntrl(15) = ztot0
80     tab_cntrl(16) = stot0
81     tab_cntrl(17) = ang0
82     tab_cntrl(18) = pa
83     tab_cntrl(19) = preff
84 guez 5
85 guez 7 ! Paramètres pour le zoom :
86 guez 3
87 guez 7 tab_cntrl(20) = clon
88     tab_cntrl(21) = clat
89     tab_cntrl(22) = grossismx
90     tab_cntrl(23) = grossismy
91 guez 5
92 guez 7 IF ( fxyhypb ) THEN
93     tab_cntrl(24) = 1.
94     tab_cntrl(25) = dzoomx
95     tab_cntrl(26) = dzoomy
96     tab_cntrl(27) = 0.
97     tab_cntrl(28) = taux
98     tab_cntrl(29) = tauy
99     ELSE
100     tab_cntrl(24) = 0.
101     tab_cntrl(25) = dzoomx
102     tab_cntrl(26) = dzoomy
103     tab_cntrl(27) = 0.
104     tab_cntrl(28) = 0.
105     tab_cntrl(29) = 0.
106     IF( ysinus ) tab_cntrl(27) = 1.
107     ENDIF
108 guez 5
109 guez 7 tab_cntrl(30) = REAL(iday_end)
110     tab_cntrl(31) = REAL(itau_dyn + itaufin)
111 guez 5
112 guez 7 call nf95_create(fichnom, NF90_CLOBBER, nid)
113     call nf95_put_att(nid, NF90_GLOBAL, "title", &
114     "Fichier de démarrage dynamique")
115 guez 5
116 guez 7 ! Definir les dimensions du fichiers:
117 guez 5
118 guez 7 call nf95_def_dim(nid, "index", length, idim_index)
119     call NF95_DEF_DIM(nid, "rlonu", iip1, idim_rlonu)
120     call NF95_DEF_DIM(nid, "rlatu", jjp1, idim_rlatu)
121     call NF95_DEF_DIM(nid, "rlonv", iip1, idim_rlonv)
122     call NF95_DEF_DIM(nid, "rlatv", jjm, idim_rlatv)
123     call NF95_DEF_DIM(nid, "sigs", llm, idim_s)
124     call NF95_DEF_DIM(nid, "sig", llmp1, idim_sig)
125     call NF95_DEF_DIM(nid, "temps", NF90_UNLIMITED, idim_tim)
126 guez 5
127 guez 7 ! Definir et enregistrer certains champs invariants:
128 guez 5
129 guez 7 call nf95_def_var(nid, "controle", NF90_FLOAT, idim_index, nvarid)
130     call nf95_put_att(nid, nvarid, "title", "Parametres de controle")
131     call NF95_ENDDEF(nid)
132     call NF95_PUT_VAR(nid, nvarid, tab_cntrl)
133 guez 5
134 guez 7 ierr = NF_REDEF (nid)
135     call nf95_def_var(nid, "rlonu", NF90_FLOAT, idim_rlonu, nvarid)
136     call nf95_put_att(nid, nvarid, "title", "Longitudes des points U")
137     call NF95_ENDDEF(nid)
138     call NF95_PUT_VAR(nid, nvarid, rlonu)
139 guez 5
140 guez 7 ierr = NF_REDEF (nid)
141     call nf95_def_var(nid, "rlatu", NF90_FLOAT, idim_rlatu, nvarid)
142     call nf95_put_att(nid, nvarid, "title", "Latitudes des points U")
143     call NF95_ENDDEF(nid)
144     call NF95_PUT_VAR(nid, nvarid, rlatu)
145 guez 5
146 guez 7 ierr = NF_REDEF (nid)
147     call nf95_def_var(nid, "rlonv", NF90_FLOAT, idim_rlonv, nvarid)
148     call nf95_put_att(nid, nvarid, "title", "Longitudes des points V")
149     call NF95_ENDDEF(nid)
150     call NF95_PUT_VAR(nid, nvarid, rlonv)
151 guez 5
152 guez 7 ierr = NF_REDEF (nid)
153     call nf95_def_var(nid, "rlatv", NF90_FLOAT, idim_rlatv, nvarid)
154     call nf95_put_att(nid, nvarid, "title", "Latitudes des points V")
155     call NF95_ENDDEF(nid)
156     call NF95_PUT_VAR(nid, nvarid, rlatv)
157 guez 5
158 guez 7 ierr = NF_REDEF (nid)
159     call nf95_def_var(nid, "nivsigs", NF90_FLOAT, idim_s, nvarid)
160     call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches s")
161     call NF95_ENDDEF(nid)
162     call NF95_PUT_VAR(nid, nvarid, nivsigs)
163 guez 5
164 guez 7 ierr = NF_REDEF (nid)
165     call nf95_def_var(nid, "nivsig", NF90_FLOAT, idim_sig, nvarid)
166     call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches sigma")
167     call NF95_ENDDEF(nid)
168     call NF95_PUT_VAR(nid, nvarid, nivsig)
169 guez 5
170 guez 7 ierr = NF_REDEF (nid)
171     call nf95_def_var(nid, "ap", NF90_FLOAT, idim_sig, nvarid)
172     call nf95_put_att(nid, nvarid, "title", "Coefficient A pour hybride")
173     call NF95_ENDDEF(nid)
174     call NF95_PUT_VAR(nid, nvarid, ap)
175 guez 5
176 guez 7 ierr = NF_REDEF (nid)
177     call nf95_def_var(nid, "bp", NF90_FLOAT, idim_sig, nvarid)
178     call nf95_put_att(nid, nvarid, "title", "Coefficient B pour hybride")
179     call NF95_ENDDEF(nid)
180     call NF95_PUT_VAR(nid, nvarid, bp)
181 guez 5
182 guez 7 ierr = NF_REDEF (nid)
183     call nf95_def_var(nid, "presnivs", NF90_FLOAT, idim_s, nvarid)
184     call NF95_ENDDEF(nid)
185     call NF95_PUT_VAR(nid, nvarid, presnivs)
186 guez 5
187 guez 7 ! Coefficients de passage cov. <-> contra. <--> naturel
188 guez 5
189 guez 7 ierr = NF_REDEF (nid)
190     dims2(1) = idim_rlonu
191     dims2(2) = idim_rlatu
192     call nf95_def_var(nid, "cu", NF90_FLOAT, dims2, nvarid)
193     call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour U")
194     call NF95_ENDDEF(nid)
195     call NF95_PUT_VAR(nid, nvarid, cu_2d)
196 guez 5
197 guez 7 ierr = NF_REDEF (nid)
198     dims2(1) = idim_rlonv
199     dims2(2) = idim_rlatv
200     call nf95_def_var(nid, "cv", NF90_FLOAT, dims2, nvarid)
201     call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour V")
202     call NF95_ENDDEF(nid)
203     call NF95_PUT_VAR(nid, nvarid, cv_2d)
204 guez 5
205 guez 7 ! Aire de chaque maille:
206 guez 5
207 guez 7 ierr = NF_REDEF (nid)
208     dims2(1) = idim_rlonv
209     dims2(2) = idim_rlatu
210     call nf95_def_var(nid, "aire", NF90_FLOAT, dims2, nvarid)
211     call nf95_put_att(nid, nvarid, "title", "Aires de chaque maille")
212     call NF95_ENDDEF(nid)
213     call NF95_PUT_VAR(nid, nvarid, aire_2d)
214 guez 5
215 guez 7 ! Geopentiel au sol:
216 guez 5
217 guez 7 ierr = NF_REDEF (nid)
218     dims2(1) = idim_rlonv
219     dims2(2) = idim_rlatu
220     call nf95_def_var(nid, "phisinit", NF90_FLOAT, dims2, nvarid)
221     call nf95_put_att(nid, nvarid, "title", "Geopotentiel au sol")
222     call NF95_ENDDEF(nid)
223     call NF95_PUT_VAR(nid, nvarid, phis)
224 guez 5
225 guez 7 ! Definir les variables pour pouvoir les enregistrer plus tard:
226 guez 3
227 guez 7 ierr = NF_REDEF (nid) ! entrer dans le mode de definition
228    
229     call nf95_def_var(nid, "temps", NF90_FLOAT, idim_tim, nvarid)
230     call nf95_put_att(nid, nvarid, "title", "Temps de simulation")
231     write(unites, 200)yyears0, mmois0, jjour0
232 guez 5 200 format('days since ', i4, '-', i2.2, '-', i2.2, ' 00:00:00')
233 guez 7 call nf95_put_att(nid, nvarid, "units", unites)
234 guez 5
235    
236 guez 7 dims4(1) = idim_rlonu
237     dims4(2) = idim_rlatu
238     dims4(3) = idim_s
239     dims4(4) = idim_tim
240     call nf95_def_var(nid, "ucov", NF90_FLOAT, dims4, nvarid)
241     call nf95_put_att(nid, nvarid, "title", "Vitesse U")
242 guez 5
243 guez 7 dims4(1) = idim_rlonv
244     dims4(2) = idim_rlatv
245     dims4(3) = idim_s
246     dims4(4) = idim_tim
247     call nf95_def_var(nid, "vcov", NF90_FLOAT, dims4, nvarid)
248     call nf95_put_att(nid, nvarid, "title", "Vitesse V")
249 guez 5
250 guez 7 dims4(1) = idim_rlonv
251     dims4(2) = idim_rlatu
252     dims4(3) = idim_s
253     dims4(4) = idim_tim
254     call nf95_def_var(nid, "teta", NF90_FLOAT, dims4, nvarid)
255     call nf95_put_att(nid, nvarid, "title", "Temperature")
256 guez 5
257 guez 7 dims4(1) = idim_rlonv
258     dims4(2) = idim_rlatu
259     dims4(3) = idim_s
260     dims4(4) = idim_tim
261     DO iq=1, nqmx
262     call nf95_def_var(nid, tname(iq), NF90_FLOAT, dims4, nvarid)
263     call nf95_put_att(nid, nvarid, "title", ttext(iq))
264     ENDDO
265 guez 5
266 guez 7 dims4(1) = idim_rlonv
267     dims4(2) = idim_rlatu
268     dims4(3) = idim_s
269     dims4(4) = idim_tim
270     call nf95_def_var(nid, "masse", NF90_FLOAT, dims4, nvarid)
271     call nf95_put_att(nid, nvarid, "title", "C est quoi ?")
272 guez 5
273 guez 7 dims3(1) = idim_rlonv
274     dims3(2) = idim_rlatu
275     dims3(3) = idim_tim
276     call nf95_def_var(nid, "ps", NF90_FLOAT, dims3, nvarid)
277     call nf95_put_att(nid, nvarid, "title", "Pression au sol")
278 guez 5
279 guez 7 ierr = NF_ENDDEF(nid) ! sortir du mode de definition
280     ierr = NF_CLOSE(nid) ! fermer le fichier
281 guez 3
282 guez 7 PRINT*, 'iim, jjm, llm, iday_end', iim, jjm, llm, iday_end
283     PRINT*, 'rad, omeg, g, cpp, kappa', rad, omeg, g, cpp, kappa
284 guez 3
285 guez 7 END SUBROUTINE dynredem0
286    
287     end module dynredem0_m
  ViewVC Help
Powered by ViewVC 1.1.21