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