[RegCNET] how to modify albedo?

[Yan Libin]

Hi, Dear Regcmers,

I want to modify the albedo where elevation is higher than 3000.
In the albedov.F file,
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.
Here is my added codes,
if(ht(np)>=3000) then
       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
end if

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

Libin Yan

2009-10-10 



Yan Libin 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20091010/a7cc373f/attachment-0001.html>