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

Annotation of /trunk/dyn3d/dynetat0.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 38 - (hide annotations)
Thu Jan 6 17:52:19 2011 UTC (13 years, 4 months ago) by guez
Original Path: trunk/libf/dyn3d/dynetat0.f90
File size: 4755 byte(s) 
Extracted ASCII art from "inigeom" into a separate text file in the
documentation.

"test_disvert" now creates a separate file for layer thicknesses.

Moved variables from module "yomcst" to module "suphec_m" because this
is where those variables are defined. Kept in "yomcst" only parameters
of Earth orbit. Gave the attribute "parameter" to some variables of
module "suphec_m".

Variables of module "yoethf" were defined in procedure "suphec". Moved
these definitions to a new procedure "yoethf" in module "yoethf_m".
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 3 ! Objet : lecture de l'état initial
14    
15 guez 38 use comconst, only: im, dtvr, jm, lllm, omeg
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     REAL, intent(out):: q(:, :, :), masse(:, :)
31 guez 10 REAL, intent(out):: ps(:) ! in Pa
32     REAL, intent(out):: phis(:, :)
33 guez 25 REAL, intent(out):: time_0
34 guez 3
35 guez 38 ! Variables
36     INTEGER iq
37     INTEGER, PARAMETER:: length = 100
38     REAL tab_cntrl(length) ! tableau des paramètres du run
39     INTEGER ierr, ncid, varid
40 guez 3
41     !-----------------------------------------------------------------------
42    
43     print *, "Call sequence information: dynetat0"
44    
45     call assert(size(vcov, 1) == (iim + 1) * jjm, "dynetat0 vcov 1")
46     call assert((/size(ucov, 1), size(teta, 1), size(q, 1), size(masse, 1), &
47 guez 7 size(ps)/) == (iim + 1) * (jjm + 1), "dynetat0 (iim + 1) * (jjm + 1)")
48     call assert(shape(phis) == (/iim + 1, jjm + 1/), "dynetat0 phis")
49 guez 3 call assert((/size(vcov, 2), size(ucov, 2), size(teta, 2), size(q, 2), &
50     size(masse, 2)/) == llm, "dynetat0 llm")
51     call assert(size(q, 3) == nqmx, "dynetat0 q 3")
52    
53     ! Fichier état initial :
54 guez 25 call nf95_open("start.nc", NF90_NOWRITE, ncid)
55 guez 3
56 guez 38 call nf95_inq_varid(ncid, "controle", varid)
57     call NF95_GET_VAR(ncid, varid, tab_cntrl)
58 guez 3
59 guez 38 im = int(tab_cntrl(1))
60     jm = int(tab_cntrl(2))
61     lllm = int(tab_cntrl(3))
62     call assert(im == iim, "dynetat0 im iim")
63     call assert(jm == jjm, "dynetat0 jm jjm")
64     call assert(lllm == llm, "dynetat0 lllm llm")
65 guez 3
66 guez 38 day_ref = int(tab_cntrl(4))
67     annee_ref = int(tab_cntrl(5))
68 guez 3
69 guez 38 omeg = tab_cntrl(7)
70     PRINT *, 'omeg = ', omeg
71    
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     call NF95_GET_VAR(ncid, varid, ucov, count_nc=(/iim + 1, jjm + 1, llm/))
120 guez 3
121 guez 38 call NF95_INQ_VARID (ncid, "vcov", varid)
122     call NF95_GET_VAR(ncid, varid, vcov, count_nc=(/iim + 1, jjm, llm/))
123 guez 3
124 guez 38 call NF95_INQ_VARID (ncid, "teta", varid)
125     call NF95_GET_VAR(ncid, varid, teta, count_nc=(/iim + 1, jjm + 1, llm/))
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 3 q(:, :, iq) = 0.
133     ELSE
134 guez 38 call NF95_GET_VAR(ncid, varid, q(:, :, iq), &
135     count_nc=(/iim + 1, jjm + 1, llm/))
136 guez 3 ENDIF
137     ENDDO
138    
139 guez 38 call NF95_INQ_VARID (ncid, "masse", varid)
140     call NF95_GET_VAR(ncid, varid, masse, count_nc=(/iim + 1, jjm + 1, llm/))
141 guez 3
142 guez 38 call NF95_INQ_VARID (ncid, "ps", varid)
143     call NF95_GET_VAR(ncid, varid, ps, count_nc=(/iim + 1, jjm + 1/))
144 guez 3
145 guez 25 call NF95_CLOSE(ncid)
146 guez 3
147     END SUBROUTINE dynetat0
148    
149     end module dynetat0_m
  ViewVC Help
Powered by ViewVC 1.1.21