1 |
guez |
3 |
! |
2 |
|
|
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/iniphysiq.F,v 1.1.1.1 2004/05/19 12:53:08 lmdzadmin Exp $ |
3 |
|
|
! |
4 |
|
|
c |
5 |
|
|
c |
6 |
|
|
SUBROUTINE iniphysiq(ngrid,nlayer, |
7 |
|
|
$ punjours, |
8 |
|
|
$ pdayref,ptimestep, |
9 |
|
|
$ plat,plon,parea,pcu,pcv, |
10 |
|
|
$ prad,pg,pr,pcpp) |
11 |
|
|
use dimens_m |
12 |
|
|
use dimphy |
13 |
|
|
use comgeomphy |
14 |
|
|
use suphec_m, only: suphec |
15 |
|
|
|
16 |
|
|
IMPLICIT NONE |
17 |
|
|
c |
18 |
|
|
c======================================================================= |
19 |
|
|
c |
20 |
|
|
c subject: |
21 |
|
|
c -------- |
22 |
|
|
c |
23 |
|
|
c Initialisation for the physical parametrisations of the LMD |
24 |
|
|
c martian atmospheric general circulation modele. |
25 |
|
|
c |
26 |
|
|
c author: Frederic Hourdin 15 / 10 /93 |
27 |
|
|
c ------- |
28 |
|
|
c |
29 |
|
|
c arguments: |
30 |
|
|
c ---------- |
31 |
|
|
c |
32 |
|
|
c input: |
33 |
|
|
c ------ |
34 |
|
|
c |
35 |
|
|
c ngrid Size of the horizontal grid. |
36 |
|
|
c All internal loops are performed on that grid. |
37 |
|
|
c nlayer Number of vertical layers. |
38 |
|
|
c pdayref Day of reference for the simulation |
39 |
|
|
c firstcall True at the first call |
40 |
|
|
c lastcall True at the last call |
41 |
|
|
c pday Number of days counted from the North. Spring |
42 |
|
|
c equinoxe. |
43 |
|
|
c |
44 |
|
|
c======================================================================= |
45 |
|
|
c |
46 |
|
|
c----------------------------------------------------------------------- |
47 |
|
|
c declarations: |
48 |
|
|
c ------------- |
49 |
|
|
|
50 |
|
|
|
51 |
|
|
REAL prad,pg,pr,pcpp,punjours |
52 |
|
|
|
53 |
|
|
INTEGER ngrid,nlayer |
54 |
|
|
REAL plat(ngrid),plon(ngrid),parea(klon),pcu(klon),pcv(klon) |
55 |
|
|
INTEGER pdayref |
56 |
|
|
|
57 |
|
|
REAL ptimestep |
58 |
|
|
|
59 |
|
|
IF (nlayer.NE.klev) THEN |
60 |
|
|
PRINT*,'STOP in inifis' |
61 |
|
|
PRINT*,'Probleme de dimensions :' |
62 |
|
|
PRINT*,'nlayer = ',nlayer |
63 |
|
|
PRINT*,'klev = ',klev |
64 |
|
|
STOP |
65 |
|
|
ENDIF |
66 |
|
|
|
67 |
|
|
IF (ngrid.NE.klon) THEN |
68 |
|
|
PRINT*,'STOP in inifis' |
69 |
|
|
PRINT*,'Probleme de dimensions :' |
70 |
|
|
PRINT*,'ngrid = ',ngrid |
71 |
|
|
PRINT*,'klon = ',klon |
72 |
|
|
STOP |
73 |
|
|
ENDIF |
74 |
|
|
|
75 |
|
|
airephy=parea |
76 |
|
|
cuphy=pcu |
77 |
|
|
cvphy=pcv |
78 |
|
|
rlond = plon |
79 |
|
|
rlatd = plat |
80 |
|
|
|
81 |
|
|
call suphec |
82 |
|
|
print*,'ATTENTION !!! TRAVAILLER SUR INIPHYSIQ' |
83 |
|
|
print*,'CONTROLE DES LATITUDES, LONGITUDES, PARAMETRES ...' |
84 |
|
|
|
85 |
|
|
END |