source: LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/condsurfs.f90 @ 5105

Last change on this file since 5105 was 5105, checked in by abarral, 8 weeks ago

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F to *.f90
Fix gradsdef.h formatting
Remove unnecessary "RETURN" at the end of functions/subroutines
File size: 6.5 KB
Line 
1SUBROUTINE condsurfs(jour, edgar, flag_dms, &
2        lmt_so2h, lmt_so2b, lmt_so2bb, lmt_so2ba, &
3        lmt_so2volc, lmt_altvolc, &
4        lmt_dmsbio, lmt_h2sbio, lmt_dms, lmt_dmsconc)
5  USE dimphy
6  USE netcdf, ONLY: nf90_close, nf90_noerr, nf90_inq_varid, nf90_open, &
7          nf90_nowrite, nf90_get_var
8  IMPLICIT none
9
10  ! Lire les conditions aux limites du modele pour la chimie.
11  ! --------------------------------------------------------
12
13  INCLUDE "dimensions.h"
14
15  REAL :: lmt_so2h(klon), lmt_so2b(klon), lmt_so2bb(klon)
16  REAL :: lmt_dmsbio(klon), lmt_h2sbio(klon), lmt_so2ba(klon)
17  REAL :: lmt_so2volc(klon), lmt_altvolc(klon)
18  REAL :: lmt_dms(klon), lmt_dmsconc(klon)
19  LOGICAL :: edgar
20  INTEGER :: flag_dms
21
22  INTEGER :: jour, i
23  INTEGER :: ierr
24  INTEGER :: nid, nvarid
25  INTEGER :: debut(2), epais(2)
26
27  IF (jour<0 .OR. jour>(360 - 1)) THEN
28    IF ((jour>(360 - 1)) .AND. (jour<=367)) THEN
29      jour = 360 - 1
30      print *, 'JE: jour changed to jour= ', jour
31    ELSE
32      PRINT*, 'Le jour demande n est pas correcte:', jour
33      CALL ABORT
34    ENDIF
35  ENDIF
36
37  ierr = nf90_open ("limitsoufre.nc", nf90_nowrite, nid)
38  if (ierr/=nf90_noerr) then
39    write(6, *)' Pb d''ouverture du fichier limitsoufre.nc'
40    write(6, *)' ierr = ', ierr
41    CALL exit(1)
42  endif
43
44  ! Tranche a lire:
45  debut(1) = 1
46  debut(2) = jour + 1
47  epais(1) = klon
48  epais(2) = 1
49
50  ierr = nf90_inq_varid (nid, "VOLC", nvarid)
51  ierr = nf90_get_var(nid, nvarid, lmt_so2volc, debut, epais)
52  IF (ierr /= nf90_noerr) THEN
53    PRINT*, 'Pb de lecture pour les sources so2 volcan'
54    CALL exit(1)
55  ENDIF
56
57  ierr = nf90_inq_varid (nid, "ALTI", nvarid)
58  ierr = nf90_get_var(nid, nvarid, lmt_altvolc, debut, epais)
59  IF (ierr /= nf90_noerr) THEN
60    PRINT*, 'Pb de lecture pour les altitudes volcan'
61    CALL exit(1)
62  ENDIF
63
64  IF (edgar) THEN   !--EDGAR w/o ship and biomass burning
65
66    ierr = nf90_inq_varid (nid, "SO2ED95L", nvarid)
67    ierr = nf90_get_var(nid, nvarid, lmt_so2b, debut, epais)
68    IF (ierr /= nf90_noerr) THEN
69      PRINT*, 'Pb de lecture pour les sources so2 edgar low'
70      CALL exit(1)
71    ENDIF
72
73    ierr = nf90_inq_varid (nid, "SO2ED95H", nvarid)
74    ierr = nf90_get_var(nid, nvarid, lmt_so2h, debut, epais)
75    IF (ierr /= nf90_noerr) THEN
76      PRINT*, 'Pb de lecture pour les sources so2 edgar high'
77      CALL exit(1)
78    ENDIF
79
80  ELSE  !--GEIA
81
82    ierr = nf90_inq_varid (nid, "SO2H", nvarid)
83    ierr = nf90_get_var(nid, nvarid, lmt_so2h, debut, epais)
84    IF (ierr /= nf90_noerr) THEN
85      PRINT*, 'Pb de lecture pour les sources so2 haut'
86      CALL exit(1)
87    ENDIF
88
89    ierr = nf90_inq_varid (nid, "SO2B", nvarid)
90    ierr = nf90_get_var(nid, nvarid, lmt_so2b, debut, epais)
91    IF (ierr /= nf90_noerr) THEN
92      PRINT*, 'Pb de lecture pour les sources so2 bas'
93      CALL exit(1)
94    ENDIF
95
96  ENDIF  !--edgar
97
98  ierr = nf90_inq_varid (nid, "SO2BB", nvarid)
99  ierr = nf90_get_var(nid, nvarid, lmt_so2bb, debut, epais)
100  IF (ierr /= nf90_noerr) THEN
101    PRINT*, 'Pb de lecture pour les sources so2 bb'
102    CALL exit(1)
103  ENDIF
104
105  ierr = nf90_inq_varid (nid, "SO2BA", nvarid)
106  ierr = nf90_get_var(nid, nvarid, lmt_so2ba, debut, epais)
107  IF (ierr /= nf90_noerr) THEN
108    PRINT*, 'Pb de lecture pour les sources so2 bateau'
109    CALL exit(1)
110  ENDIF
111
112  ierr = nf90_inq_varid (nid, "DMSB", nvarid)
113  ierr = nf90_get_var(nid, nvarid, lmt_dmsbio, debut, epais)
114  IF (ierr /= nf90_noerr) THEN
115    PRINT*, 'Pb de lecture pour les sources dms bio'
116    CALL exit(1)
117  ENDIF
118
119  ierr = nf90_inq_varid (nid, "H2SB", nvarid)
120  ierr = nf90_get_var(nid, nvarid, lmt_h2sbio, debut, epais)
121  IF (ierr /= nf90_noerr) THEN
122    PRINT*, 'Pb de lecture pour les sources h2s bio'
123    CALL exit(1)
124  ENDIF
125
126  IF (flag_dms==1) THEN
127
128    ierr = nf90_inq_varid (nid, "DMSL", nvarid)
129    ierr = nf90_get_var(nid, nvarid, lmt_dms, debut, epais)
130    IF (ierr /= nf90_noerr) THEN
131      PRINT*, 'Pb de lecture pour les sources dms liss'
132      CALL exit(1)
133    ENDIF
134
135  ELSEIF (flag_dms==2) THEN
136
137    ierr = nf90_inq_varid (nid, "DMSW", nvarid)
138    ierr = nf90_get_var(nid, nvarid, lmt_dms, debut, epais)
139    IF (ierr /= nf90_noerr) THEN
140      PRINT*, 'Pb de lecture pour les sources dms wann'
141      CALL exit(1)
142    ENDIF
143
144  ELSEIF (flag_dms==3) THEN
145
146    ierr = nf90_inq_varid (nid, "DMSC1", nvarid)
147    ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
148    IF (ierr /= nf90_noerr) THEN
149      PRINT*, 'Pb de lecture pour les sources dmsconc old'
150      CALL exit(1)
151    ENDIF
152
153  ELSEIF (flag_dms==4) THEN
154
155    ierr = nf90_inq_varid (nid, "DMSC2", nvarid)
156    ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
157    IF (ierr /= nf90_noerr) THEN
158      PRINT*, 'Pb de lecture pour les sources dms conc 2'
159      CALL exit(1)
160    ENDIF
161
162  ELSEIF (flag_dms==5) THEN
163
164    ierr = nf90_inq_varid (nid, "DMSC3", nvarid)
165    ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
166    IF (ierr /= nf90_noerr) THEN
167      PRINT*, 'Pb de lecture pour les sources dms conc 3'
168      CALL exit(1)
169    ENDIF
170
171  ELSEIF (flag_dms==6) THEN
172
173    ierr = nf90_inq_varid (nid, "DMSC4", nvarid)
174    ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
175    IF (ierr /= nf90_noerr) THEN
176      PRINT*, 'Pb de lecture pour les sources dms conc 4'
177      CALL exit(1)
178    ENDIF
179
180  ELSEIF (flag_dms==7) THEN
181
182    ierr = nf90_inq_varid (nid, "DMSC5", nvarid)
183    ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
184    IF (ierr /= nf90_noerr) THEN
185      PRINT*, 'Pb de lecture pour les sources dms conc 5'
186      CALL exit(1)
187    ENDIF
188
189  ELSEIF (flag_dms==8) THEN
190
191    ierr = nf90_inq_varid (nid, "DMSC6", nvarid)
192    ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
193    IF (ierr /= nf90_noerr) THEN
194      PRINT*, 'Pb de lecture pour les sources dms conc 6'
195      CALL exit(1)
196    ENDIF
197
198  ELSEIF (flag_dms==9) THEN
199
200    ierr = nf90_inq_varid (nid, "DMSC7", nvarid)
201    ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
202    IF (ierr /= nf90_noerr) THEN
203      PRINT*, 'Pb de lecture pour les sources dms conc 7'
204      CALL exit(1)
205    ENDIF
206
207  ELSEIF (flag_dms==10) THEN
208
209    ierr = nf90_inq_varid (nid, "DMSC8", nvarid)
210    ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
211    IF (ierr /= nf90_noerr) THEN
212      PRINT*, 'Pb de lecture pour les sources dms conc 8'
213      CALL exit(1)
214    ENDIF
215
216  ELSE
217
218    PRINT *, 'choix non possible pour flag_dms'
219    STOP
220
221  ENDIF
222
223  ierr = nf90_close(nid)
224
225  IF (flag_dms<=2) THEN
226    DO i = 1, klon
227      lmt_dmsconc(i) = 0.0
228    ENDDO
229  ELSE
230    DO i = 1, klon
231      lmt_dms(i) = 0.0
232    ENDDO
233  ENDIF
234
235  PRINT*, 'Sources SOUFRE lues pour jour: ', jour
236
237
238END SUBROUTINE condsurfs
Note: See TracBrowser for help on using the repository browser.