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

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

Parent Directory Parent Directory | Revision Log Revision Log


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

  ViewVC Help
Powered by ViewVC 1.1.21