source: LMDZ4/branches/LMDZ4-dev-20091210/libf/phylmd/iniphysiq.F @ 5170

Last change on this file since 5170 was 879, checked in by Laurent Fairhead, 17 years ago

Suite de la bascule vers une physique avec thermiques, nouvelle convection, poche froide ...
LF

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