²é¿´: 2545  |  »Ø¸´: 9
±¾Ìû²úÉú 1 ¸ö ³ÌÐòÇ¿Ìû £¬µã»÷ÕâÀï½øÐв鿴
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

li_clifff

Òø³æ (ÕýʽдÊÖ)

[ÇóÖú] Çë½ÌÒ»¸öfortranС³ÌÐò±àÒë³ö´íµÄÎÊÌ⣬лл

ÈçϳÌÐò±àÒë²»³É¹¦£¬²»ÖªµÀÎÊÌâ³öÔÚÄÄÀmesh.datÎļþ¸ñʽÊǵÚÒ»ÐÐÊÇnumnp, numelµÄÊýÖµ£¬µÚ¶þÐе½×îºóÊÇmeshµÄ¸÷¸öËıßÐεľßÌå×ø±ê¡£´íÎóÐÅÏ¢ÊÇ£ºError: A specification expression object must be a dummy argument, a COMMON block object, or an object accessible through host or use association   [NUMNP]  integer*4 numnp, numel   
»¹ÓÐG:\Program Files\Microsoft Visual Studio\MyProjects\mesh\pmesh.for(22) : Error: An automatic object is invalid in a main program.   [X] dimension x(2,numnp),ia(numnp)£¬Çë¶à¶àÖ¸½Ì£¬Ð»Ð»                        


  PROGRAM pmesh
c
c-----------------------------------------------------------------
c     purpose : to eliminate the string from existed mesh data
c               warning : the mesh has to be made out of quadrangles only
c     in :c
c     out :
c             x     : array containing node coordinates
c             ia    : nodes reference number
c             ien   : list of element nodes
c             mat   : material ref. numbers for all elements
c             in    : element sides ref. number
cc-----------------------------------------------------------------------
c
      implicit real*8 (a-h,o-z)
      real xreal,yreal
      integer*4 numnp, numel                           
      dimension x(2,numnp),ia(numnp)
      dimension ien(4,numel),mat(numel),in(4,numel)
   
      open(57,file='pmesh1.dat')
      open(58,file='pmesh2.dat')
      
         call sizemesh(numnp, numel)
         open(21,file='mesh.dat')
          read(21,*)
          do 300 i=1,numnp                   
          read(21,2100)idum,xreal,yreal,ia(i)
              x(1,i) = xreal
              x(2,i) = yreal
          write(57,*)idum,xreal,yreal,ia(i)
300    continue
          do 500 n=1,numel
          read(21,2300)ndum,(ien(i,n),i=1,4),(in(i,n),i=1,4)
          write(58,*)ndum,(ien(i,n),i=1,4),(in(i,n),i=1,4)
500    continue
          close(21)

2100  format(i5,2x,e15.5,2x,e12.5,2x,i2)
2300  format(5(i5,2x),4(i2,2x))
          end
c**************************************************************


     subroutine sizemesh(numnp,numel)

        integer numnp, numel

      open(21,file='mesh.dat')
      read(21,*) numnp,numel
     close(21)
     return
     end
»Ø¸´´ËÂ¥
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

li_clifff

Òø³æ (ÕýʽдÊÖ)

ÎÒ·¢ÏÖÔçÆÚf77³ÌÐò£¬¶ÔÊý×éûÓж¯Ì¬¸³Öµ£¬¶¼ÊÇÏȶ¨ÒåÒ»¸ö´óµÄÊý×é±ÈÈç
x£¨nnpmax£©,nnpmaxÔÚparameter.incÎļþÀ¶øÊµ¼ÊÉÏ×Ó³ÌÐòÀïµÄ´óСnnpÊDZ仯µÄ£¬µ«Ö»Òª²»³¬¹ýnnpmax¾ÍÐУ¬ÕâÑùµÄ·½Ê½Ò²²»´í¡£
9Â¥2011-09-28 11:46:33
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 10 ¸ö»Ø´ð

snoopyzhao

