- Timestamp:
- May 13, 2015, 2:58:47 PM (9 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- d486c43
- Parents:
- 78e62dc
- Location:
- src
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
src/readreleases.f90
r78e62dc r0f20c31 376 376 if (lroot) then 377 377 write (*,*) 'In-cloud scavenging: ON' 378 write (*,*) 'In-cloud scavenging coefficients: ',&379 &weta_in(i),wetb_in(i),i !,wetc_in(i), wetd_in(i),i378 ! write (*,*) 'In-cloud scavenging coefficients: ',& 379 ! &weta_in(i),wetb_in(i),i !,wetc_in(i), wetd_in(i),i 380 380 end if 381 381 else -
src/readspecies.f90
r78e62dc r0f20c31 262 262 else !is gas 263 263 if (lroot) then 264 write(*, '(a,f5.2)') ' Gas below-cloud scavenging parameter A ', weta(pos_spec)264 write(*,*) ' Gas below-cloud scavenging parameter A ', weta(pos_spec) 265 265 write(*,'(a,f5.2)') ' Gas below-cloud scavenging parameter B ', wetb(pos_spec) 266 266 write(*,*) ' Gas in-cloud scavenging uses default values as in Hertel et al 1995' -
src/verttransform.f90
r8a65cb0 r0f20c31 510 510 511 511 if (readclouds) write(*,*) 'using cloud water from ECMWF' 512 if (.not.readclouds) write(*,*) 'using cloud water from parameterization'512 ! if (.not.readclouds) write(*,*) 'using cloud water from parameterization' 513 513 514 514 rcw(:,:)=0 -
src/wetdepo.f90
r5f9d14a r0f20c31 264 264 if (clouds_v.ge.4) then !below cloud 265 265 266 if (weta(ks).gt.0. .or. wetb(ks).gt.0 ) then !if positive below-cloud parameters given in SPECIES file (either A or B)266 if (weta(ks).gt.0. .or. wetb(ks).gt.0.) then !if positive below-cloud parameters given in SPECIES file (either A or B) 267 267 blc_count=blc_count+1 268 268 269 269 270 !GAS … … 279 280 280 281 !ZHG snow or rain removal is applied based on the temperature. 281 if (act_temp .ge. 273 ) then !Rain282 if (act_temp .ge. 273 .and. weta(ks).gt.0.) then !Rain 282 283 283 284 !Particle RAIN scavenging coefficient based on Laakso et al 2003, the below-cloud scavenging (rain efficienty) parameter A (=weta) from SPECIES file … … 285 286 (log10(dquer_m))**(-2))+ (bclr(5)*(log10(dquer_m))**(-1))+ bclr(6)* (prec(1))**(0.5)) 286 287 287 elseif (act_temp .lt. 273 ) then !snow288 elseif (act_temp .lt. 273 .and. wetb(ks).gt.0.) then !snow 288 289 289 290 !Particle SNOW scavenging coefficient based on Kyro et al 2009, the below-cloud scavenging (Snow efficiency) parameter B (=wetb) from SPECIES file … … 305 306 if (clouds_v.lt.4) then !in-cloud 306 307 307 inc_count=inc_count+1 308 ! NIK 13 may 2015: only do incloud if positive in-cloud scavenging parameters are given in species file 309 if (weta_in(ks).gt.0. .or. wetb_in(ks).gt.0.) then !if positive in-cloud parameters given in SPECIES file (either Ai or Bi) 310 311 ! if negative coefficients (turned off) set to zero for use in equation 312 if (weta_in(ks).lt.0.) weta_in(ks)=0. 313 if (wetb_in(ks).lt.0.) wetb_in(ks)=0. 314 315 inc_count=inc_count+1 308 316 309 317 !ZHG liquid water parameterization (CLWC+CIWC) 310 if (readclouds) then !get cloud water clwc & ciwc units Kg/Kg 311 cl=clwc(ix,jy,hz,n)+ciwc(ix,jy,hz,n) 312 else !parameterize cloudwater 313 cl=2E-7*prec(1)**0.36 314 endif 315 316 ! AEROSOL 317 if (dquer(ks).gt. 0.) then ! is particle 318 if (act_temp .le. 253) then 319 liq_frac=0 320 else if (act_temp .ge. 273) then 321 liq_frac=1 322 else 323 liq_frac =((act_temp-273)/(273-253))**2 318 if (readclouds) then !get cloud water clwc & ciwc units Kg/Kg 319 cl=clwc(ix,jy,hz,n)+ciwc(ix,jy,hz,n) 320 else !parameterize cloudwater 321 cl=2E-7*prec(1)**0.36 324 322 endif 323 324 if (act_temp .le. 253) then 325 liq_frac=0 326 else if (act_temp .ge. 273) then 327 liq_frac=1 328 else 329 liq_frac =((act_temp-273)/(273-253))**2 330 endif 325 331 326 332 ! ZHG calculate the activated fraction based on the In-cloud scavenging parameters Ai (=weta_in) and Bi (=wetb_in) from SPECIES file … … 328 334 ! The default values are 0.9 for CCN and 0.1 IN 329 335 ! This parameterization is based on Verheggen et al. (2007) & Cozich et al. (2006) 330 frac_act = liq_frac*weta_in(ks) +(1-liq_frac)*wetb_in(ks)336 frac_act = liq_frac*weta_in(ks) +(1-liq_frac)*wetb_in(ks) 331 337 332 338 !ZHG Use the activated fraction and the liqid water to calculate the washout 333 S_i= frac_act/cl 339 340 ! AEROSOL 341 if (dquer(ks).gt. 0.) then ! is particle 342 343 S_i= frac_act/cl 334 344 335 345 ! GAS 336 else ! is gas 337 cle=(1-cl)/(henry(ks)*(r_air/3500.)*act_temp)+cl 338 S_i=1/cle 339 endif 346 else ! is gas 347 348 cle=(1-cl)/(henry(ks)*(r_air/3500.)*act_temp)+cl 349 S_i=frac_act/cle 350 351 endif ! gas or particle 340 352 341 353 ! scavenging coefficient based on Hertel et al 1995 - using the S_i for either gas or aerosol 342 wetscav=S_i*(prec(1)/3.6E6)/clouds_h354 wetscav=S_i*(prec(1)/3.6E6)/clouds_h 343 355 344 356 ! write(*,*) 'in-cloud, act_temp=',act_temp,',prec=',prec(1),',wetscav=',wetscav,',jpart=',jpart,',clouds_h=,', & 345 357 ! clouds_h,',cl=',cl, 'diff to old scheme=', cl-2E-7*prec(1)**0.36 346 358 endif ! positive in-cloud scavenging parameters given in Species file 347 359 endif !incloud 348 360
Note: See TracChangeset
for help on using the changeset viewer.