1 |
! |
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 |