/[lmdze]/trunk/Sources/dyn3d/dynetat0.f
ViewVC logotype

Annotation of /trunk/Sources/dyn3d/dynetat0.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 139 - (hide annotations)
Tue May 26 17:46:03 2015 UTC (8 years, 11 months ago) by guez
File size: 7166 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 dynetat0_m
2    
3 guez 139 use dimens_m, only: iim, jjm
4    
5 guez 3 IMPLICIT NONE
6    
7 guez 139 private iim, jjm
8    
9 guez 129 INTEGER day_ini
10 guez 130 ! day number at the beginning of the run, based at value 1 on
11     ! January 1st of annee_ref
12 guez 25
13 guez 139 integer:: day_ref = 1 ! jour de l'ann\'ee de l'\'etat initial
14     ! (= 350 si 20 d\'ecembre par exemple)
15 guez 129
16     integer:: annee_ref = 1998 ! Annee de l'etat initial (avec 4 chiffres)
17    
18 guez 139 REAL clon ! longitude of the center of the zoom, in rad
19     real clat ! latitude of the center of the zoom, in rad
20    
21     real grossismx, grossismy
22     ! facteurs de grossissement du zoom, selon la longitude et la latitude
23     ! = 2 si 2 fois, = 3 si 3 fois, etc.
24    
25     real dzoomx, dzoomy
26     ! extensions en longitude et latitude de la zone du zoom (fractions
27     ! de la zone totale)
28    
29     real taux, tauy
30     ! raideur de la transition de l'int\'erieur à l'ext\'erieur du zoom
31    
32     real rlatu(jjm + 1)
33     ! (latitudes of points of the "scalar" and "u" grid, in rad)
34    
35     real rlatv(jjm)
36     ! (latitudes of points of the "v" grid, in rad, in decreasing order)
37    
38     real rlonu(iim + 1) ! longitudes of points of the "u" grid, in rad
39    
40     real rlonv(iim + 1)
41     ! (longitudes of points of the "scalar" and "v" grid, in rad)
42    
43     real xprimu(iim + 1), xprimv(iim + 1)
44     ! xprimu et xprimv sont respectivement les valeurs de dx / dX aux
45     ! points u et v.
46    
47     REAL xprimm025(iim + 1), xprimp025(iim + 1)
48     REAL rlatu1(jjm), rlatu2(jjm), yprimu1(jjm), yprimu2(jjm)
49    
50     save
51    
52 guez 3 contains
53    
54 guez 128 SUBROUTINE dynetat0(vcov, ucov, teta, q, masse, ps, phis)
55 guez 3
56 guez 38 ! From dynetat0.F, version 1.2, 2004/06/22 11:45:30
57     ! Authors: P. Le Van, L. Fairhead
58 guez 39 ! This procedure reads the initial state of the atmosphere.
59 guez 3
60 guez 79 use comconst, only: dtvr
61 guez 129 use conf_gcm_m, only: raz_date
62 guez 38 use dimens_m, only: iim, jjm, llm, nqmx
63 guez 67 use disvert_m, only: pa
64 guez 3 use ener, only: etot0, ang0, ptot0, stot0, ztot0
65 guez 18 use iniadvtrac_m, only: tname
66 guez 67 use netcdf, only: NF90_NOWRITE, NF90_NOERR
67 guez 44 use netcdf95, only: NF95_GET_VAR, nf95_open, nf95_inq_varid, NF95_CLOSE, &
68     NF95_Gw_VAR
69 guez 36 use nr_util, only: assert
70 guez 129 use temps, only: itau_dyn
71     use unit_nml_m, only: unit_nml
72 guez 3
73 guez 55 REAL, intent(out):: vcov(: , :, :) ! (iim + 1, jjm, llm)
74     REAL, intent(out):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm)
75 guez 43 REAL, intent(out):: teta(:, :, :) ! (iim + 1, jjm + 1, llm)
76 guez 40 REAL, intent(out):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx)
77 guez 55 REAL, intent(out):: masse(:, :, :) ! (iim + 1, jjm + 1, llm)
78 guez 39 REAL, intent(out):: ps(:, :) ! (iim + 1, jjm + 1) in Pa
79     REAL, intent(out):: phis(:, :) ! (iim + 1, jjm + 1)
80 guez 3
81 guez 44 ! Local variables:
82 guez 38 INTEGER iq
83 guez 139 REAL, pointer:: tab_cntrl(:) ! tableau des param\`etres du run
84 guez 38 INTEGER ierr, ncid, varid
85 guez 3
86 guez 129 namelist /dynetat0_nml/ day_ref, annee_ref
87    
88 guez 3 !-----------------------------------------------------------------------
89    
90     print *, "Call sequence information: dynetat0"
91    
92 guez 55 call assert((/size(ucov, 1), size(vcov, 1), size(masse, 1), size(ps, 1), &
93     size(phis, 1), size(q, 1), size(teta, 1)/) == iim + 1, "dynetat0 iim")
94     call assert((/size(ucov, 2), size(vcov, 2) + 1, size(masse, 2), &
95     size(ps, 2), size(phis, 2), size(q, 2), size(teta, 2)/) == jjm + 1, &
96     "dynetat0 jjm")
97     call assert((/size(vcov, 3), size(ucov, 3), size(teta, 3), size(q, 3), &
98     size(masse, 3)/) == llm, "dynetat0 llm")
99 guez 40 call assert(size(q, 4) == nqmx, "dynetat0 q nqmx")
100 guez 3
101 guez 139 ! Fichier \'etat initial :
102 guez 25 call nf95_open("start.nc", NF90_NOWRITE, ncid)
103 guez 3
104 guez 38 call nf95_inq_varid(ncid, "controle", varid)
105 guez 44 call NF95_Gw_VAR(ncid, varid, tab_cntrl)
106 guez 3
107 guez 79 call assert(int(tab_cntrl(1)) == iim, "dynetat0 tab_cntrl iim")
108     call assert(int(tab_cntrl(2)) == jjm, "dynetat0 tab_cntrl jjm")
109     call assert(int(tab_cntrl(3)) == llm, "dynetat0 tab_cntrl llm")
110 guez 3
111 guez 79 IF (dtvr /= tab_cntrl(12)) THEN
112     print *, 'Warning: the time steps from day_step and "start.nc" ' // &
113 guez 129 'are different.'
114 guez 79 print *, 'dtvr from day_step: ', dtvr
115     print *, 'dtvr from "start.nc": ', tab_cntrl(12)
116     print *, 'Using the value from day_step.'
117     ENDIF
118    
119 guez 38 etot0 = tab_cntrl(13)
120     ptot0 = tab_cntrl(14)
121     ztot0 = tab_cntrl(15)
122     stot0 = tab_cntrl(16)
123     ang0 = tab_cntrl(17)
124     pa = tab_cntrl(18)
125 guez 113
126 guez 38 clon = tab_cntrl(20)
127     clat = tab_cntrl(21)
128     grossismx = tab_cntrl(22)
129     grossismy = tab_cntrl(23)
130 guez 113 dzoomx = tab_cntrl(25)
131     dzoomy = tab_cntrl(26)
132     taux = tab_cntrl(28)
133     tauy = tab_cntrl(29)
134    
135 guez 129 print *, "Enter namelist 'dynetat0_nml'."
136     read(unit=*, nml=dynetat0_nml)
137     write(unit_nml, nml=dynetat0_nml)
138 guez 3
139 guez 129 if (raz_date) then
140 guez 139 print *, 'Resetting the date, using the namelist.'
141 guez 129 day_ini = day_ref
142     itau_dyn = 0
143     else
144     day_ref = tab_cntrl(4)
145     annee_ref = tab_cntrl(5)
146     itau_dyn = tab_cntrl(31)
147     day_ini = tab_cntrl(30)
148     end if
149    
150     print *, "day_ini = ", day_ini
151    
152     deallocate(tab_cntrl) ! pointer
153    
154 guez 38 call NF95_INQ_VARID (ncid, "rlonu", varid)
155     call NF95_GET_VAR(ncid, varid, rlonu)
156 guez 3
157 guez 38 call NF95_INQ_VARID (ncid, "rlatu", varid)
158     call NF95_GET_VAR(ncid, varid, rlatu)
159 guez 3
160 guez 38 call NF95_INQ_VARID (ncid, "rlonv", varid)
161     call NF95_GET_VAR(ncid, varid, rlonv)
162 guez 3
163 guez 38 call NF95_INQ_VARID (ncid, "rlatv", varid)
164     call NF95_GET_VAR(ncid, varid, rlatv)
165 guez 3
166 guez 139 CALL nf95_inq_varid(ncid, 'xprimu', varid)
167     CALL nf95_get_var(ncid, varid, xprimu)
168 guez 3
169 guez 139 CALL nf95_inq_varid(ncid, 'xprimv', varid)
170     CALL nf95_get_var(ncid, varid, xprimv)
171 guez 3
172 guez 139 CALL nf95_inq_varid(ncid, 'xprimm025', varid)
173     CALL nf95_get_var(ncid, varid, xprimm025)
174 guez 3
175 guez 139 CALL nf95_inq_varid(ncid, 'xprimp025', varid)
176     CALL nf95_get_var(ncid, varid, xprimp025)
177    
178     call NF95_INQ_VARID (ncid, "rlatu1", varid)
179     call NF95_GET_VAR(ncid, varid, rlatu1)
180    
181     call NF95_INQ_VARID (ncid, "rlatu2", varid)
182     call NF95_GET_VAR(ncid, varid, rlatu2)
183    
184     CALL nf95_inq_varid(ncid, 'yprimu1', varid)
185     CALL nf95_get_var(ncid, varid, yprimu1)
186    
187     CALL nf95_inq_varid(ncid, 'yprimu2', varid)
188     CALL nf95_get_var(ncid, varid, yprimu2)
189    
190 guez 38 call NF95_INQ_VARID (ncid, "phisinit", varid)
191     call NF95_GET_VAR(ncid, varid, phis)
192 guez 3
193 guez 38 call NF95_INQ_VARID (ncid, "ucov", varid)
194 guez 55 call NF95_GET_VAR(ncid, varid, ucov)
195 guez 3
196 guez 38 call NF95_INQ_VARID (ncid, "vcov", varid)
197 guez 55 call NF95_GET_VAR(ncid, varid, vcov)
198 guez 3
199 guez 38 call NF95_INQ_VARID (ncid, "teta", varid)
200 guez 43 call NF95_GET_VAR(ncid, varid, teta)
201 guez 3
202     DO iq = 1, nqmx
203 guez 38 call NF95_INQ_VARID(ncid, tname(iq), varid, ierr)
204     IF (ierr /= NF90_NOERR) THEN
205     PRINT *, 'dynetat0: "' // tname(iq) // '" not found, ' // &
206     "setting it to zero..."
207 guez 40 q(:, :, :, iq) = 0.
208 guez 3 ELSE
209 guez 40 call NF95_GET_VAR(ncid, varid, q(:, :, :, iq))
210 guez 3 ENDIF
211     ENDDO
212    
213 guez 38 call NF95_INQ_VARID (ncid, "masse", varid)
214 guez 55 call NF95_GET_VAR(ncid, varid, masse)
215 guez 3
216 guez 38 call NF95_INQ_VARID (ncid, "ps", varid)
217 guez 39 call NF95_GET_VAR(ncid, varid, ps)
218 guez 85 ! Check that there is a single value at each pole:
219     call assert(ps(1, 1) == ps(2:, 1), "dynetat0 ps north pole")
220     call assert(ps(1, jjm + 1) == ps(2:, jjm + 1), "dynetat0 ps south pole")
221 guez 3
222 guez 25 call NF95_CLOSE(ncid)
223 guez 3
224     END SUBROUTINE dynetat0
225    
226     end module dynetat0_m
  ViewVC Help
Powered by ViewVC 1.1.21