/[lmdze]/trunk/libf/bibio/initdynav.f90
ViewVC logotype

Annotation of /trunk/libf/bibio/initdynav.f90

Parent Directory Parent Directory | Revision Log Revision Log

Revision 15 - (hide annotations)
Fri Aug 1 15:24:12 2008 UTC (15 years, 9 months ago) by guez
File size: 4256 byte(s) 
-- Minor modification of input/output:

Added variable "Sigma_O3_Royer" to "histday.nc". "ecrit_day" is not
modified in "physiq". Removed variables "pyu1", "pyv1", "ftsol1",
"ftsol2", "ftsol3", "ftsol4", "psrf1", "psrf2", "psrf3", "psrf4"
"mfu", "mfd", "en_u", "en_d", "de_d", "de_u", "coefh" from
"histrac.nc".

Variable "raz_date" of module "conf_gcm_m" has logical type instead of
integer type.

-- Should not change any result at run time:

Modified calls to "IOIPSL_Lionel" procedures because the interfaces of
these procedures have been simplified.

Changed name of variable in module "start_init_orog_m": "masque" to
"mask".

Created a module containing procedure "phyredem".

Removed arguments "punjours", "pdayref" and "ptimestep" of procedure
"iniphysiq".

Renamed procedure "gr_phy_write" to "gr_phy_write_2d". Created
procedure "gr_phy_write_3d".

Removed procedures "ini_undefstd", "moy_undefSTD", "calcul_STDlev",
"calcul_divers".
1 guez 3 module initdynav_m
2    
3     ! This module is clean: no C preprocessor directive, no include line
4    
5     implicit none
6    
7     contains
8    
9     subroutine initdynav(day0, anne0, tstep, nq, fileid, infile, t_ops, t_wrt)
10    
11     ! From initdynav.F,v 1.1.1.1 2004/05/19 12:53:05
12    
13     USE IOIPSL, only: ymds2ju, histbeg_totreg, histvert, histdef, histend
14     use dimens_m
15     use paramet_m
16     use comconst, only: pi
17     use comvert, only: nivsigs
18     use logic
19     use comgeom
20     use serre
21     use temps
22     use ener
23     use advtrac_m, only: ttext
24    
25     ! Routine d'initialisation des ecritures des fichiers histoires LMDZ
26     ! au format IOIPSL. Initialisation du fichier histoire moyenne.
27    
28     ! Appels succesifs des routines: histbeg
29     ! histhori
30     ! histver
31     ! histdef
32     ! histend
33    
34     ! Entree:
35     ! infile: nom du fichier histoire a creer
36     ! day0,anne0: date de reference
37     ! tstep : frequence d'ecriture
38     ! t_ops: frequence de l'operation pour IOIPSL
39     ! t_wrt: frequence d'ecriture sur le fichier
40     ! nq: nombre de traceurs
41    
42     ! Sortie:
43     ! fileid: ID du fichier netcdf cree
44    
45     ! L. Fairhead, LMD, 03/99
46    
47     ! Arguments
48     character(len=*) infile
49     integer day0, anne0
50     real, intent(in):: tstep, t_ops, t_wrt
51     integer fileid
52     integer nq
53     integer thoriid, zvertiid
54    
55     ! Variables locales
56    
57     integer tau0
58     real zjulian
59     integer iq
60     real rlong(iip1,jjp1), rlat(iip1,jjp1)
61     integer ii,jj
62     integer zan, dayref
63    
64     !----------------------------------------------------
65    
66     ! Initialisations
67    
68     pi = 4. * atan (1.)
69     !
70     ! Appel a histbeg: creation du fichier netcdf et initialisations diverses
71     !
72    
73     zan = anne0
74     dayref = day0
75     CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
76     tau0 = itau_dyn
77    
78     do jj = 1, jjp1
79     do ii = 1, iip1
80     rlong(ii,jj) = rlonv(ii) * 180. / pi
81     rlat(ii,jj) = rlatu(jj) * 180. / pi
82     enddo
83     enddo
84    
85 guez 15 call histbeg_totreg(infile, rlong(:,1), rlat(1,:), &
86 guez 3 1, iip1, 1, jjp1, &
87     tau0, zjulian, tstep, thoriid, fileid)
88    
89     !
90     ! Appel a histvert pour la grille verticale
91     !
92     call histvert(fileid, 'sigss', 'Niveaux sigma','Pa', &
93     llm, nivsigs, zvertiid)
94     !
95     ! Appels a histdef pour la definition des variables a sauvegarder
96     !
97     ! Vents U
98     !
99     write(6,*)'inithistave',tstep
100     call histdef(fileid, 'u', 'vents u scalaires moyennes', &
101     'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
102 guez 15 'ave(X)', t_ops, t_wrt)
103 guez 3
104     !
105     ! Vents V
106     !
107     call histdef(fileid, 'v', 'vents v scalaires moyennes', &
108     'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
109 guez 15 'ave(X)', t_ops, t_wrt)
110 guez 3
111     !
112     ! Temperature
113     !
114     call histdef(fileid, 'temp', 'temperature moyennee', 'K', &
115     iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
116 guez 15 'ave(X)', t_ops, t_wrt)
117 guez 3 !
118     ! Temperature potentielle
119     !
120     call histdef(fileid, 'theta', 'temperature potentielle', 'K', &
121     iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
122 guez 15 'ave(X)', t_ops, t_wrt)
123 guez 3
124    
125     !
126     ! Geopotentiel
127     !
128     call histdef(fileid, 'phi', 'geopotentiel moyenne', '-', &
129     iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
130 guez 15 'ave(X)', t_ops, t_wrt)
131 guez 3 !
132     ! Traceurs
133     !
134     DO iq=1,nq
135     call histdef(fileid, ttext(iq), ttext(iq), '-', &
136     iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
137 guez 15 'ave(X)', t_ops, t_wrt)
138 guez 3 enddo
139     !
140     ! Masse
141     !
142     call histdef(fileid, 'masse', 'masse', 'kg', &
143     iip1, jjp1, thoriid, 1, 1, 1, -99, &
144 guez 15 'ave(X)', t_ops, t_wrt)
145 guez 3 !
146     ! Pression au sol
147     !
148     call histdef(fileid, 'ps', 'pression naturelle au sol', 'Pa', &
149     iip1, jjp1, thoriid, 1, 1, 1, -99, &
150 guez 15 'ave(X)', t_ops, t_wrt)
151 guez 3 !
152     ! Pression au sol
153     !
154     call histdef(fileid, 'phis', 'geopotentiel au sol', '-', &
155     iip1, jjp1, thoriid, 1, 1, 1, -99, &
156 guez 15 'ave(X)', t_ops, t_wrt)
157 guez 3 !
158     ! Fin
159     !
160     call histend(fileid)
161    
162     end subroutine initdynav
163    
164     end module initdynav_m
  ViewVC Help
Powered by ViewVC 1.1.21