I want to modify the albedo where elevation is
higher than 3000.
I wanna add some codes to reach the goal, but I am
not sure how to cite ht variable.
This programme include main file,so ht is included
in.
albvs(np) = albvs(np)/dble(NNSG)*0.95
albvl(np) = albvl(np)/dble(NNSG)*0.95
aldirs(np) = aldirs(np)/dble(NNSG)*0.95
aldirl(np) = aldirl(np)/dble(NNSG)*0.95
aldifs(np) = aldifs(np)/dble(NNSG)*0.95
aldifl(np) = aldifl(np)/dble(NNSG)*0.95
I want to know how to cite 'ht' in the right
way.
Please give me a hand.
Thank you in advance.
c
c ***************not part of albedo in the ccm ****************
c
aldirs_s(n)=(1.-veg1d(n,np))*albgs +veg1d(n,np)*albs
aldirl_s(n)=(1.-veg1d(n,np))*albgl +veg1d(n,np)*albl
aldifs_s(n)=(1.-veg1d(n,np))*albgsd+veg1d(n,np)*albsd
aldifl_s(n)=(1.-veg1d(n,np))*albgld+veg1d(n,np)*albld
enddo
albvs(np) = albvs_s(1)
albvl(np) = albvl_s(1)
aldirs(np) = aldirs_s(1)
aldirl(np) = aldirl_s(1)
aldifs(np) = aldifs_s(1)
aldifl(np) = aldifl_s(1)
do n=2,NNSG
albvs(np) = albvs(np) +albvs_s(n)
albvl(np) = albvl(np) +albvl_s(n)
aldirs(np) = aldirs(np)+aldirs_s(n)
aldirl(np) = aldirl(np)+aldirl_s(n)
aldifs(np) = aldifs(np)+aldifs_s(n)
aldifl(np) = aldifl(np)+aldifl_s(n)
enddo
albvs(np) = albvs(np)/dble(NNSG)
albvl(np) = albvl(np)/dble(NNSG)
aldirs(np) = aldirs(np)/dble(NNSG)
aldirl(np) = aldirl(np)/dble(NNSG)
241,0-1 96%
aldifl(np) = aldifl_s(1)
do n=2,NNSG
albvs(np) = albvs(np) +albvs_s(n)
albvl(np) = albvl(np) +albvl_s(n)
aldirs(np) = aldirs(np)+aldirs_s(n)
aldirl(np) = aldirl(np)+aldirl_s(n)
aldifs(np) = aldifs(np)+aldifs_s(n)
aldifl(np) = aldifl(np)+aldifl_s(n)
enddo
albvs(np) = albvs(np)/dble(NNSG)
albvl(np) = albvl(np)/dble(NNSG)
aldirs(np) = aldirs(np)/dble(NNSG)
aldirl(np) = aldirl(np)/dble(NNSG)
aldifs(np) = aldifs(np)/dble(NNSG)
aldifl(np) = aldifl(np)/dble(NNSG)
c ****** fsw1d(np),sabveg(np),solis(np) computed in colrad
90 continue !end np=np1,npts loop
return
end