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

Annotation of /trunk/Sources/bibio/initfluxsto.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: 6211 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 31 SUBROUTINE initfluxsto(tstep, t_ops, t_wrt, nq, fileid, filevid, filedid)
2 guez 3
3 guez 31 ! From bibio/initfluxsto.F, v 1.1.1.1 2004/05/19 12:53:05
4 guez 3
5 guez 31 ! Routine d'initialisation des ecritures des fichiers histoires LMDZ
6     ! au format IOIPSL
7     ! Appels succesifs des routines: histbeg
8     ! histhori
9     ! histver
10     ! histdef
11     ! histend
12 guez 3
13 guez 31 ! Entree:
14     ! day0, anne0: date de reference
15     ! tstep: duree du pas de temps en seconde
16     ! t_ops: frequence de l'operation pour IOIPSL
17     ! t_wrt: frequence d'ecriture sur le fichier
18     ! nq: nombre de traceurs
19 guez 3
20 guez 31 ! Sortie:
21     ! fileid: ID du fichier netcdf cree
22     ! filevid:ID du fichier netcdf pour la grille v
23 guez 3
24 guez 31 ! L. Fairhead, 03/99
25 guez 3
26 guez 92 USE comconst
27     use conf_gcm_m
28     USE dimens_m
29     USE disvert_m
30 guez 139 use dynetat0_m, only: day_ref, annee_ref, rlonu, rlatu, rlonv, rlatv
31 guez 61 USE histbeg_totreg_m, ONLY : histbeg_totreg
32     USE histdef_m, ONLY : histdef
33     USE histend_m, ONLY : histend
34     use histhori_regular_m, only: histhori_regular
35     use histsync_m, only: histsync
36     USE histvert_m, ONLY : histvert
37 guez 92 USE nr_util, ONLY : pi
38 guez 31 USE paramet_m
39 guez 129 USE temps, ONLY : itau_dyn
40 guez 92 use ymds2ju_m, only: ymds2ju
41 guez 3
42 guez 31 IMPLICIT NONE
43 guez 3
44 guez 31 ! Arguments
45     INTEGER itau
46     REAL, INTENT (IN) :: tstep
47     REAL t_ops, t_wrt
48     INTEGER fileid, filevid, filedid
49     INTEGER nq, ndex(1)
50     REAL nivd(1)
51 guez 3
52 guez 31 ! Variables locales
53     REAL zjulian
54 guez 108 CHARACTER(len=3) str
55     CHARACTER(len=10) ctrac
56 guez 31 INTEGER iq
57     REAL rlong(iip1, jjp1), rlat(iip1, jjp1)
58     INTEGER uhoriid, vhoriid, thoriid, zvertiid, dhoriid, dvertiid
59 guez 83 INTEGER ii, jj, l
60 guez 31 LOGICAL ok_sync
61 guez 3
62 guez 31 !---------------------------------------------------------
63 guez 3
64 guez 31 ! Initialisations
65     str = 'q '
66     ctrac = 'traceur '
67     ok_sync = .TRUE.
68 guez 3
69 guez 31 ! Appel a histbeg: creation du fichier netcdf et initialisations diverses
70 guez 3
71 guez 129 CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
72 guez 3
73 guez 31 DO jj = 1, jjp1
74     DO ii = 1, iip1
75     rlong(ii, jj) = rlonu(ii)*180./pi
76     rlat(ii, jj) = rlatu(jj)*180./pi
77     END DO
78     END DO
79 guez 3
80 guez 31 CALL histbeg_totreg('fluxstoke', rlong(:, 1), rlat(1, :), 1, iip1, 1, jjp1, &
81     itau_dyn, zjulian, tstep, uhoriid, fileid)
82    
83     ! Creation du fichier histoire pour la grille en V (oblige pour l'instant,
84     ! IOIPSL ne permet pas de grilles avec des nombres de point differents dans
85     ! un meme fichier)
86    
87     DO jj = 1, jjm
88     DO ii = 1, iip1
89     rlong(ii, jj) = rlonv(ii)*180./pi
90     rlat(ii, jj) = rlatv(jj)*180./pi
91     END DO
92     END DO
93    
94     CALL histbeg_totreg('fluxstokev.nc', rlong(:, 1), rlat(1, :jjm), 1, iip1, &
95     1, jjm, itau_dyn, zjulian, tstep, vhoriid, filevid)
96    
97     CALL histbeg_totreg('defstoke.nc', (/1./), (/1./), 1, 1, 1, 1, itau_dyn, &
98     zjulian, tstep, dhoriid, filedid)
99    
100     ! Appel a histhori pour rajouter les autres grilles horizontales
101    
102     DO jj = 1, jjp1
103     DO ii = 1, iip1
104     rlong(ii, jj) = rlonv(ii)*180./pi
105     rlat(ii, jj) = rlatu(jj)*180./pi
106     END DO
107     END DO
108    
109     CALL histhori_regular(fileid, iip1, rlong, jjp1, rlat, 'scalar', &
110     'Grille points scalaires', thoriid)
111    
112     ! Appel a histvert pour la grille verticale
113    
114 guez 83 CALL histvert(fileid, 'sig_s', 'Niveaux sigma', 'sigma_level', &
115     (/(real(l), l = 1, llm)/), zvertiid)
116 guez 31 ! Pour le fichier V
117 guez 83 CALL histvert(filevid, 'sig_s', 'Niveaux sigma', 'sigma_level', &
118     (/(real(l), l = 1, llm)/), zvertiid)
119 guez 31 ! pour le fichier def
120     nivd(1) = 1
121 guez 67 CALL histvert(filedid, 'sig_s', 'Niveaux sigma', 'sigma_level', nivd, &
122 guez 31 dvertiid)
123    
124     ! Appels a histdef pour la definition des variables a sauvegarder
125     CALL histdef(fileid, 'phis', 'Surface geop. height', '-', iip1, jjp1, &
126     thoriid, 1, 1, 1, -99, 'once', t_ops, t_wrt)
127     CALL histdef(fileid, 'aire', 'Grid area', '-', iip1, jjp1, thoriid, 1, 1, &
128     1, -99, 'once', t_ops, t_wrt)
129     CALL histdef(filedid, 'dtvr', 'tps dyn', 's', 1, 1, dhoriid, 1, 1, 1, -99, &
130     'once', t_ops, t_wrt)
131     CALL histdef(filedid, 'istdyn', 'tps stock', 's', 1, 1, dhoriid, 1, 1, 1, &
132     -99, 'once', t_ops, t_wrt)
133     CALL histdef(filedid, 'istphy', 'tps stock phy', 's', 1, 1, dhoriid, 1, 1, &
134     1, -99, 'once', t_ops, t_wrt)
135     CALL histdef(fileid, 'masse', 'Masse', 'kg', iip1, jjp1, thoriid, llm, 1, &
136     llm, zvertiid, 'inst(X)', t_ops, t_wrt)
137     CALL histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s', iip1, jjp1, &
138     uhoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
139     CALL histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s', iip1, jjm, &
140     vhoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
141     CALL histdef(fileid, 'w', 'flx de masse vert', 'kg m/s', iip1, jjp1, &
142     thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
143     CALL histdef(fileid, 'teta', 'temperature potentielle', '-', iip1, jjp1, &
144     thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
145     CALL histdef(fileid, 'phi', 'geopotentiel instantane', '-', iip1, jjp1, &
146     thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
147    
148     CALL histend(fileid)
149     CALL histend(filevid)
150     CALL histend(filedid)
151     IF (ok_sync) THEN
152     CALL histsync(fileid)
153     CALL histsync(filevid)
154     CALL histsync(filedid)
155     END IF
156    
157     END SUBROUTINE initfluxsto
  ViewVC Help
Powered by ViewVC 1.1.21