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

Annotation of /trunk/dyn3d/dynetat0.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 40 - (hide annotations)
Tue Feb 22 13:49:36 2011 UTC (13 years, 3 months ago) by guez
Original Path: trunk/libf/dyn3d/dynetat0.f90
File size: 4860 byte(s) 
"alpha" useless, always 0, in "exner_hyb".
1 guez 3 module dynetat0_m
2    
3     IMPLICIT NONE
4    
5 guez 25 INTEGER day_ini
6    
7 guez 3 contains
8    
9 guez 25 SUBROUTINE dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)
10 guez 3
11 guez 38 ! From dynetat0.F, version 1.2, 2004/06/22 11:45:30
12     ! Authors: P. Le Van, L. Fairhead
13 guez 39 ! This procedure reads the initial state of the atmosphere.
14 guez 3
15 guez 39 use comconst, only: im, dtvr, jm, lllm
16 guez 3 use comvert, only: pa
17     use comgeom, only: rlonu, rlatu, rlonv, rlatv, cu_2d, cv_2d, aire_2d
18 guez 38 use dimens_m, only: iim, jjm, llm, nqmx
19 guez 3 use ener, only: etot0, ang0, ptot0, stot0, ztot0
20 guez 18 use iniadvtrac_m, only: tname
21 guez 38 use logic, only: fxyhypb, ysinus
22     use serre, only: clon, clat, grossismy, grossismx
23     use netcdf95, only: NF95_GET_VAR, nf95_open, nf95_inq_varid, NF95_CLOSE
24     use netcdf, only: NF90_NOWRITE, NF90_NOERR
25 guez 36 use nr_util, only: assert
26 guez 38 use temps, only: day_ref, itau_dyn, annee_ref
27 guez 3
28 guez 38 ! Arguments:
29 guez 3 REAL, intent(out):: vcov(: , :), ucov(:, :), teta(:, :)
30 guez 40 REAL, intent(out):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx)
31     REAL, intent(out):: masse(:, :)
32 guez 39 REAL, intent(out):: ps(:, :) ! (iim + 1, jjm + 1) in Pa
33     REAL, intent(out):: phis(:, :) ! (iim + 1, jjm + 1)
34 guez 25 REAL, intent(out):: time_0
35 guez 3
36 guez 38 ! Variables
37     INTEGER iq
38     INTEGER, PARAMETER:: length = 100
39     REAL tab_cntrl(length) ! tableau des paramètres du run
40     INTEGER ierr, ncid, varid
41 guez 3
42     !-----------------------------------------------------------------------
43    
44     print *, "Call sequence information: dynetat0"
45    
46     call assert(size(vcov, 1) == (iim + 1) * jjm, "dynetat0 vcov 1")
47 guez 40 call assert((/size(ucov, 1), size(teta, 1), size(masse, 1)/) &
48 guez 39 == (iim + 1) * (jjm + 1), "dynetat0 (iim + 1) * (jjm + 1)")
49 guez 40 call assert((/size(ps, 1), size(phis, 1), size(q, 1)/) == iim + 1, &
50     "dynetat0 iim")
51     call assert((/size(ps, 2), size(phis, 2), size(q, 2)/) == jjm + 1, &
52     "dynetat0 jjm")
53     call assert((/size(vcov, 2), size(ucov, 2), size(teta, 2), size(q, 3), &
54 guez 3 size(masse, 2)/) == llm, "dynetat0 llm")
55 guez 40 call assert(size(q, 4) == nqmx, "dynetat0 q nqmx")
56 guez 3
57     ! Fichier état initial :
58 guez 25 call nf95_open("start.nc", NF90_NOWRITE, ncid)
59 guez 3
60 guez 38 call nf95_inq_varid(ncid, "controle", varid)
61     call NF95_GET_VAR(ncid, varid, tab_cntrl)
62 guez 3
63 guez 38 im = int(tab_cntrl(1))
64     jm = int(tab_cntrl(2))
65     lllm = int(tab_cntrl(3))
66     call assert(im == iim, "dynetat0 im iim")
67     call assert(jm == jjm, "dynetat0 jm jjm")
68     call assert(lllm == llm, "dynetat0 lllm llm")
69 guez 3
70 guez 38 day_ref = int(tab_cntrl(4))
71     annee_ref = int(tab_cntrl(5))
72     dtvr = tab_cntrl(12)
73     etot0 = tab_cntrl(13)
74     ptot0 = tab_cntrl(14)
75     ztot0 = tab_cntrl(15)
76     stot0 = tab_cntrl(16)
77     ang0 = tab_cntrl(17)
78     pa = tab_cntrl(18)
79     clon = tab_cntrl(20)
80     clat = tab_cntrl(21)
81     grossismx = tab_cntrl(22)
82     grossismy = tab_cntrl(23)
83     fxyhypb = tab_cntrl(24) == 1.
84     if (.not. fxyhypb) ysinus = tab_cntrl(27) == 1.
85 guez 3 itau_dyn = tab_cntrl(31)
86    
87 guez 38 call NF95_INQ_VARID (ncid, "rlonu", varid)
88     call NF95_GET_VAR(ncid, varid, rlonu)
89 guez 3
90 guez 38 call NF95_INQ_VARID (ncid, "rlatu", varid)
91     call NF95_GET_VAR(ncid, varid, rlatu)
92 guez 3
93 guez 38 call NF95_INQ_VARID (ncid, "rlonv", varid)
94     call NF95_GET_VAR(ncid, varid, rlonv)
95 guez 3
96 guez 38 call NF95_INQ_VARID (ncid, "rlatv", varid)
97     call NF95_GET_VAR(ncid, varid, rlatv)
98 guez 3
99 guez 38 call NF95_INQ_VARID (ncid, "cu", varid)
100     call NF95_GET_VAR(ncid, varid, cu_2d)
101 guez 3
102 guez 38 call NF95_INQ_VARID (ncid, "cv", varid)
103     call NF95_GET_VAR(ncid, varid, cv_2d)
104 guez 3
105 guez 38 call NF95_INQ_VARID (ncid, "aire", varid)
106     call NF95_GET_VAR(ncid, varid, aire_2d)
107 guez 3
108 guez 38 call NF95_INQ_VARID (ncid, "phisinit", varid)
109     call NF95_GET_VAR(ncid, varid, phis)
110 guez 3
111 guez 38 call NF95_INQ_VARID (ncid, "temps", varid)
112     call NF95_GET_VAR(ncid, varid, time_0)
113 guez 3
114 guez 38 day_ini = tab_cntrl(30) + INT(time_0)
115     time_0 = time_0 - INT(time_0)
116     ! {0 <= time0 < 1}
117 guez 3
118 guez 38 call NF95_INQ_VARID (ncid, "ucov", varid)
119 guez 40 call NF95_GET_VAR(ncid, varid, ucov, count_nc=(/iim + 1, jjm + 1, llm, 1/))
120 guez 3
121 guez 38 call NF95_INQ_VARID (ncid, "vcov", varid)
122 guez 40 call NF95_GET_VAR(ncid, varid, vcov, count_nc=(/iim + 1, jjm, llm, 1/))
123 guez 3
124 guez 38 call NF95_INQ_VARID (ncid, "teta", varid)
125 guez 40 call NF95_GET_VAR(ncid, varid, teta, count_nc=(/iim + 1, jjm + 1, llm, 1/))
126 guez 3
127     DO iq = 1, nqmx
128 guez 38 call NF95_INQ_VARID(ncid, tname(iq), varid, ierr)
129     IF (ierr /= NF90_NOERR) THEN
130     PRINT *, 'dynetat0: "' // tname(iq) // '" not found, ' // &
131     "setting it to zero..."
132 guez 40 q(:, :, :, iq) = 0.
133 guez 3 ELSE
134 guez 40 call NF95_GET_VAR(ncid, varid, q(:, :, :, iq))
135 guez 3 ENDIF
136     ENDDO
137    
138 guez 38 call NF95_INQ_VARID (ncid, "masse", varid)
139 guez 40 call NF95_GET_VAR(ncid, varid, masse, count_nc=(/iim + 1, jjm + 1, llm, 1/))
140 guez 3
141 guez 38 call NF95_INQ_VARID (ncid, "ps", varid)
142 guez 39 call NF95_GET_VAR(ncid, varid, ps)
143 guez 3
144 guez 25 call NF95_CLOSE(ncid)
145 guez 3
146     END SUBROUTINE dynetat0
147    
148     end module dynetat0_m
  ViewVC Help
Powered by ViewVC 1.1.21