+ if (lang.gt.0 .and. .not.surfarea) then
+ do i=nnt,nct-1
+ mnum=(molnum(i))
+ stdforcp(i)=stdfp(mnum)*dsqrt(gamp(mnum))
+! write(iout,*) "stdforcp=",stdforcp(i),itype(i,mnum),i
+ enddo
+ do i=nnt,nct
+ mnum=molnum(i)
+ if (itype(i,mnum).ne.ntyp1) stdforcsc(i)=stdfsc(iabs(itype(i,mnum)),mnum)&
+ *dsqrt(gamsc(iabs(itype(i,mnum)),mnum))
+! write(iout,*) "stdforcsc=",stdforcsc(i),itype(i,mnum),i
+ enddo
+ endif