Changeset 5075 for LMDZ6/trunk/libf/phylmd/Dust/condsurfs_new.F
- Timestamp:
- Jul 19, 2024, 10:05:57 AM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/Dust/condsurfs_new.F
r4593 r5075 9 9 USE mod_phys_lmdz_para 10 10 USE dimphy 11 USE lmdz_netcdf, ONLY: nf90_get_var,nf_inq_varid,nf_close,nf_noerr,nf_open,nf_nowrite 11 12 IMPLICIT none 12 13 c … … 15 16 c 16 17 INCLUDE "dimensions.h" 17 INCLUDE "netcdf.inc"18 18 c 19 19 REAL lmt_so2b(klon), lmt_so2h(klon), lmt_so2nff(klon) … … 40 40 INTEGER debut(2),epais(2) 41 41 c 42 IF (jour .LT.0 .OR. jour.GT.(366-1)) THEN42 IF (jour<0 .OR. jour>(366-1)) THEN 43 43 PRINT*,'Le jour demande n est pas correcte:', jour 44 44 print *,'JE: FORCED TO CONTINUE (emissions have … … 62 62 ! 63 63 ierr = NF_OPEN ("sulphur_emissions_antro.nc", NF_NOWRITE, nid) 64 if (ierr .ne.NF_NOERR) then64 if (ierr/=NF_NOERR) then 65 65 write(6,*)' Pb d''ouverture du fichier sulphur_emissions_antro' 66 66 write(6,*)' ierr = ', ierr … … 72 72 ! 73 73 ierr = NF_INQ_VARID (nid, "SO2FF_LOW", nvarid) 74 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,lmt_so2b_glo)75 IF (ierr .NE.NF_NOERR) THEN74 ierr = nf90_get_var(nid, nvarid, lmt_so2b_glo, debut, epais) 75 IF (ierr /= NF_NOERR) THEN 76 76 PRINT*, 'Pb de lecture pour les sources so2 low' 77 77 print *,'JE klon, jour, debut ,epais ',klon_glo,jour,debut,epais … … 84 84 ! 85 85 ierr = NF_INQ_VARID (nid, "SO2FF_HIGH", nvarid) 86 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,lmt_so2h_glo)87 IF (ierr .NE.NF_NOERR) THEN86 ierr = nf90_get_var(nid, nvarid, lmt_so2h_glo, debut, epais) 87 IF (ierr /= NF_NOERR) THEN 88 88 PRINT*, 'Pb de lecture pour les sources so2 high' 89 89 CALL exit(1) … … 93 93 ! 94 94 ierr = NF_INQ_VARID (nid, "SO2BBH", nvarid) 95 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, 96 . epais, lmt_so2bb_h_glo) 97 IF (ierr .NE. NF_NOERR) THEN 95 ierr = nf90_get_var(nid, nvarid, lmt_so2bb_h_glo, debut, epais) 96 IF (ierr /= NF_NOERR) THEN 98 97 PRINT*, 'Pb de lecture pour les sources so2 BB high' 99 98 CALL exit(1) … … 103 102 ! 104 103 ierr = NF_INQ_VARID (nid, "SO2BBL", nvarid) 105 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, 106 . epais, lmt_so2bb_l_glo) 107 IF (ierr .NE. NF_NOERR) THEN 104 ierr = nf90_get_var(nid, nvarid, lmt_so2bb_l_glo, debut, epais) 105 IF (ierr /= NF_NOERR) THEN 108 106 PRINT*, 'Pb de lecture pour les sources so2 BB low' 109 107 CALL exit(1) … … 113 111 ! 114 112 ierr = NF_INQ_VARID (nid, "SO2BA", nvarid) 115 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_so2ba_glo)116 IF (ierr .NE.NF_NOERR) THEN113 ierr = nf90_get_var(nid, nvarid, lmt_so2ba_glo, debut, epais) 114 IF (ierr /= NF_NOERR) THEN 117 115 PRINT*, 'Pb de lecture pour les sources so2 ship' 118 116 CALL exit(1) … … 122 120 ! 123 121 ierr = NF_INQ_VARID (nid, "SO2NFF", nvarid) 124 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 125 . lmt_so2nff_glo) 126 IF (ierr .NE. NF_NOERR) THEN 122 ierr = nf90_get_var(nid, nvarid, lmt_so2nff_glo, debut, epais) 123 IF (ierr /= NF_NOERR) THEN 127 124 PRINT*, 'Pb de lecture pour les sources so2 non FF' 128 125 CALL exit(1) … … 135 132 !======================================================================= 136 133 ierr = NF_OPEN ("sulphur_emissions_nat.nc", NF_NOWRITE, nid) 137 if (ierr .ne.NF_NOERR) then134 if (ierr/=NF_NOERR) then 138 135 write(6,*)' Pb d''ouverture du fichier sulphur_emissions_nat' 139 136 write(6,*)' ierr = ', ierr … … 144 141 c 145 142 ierr = NF_INQ_VARID (nid, "DMSB", nvarid) 146 ierr = NF_GET_VARA_DOUBLE (nid, nvarid,debut,epais,lmt_dmsbio_glo)147 IF (ierr .NE.NF_NOERR) THEN143 ierr = nf90_get_var(nid, nvarid, lmt_dmsbio_glo, debut, epais) 144 IF (ierr /= NF_NOERR) THEN 148 145 PRINT*, 'Pb de lecture pour les sources dms bio' 149 146 CALL exit(1) … … 153 150 c 154 151 ierr = NF_INQ_VARID (nid, "H2SB", nvarid) 155 ierr = NF_GET_VARA_DOUBLE (nid, nvarid,debut,epais,lmt_h2sbio_glo)156 IF (ierr .NE.NF_NOERR) THEN152 ierr = nf90_get_var(nid, nvarid, lmt_h2sbio_glo, debut, epais) 153 IF (ierr /= NF_NOERR) THEN 157 154 PRINT*, 'Pb de lecture pour les sources h2s bio' 158 155 CALL exit(1) … … 161 158 c Ocean surface concentration of dms (emissions are computed later) 162 159 c 163 IF (flag_dms .EQ.4) THEN160 IF (flag_dms==4) THEN 164 161 c 165 162 ierr = NF_INQ_VARID (nid, "DMSC2", nvarid) 166 ierr = NF_GET_VARA_DOUBLE (nid,nvarid,debut,epais,lmt_dmsconc_glo)167 IF (ierr .NE.NF_NOERR) THEN163 ierr = nf90_get_var(nid, nvarid, lmt_dmsconc_glo, debut, epais) 164 IF (ierr /= NF_NOERR) THEN 168 165 PRINT*, 'Pb de lecture pour les sources dms conc 2' 169 166 CALL exit(1) … … 190 187 print *,' Jour = ',jour 191 188 ierr = NF_OPEN ("sulphur_emissions_volc.nc", NF_NOWRITE, nid) 192 if (ierr .ne.NF_NOERR) then189 if (ierr/=NF_NOERR) then 193 190 write(6,*)' Pb d''ouverture du fichier sulphur_emissions_volc' 194 191 write(6,*)' ierr = ', ierr … … 200 197 ! ierr = NF_INQ_VARID (nid, "VOLC", nvarid) 201 198 ierr = NF_INQ_VARID (nid, "flx_volc_cont", nvarid) 202 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 203 . lmt_so2volc_cont_glo) 204 IF (ierr .NE. NF_NOERR) THEN 199 ierr = nf90_get_var(nid, nvarid, lmt_so2volc_cont_glo, debut, epais) 200 IF (ierr /= NF_NOERR) THEN 205 201 PRINT*, 'Pb de lecture pour les sources so2 volcan (cont)' 206 202 CALL exit(1) … … 214 210 ! ierr = NF_INQ_VARID (nid, "ALTI", nvarid) 215 211 ierr = NF_INQ_VARID (nid, "flx_volc_altcont", nvarid) 216 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 217 . lmt_altvolc_cont_glo) 218 IF (ierr .NE. NF_NOERR) THEN 212 ierr = nf90_get_var(nid, nvarid, lmt_altvolc_cont_glo, debut, epais) 213 IF (ierr /= NF_NOERR) THEN 219 214 PRINT*, 'Pb de lecture pour les altitudes volcan (cont)' 220 215 CALL exit(1) … … 224 219 c 225 220 ierr = NF_INQ_VARID (nid, "flx_volc_expl", nvarid) 226 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 227 . lmt_so2volc_expl_glo) 228 IF (ierr .NE. NF_NOERR) THEN 221 ierr = nf90_get_var(nid, nvarid, lmt_so2volc_expl_glo, debut, epais) 222 IF (ierr /= NF_NOERR) THEN 229 223 PRINT*, 'Pb de lecture pour les sources so2 volcan (expl)' 230 224 CALL exit(1) … … 237 231 c 238 232 ierr = NF_INQ_VARID (nid, "flx_volc_altexpl", nvarid) 239 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 240 . lmt_altvolc_expl_glo) 241 IF (ierr .NE. NF_NOERR) THEN 233 ierr = nf90_get_var(nid, nvarid, lmt_altvolc_expl_glo, debut, epais) 234 IF (ierr /= NF_NOERR) THEN 242 235 PRINT*, 'Pb de lecture pour les altitudes volcan' 243 236 CALL exit(1)
Note: See TracChangeset
for help on using the changeset viewer.