ÖÁ×ðľ³æ (Ö°Òµ×÷¼Ò)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï
ÓàÔó³É(½ð±Ò+2): лл²ÎÓëÓ¦Öú£¡ 2011-08-24 16:26:52
ben_ladeng: ר¼Ò¿¼ºË´æµµ 2011-09-28 11:39:50
ÄãÖ÷³ÌÖÐµÄ numnp ºÍ numel ±ØÐ븳ֵ£¬·ñÔò±ØÐëʹÓö¯Ì¬Êý×飬ÄÇôÄãµÄ x, ia, ien, mat, in µÈÊý×é¾ÍÐèÒªÓà allocatable Ö®ÀàµÄ¹Ø¼ü´ÊÖØÐÂÉùÃ÷£¬È»ºóÔÚµ÷ÓÃ×Ó³ÌÐò sizemesh Ö®ºó£¬Óà allocated À´·ÖÅäÄÚ´æ¡£
2Â¥2011-08-24 11:13:08
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

snoopyzhao

ÖÁ×ðľ³æ (Ö°Òµ×÷¼Ò)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
li_clifff(½ð±Ò+20): ȷʵÊÇÕâÑù£¬¶¯Ì¬Êý×飬ѧϰÁË£¬Ð»Ð» 2011-08-24 11:36:10
ÓàÔó³É(½ð±Ò+3): ÐÁ¿àÁË£¡ 2011-08-24 16:27:14
dubo(½ð±Ò+5, ³ÌÐòÇ¿Ìû+1): ¶àлӦÖí 2011-08-27 16:42:39
ben_ladeng: ר¼Ò¿¼ºË´æµµ 2011-09-28 11:40:00
ÉÏÃæ´íÁËÒ»µã¶ù£¬²»ÊÇÓà allocated£¬¶øÊÇÓà allocate À´·ÖÅäÄÚ´æ¡£ÏÂÃæÊÇÐ޸ĺóµÄ³ÌÐò£º
CODE:
      PROGRAM pmesh
c
c-----------------------------------------------------------------
c     purpose : to eliminate the string from existed mesh data
c               warning : the mesh has to be made out of quadrangles only
c     in :c
c     out :
c             x     : array containing node coordinates
c             ia    : nodes reference number
c             ien   : list of element nodes
c             mat   : material ref. numbers for all elements
c             in    : element sides ref. number
cc-----------------------------------------------------------------------
c
      implicit real*8 (a-h,o-z)
      real xreal,yreal
      integer*4 numnp, numel                           
      real, allocatable :: x(:,:)
      integer*4, allocatable :: ien(:,:), ia(:), mat(:), in(:,:)
!     dimension x(2,numnp),ia(numnp)
!     dimension ien(4,numel),mat(numel),in(4,numel)
   
      open(57,file='pmesh1.dat')
      open(58,file='pmesh2.dat')
      
         call sizemesh(numnp, numel)
         allocate(x(2,numnp),ia(numnp),ien(4,numel),mat(numel),
     +             in(4,numel))
         open(21,file='mesh.dat')
          read(21,*)
          do 300 i=1,numnp                  
          read(21,2100)idum,xreal,yreal,ia(i)
              x(1,i) = xreal
              x(2,i) = yreal
          write(57,*)idum,xreal,yreal,ia(i)
300    continue
          do 500 n=1,numel
          read(21,2300)ndum,(ien(i,n),i=1,4),(in(i,n),i=1,4)
          write(58,*)ndum,(ien(i,n),i=1,4),(in(i,n),i=1,4)
500    continue
          close(21)

2100  format(i5,2x,e15.5,2x,e12.5,2x,i2)
2300  format(5(i5,2x),4(i2,2x))
          end
c**************************************************************


      subroutine sizemesh(numnp,numel)

        integer numnp, numel

      open(21,file='mesh.dat')
      read(21,*) numnp,numel
      close(21)
      return
      end

3Â¥2011-08-24 11:26:59
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

li_clifff

Òø³æ (ÕýʽдÊÖ)

