/[lmdze]/trunk/libf/phylmd/iniphysiq.f
ViewVC logotype

Contents of /trunk/libf/phylmd/iniphysiq.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3 - (show annotations)
Wed Feb 27 13:16:39 2008 UTC (16 years, 2 months ago) by guez
File size: 2208 byte(s)
Initial import
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

  ViewVC Help
Powered by ViewVC 1.1.21