#include "ctrparam.h" subroutine rdata4clm(var) character *3 var #include "CLM.COM" dimension ps4clmc(lsmlon,lsmlat),pcpl4clmc(lsmlon,lsmlat), & pcpc4clmc(lsmlon,lsmlat),tpr4clmc(lsmlon,lsmlat), & tsl4clmc(lsmlon,lsmlat), & qs4clmc(lsmlon,lsmlat),ws4clmc(lsmlon,lsmlat) & ,us4clmc(lsmlon,lsmlat),vs4clmc(lsmlon,lsmlat), & dsw4clmc(lsmlon,lsmlat), & dlw4clmc(lsmlon,lsmlat),pco24clmc(lsmlon,lsmlat) & ,swinr4clmc(lsmlon,lsmlat),swvis4clmc(lsmlon,lsmlat) read (945),ps4clm,pcpl4clm, & pcpc4clm,tpr4clm, & tsl4clm, & qs4clm,ws4clm & ,us4clm,vs4clm, & dsw4clm, & dlw4clm,pco24clm & ,swinr4clm,swvis4clm if(var.ne.'ALL')then read (935),ps4clmc,pcpl4clmc, & pcpc4clmc,tpr4clmc, & tsl4clmc, & qs4clmc,ws4clmc & ,us4clmc,vs4clmc, & dsw4clmc, & dlw4clmc,pco24clmc & ,swinr4clmc,swvis4clmc endif if(var.eq.'TEM')then ! print *,'Temperature from control run' tsl4clm = tsl4clmc dsw4clm = dsw4clmc dlw4clm = dlw4clmc endif if(var.eq.'PRE'.or.var.eq.'PC2')then ! print *,'Precip from control run' do j=1,jm0 pcpl4clm(j) = pcpl4clmc(j) pcpc4clm(j) = pcpc4clmc(j) tpr4clm(j) = tpr4clmc(j) enddo endif return end