1 | |
---|
2 | ! $Id: inifis.F90 1992 2014-03-05 13:19:12Z musat $ |
---|
3 | |
---|
4 | SUBROUTINE inifis(ngrid, nlayer, punjours, pdayref, ptimestep, plat, plon, & |
---|
5 | parea, prad, pg, pr, pcpp) |
---|
6 | USE dimphy |
---|
7 | IMPLICIT NONE |
---|
8 | |
---|
9 | ! ======================================================================= |
---|
10 | |
---|
11 | ! subject: |
---|
12 | ! -------- |
---|
13 | |
---|
14 | ! Initialisation for the physical parametrisations of the LMD |
---|
15 | ! martian atmospheric general circulation modele. |
---|
16 | |
---|
17 | ! author: Frederic Hourdin 15 / 10 /93 |
---|
18 | ! ------- |
---|
19 | |
---|
20 | ! arguments: |
---|
21 | ! ---------- |
---|
22 | |
---|
23 | ! input: |
---|
24 | ! ------ |
---|
25 | |
---|
26 | ! ngrid Size of the horizontal grid. |
---|
27 | ! All internal loops are performed on that grid. |
---|
28 | ! nlayer Number of vertical layers. |
---|
29 | ! pdayref Day of reference for the simulation |
---|
30 | ! firstcall True at the first call |
---|
31 | ! lastcall True at the last call |
---|
32 | ! pday Number of days counted from the North. Spring |
---|
33 | ! equinoxe. |
---|
34 | |
---|
35 | ! ======================================================================= |
---|
36 | |
---|
37 | ! ----------------------------------------------------------------------- |
---|
38 | ! declarations: |
---|
39 | ! ------------- |
---|
40 | |
---|
41 | ! ym#include "dimensions.h" |
---|
42 | ! ym#include "dimphy.h" |
---|
43 | |
---|
44 | include 'iniprint.h' |
---|
45 | REAL prad, pg, pr, pcpp, punjours |
---|
46 | |
---|
47 | INTEGER ngrid, nlayer |
---|
48 | REAL plat(ngrid), plon(ngrid), parea(klon) |
---|
49 | INTEGER pdayref |
---|
50 | |
---|
51 | REAL ptimestep |
---|
52 | CHARACTER (LEN=20) :: modname = 'inifis' |
---|
53 | CHARACTER (LEN=80) :: abort_message |
---|
54 | |
---|
55 | |
---|
56 | IF (nlayer/=klev) THEN |
---|
57 | PRINT *, 'STOP in inifis' |
---|
58 | PRINT *, 'Probleme de dimensions :' |
---|
59 | PRINT *, 'nlayer = ', nlayer |
---|
60 | PRINT *, 'klev = ', klev |
---|
61 | abort_message = '' |
---|
62 | CALL abort_gcm(modname, abort_message, 1) |
---|
63 | END IF |
---|
64 | |
---|
65 | IF (ngrid/=klon) THEN |
---|
66 | PRINT *, 'STOP in inifis' |
---|
67 | PRINT *, 'Probleme de dimensions :' |
---|
68 | PRINT *, 'ngrid = ', ngrid |
---|
69 | PRINT *, 'klon = ', klon |
---|
70 | abort_message = '' |
---|
71 | CALL abort_gcm(modname, abort_message, 1) |
---|
72 | END IF |
---|
73 | |
---|
74 | RETURN |
---|
75 | abort_message = 'Cette version demande les fichier rnatur.dat & |
---|
76 | & & |
---|
77 | & et surf.def' |
---|
78 | CALL abort_gcm(modname, abort_message, 1) |
---|
79 | |
---|
80 | END SUBROUTINE inifis |
---|