source: LMDZ4/branches/V3_test/libf/phylmd/iniphysiq.F @ 1198

Last change on this file since 1198 was 704, checked in by Laurent Fairhead, 18 years ago

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.5 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      USE dimphy,klon2=>klon,klon=>klon2
12      USE comgeomphy
13      IMPLICIT NONE
14c
15c=======================================================================
16c
17c   subject:
18c   --------
19c
20c   Initialisation for the physical parametrisations of the LMD
21c   martian atmospheric general circulation modele.
22c
23c   author: Frederic Hourdin 15 / 10 /93
24c   -------
25c
26c   arguments:
27c   ----------
28c
29c   input:
30c   ------
31c
32c    ngrid                 Size of the horizontal grid.
33c                          All internal loops are performed on that grid.
34c    nlayer                Number of vertical layers.
35c    pdayref               Day of reference for the simulation
36c    firstcall             True at the first call
37c    lastcall              True at the last call
38c    pday                  Number of days counted from the North. Spring
39c                          equinoxe.
40c
41c=======================================================================
42c
43c-----------------------------------------------------------------------
44c   declarations:
45c   -------------
46 
47cym#include "dimensions.h"
48cym#include "dimphy.h"
49cym#include "comgeomphy.h"
50#include "YOMCST.h"
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      INTEGER :: ibegin,iend
57 
58      REAL ptimestep
59 
60      IF (nlayer.NE.klev) THEN
61         PRINT*,'STOP in inifis'
62         PRINT*,'Probleme de dimensions :'
63         PRINT*,'nlayer     = ',nlayer
64         PRINT*,'klev   = ',klev
65         STOP
66      ENDIF
67
68      IF (ngrid.NE.klon) THEN
69         PRINT*,'STOP in inifis'
70         PRINT*,'Probleme de dimensions :'
71         PRINT*,'ngrid     = ',ngrid
72         PRINT*,'klon   = ',klon
73         STOP
74      ENDIF
75c$OMP PARALLEL PRIVATE(ibegin,iend)
76c$OMP+         SHARED(parea,pcu,pcv,plon,plat)
77      ibegin=klon_begin+klon_omp_begin(omp_rank)-1
78      iend=ibegin+klon_omp-1
79
80      airephy(1:klon_omp)=parea(ibegin:iend)
81      cuphy(1:klon_omp)=pcu(ibegin:iend)
82      cvphy(1:klon_omp)=pcv(ibegin:iend)
83      rlond(1:klon_omp) = plon(ibegin:iend)
84      rlatd(1:klon_omp) = plat(ibegin:iend)
85
86      call suphec
87
88c$OMP END PARALLEL
89
90      print*,'ATTENTION !!! TRAVAILLER SUR INIPHYSIQ'
91      print*,'CONTROLE DES LATITUDES, LONGITUDES, PARAMETRES ...'
92
93
94      RETURN
959999  STOP'Cette version demande les fichier rnatur.dat et surf.def'
96      END
Note: See TracBrowser for help on using the repository browser.