Changeset 5075 for LMDZ6/trunk/libf/phylmd/Dust
- Timestamp:
- Jul 19, 2024, 10:05:57 AM (12 months ago)
- Location:
- LMDZ6/trunk/libf/phylmd/Dust
- Files:
-
- 9 edited
-
condsurfc.F (modified) (2 diffs)
-
condsurfc_new.F (modified) (15 diffs)
-
condsurfs.F (modified) (2 diffs)
-
condsurfs_new.F (modified) (19 diffs)
-
phys_output_write_spl_mod.F90 (modified) (1 diff)
-
phytracr_spl_mod.F90 (modified) (20 diffs)
-
read_dust.F (modified) (2 diffs)
-
read_surface.F90 (modified) (1 diff)
-
read_vent.F (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/Dust/condsurfc.F
r5073 r5075 4 4 . lmt_omnat) 5 5 USE dimphy 6 USE lmdz_netcdf, ONLY: nf_close,nf_noerr,nf_inq_varid,nf_open,nf_nowrite,nf90_get_var 6 7 IMPLICIT none 7 8 ! … … 10 11 ! 11 12 INCLUDE "dimensions.h" 12 INCLUDE "netcdf.inc" 13 13 14 14 REAL lmt_bcff(klon), lmt_bcbb(klon),lmt_bc_penner(klon) 15 15 REAL lmt_omff(klon), lmt_ombb(klon) -
LMDZ6/trunk/libf/phylmd/Dust/condsurfc_new.F
r4593 r5075 6 6 USE mod_phys_lmdz_para 7 7 USE dimphy 8 USE lmdz_netcdf, ONLY:nf90_get_var,nf_close,nf_noerr,nf_inq_varid,nf_open,nf_nowrite 8 9 IMPLICIT none 9 10 c … … 12 13 c 13 14 INCLUDE "dimensions.h" 14 INCLUDE "netcdf.inc" 15 15 16 16 REAL lmt_bcff(klon), lmt_bcnff(klon), lmt_bcba(klon) 17 17 REAL lmt_omff(klon), lmt_omnff(klon), lmt_ombb(klon) … … 36 36 c 37 37 ! IF (jour.LT.0 .OR. jour.GT.(366-1)) THEN 38 IF (jour .LT.0 .OR. jour.GT.366) THEN38 IF (jour<0 .OR. jour>366) THEN 39 39 PRINT*,'Le jour demande n est pas correcte:', jour 40 40 print *,'JE: FORCED TO CONTINUE (emissions have … … 58 58 ! 59 59 ierr = NF_OPEN ("carbon_emissions.nc", NF_NOWRITE, nid1) 60 if (ierr .ne.NF_NOERR) then60 if (ierr/=NF_NOERR) then 61 61 write(6,*)' Pb d''ouverture du fichier limitbc.nc' 62 62 write(6,*)' ierr = ', ierr … … 67 67 ! 68 68 ierr = NF_INQ_VARID (nid1, "BCFF", nvarid) 69 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 70 . lmt_bcff_glo) 71 IF (ierr .NE. NF_NOERR) THEN 69 ierr = nf90_get_var (nid1, nvarid, lmt_bcff_glo, debut, epais) 70 IF (ierr /= NF_NOERR) THEN 72 71 PRINT*, 'Pb de lecture pour les sources BC' 73 72 CALL exit(1) … … 79 78 ! 80 79 ierr = NF_INQ_VARID (nid1, "BCNFF", nvarid) 81 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 82 . lmt_bcnff_glo) 83 IF (ierr .NE. NF_NOERR) THEN 80 ierr = nf90_get_var (nid1, nvarid, lmt_bcnff_glo, debut, epais) 81 IF (ierr /= NF_NOERR) THEN 84 82 PRINT*, 'Pb de lecture pour les sources BC' 85 83 CALL exit(1) … … 89 87 ! 90 88 ierr = NF_INQ_VARID (nid1, "BCBBL", nvarid) 91 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 92 . lmt_bcbbl_glo) 93 IF (ierr .NE. NF_NOERR) THEN 89 ierr = nf90_get_var (nid1, nvarid, lmt_bcbbl_glo, debut, epais) 90 IF (ierr /= NF_NOERR) THEN 94 91 PRINT*, 'Pb de lecture pour les sources BC low' 95 92 CALL exit(1) … … 99 96 ! 100 97 ierr = NF_INQ_VARID (nid1, "BCBBH", nvarid) 101 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 102 . lmt_bcbbh_glo) 103 IF (ierr .NE. NF_NOERR) THEN 98 ierr = nf90_get_var (nid1, nvarid, lmt_bcbbh_glo, debut, epais) 99 IF (ierr /= NF_NOERR) THEN 104 100 PRINT*, 'Pb de lecture pour les sources BC high' 105 101 CALL exit(1) … … 109 105 ! 110 106 ierr = NF_INQ_VARID (nid1, "BCBA", nvarid) 111 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 112 . lmt_bcba_glo) 113 IF (ierr .NE. NF_NOERR) THEN 107 ierr = nf90_get_var (nid1, nvarid, lmt_bcba_glo, debut, epais) 108 IF (ierr /= NF_NOERR) THEN 114 109 PRINT*, 'Pb de lecture pour les sources BC' 115 110 CALL exit(1) … … 125 120 ! 126 121 ierr = NF_INQ_VARID (nid1, "OMFF", nvarid) 127 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 128 . lmt_omff_glo) 129 IF (ierr .NE. NF_NOERR) THEN 122 ierr = nf90_get_var (nid1, nvarid, lmt_omff_glo, debut, epais) 123 IF (ierr /= NF_NOERR) THEN 130 124 PRINT*, 'Pb de lecture pour les sources OM' 131 125 CALL exit(1) … … 135 129 ! 136 130 ierr = NF_INQ_VARID (nid1, "OMNFF", nvarid) 137 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 138 . lmt_omnff_glo) 139 IF (ierr .NE. NF_NOERR) THEN 131 ierr = nf90_get_var (nid1, nvarid, lmt_omnff_glo, debut, epais) 132 IF (ierr /= NF_NOERR) THEN 140 133 PRINT*, 'Pb de lecture pour les sources OM' 141 134 CALL exit(1) … … 145 138 ! 146 139 ierr = NF_INQ_VARID (nid1, "OMBBL", nvarid) 147 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 148 . lmt_ombbl_glo) 149 IF (ierr .NE. NF_NOERR) THEN 140 ierr = nf90_get_var (nid1, nvarid, lmt_ombbl_glo, debut, epais) 141 IF (ierr /= NF_NOERR) THEN 150 142 PRINT*, 'Pb de lecture pour les sources OM low' 151 143 CALL exit(1) … … 155 147 ! 156 148 ierr = NF_INQ_VARID (nid1, "OMBBH", nvarid) 157 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 158 . lmt_ombbh_glo) 159 IF (ierr .NE. NF_NOERR) THEN 149 ierr = nf90_get_var (nid1, nvarid, lmt_ombbh_glo, debut, epais) 150 IF (ierr /= NF_NOERR) THEN 160 151 PRINT*, 'Pb de lecture pour les sources OM high' 161 152 CALL exit(1) … … 165 156 ! 166 157 ierr = NF_INQ_VARID (nid1, "OMBA", nvarid) 167 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 168 . lmt_omba_glo) 169 IF (ierr .NE. NF_NOERR) THEN 158 ierr = nf90_get_var (nid1, nvarid, lmt_omba_glo, debut, epais) 159 IF (ierr /= NF_NOERR) THEN 170 160 PRINT*, 'Pb de lecture pour les sources OM ship' 171 161 CALL exit(1) … … 175 165 ! 176 166 ierr = NF_INQ_VARID (nid1, "TERP", nvarid) 177 ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, 178 . lmt_terp_glo) 179 IF (ierr .NE. NF_NOERR) THEN 167 ierr = nf90_get_var (nid1, nvarid, lmt_terp_glo, debut, epais) 168 IF (ierr /= NF_NOERR) THEN 180 169 PRINT*, 'Pb de lecture pour les sources Terpene' 181 170 CALL exit(1) -
LMDZ6/trunk/libf/phylmd/Dust/condsurfs.F
r5073 r5075 4 4 . lmt_dmsbio, lmt_h2sbio, lmt_dms, lmt_dmsconc) 5 5 USE dimphy 6 USE lmdz_netcdf, ONLY:nf_close,nf_noerr,nf_inq_varid,nf_open,nf_nowrite,nf90_get_var 6 7 IMPLICIT none 7 8 c … … 10 11 c 11 12 INCLUDE "dimensions.h" 12 INCLUDE "netcdf.inc"13 13 c 14 14 REAL lmt_so2h(klon), lmt_so2b(klon), lmt_so2bb(klon) -
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) -
LMDZ6/trunk/libf/phylmd/Dust/phys_output_write_spl_mod.F90
r5024 r5075 390 390 USE ioipsl, ONLY: histend, histsync 391 391 USE iophy, ONLY: set_itau_iophy, histwrite_phy 392 USE netcdf, ONLY: nf90_fill_real392 USE lmdz_netcdf, ONLY: nf90_fill_real 393 393 ! ug Pour les sorties XIOS 394 394 USE lmdz_xios, ONLY: xios_update_calendar, using_xios -
LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90
r4618 r5075 1441 1441 call abort_gcm('phytracr_mod', 'pb in ok_chimdust 1 SCDU',1) 1442 1442 endif 1443 if ( (id_codu .le. 0) .or. ( id_fine.le.0) ) then1443 if ( (id_codu <= 0) .or. ( id_fine<=0) ) then 1444 1444 call abort_gcm('phytracr_mod', 'pb in ok_chimdust 1',1) 1445 1445 endif … … 2437 2437 ENDDO 2438 2438 ENDDO 2439 IF (iflag_conv .EQ.2) THEN2439 IF (iflag_conv==2) THEN 2440 2440 ! Tiedke 2441 2441 CALL cltrac_spl(pdtphys,coefh,yu1,yv1,t_seri,tmp_var, & 2442 2442 aux_var2,paprs,pplay,aux_var3) 2443 2443 2444 ELSE IF (iflag_conv .GE.3) THEN2444 ELSE IF (iflag_conv>=3) THEN 2445 2445 !KE 2446 2446 CALL cltrac(pdtphys, coefh,t_seri,tmp_var,aux_var2,paprs,pplay, & … … 2494 2494 2495 2495 2496 IF (iflag_conv .GE.3) THEN2496 IF (iflag_conv>=3) THEN 2497 2497 2498 2498 IF (logitime) THEN … … 2786 2786 2787 2787 2788 IF (iflag_conv .EQ.2) THEN2788 IF (iflag_conv==2) THEN 2789 2789 2790 2790 IF (logitime) THEN … … 2839 2839 print *,'iflag_conv bef incloud',iflag_conv 2840 2840 2841 IF (iflag_conv .EQ.2) THEN2841 IF (iflag_conv==2) THEN 2842 2842 ! Tiedke 2843 2843 CALL incloud_scav(.false.,qmin,qmax,masse,henry,kk,prfl, & … … 2877 2877 ! . his_dhbclsc,his_dhbccon,tr_seri) 2878 2878 2879 IF (iflag_conv .EQ.2) THEN2879 IF (iflag_conv==2) THEN 2880 2880 ! Tiedke 2881 2881 … … 2991 2991 ! . dtrconv,tr_seri) 2992 2992 ! ------------------------------------------------------------- 2993 IF (iflag_conv .EQ.2) THEN2993 IF (iflag_conv==2) THEN 2994 2994 ! Tiedke 2995 2995 CALL trconvect(pplay,t_seri,pdtphys,pmfu,pmfd,pen_u,pde_u, & … … 3000 3000 ENDDO 3001 3001 3002 ELSE IF (iflag_conv .GE.3) THEN3002 ELSE IF (iflag_conv>=3) THEN 3003 3003 ! KE 3004 3004 print *,'JE: KE in phytracr_spl' … … 3164 3164 3165 3165 3166 IF (iflag_conv .GE.3) THEN3166 IF (iflag_conv>=3) THEN 3167 3167 IF (logitime) THEN 3168 3168 CALL SYSTEM_CLOCK(COUNT=clock_start) … … 3195 3195 ql_incl = ql_incloud_ref 3196 3196 ! choix du lessivage 3197 IF (iflag_lscav .EQ. 3 .OR. iflag_lscav .EQ.4) THEN3197 IF (iflag_lscav == 3 .OR. iflag_lscav == 4) THEN 3198 3198 !IF (.false.) THEN ! test #DFB (Binta) sans lsc_scav_spl 3199 3199 print *,'JE iflag_lscav',iflag_lscav … … 3362 3362 CALL satellite_out_spla(jD_cur,jH_cur,pdtphys,rlat,rlon, & 3363 3363 masque_aqua_cur, masque_terra_cur ) 3364 IF (jH_cur-pdtphys/86400. .LT.0.) THEN3364 IF (jH_cur-pdtphys/86400. < 0.) THEN 3365 3365 !new utc day: put in 0 everything 3366 3366 !JE20150518<< … … 3470 3470 ENDDO 3471 3471 3472 IF (jH_cur+pdtphys/86400. .GE. 1.) THEN3472 IF (jH_cur+pdtphys/86400. >= 1.) THEN 3473 3473 ! print *,'last step of the day' 3474 3474 DO i=1,klon 3475 IF (masque_aqua(i) .GT.0) THEN3475 IF (masque_aqua(i)> 0) THEN 3476 3476 aod550_aqua(i)=aod550_aqua(i)/masque_aqua(i) 3477 3477 aod670_aqua(i)=aod670_aqua(i)/masque_aqua(i) … … 3506 3506 aod865_dustsco_aqua(i)= -999. 3507 3507 ENDIF 3508 IF (masque_terra(i) .GT.0) THEN3508 IF (masque_terra(i)> 0) THEN 3509 3509 aod550_terra(i)=aod550_terra(i)/masque_terra(i) 3510 3510 aod670_terra(i)=aod670_terra(i)/masque_terra(i) … … 3635 3635 fluxss(:)=0.0 3636 3636 DO i=1, klon 3637 IF (iregion_ind(i) .GT.0) THEN ! LAND3637 IF (iregion_ind(i)>0) THEN ! LAND 3638 3638 ! SULFUR EMISSIONS 3639 3639 fluxh2sff(i)= (lmt_so2ff_l(i)+lmt_so2ff_h(i))*frach2sofso2* & … … 3656 3656 fluxff(i)=fluxbcff(i)+fluxomff(i) 3657 3657 ENDIF 3658 IF (iregion_bb(i) .GT.0) THEN ! LAND3658 IF (iregion_bb(i)>0) THEN ! LAND 3659 3659 ! SULFUR EMISSIONS 3660 3660 fluxso2bb(i) =scale_param_bb(iregion_bb(i)) * fracso2emis * & … … 4515 4515 ENDIF 4516 4516 4517 IF (test_sca .EQ.0 ) THEN4517 IF (test_sca == 0 ) THEN 4518 4518 ! READ file!! 4519 4519 call read_scalenc(filescaleparams,paramname_ind, & … … 4556 4556 4557 4557 jH_sca=jH_sca+pdtphys/(24.*3600.) 4558 IF (jH_sca .GT.(sca_resol)/24.) THEN4558 IF (jH_sca>(sca_resol)/24.) THEN 4559 4559 test_sca=0 4560 4560 jH_sca=jH_ini … … 4568 4568 USE mod_grid_phy_lmdz 4569 4569 USE mod_phys_lmdz_para 4570 USE lmdz_netcdf, ONLY:nf_open,nf_close,nf_inq_varid,nf_nowrite,nf_noerr,nf90_get_var 4570 4571 IMPLICIT NONE 4571 4572 include "netcdf.inc"4573 4572 4574 4573 CHARACTER*800 filescaleparams … … 4589 4588 !nci=NCOPN(trim(adjustl(filescaleparams)),NCNOWRIT,rcode) 4590 4589 ierr = NF_OPEN (trim(adjustl(filescaleparams)),NF_NOWRITE, nid) 4591 if (ierr .EQ.NF_NOERR) THEN4590 if (ierr == NF_NOERR) THEN 4592 4591 debutread=step_sca 4593 4592 countread=1 … … 4598 4597 print *,varname 4599 4598 ierr = NF_INQ_VARID (nid,trim(adjustl(varname)), nvarid) 4600 ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debutread, & 4601 countread, auxreal) 4602 IF (ierr .NE. NF_NOERR) THEN 4599 ierr = nf90_get_var (nid, nvarid, auxreal, debutread, countread) 4600 IF (ierr /= NF_NOERR) THEN 4603 4601 PRINT*, 'Pb de lecture pour modvalues' 4604 4602 print *,'JE scale_var, step_sca',trim(adjustl(varname)),step_sca -
LMDZ6/trunk/libf/phylmd/Dust/read_dust.F
r5073 r5075 3 3 USE mod_grid_phy_lmdz 4 4 USE mod_phys_lmdz_para 5 USE lmdz_netcdf, ONLY:nf90_get_var 5 6 IMPLICIT NONE 6 7 c 7 8 INCLUDE "dimensions.h" 8 9 INCLUDE "paramet.h" 9 INCLUDE "netcdf.inc"10 10 c 11 11 INTEGER step, nbjour … … 45 45 c 46 46 start(3)=step 47 c 48 ! status=NF_GET_VARA_DOUBLE(ncid1,varid1,start,count,dust_nc) 47 49 48 status=nf90_get_var(ncid1,varid1,dust_nc_glo,start,count) 50 c 49 51 50 ! call correctbid(iim,jjp1,dust_nc) 52 51 call correctbid(nbp_lon,nbp_lat,dust_nc_glo) -
LMDZ6/trunk/libf/phylmd/Dust/read_surface.F90
r5073 r5075 10 10 USE mod_phys_lmdz_para 11 11 USE iophy 12 ! USE netcdf 12 USE lmdz_netcdf, ONLY:nf_inq_varid,nf_noerr,nf90_get_var 13 13 IMPLICIT NONE 14 14 15 INCLUDE "netcdf.inc"16 15 INCLUDE "dimensions.h" 17 16 INCLUDE "paramet.h" -
LMDZ6/trunk/libf/phylmd/Dust/read_vent.F
r5073 r5075 3 3 USE mod_grid_phy_lmdz 4 4 USE mod_phys_lmdz_para 5 USE lmdz_netcdf, ONLY: nf90_get_var 5 6 ! USE write_field_phy 6 7 IMPLICIT NONE … … 8 9 c INCLUDE "dimphy.h" 9 10 INCLUDE "paramet.h" 10 INCLUDE "netcdf.inc"11 11 c 12 12 INTEGER step, nbjour
Note: See TracChangeset
for help on using the changeset viewer.
