source: trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F90 @ 937

Last change on this file since 937 was 892, checked in by slebonnois, 12 years ago

SL: Important commit ! Adaptation of Venus physics to parallel computation / template for arch on the LMD servers using ifort / documentation for 1D column physics and for parallel computations
File size: 8.2 KB
Line 
1!
2! $Id $
3!
4subroutine phyetat0(fichnom)
5! Load initial state for the physics
6! and do some resulting initializations
7
8      USE dimphy
9      USE mod_grid_phy_lmdz
10      USE mod_phys_lmdz_para
11      USE iophy
12      USE phys_state_var_mod
13      USE iostart
14      USE infotrac
15      USE comgeomphy,  only: rlatd,rlond
16      USE control_mod, only: raz_date
17
18implicit none
19!======================================================================
20! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
21! Objet: Lecture de l'etat initial pour la physique
22!======================================================================
23#include "dimensions.h"
24#include "netcdf.inc"
25#include "dimsoil.h"
26#include "clesphys.h"
27#include "tabcontrol.h"
28#include "temps.h"
29!======================================================================
30
31character(len=*),intent(in) :: fichnom ! input file name
32REAL    :: xmin, xmax
33LOGICAL :: found
34REAL    :: tab_cntrl(length)
35integer :: i,isoil
36CHARACTER(len=2) :: str2
37
38! les variables globales lues dans le fichier restart
39
40! open physics initial state file:
41call open_startphy(fichnom)
42
43!
44! Lecture des parametres de controle:
45!
46      CALL get_var("controle",tab_cntrl,found)
47      IF (.not.found) THEN
48         PRINT*, 'phyetat0: Le champ <controle> est absent'
49         CALL abort
50      ENDIF
51       
52      DO i = 1, length
53           tabcntr0( i ) = tab_cntrl( i )
54      ENDDO
55
56
57      dtime        = tab_cntrl(1)
58      radpas       = tab_cntrl(2)
59
60      itau_phy = tab_cntrl(15)
61
62! Attention si raz_date est active :
63! il faut remettre a zero itau_phy apres phyetat0 !
64      IF (raz_date.eq.1) THEN
65        itau_phy=0
66      ENDIF
67
68! read latitudes
69call get_field("latitude",rlatd,found)
70      IF (.not.found) THEN
71         PRINT*, 'phyetat0: Le champ <latitude> est absent'
72         CALL abort
73      ENDIF
74
75! read longitudes
76call get_field("longitude",rlond,found)
77      IF (.not.found) THEN
78         PRINT*, 'phyetat0: Le champ <longitude> est absent'
79         CALL abort
80      ENDIF
81
82! read in other variables here ...
83
84! Lecture des temperatures du sol:
85
86       CALL get_field("TS",ftsol(:),found)
87      IF (.not.found) THEN
88         PRINT*, 'phyetat0: Le champ <TS> est absent'
89         PRINT*, "phyetat0: Lecture echouee pour <TS>"
90         CALL abort
91      ELSE
92         PRINT*, 'phyetat0: Le champ <TS> est present'
93         xmin = 1.0E+20
94         xmax = -1.0E+20
95         DO i = 1, klon
96            xmin = MIN(ftsol(i),xmin)
97            xmax = MAX(ftsol(i),xmax)
98         ENDDO
99         PRINT*,'Temperature du sol <TS>', xmin, xmax
100      ENDIF
101
102
103! Lecture des temperatures du sol profond:
104
105      DO isoil=1, nsoilmx
106      IF (isoil.GT.99) THEN
107         PRINT*, "Trop de couches"
108         CALL abort
109      ENDIF
110      WRITE(str2,'(i2.2)') isoil
111      CALL get_field('Tsoil'//str2,ftsoil(:,isoil),found)
112      IF (.not.found) THEN
113         PRINT*, "phyetat0: Le champ <Tsoil"//str2//"> est absent"
114         PRINT*, "          Il prend donc la valeur de surface"
115         DO i=1, klon
116             ftsoil(i,isoil)=ftsol(i)
117         ENDDO
118      ENDIF
119      ENDDO
120
121! Lecture de albedo au sol:
122
123      CALL get_field("ALBE", falbe,found)
124      IF (.not.found) THEN
125         PRINT*, 'phyetat0: Le champ <ALBE> est absent'
126         PRINT*, "phyetat0: Lecture echouee pour <ALBE>"
127         CALL abort
128      ELSE
129         xmin = 1.0E+20
130         xmax = -1.0E+20
131         DO i = 1, klon
132            xmin = MIN(falbe(i),xmin)
133            xmax = MAX(falbe(i),xmax)
134         ENDDO
135         PRINT*,'Albedo du sol <ALBE>', xmin, xmax
136      ENDIF
137
138! Lecture rayonnement solaire au sol:
139
140      CALL get_field("solsw",solsw,found)
141      IF (.not.found) THEN
142         PRINT*, 'phyetat0: Le champ <solsw> est absent'
143         PRINT*, 'mis a zero'
144         solsw = 0.
145      ENDIF
146      xmin = 1.0E+20
147      xmax = -1.0E+20
148      DO i = 1, klon
149         xmin = MIN(solsw(i),xmin)
150         xmax = MAX(solsw(i),xmax)
151      ENDDO
152      PRINT*,'Rayonnement solaire au sol solsw:', xmin, xmax
153
154! Lecture rayonnement IF au sol:
155
156      CALL get_field("sollw",sollw,found)
157      IF (.not.found) THEN
158         PRINT*, 'phyetat0: Le champ <sollw> est absent'
159         PRINT*, 'mis a zero'
160         sollw = 0.
161      ENDIF
162      xmin = 1.0E+20
163      xmax = -1.0E+20
164      DO i = 1, klon
165         xmin = MIN(sollw(i),xmin)
166         xmax = MAX(sollw(i),xmax)
167      ENDDO
168      PRINT*,'Rayonnement IF au sol sollw:', xmin, xmax
169
170! Lecture derive des flux:
171
172      CALL get_field("fder",fder,found)
173      IF (.not.found) THEN
174         PRINT*, 'phyetat0: Le champ <fder> est absent'
175         PRINT*, 'mis a zero'
176         fder = 0.
177      ENDIF
178      xmin = 1.0E+20
179      xmax = -1.0E+20
180      DO i = 1, klon
181         xmin = MIN(fder(i),xmin)
182         xmax = MAX(fder(i),xmax)
183      ENDDO
184      PRINT*,'Derive des flux fder:', xmin, xmax
185
186! Lecture du rayonnement net au sol:
187
188      CALL get_field("RADS",radsol,found)
189      IF (.not.found) THEN
190         PRINT*, 'phyetat0: Le champ <RADS> est absent'
191         CALL abort
192      ENDIF
193      xmin = 1.0E+20
194      xmax = -1.0E+20
195      DO i = 1, klon
196         xmin = MIN(radsol(i),xmin)
197         xmax = MAX(radsol(i),xmax)
198      ENDDO
199      PRINT*,'Rayonnement net au sol radsol:', xmin, xmax
200
201! Lecture de l'orographie sous-maille si ok_orodr:
202
203      if(ok_orodr) then
204     
205      CALL get_field("ZMEA",zmea,found)
206      IF (.not.found) THEN
207         PRINT*, 'phyetat0: Le champ <ZMEA> est absent'
208         CALL abort
209      ENDIF
210      xmin = 1.0E+20
211      xmax = -1.0E+20
212      DO i = 1, klon
213         xmin = MIN(zmea(i),xmin)
214         xmax = MAX(zmea(i),xmax)
215      ENDDO
216      PRINT*,'OROGRAPHIE SOUS-MAILLE zmea:', xmin, xmax
217
218      CALL get_field("ZSTD",zstd,found)
219      IF (.not.found) THEN
220         PRINT*, 'phyetat0: Le champ <ZSTD> est absent'
221         CALL abort
222      ENDIF
223      xmin = 1.0E+20
224      xmax = -1.0E+20
225      DO i = 1, klon
226         xmin = MIN(zstd(i),xmin)
227         xmax = MAX(zstd(i),xmax)
228      ENDDO
229      PRINT*,'OROGRAPHIE SOUS-MAILLE zstd:', xmin, xmax
230
231      CALL get_field("ZSIG",zsig,found)
232      IF (.not.found) THEN
233         PRINT*, 'phyetat0: Le champ <ZSIG> est absent'
234         CALL abort
235      ENDIF
236      xmin = 1.0E+20
237      xmax = -1.0E+20
238      DO i = 1, klon
239         xmin = MIN(zsig(i),xmin)
240         xmax = MAX(zsig(i),xmax)
241      ENDDO
242      PRINT*,'OROGRAPHIE SOUS-MAILLE zsig:', xmin, xmax
243
244      CALL get_field("ZGAM",zgam,found)
245      IF (.not.found) THEN
246         PRINT*, 'phyetat0: Le champ <ZGAM> est absent'
247         CALL abort
248      ENDIF
249      xmin = 1.0E+20
250      xmax = -1.0E+20
251      DO i = 1, klon
252         xmin = MIN(zgam(i),xmin)
253         xmax = MAX(zgam(i),xmax)
254      ENDDO
255      PRINT*,'OROGRAPHIE SOUS-MAILLE zgam:', xmin, xmax
256
257      CALL get_field("ZTHE",zthe,found)
258      IF (.not.found) THEN
259         PRINT*, 'phyetat0: Le champ <ZTHE> est absent'
260         CALL abort
261      ENDIF
262      xmin = 1.0E+20
263      xmax = -1.0E+20
264      DO i = 1, klon
265         xmin = MIN(zthe(i),xmin)
266         xmax = MAX(zthe(i),xmax)
267      ENDDO
268      PRINT*,'OROGRAPHIE SOUS-MAILLE zthe:', xmin, xmax
269
270      CALL get_field("ZPIC",zpic,found)
271      IF (.not.found) THEN
272         PRINT*, 'phyetat0: Le champ <ZPIC> est absent'
273         CALL abort
274      ENDIF
275      xmin = 1.0E+20
276      xmax = -1.0E+20
277      DO i = 1, klon
278         xmin = MIN(zpic(i),xmin)
279         xmax = MAX(zpic(i),xmax)
280      ENDDO
281      PRINT*,'OROGRAPHIE SOUS-MAILLE zpic:', xmin, xmax
282
283      CALL get_field("ZVAL",zval,found)
284      IF (.not.found) THEN
285         PRINT*, 'phyetat0: Le champ <ZVAL> est absent'
286         CALL abort
287      ENDIF
288      xmin = 1.0E+20
289      xmax = -1.0E+20
290      DO i = 1, klon
291         xmin = MIN(zval(i),xmin)
292         xmax = MAX(zval(i),xmax)
293      ENDDO
294      PRINT*,'OROGRAPHIE SOUS-MAILLE zval:', xmin, xmax
295
296      else
297         zmea = 0.
298         zstd = 0.
299         zsig = 0.
300         zgam = 0.
301         zthe = 0.
302         zpic = 0.
303         zval = 0.
304
305      endif   ! fin test sur ok_orodr
306
307! Lecture de TANCIEN:
308
309      ancien_ok = .TRUE.
310
311      CALL get_field("TANCIEN",t_ancien,found)
312      IF (.not.found) THEN
313         PRINT*, "phyetat0: Le champ <TANCIEN> est absent"
314         PRINT*, "Depart legerement fausse. Mais je continue"
315         ancien_ok = .FALSE.
316      ENDIF
317
318! close file
319call close_startphy
320
321! do some more initializations
322call init_iophy_new(rlatd,rlond)
323
324end subroutine phyetat0
Note: See TracBrowser for help on using the repository browser.