Changes in src/get_wetscav.f90 [e9e0f06:79e0349] in flexpart.git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/get_wetscav.f90
re9e0f06 r79e0349 151 151 if (height(il).gt.ztra1(jpart)) then 152 152 hz=il-1 153 ! goto 26 153 154 exit 154 155 endif 155 156 end do 157 !26 continue 158 156 159 157 160 if (ngrid.eq.0) then … … 200 203 201 204 202 !ZHG oct 2014 : Calculated for 1) both 2) lsp 3) convp - 2 and 3 not used removed by SE205 !ZHG oct 2014 : Calculated for 1) both 2) lsp 3) convp 203 206 ! Tentatively differentiate the grfraction for lsp and convp for treating differently the two forms 204 207 ! for now they are treated the same 205 208 grfraction(1)=max(0.05,cc*(lsp*lfr(i)+convp*cfr(j))/(lsp+convp)) 209 grfraction(2)=max(0.05,cc*(lfr(i))) 210 grfraction(3)=max(0.05,cc*(cfr(j))) 211 206 212 207 213 ! 2) Computation of precipitation rate in sub-grid cell 208 214 !****************************************************** 209 215 prec(1)=(lsp+convp)/grfraction(1) 216 prec(2)=(lsp)/grfraction(2) 217 prec(3)=(convp)/grfraction(3) 218 210 219 211 220 ! 3) Computation of scavenging coefficients for all species 212 221 ! Computation of wet deposition 213 222 !********************************************************** 223 214 224 215 225 if (ngrid.gt.0) then … … 227 237 !****************************************************************** 228 238 if ((dquer(ks).le.0.).and.(weta_gas(ks).gt.0..or.wetb_gas(ks).gt.0.)) then 239 ! if (weta(ks).gt.0. .or. wetb(ks).gt.0.) then 229 240 blc_count(ks)=blc_count(ks)+1 230 241 wetscav=weta_gas(ks)*prec(1)**wetb_gas(ks) … … 260 271 endif 261 272 273 ! write(*,*) 'bl-cloud, act_temp=',act_temp, ',prec=',prec(1),',wetscav=', wetscav, ', jpart=',jpart 274 262 275 endif ! gas or particle 263 276 ! endif ! positive below-cloud scavenging parameters given in Species file … … 272 285 if ((ccn_aero(ks).gt.0. .or. in_aero(ks).gt.0.).or.(henry(ks).gt.0.and.dquer(ks).le.0)) then 273 286 inc_count(ks)=inc_count(ks)+1 287 ! write(*,*) 'Incloud: ',inc_count 274 288 ! if negative coefficients (turned off) set to zero for use in equation 275 289 if (ccn_aero(ks).lt.0.) ccn_aero(ks)=0. … … 286 300 ! sec test 287 301 ! cl=1E6*1E-7*prec(1)**0.3 !Sec GFS new 288 cl=1E6*2E-7*prec(1)**0.36 !Sec ECMWF new , is also suitable for GFS302 cl=1E6*2E-7*prec(1)**0.36 !Sec ECMWF new 289 303 ! cl=2E-7*prec(1)**0.36 !Andreas 290 304 ! cl=1.6E-6*prec(1)**0.36 !Henrik … … 308 322 if (dquer(ks).gt.0.) then 309 323 S_i= frac_act/cl 324 ! write(*,*) 'Si: ',S_i 325 310 326 ! GAS 311 327 !**** 312 328 else 329 313 330 cle=(1-cl)/(henry(ks)*(r_air/3500.)*act_temp)+cl 331 !REPLACE to switch old/ new scheme 332 ! S_i=frac_act/cle 314 333 S_i=1/cle 315 334 endif ! gas or particle 316 335 317 336 ! scavenging coefficient based on Hertel et al 1995 - using the S_i for either gas or aerosol 318 !SEC wetscav fix, the cloud height is no longer needed, it gives wrong results 337 !OLD 338 if ((readclouds.and.ngrid.eq.0).or.(readclouds_this_nest.and.ngrid.gt.0)) then 319 339 wetscav=incloud_ratio*S_i*(prec(1)/3.6E6) 340 else 341 !SEC wetscav fix 342 wetscav=incloud_ratio*S_i*(prec(1)/3.6E6) 343 ! wetscav=incloud_ratio*S_i*(prec(1)/3.6E6)/clouds_h 344 endif 320 345 endif ! positive in-cloud scavenging parameters given in Species file 321 346 endif !incloud
Note: See TracChangeset
for help on using the changeset viewer.