¡ï
dubo(½ð±Ò+1): »¶Ó­ÌÖÂÛ 2011-08-27 16:41:30
µ«ÊÇÎÒÕæ¿´µ½ÓеÄf77³ÌÐòÖ÷³ÌÐò²»ÓÃallocatable£¬Ö±½Óµ÷ÓÃ×Ó³ÌÐòÉú³ÉµÄnumnp, numel¾Í¿ÉÒÔ£¬¾ÍÊǰ´ÎÒÄÇÑù¿ªÊ¼Ð´µÄ£¬²»ÖªµÀ×Ô¼º°´ÄÇÑùµÄ·½Ê½ÔõôÐв»Í¨ÁË
4Â¥2011-08-24 12:53:08
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Ò»Ö¾Ô¸¼ª´ó»¯Ñ§322Çóµ÷¼Á +3 17501029541 2026-03-23 4/200 2026-03-23 23:47 by Txy@872106
[¿¼ÑÐ] 341Çóµ÷¼Á(Ò»Ö¾Ô¸ºþÄÏ´óѧ070300) +5 ·¬ÇÑÍ·--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[¿¼ÑÐ] 279·ÖÇóµ÷¼Á Ò»Ö¾Ô¸211 +17 chaojifeixia 2026-03-19 19/950 2026-03-23 23:26 by ´ô´ôʦ½ã
[¿¼ÑÐ] 0854 ¿¼Ñе÷¼Á ÕÐÉúÁË£¡AI ·½Ïò +4 pk3725069 2026-03-19 16/800 2026-03-23 23:09 by Íô£¡£¿£¡
[¿¼ÑÐ] 269Çóµ÷¼Á +4 ÎÒÏë¶ÁÑÐ11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[¿¼ÑÐ] Çóµ÷¼ÁÒ»Ö¾Ô¸Î人Àí¹¤´óѧ²ÄÁϹ¤³Ì£¨085601£© +3 WW.' 2026-03-23 5/250 2026-03-23 17:18 by ·ãÒíljj
[¿¼ÑÐ] ÇóÀÏʦÊÕÎÒ +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[¿¼ÑÐ] ÉúÎïѧ071000 329·ÖÇóµ÷¼Á +5 ÎÒ°®ÉúÎïÉúÎﰮΠ2026-03-17 5/250 2026-03-22 16:42 by tcx007
[¿¼ÑÐ] 303Çóµ÷¼Á +5 °²ÒäÁé 2026-03-22 6/300 2026-03-22 12:46 by ËØÑÕÇã³Ç1988
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[¿¼ÑÐ] ²ÄÁÏѧ˶301·ÖÇóµ÷¼Á +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[¿¼ÑÐ] ¿¼Ñе÷¼Á +3 ºôºô£¿~+123456 2026-03-21 3/150 2026-03-21 20:04 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 307Çóµ÷¼Á +3 ÓàÒâÇä 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[¿¼ÑÐ] 22 350 ±¾¿Æ985Çóµ÷¼Á£¬ÇóÀϵÇÊÕÁô +3 ÀîéóÄÐ003 2026-03-20 3/150 2026-03-21 13:28 by ²«»÷518
[¿¼ÑÐ] Ò»Ö¾Ô¸Öк£Ñó²ÄÁϹ¤³Ìר˶330·ÖÇóµ÷¼Á +8 С²Ä»¯±¾¿Æ 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[¿¼ÑÐ] Çóµ÷¼Á£¬Ò»Ö¾Ô¸:ÄϾ©º½¿Õº½Ìì´óѧ´óѧ £¬080500²ÄÁÏ¿ÆÑ§Ó빤³Ìѧ˶£¬×Ü·Ö289·Ö +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[¿¼ÑÐ] 353Çóµ÷¼Á +3 À­¹³²»Ðí±ä 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[¿¼ÑÐ] ²ÄÁÏѧ˶318Çóµ÷¼Á +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[¿¼ÑÐ] ÉúÎïѧµ÷¼ÁÕÐÈË£¡£¡£¡ +3 ɽº£Ììá° 2026-03-17 4/200 2026-03-19 21:34 by ÔõôÊÍ»³
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +3 ÄÝÄÝninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û