Changeset 5084 for LMDZ6/trunk/libf/phylmd/Dust/condsurfs_new.F
- Timestamp:
- Jul 19, 2024, 6:40:44 PM (4 months ago)
- File:
-
- 1 edited
-
LMDZ6/trunk/libf/phylmd/Dust/condsurfs_new.F (modified) (19 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/Dust/condsurfs_new.F
r5075 r5084 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_nowrite12 11 IMPLICIT none 13 12 c … … 16 15 c 17 16 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 <0 .OR. jour>(366-1)) THEN42 IF (jour.LT.0 .OR. jour.GT.(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 /=NF_NOERR) then64 if (ierr.ne.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 = nf90_get_var(nid, nvarid, lmt_so2b_glo, debut, epais)75 IF (ierr /=NF_NOERR) THEN74 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,lmt_so2b_glo) 75 IF (ierr .NE. 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 = nf90_get_var(nid, nvarid, lmt_so2h_glo, debut, epais)87 IF (ierr /=NF_NOERR) THEN86 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,lmt_so2h_glo) 87 IF (ierr .NE. 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 = nf90_get_var(nid, nvarid, lmt_so2bb_h_glo, debut, epais) 96 IF (ierr /= NF_NOERR) THEN 95 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, 96 . epais, lmt_so2bb_h_glo) 97 IF (ierr .NE. NF_NOERR) THEN 97 98 PRINT*, 'Pb de lecture pour les sources so2 BB high' 98 99 CALL exit(1) … … 102 103 ! 103 104 ierr = NF_INQ_VARID (nid, "SO2BBL", nvarid) 104 ierr = nf90_get_var(nid, nvarid, lmt_so2bb_l_glo, debut, epais) 105 IF (ierr /= NF_NOERR) THEN 105 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, 106 . epais, lmt_so2bb_l_glo) 107 IF (ierr .NE. NF_NOERR) THEN 106 108 PRINT*, 'Pb de lecture pour les sources so2 BB low' 107 109 CALL exit(1) … … 111 113 ! 112 114 ierr = NF_INQ_VARID (nid, "SO2BA", nvarid) 113 ierr = nf90_get_var(nid, nvarid, lmt_so2ba_glo, debut, epais)114 IF (ierr /=NF_NOERR) THEN115 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_so2ba_glo) 116 IF (ierr .NE. NF_NOERR) THEN 115 117 PRINT*, 'Pb de lecture pour les sources so2 ship' 116 118 CALL exit(1) … … 120 122 ! 121 123 ierr = NF_INQ_VARID (nid, "SO2NFF", nvarid) 122 ierr = nf90_get_var(nid, nvarid, lmt_so2nff_glo, debut, epais) 123 IF (ierr /= NF_NOERR) THEN 124 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 125 . lmt_so2nff_glo) 126 IF (ierr .NE. NF_NOERR) THEN 124 127 PRINT*, 'Pb de lecture pour les sources so2 non FF' 125 128 CALL exit(1) … … 132 135 !======================================================================= 133 136 ierr = NF_OPEN ("sulphur_emissions_nat.nc", NF_NOWRITE, nid) 134 if (ierr /=NF_NOERR) then137 if (ierr.ne.NF_NOERR) then 135 138 write(6,*)' Pb d''ouverture du fichier sulphur_emissions_nat' 136 139 write(6,*)' ierr = ', ierr … … 141 144 c 142 145 ierr = NF_INQ_VARID (nid, "DMSB", nvarid) 143 ierr = nf90_get_var(nid, nvarid, lmt_dmsbio_glo, debut, epais)144 IF (ierr /=NF_NOERR) THEN146 ierr = NF_GET_VARA_DOUBLE (nid, nvarid,debut,epais,lmt_dmsbio_glo) 147 IF (ierr .NE. NF_NOERR) THEN 145 148 PRINT*, 'Pb de lecture pour les sources dms bio' 146 149 CALL exit(1) … … 150 153 c 151 154 ierr = NF_INQ_VARID (nid, "H2SB", nvarid) 152 ierr = nf90_get_var(nid, nvarid, lmt_h2sbio_glo, debut, epais)153 IF (ierr /=NF_NOERR) THEN155 ierr = NF_GET_VARA_DOUBLE (nid, nvarid,debut,epais,lmt_h2sbio_glo) 156 IF (ierr .NE. NF_NOERR) THEN 154 157 PRINT*, 'Pb de lecture pour les sources h2s bio' 155 158 CALL exit(1) … … 158 161 c Ocean surface concentration of dms (emissions are computed later) 159 162 c 160 IF (flag_dms ==4) THEN163 IF (flag_dms.EQ.4) THEN 161 164 c 162 165 ierr = NF_INQ_VARID (nid, "DMSC2", nvarid) 163 ierr = nf90_get_var(nid, nvarid, lmt_dmsconc_glo, debut, epais)164 IF (ierr /=NF_NOERR) THEN166 ierr = NF_GET_VARA_DOUBLE (nid,nvarid,debut,epais,lmt_dmsconc_glo) 167 IF (ierr .NE. NF_NOERR) THEN 165 168 PRINT*, 'Pb de lecture pour les sources dms conc 2' 166 169 CALL exit(1) … … 187 190 print *,' Jour = ',jour 188 191 ierr = NF_OPEN ("sulphur_emissions_volc.nc", NF_NOWRITE, nid) 189 if (ierr /=NF_NOERR) then192 if (ierr.ne.NF_NOERR) then 190 193 write(6,*)' Pb d''ouverture du fichier sulphur_emissions_volc' 191 194 write(6,*)' ierr = ', ierr … … 197 200 ! ierr = NF_INQ_VARID (nid, "VOLC", nvarid) 198 201 ierr = NF_INQ_VARID (nid, "flx_volc_cont", nvarid) 199 ierr = nf90_get_var(nid, nvarid, lmt_so2volc_cont_glo, debut, epais) 200 IF (ierr /= NF_NOERR) THEN 202 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 203 . lmt_so2volc_cont_glo) 204 IF (ierr .NE. NF_NOERR) THEN 201 205 PRINT*, 'Pb de lecture pour les sources so2 volcan (cont)' 202 206 CALL exit(1) … … 210 214 ! ierr = NF_INQ_VARID (nid, "ALTI", nvarid) 211 215 ierr = NF_INQ_VARID (nid, "flx_volc_altcont", nvarid) 212 ierr = nf90_get_var(nid, nvarid, lmt_altvolc_cont_glo, debut, epais) 213 IF (ierr /= NF_NOERR) THEN 216 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 217 . lmt_altvolc_cont_glo) 218 IF (ierr .NE. NF_NOERR) THEN 214 219 PRINT*, 'Pb de lecture pour les altitudes volcan (cont)' 215 220 CALL exit(1) … … 219 224 c 220 225 ierr = NF_INQ_VARID (nid, "flx_volc_expl", nvarid) 221 ierr = nf90_get_var(nid, nvarid, lmt_so2volc_expl_glo, debut, epais) 222 IF (ierr /= NF_NOERR) THEN 226 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 227 . lmt_so2volc_expl_glo) 228 IF (ierr .NE. NF_NOERR) THEN 223 229 PRINT*, 'Pb de lecture pour les sources so2 volcan (expl)' 224 230 CALL exit(1) … … 231 237 c 232 238 ierr = NF_INQ_VARID (nid, "flx_volc_altexpl", nvarid) 233 ierr = nf90_get_var(nid, nvarid, lmt_altvolc_expl_glo, debut, epais) 234 IF (ierr /= NF_NOERR) THEN 239 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, 240 . lmt_altvolc_expl_glo) 241 IF (ierr .NE. NF_NOERR) THEN 235 242 PRINT*, 'Pb de lecture pour les altitudes volcan' 236 243 CALL exit(1)
Note: See TracChangeset
for help on using the changeset viewer.
