User Tools

Site Tools


rodos:alsmc_quickpro_adm_atstep.f_2008.05

All patch corrections in ALSMC Prognosis Atprogn.f are also carried out in Quickpro Atstep.f

Additional Correction of Depletion Velocities in Quickpro:

Parallel depletion to soil was the reason for depletion being 2 or 3 times too strong

Before Patch 2008.05:

TOTAL RESISTANCES for GRAS,AGRI,URBAN,SOIL,LAWN: it is assumed that RESATM of GRAS, SOIL, and LAWN are the same because of equal z0=0.05m. All atmospheric resistances RESAT are the resistances of the surrounding area roughness type IRAU=1 or 2, i.e. MOL or KAJUE:

For NN=2,3,4: (elem.Iodine, org.Iodine, Aerosols)
 
RSOILA(NN)=RSOIL(NN)+RESAT
RLAWNA(NN)=RLAWN(NN)+RESAT

TOTAL DEPOSITION VELOCITIES for AGRI,URBAN,LAWN: AGRICULTURAL or URBAN surface is used for DEPLETION of the cloud, LAWN surface is used for REFERENCE DRY DEPOSITION. The parallel resistance RSOILA considers depo to soil.

IF(IRAU.EQ.1) THEN
  RAGRIA(NN)=RAGRI(NN)+RESAT
  VDEPL(NN)=1./RAGRIA(NN) + 1./RSOILA(NN)
END IF
 
IF(IRAU.EQ.2) THEN
  RURBAA(NN)=RURBA(NN)+RESAT
  VDEPL(NN)=1./RURBAA(NN)
END IF
 
VLAWN(NN)=1./RLAWNA(NN) + 1./RSOILA(NN)

Now:

For NN=2,3,4:
 
RLAWNA(NN)=RLAWN(NN)+RESAT

The parallel resistance RSOILA is omitted now:

IF(IRAU.EQ.1) THEN
  RAGRIA(NN)=RAGRI(NN)+RESAT
  VDEPL(NN)=1./RAGRIA(NN)
END IF   
 
IF(IRAU.EQ.2) THEN
  RURBAA(NN)=RURBA(NN)+RESAT
  VDEPL(NN)=1./RURBAA(NN)
END IF
 
VLAWN(NN)=1./RLAWNA(NN)
rodos/alsmc_quickpro_adm_atstep.f_2008.05.txt · Last modified: 2015/03/24 15:32 (external edit)