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

Annotation of /trunk/dyn3d/dynetat0.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 130 - (hide annotations)
Tue Feb 24 15:43:51 2015 UTC (9 years, 3 months ago) by guez
File size: 5756 byte(s) 
The information in argument rdayvrai of calfis was redundant with the
information in argument time. Furthermore, in the physics part of gcm,
we need separately the day number (an integer) and the time of
day. So, replaced real argument rdayvrai of calfis containing elapsed
time by integer argument dayvrai containing day number. Corresponding
change in leapfrog. In procedure physiq, replaced real argument
rdayvrai by integer argument dayvrai. In procedures readsulfate and
readsulfate_preind, replaced real argument r_day by arguments dayvrai
and time.

In procedure alboc, replaced real argument rjour by integer argument
jour. alboc was always called by interfsurf_hq with actual argument
real(jour), and the meaning of the dummy argument in alboc seems to be
that it should be an integer.

In procedure leapfrog, local variable time could not be > 1. Removed
test.

In physiq, replaced nint(rdayvrai) by dayvrai. This changes the
results since julien now changes at 0 h instead of 12 h. This follows
LMDZ, where the argument of ozonecm is days_elapsed+1.
1 guez 3 module dynetat0_m
2    
3     IMPLICIT NONE
4    
5 guez 129 INTEGER day_ini
6 guez 130 ! day number at the beginning of the run, based at value 1 on
7     ! January 1st of annee_ref
8 guez 25
9 guez 129 integer:: day_ref = 1 ! jour de l'année de l'état initial
10     ! (= 350 si 20 décembre par exemple)
11    
12     integer:: annee_ref = 1998 ! Annee de l'etat initial (avec 4 chiffres)
13    
14 guez 3 contains
15    
16 guez 128 SUBROUTINE dynetat0(vcov, ucov, teta, q, masse, ps, phis)
17 guez 3
18 guez 38 ! From dynetat0.F, version 1.2, 2004/06/22 11:45:30
19     ! Authors: P. Le Van, L. Fairhead
20 guez 39 ! This procedure reads the initial state of the atmosphere.
21 guez 3
22 guez 79 use comconst, only: dtvr
23 guez 3 use comgeom, only: rlonu, rlatu, rlonv, rlatv, cu_2d, cv_2d, aire_2d
24 guez 129 use conf_gcm_m, only: raz_date
25 guez 38 use dimens_m, only: iim, jjm, llm, nqmx
26 guez 67 use disvert_m, only: pa
27 guez 3 use ener, only: etot0, ang0, ptot0, stot0, ztot0
28 guez 18 use iniadvtrac_m, only: tname
29 guez 67 use netcdf, only: NF90_NOWRITE, NF90_NOERR
30 guez 44 use netcdf95, only: NF95_GET_VAR, nf95_open, nf95_inq_varid, NF95_CLOSE, &
31     NF95_Gw_VAR
32 guez 36 use nr_util, only: assert
33 guez 113 use serre, only: clon, clat, grossismy, grossismx, dzoomx, dzoomy, taux, &
34     tauy
35 guez 129 use temps, only: itau_dyn
36     use unit_nml_m, only: unit_nml
37 guez 3
38 guez 55 REAL, intent(out):: vcov(: , :, :) ! (iim + 1, jjm, llm)
39     REAL, intent(out):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm)
40 guez 43 REAL, intent(out):: teta(:, :, :) ! (iim + 1, jjm + 1, llm)
41 guez 40 REAL, intent(out):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx)
42 guez 55 REAL, intent(out):: masse(:, :, :) ! (iim + 1, jjm + 1, llm)
43 guez 39 REAL, intent(out):: ps(:, :) ! (iim + 1, jjm + 1) in Pa
44     REAL, intent(out):: phis(:, :) ! (iim + 1, jjm + 1)
45 guez 3
46 guez 44 ! Local variables:
47 guez 38 INTEGER iq
48 guez 44 REAL, pointer:: tab_cntrl(:) ! tableau des paramètres du run
49 guez 38 INTEGER ierr, ncid, varid
50 guez 3
51 guez 129 namelist /dynetat0_nml/ day_ref, annee_ref
52    
53 guez 3 !-----------------------------------------------------------------------
54    
55     print *, "Call sequence information: dynetat0"
56    
57 guez 55 call assert((/size(ucov, 1), size(vcov, 1), size(masse, 1), size(ps, 1), &
58     size(phis, 1), size(q, 1), size(teta, 1)/) == iim + 1, "dynetat0 iim")
59     call assert((/size(ucov, 2), size(vcov, 2) + 1, size(masse, 2), &
60     size(ps, 2), size(phis, 2), size(q, 2), size(teta, 2)/) == jjm + 1, &
61     "dynetat0 jjm")
62     call assert((/size(vcov, 3), size(ucov, 3), size(teta, 3), size(q, 3), &
63     size(masse, 3)/) == llm, "dynetat0 llm")
64 guez 40 call assert(size(q, 4) == nqmx, "dynetat0 q nqmx")
65 guez 3
66     ! Fichier état initial :
67 guez 25 call nf95_open("start.nc", NF90_NOWRITE, ncid)
68 guez 3
69 guez 38 call nf95_inq_varid(ncid, "controle", varid)
70 guez 44 call NF95_Gw_VAR(ncid, varid, tab_cntrl)
71 guez 3
72 guez 79 call assert(int(tab_cntrl(1)) == iim, "dynetat0 tab_cntrl iim")
73     call assert(int(tab_cntrl(2)) == jjm, "dynetat0 tab_cntrl jjm")
74     call assert(int(tab_cntrl(3)) == llm, "dynetat0 tab_cntrl llm")
75 guez 3
76 guez 79 IF (dtvr /= tab_cntrl(12)) THEN
77     print *, 'Warning: the time steps from day_step and "start.nc" ' // &
78 guez 129 'are different.'
79 guez 79 print *, 'dtvr from day_step: ', dtvr
80     print *, 'dtvr from "start.nc": ', tab_cntrl(12)
81     print *, 'Using the value from day_step.'
82     ENDIF
83    
84 guez 38 etot0 = tab_cntrl(13)
85     ptot0 = tab_cntrl(14)
86     ztot0 = tab_cntrl(15)
87     stot0 = tab_cntrl(16)
88     ang0 = tab_cntrl(17)
89     pa = tab_cntrl(18)
90 guez 113
91 guez 38 clon = tab_cntrl(20)
92     clat = tab_cntrl(21)
93     grossismx = tab_cntrl(22)
94     grossismy = tab_cntrl(23)
95 guez 113 dzoomx = tab_cntrl(25)
96     dzoomy = tab_cntrl(26)
97     taux = tab_cntrl(28)
98     tauy = tab_cntrl(29)
99    
100 guez 129 print *, "Enter namelist 'dynetat0_nml'."
101     read(unit=*, nml=dynetat0_nml)
102     write(unit_nml, nml=dynetat0_nml)
103 guez 3
104 guez 129 if (raz_date) then
105     print *, 'On réinitialise à la date lue dans la namelist.'
106     day_ini = day_ref
107     itau_dyn = 0
108     else
109     day_ref = tab_cntrl(4)
110     annee_ref = tab_cntrl(5)
111     itau_dyn = tab_cntrl(31)
112     day_ini = tab_cntrl(30)
113     end if
114    
115     print *, "day_ini = ", day_ini
116    
117     deallocate(tab_cntrl) ! pointer
118    
119 guez 38 call NF95_INQ_VARID (ncid, "rlonu", varid)
120     call NF95_GET_VAR(ncid, varid, rlonu)
121 guez 3
122 guez 38 call NF95_INQ_VARID (ncid, "rlatu", varid)
123     call NF95_GET_VAR(ncid, varid, rlatu)
124 guez 3
125 guez 38 call NF95_INQ_VARID (ncid, "rlonv", varid)
126     call NF95_GET_VAR(ncid, varid, rlonv)
127 guez 3
128 guez 38 call NF95_INQ_VARID (ncid, "rlatv", varid)
129     call NF95_GET_VAR(ncid, varid, rlatv)
130 guez 3
131 guez 38 call NF95_INQ_VARID (ncid, "cu", varid)
132     call NF95_GET_VAR(ncid, varid, cu_2d)
133 guez 3
134 guez 38 call NF95_INQ_VARID (ncid, "cv", varid)
135     call NF95_GET_VAR(ncid, varid, cv_2d)
136 guez 3
137 guez 38 call NF95_INQ_VARID (ncid, "aire", varid)
138     call NF95_GET_VAR(ncid, varid, aire_2d)
139 guez 3
140 guez 38 call NF95_INQ_VARID (ncid, "phisinit", varid)
141     call NF95_GET_VAR(ncid, varid, phis)
142 guez 3
143 guez 38 call NF95_INQ_VARID (ncid, "ucov", varid)
144 guez 55 call NF95_GET_VAR(ncid, varid, ucov)
145 guez 3
146 guez 38 call NF95_INQ_VARID (ncid, "vcov", varid)
147 guez 55 call NF95_GET_VAR(ncid, varid, vcov)
148 guez 3
149 guez 38 call NF95_INQ_VARID (ncid, "teta", varid)
150 guez 43 call NF95_GET_VAR(ncid, varid, teta)
151 guez 3
152     DO iq = 1, nqmx
153 guez 38 call NF95_INQ_VARID(ncid, tname(iq), varid, ierr)
154     IF (ierr /= NF90_NOERR) THEN
155     PRINT *, 'dynetat0: "' // tname(iq) // '" not found, ' // &
156     "setting it to zero..."
157 guez 40 q(:, :, :, iq) = 0.
158 guez 3 ELSE
159 guez 40 call NF95_GET_VAR(ncid, varid, q(:, :, :, iq))
160 guez 3 ENDIF
161     ENDDO
162    
163 guez 38 call NF95_INQ_VARID (ncid, "masse", varid)
164 guez 55 call NF95_GET_VAR(ncid, varid, masse)
165 guez 3
166 guez 38 call NF95_INQ_VARID (ncid, "ps", varid)
167 guez 39 call NF95_GET_VAR(ncid, varid, ps)
168 guez 85 ! Check that there is a single value at each pole:
169     call assert(ps(1, 1) == ps(2:, 1), "dynetat0 ps north pole")
170     call assert(ps(1, jjm + 1) == ps(2:, jjm + 1), "dynetat0 ps south pole")
171 guez 3
172 guez 25 call NF95_CLOSE(ncid)
173 guez 3
174     END SUBROUTINE dynetat0
175    
176     end module dynetat0_m
  ViewVC Help
Powered by ViewVC 1.1.21