Search:

Login

Help/Guide

About Trac

Preferences

Wiki

Timeline

Roadmap

Browse Source

View Tickets

Search

Context Navigation

← Previous Revision

Latest Revision

Next Revision →

Blame

Revision Log

source: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/iniphysiq.F @ 3400

View diff against:

View revision:

Last change on this file since 3400 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 | !
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 dimphy
12 | USE mod_grid_phy_lmdz
13 | USE mod_phys_lmdz_para
14 | USE comgeomphy
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 | cym#include "dimensions.h"
51 | cym#include "dimphy.h"
52 | cym#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
79 | c$OMP PARALLEL PRIVATE(ibegin,iend)
80 | c$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 |
92 | c$OMP END PARALLEL
93 |
94 | print*,'ATTENTION !!! TRAVAILLER SUR INIPHYSIQ'
95 | print*,'CONTROLE DES LATITUDES, LONGITUDES, PARAMETRES ...'
96 |
97 | RETURN
98 | 9999  STOP'Cette version demande les fichier rnatur.dat et surf.def'
99 | END

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format

Powered by Trac 1.2
By Edgewall Software.

Visit the Trac open source project at
http://trac.edgewall.org/