/[lmdze]/trunk/Sources/bibio/inithist.f
ViewVC logotype

Annotation of /trunk/Sources/bibio/inithist.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 139 - (hide annotations)
Tue May 26 17:46:03 2015 UTC (9 years ago) by guez
File size: 4109 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 guez 3 module inithist_m
2    
3     implicit none
4    
5     contains
6    
7 guez 129 subroutine inithist(tstep, nq, t_ops, t_wrt)
8 guez 3
9 guez 39 ! From inithist.F, version 1.1.1.1 2004/05/19 12:53:05
10 guez 61 ! L. Fairhead, LMD, 03/99
11 guez 3
12 guez 61 ! Routine d'initialisation des écritures des fichiers histoires
13     ! LMDZ au format IOIPSL.
14 guez 3
15 guez 61 USE com_io_dyn, ONLY: histid, histuid, histvid
16 guez 92 USE dimens_m, ONLY: jjm, llm
17     USE disvert_m, ONLY: presnivs
18 guez 139 use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv
19 guez 61 USE histbeg_totreg_m, ONLY : histbeg_totreg
20     USE histdef_m, ONLY : histdef
21     USE histend_m, ONLY : histend
22     USE histvert_m, ONLY : histvert
23 guez 92 USE iniadvtrac_m, ONLY: ttext
24     USE nr_util, ONLY: pi
25 guez 61 USE paramet_m, ONLY: iip1, jjp1
26     USE temps, ONLY: itau_dyn
27 guez 92 USE ymds2ju_m, ONLY: ymds2ju
28 guez 3
29 guez 61 real, intent(in):: tstep ! durée du pas de temps en secondes
30     integer, intent(in):: nq ! nombre de traceurs
31     real, intent(in):: t_ops ! fréquence de l'opération pour IOIPSL
32     real, intent(in):: t_wrt ! fréquence d'écriture sur le fichier
33 guez 3
34 guez 61 ! Variables locales:
35 guez 3 real zjulian
36     integer iq
37 guez 39 real rlong(iip1, jjp1), rlat(iip1, jjp1)
38 guez 3 integer uhoriid, vhoriid, thoriid, zvertiid
39 guez 39 integer ii, jj
40 guez 3
41     !-----------------------------------------------------------------------
42    
43 guez 129 CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
44 guez 3
45     do jj = 1, jjp1
46     do ii = 1, iip1
47 guez 39 rlong(ii, jj) = rlonu(ii) * 180. / pi
48     rlat(ii, jj) = rlatu(jj) * 180. / pi
49 guez 3 enddo
50     enddo
51    
52 guez 56 call histbeg_totreg("dyn_histu.nc", rlong(:,1), rlat(1,:), 1, iip1, 1, &
53     jjp1, itau_dyn, zjulian, tstep, uhoriid, histuid)
54 guez 3
55     do jj = 1, jjm
56     do ii = 1, iip1
57 guez 39 rlong(ii, jj) = rlonv(ii) * 180. / pi
58     rlat(ii, jj) = rlatv(jj) * 180. / pi
59 guez 3 enddo
60     enddo
61    
62 guez 56 ! Creation du fichier histoire pour la grille en V (oblige pour l'instant,
63     ! IOIPSL ne permet pas de grilles avec des nombres de point differents dans
64     ! un meme fichier)
65    
66 guez 61 call histbeg_totreg('dyn_histv.nc', rlong(:, 1), rlat(1, :jjm), 1, iip1, &
67     1, jjm, itau_dyn, zjulian, tstep, vhoriid, histvid)
68    
69 guez 3 do jj = 1, jjp1
70     do ii = 1, iip1
71 guez 39 rlong(ii, jj) = rlonv(ii) * 180. / pi
72     rlat(ii, jj) = rlatu(jj) * 180. / pi
73 guez 3 enddo
74     enddo
75    
76 guez 56 call histbeg_totreg("dyn_hist.nc", rlong(:, 1), rlat(1, :), 1, iip1, 1, &
77     jjp1, itau_dyn, zjulian, tstep, thoriid, histid)
78 guez 61
79 guez 39 ! Appel a histvert pour la grille verticale
80 guez 61
81 guez 67 call histvert(histid, 'presnivs', 'Niveaux pression','mb', presnivs/100., &
82     zvertiid,'down')
83     call histvert(histvid, 'presnivs', 'Niveaux pression','mb', &
84 guez 61 presnivs/100., zvertiid,'down')
85 guez 67 call histvert(histuid, 'presnivs', 'Niveaux pression','mb', &
86 guez 61 presnivs/100., zvertiid,'down')
87    
88 guez 39 ! Appels a histdef pour la definition des variables a sauvegarder
89 guez 3
90 guez 61 call histdef(histuid, 'u', 'vent u', 'm/s', iip1, jjp1, uhoriid, llm, 1, &
91     llm, zvertiid, 'inst(X)', t_ops, t_wrt)
92     call histdef(histvid, 'v', 'vent v', 'm/s', iip1, jjm, vhoriid, llm, 1, &
93     llm, zvertiid, 'inst(X)', t_ops, t_wrt)
94     call histdef(histid, 'teta', 'temperature potentielle', '-', iip1, jjp1, &
95     thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
96     call histdef(histid, 'phi', 'geopotentiel', '-', iip1, jjp1, thoriid, &
97     llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
98    
99 guez 39 ! Traceurs
100     DO iq=1, nq
101 guez 61 call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &
102     llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
103 guez 3 enddo
104 guez 61
105     call histdef(histid, 'masse', 'masse', 'kg', iip1, jjp1, thoriid, llm, 1, &
106     llm, zvertiid, 'inst(X)', t_ops, t_wrt)
107     call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', iip1, jjp1, &
108     thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
109     call histdef(histid, 'phis', 'geopotentiel au sol', '-', iip1, jjp1, &
110     thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
111    
112 guez 56 call histend(histid)
113     call histend(histuid)
114     call histend(histvid)
115 guez 3
116     end subroutine inithist
117    
118     end module inithist_m
  ViewVC Help
Powered by ViewVC 1.1.21