²é¿´: 363  |  »Ø¸´: 0

lxz73

Í­³æ (СÓÐÃûÆø)

[ÇóÖú] ¸ßÊÖ°ï°ïæ°¡

Õâ¸ö³ÌÐòÓ¦¸ÃÊǼÆËãÃܶȷֲ¼µÄ µ«ÊÇ¿´²»¶®£¬Çë½Ì¸ßÊÖ
¶øÇÒ±àÒëµÄʱºò»¹Óо¯¸æ£¡
c program calculate density profile

c  density profiles : (1) polymer; (2) nano; (3) solvent; (4) total
        parameter (NdpdMax=1400000)
        parameter (Ndenpro=500)
        parameter (nGn=15)

         implicit double precision (a-h,o-z)
        character*25 file_in,file_out

       real*8 rdpd(3,NdpdMax),vdpd(3,NdpdMax),rdpdcont(3,NdpdMax)
         real*8 denpro(4,Ndenpro)
         real*8 dz
         real*8 s_size(3),cell(3),sd,rho,D0
        Integer*4 species(NdpdMax)

        Integer*4 n_dpd_polymer,n_dpd_ps,n_dpd_nano
        integer*4 mdenpro
        Real*8 Ene(-1:3,-1:3),elsc,d_els,rm,T_solv,r_eq

                open(2, file='density.in')
              read(2,*)file_in, file_out, dz
            close(2)

            open(40, file=file_in)
                read(40,*)itt,ndpd,npoly,lchain,lchainA,n_dpd_ps
     &                ,n_dpd_nano
            read(40,*)T_solv
            read(40,*)ENE(-1,-1),ENE(1,1),ENE(-1,1),ENE(0,-1)
     &        , ENE(0,1),ENE(0,0),ENE(2,2),ENE(1,2)
          read(40,*)ENE(-1,2), ENE(0,2),ENE(3,3),ENE(2,3)
     &        ,  ENE(1,3),ENE(-1,3),ENE(0,3)
            read(40,*)s_size,D0

                D0=D0+2.75
                dv=s_size(1)*s_size(2)*dz

                n_dpd_polymer=npoly*lchain

                mdenpro=D0/dz+0.5
            do jn=1,mdenpro
                        denpro(1,jn)=0.0
                        denpro(2,jn)=0.0
                        denpro(3,jn)=0.0
                        denpro(4,jn)=0.0
            enddo
         
                do i=1,n_dpd_ps
            read(40,*)rdpd(1,i),rdpd(2,i),rdpd(3,i)
     >           ,vdpd(1,i),vdpd(2,i),vdpd(3,i)
     &           ,rdpdcont(1,i),rdpdcont(2,i),rdpdcont(3,i)
     &           ,species(i)

           m=rdpd(3,i)/dz+1

          denpro(4,m)=denpro(4,m)+1


         if(i.le.n_dpd_polymer) denpro(1,m)=denpro(1,m)+1
         if(i.gt.n_dpd_polymer.and.i.le.(n_dpd_polymer+n_dpd_nano)) then
           denpro(2,m)=denpro(2,m)+1
         endif
         if(i.gt.(n_dpd_polymer+n_dpd_nano).and.i.le.n_dpd_ps) then
           denpro(3,m)=denpro(3,m)+1
         endif

         enddo
                       
           open(13, file=file_out)
        do k=1,mdenpro
             write(13, '(f10.3,4f18.7)')(k-0.5)*dz,
     &        denpro(1,k)/dv,
     &            denpro(2,k)/dv,
     &        denpro(3,k)/dv,
     &        denpro(4,k)/dv
        enddo



       stop
        end
»Ø¸´´ËÂ¥
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÖÇÄÜ»úÆ÷ÈË

Robot (super robot)

ÎÒÃǶ¼°®Ð¡Ä¾³æ

ÕÒµ½Ò»Ð©Ïà¹ØµÄ¾«»ªÌû×Ó£¬Ï£ÍûÓÐÓÃŶ~

