source: LMDZ4/branches/pre_V3/libf/phylmd/iniphysiq.F @ 5428

Last change on this file since 5428 was 701, checked in by (none), 19 years ago

This commit was manufactured by cvs2svn to create branch 'pre_V3'.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.1 KB
Line 
1!
2! $Header$
3!
4c
5c
6      SUBROUTINE iniphysiq(ngrid,nlayer,
7     $           punjours,
8     $           pdayref,ptimestep,
9     $           plat,plon,parea,pcu,pcv,
10     $           prad,pg,pr,pcpp)
11      IMPLICIT NONE
12c
13c=======================================================================
14c
15c   subject:
16c   --------
17c
18c   Initialisation for the physical parametrisations of the LMD
19c   martian atmospheric general circulation modele.
20c
21c   author: Frederic Hourdin 15 / 10 /93
22c   -------
23c
24c   arguments:
25c   ----------
26c
27c   input:
28c   ------
29c
30c    ngrid                 Size of the horizontal grid.
31c                          All internal loops are performed on that grid.
32c    nlayer                Number of vertical layers.
33c    pdayref               Day of reference for the simulation
34c    firstcall             True at the first call
35c    lastcall              True at the last call
36c    pday                  Number of days counted from the North. Spring
37c                          equinoxe.
38c
39c=======================================================================
40c
41c-----------------------------------------------------------------------
42c   declarations:
43c   -------------
44 
45#include "dimensions.h"
46#include "dimphy.h"
47#include "comgeomphy.h"
48
49      REAL prad,pg,pr,pcpp,punjours
50 
51      INTEGER ngrid,nlayer
52      REAL plat(ngrid),plon(ngrid),parea(klon),pcu(klon),pcv(klon)
53      INTEGER pdayref
54 
55      REAL ptimestep
56 
57      IF (nlayer.NE.klev) THEN
58         PRINT*,'STOP in inifis'
59         PRINT*,'Probleme de dimensions :'
60         PRINT*,'nlayer     = ',nlayer
61         PRINT*,'klev   = ',klev
62         STOP
63      ENDIF
64
65      IF (ngrid.NE.klon) THEN
66         PRINT*,'STOP in inifis'
67         PRINT*,'Probleme de dimensions :'
68         PRINT*,'ngrid     = ',ngrid
69         PRINT*,'klon   = ',klon
70         STOP
71      ENDIF
72
73      airephy=parea
74      cuphy=pcu
75      cvphy=pcv
76      rlond = plon
77      rlatd = plat
78
79      call suphec
80      print*,'ATTENTION !!! TRAVAILLER SUR INIPHYSIQ'
81      print*,'CONTROLE DES LATITUDES, LONGITUDES, PARAMETRES ...'
82
83
84      RETURN
859999  STOP'Cette version demande les fichier rnatur.dat et surf.def'
86      END
Note: See TracBrowser for help on using the repository browser.