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