¿ÆÑдÓСľ³æ¿ªÊ¼£¬ÈËÈËΪÎÒ£¬ÎÒΪÈËÈË
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ lxz73 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[»ù½ðÉêÇë] ½ñÌìÎñί»á¿ªÍêÁË£¬Ã÷Ìì³ö½á¹ûÂð +14 angus9576 2026-08-25 17/850 2026-08-25 20:24 by huagongfeihu
[»ù½ðÉêÇë] ÓÐûÓдóÉñ°ïÎÒ¿´¿´»ù½ð´úÂë 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[»ù½ðÉêÇë] ·Å°ñǰµÄ²»µ­¶¨ 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[»ù½ðÉêÇë] ÔÚ¼á±ù»¹¸Ç×ű±º£µÄʱºò£¬ÎÒ¿´µ½ÁËÅ­·ÅµÄ÷»¨¡£ +6 ziyangfang 2026-08-25 8/400 2026-08-25 16:28 by ±¦±´µÄÖí
[»ù½ðÉêÇë] ½ñÈÕ²»·Å°ñ£¿Íø´«¹ú×ÔȻԤ¼Æ 8 Ô 27 Èտɲé½á¹û +17 ҽѧÀÏÄк¢ 2026-08-20 22/1100 2026-08-25 15:36 by ҽѧÀÏÄк¢
[»ù½ðÉêÇë] ÈËÆø²»ÐÐÁË +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by ¹Â¶ÀµÄÓ¢ÐÛ6
[½Ìʦ֮¼Ò] µ¼Ê¦Í²ۣºÎÒÔõô̯ÉÏÁËÕâô¸ö¼«Æ·Ñо¿Éú£¡ +3 ËÕ¶«ÆÂ¶þÊÀ 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[»ù½ðÉêÇë] 2026ÄêµÄ¹ú¼ÒÉç¿Æ»ù½ðÏîĿͨѶÆÀÉóµÄйæÔòÓëж¯Ïò¡¢ÐÂÌôÕ½ +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[»ù½ðÉêÇë] ÎÒÃæÉÏÍêµ°ÁË +13 ÇÒÌý»¢Ð¥ 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[»ù½ðÉêÇë] ·¶½øÖоÙÒ»ÎĵÄÖÐÐÄ˼Ïë +7 Ñ׻ƹóëÐ 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[»ù½ðÉêÇë] ÄÜ·ñÍ˳ö²ÎÓëµÄÃæÉÏÏîÄ¿½â³ýÏÞÏî +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by ¼ÒÓëÔ¶·½
[»ù½ðÉêÇë] filecode£¬4¸öjtjcÁË +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by ¹þ¹þ¸ò£¿
[»ù½ðÉêÇë] ÈÃÎÒÖÐÒ»¸öÃæÉϰɣ¡ +13 ´óƼ1987 2026-08-20 16/800 2026-08-24 10:23 by ̫ɵÁË
[»ù½ðÉêÇë] ¿ÆÑй¶ùÌ«ÄÑÁË +18 ÎÒ4´ó°×²Ë 2026-08-20 19/950 2026-08-24 09:47 by ¿­¶÷¹ãÊ¢´ß»¯·ÖÎ
[½Ìʦ֮¼Ò] Ìø²ÛºóÔÚÑÐÏîÄ¿Ôõô°ì£¿ +5 ¼òµ¥»¯xn 2026-08-22 10/500 2026-08-23 12:38 by ¼òµ¥»¯xn
[»ù½ðÉêÇë] ¿´À´½ñÌì²»»á·Å°ñÁË£¿ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[»ù½ðÉêÇë] ʱ¼ä´ÁÓÖ±äÁË +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by ×Ïɼ´¼
[»ù½ðÉêÇë] Ó¦¸ÃÊÇÏÂÖÜÈý26ÈÕ¹«²¼Á˰ɣ¿ +4 ¹þ¹þ¸ò£¿ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[»ù½ðÉêÇë] ½ñÌì·Å°ñûϷÁ衃 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[»ù½ðÉêÇë] »ù½ð°¡»ù½ð +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û