²é¿´: 3526  |  »Ø¸´: 11
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

ºØÒÇ

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

[ÇóÖú] fortran±àÒë read³ö´í Íê±Ï

´íÎóÌáʾ
POSITION.F(31): error #6358: Constants and expressions are invalid in read-only I/O lists.   [MSP]
      read(15,*) (msp(i),i=1,10)
Ç°ÃæmspÓж¨Òå
character msp(10)
  read(15,*) (msp(i),i=1,10)
Óï·¨ÓÐÎÊÌâÂð£¿

Ò»ÏÂÊÇÈ«²¿³ÌÐò
      program getpos
ccccccccc  12 is SPPARKS OUTPUT  cccccccc
ccccccccc 15 is the molecure  ccccccccc
ccccccccc  13 is OUTPUT   ccccccccccc
      integer mnum nsp
      real coord(2000,3)
      real species(10,3)
      integer nsp(10)
      real x y z
      character tile cellp cella cellb cellz msp(10)  
      open(12,file='POSCAR',status='old')
      open(15,file='POSCAR',status='old')
      open(13,file='CONTCAR',status='new')
ccccccccccccccccccc sys param a b c cccccccccccccc
      read(12,*) tile
      write(13,*) tile
      read(12,*) cellp
      write(13,*) cellp
      read(12,*) cella ;
      write(13,*) cella
      read(12,*) cellb ;
      write(13,*) cellb
      read(12,*) cellc ;
      write(13,*) cellc
      read(15,*);
      read(15,*);
      read(15,*);
      read(15,*);
      read(15,*);
cccccccccccccccccccc atom char ccccccccccccccccccccccc
      read(15,*) (msp(i),i=1,10)
      write(13,*) (msp(i),i=1,10)
      read(12,*);
cccccccccccccccccccc atom num  cccccccccccccccccccccccc
      read(12,*) mnum
      read(15,*) (nsp(i),i=1,10)
      DO i=1, 10
         N_anum=mnum*nsp(i)
         write(13,N_anum)
      END DO
      read(12,*)
      read(15,*)
      write(13,*) 'R'
cccccccccccccccc  read molecue coord cccccccccccc
      DO i=1, mnum
         read(12,*) (coord(i,j),j=1,13)
      ENDDO
      close(12)
cc
      DO j=1, 10                !!! species cyc !!!!!!
            DO i=1, nsp(j)          !!! the j atom cyc !!!!!!
                  read(15,*) x y z  !!!!!!!!!!!!!!!!!!!!!!!! molecure coord !!!!!!!!!!!!!
                      DO k=1, mnum
                                nx=x+coord(k,1)
                                ny=y+coord(k,2)
                                nz=z+coord(k,3)
                                write(13,*) nx ,ny ,nz
                      ENDDO
            ENDDO
      ENDDO
          close(15)
          close(13)
      END

[ Last edited by ºØÒÇ on 2011-8-23 at 14:53 ]
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

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

snoopyzhao

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

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

¡ï
ºØÒÇ(½ð±Ò+5): ¶àлָ½Ì£¡ 2011-08-23 15:30:10
ÓàÔó³É(½ð±Ò+1): ллӦÖúµ½µ×£¬ºÇºÇ£¡ 2011-08-24 16:23:08
ÒýÓûØÌû:
8Â¥: Originally posted by ºØÒÇ at 2011-08-23 15:14:08:
ÎÒÓÃifort -O2 ÒѾ­Í¨¹ý±àÒëÁË
Äõ½×Ô¼ºÏëÒªµÄÎļþÁË ¾Í²»¹ÜÄÇô¶àÁË ºÇºÇ

ÄãÔËÐÐһϿ´¿´µÃµ½µÄ½á¹ûÊDz»ÊÇÄãÏëÒªµÄ£¬ºÇºÇ¡­¡­

ÔÚ Fortran ÖУ¬Ïñ real x y z ÕâÑùµÄÓï¾ä£¬Äã±¾À´ÉùÃûÁËÈý¸ö±äÁ¿£¬µ«±àÒëÖ®ºóÊÇÒ»¸ö±äÁ¿ xyz¡£Äã»òÐí»á˵£¬ÎªÊ²Ã´Ã»Óб¨´í£¬ÄãÊÔ×ÅÔÚ program ÄǾäÖ®ºó¼ÓÒ»¾ä implicit none ¿´¿´¡­¡­
9Â¥2011-08-23 15:16:59
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 12 ¸ö»Ø´ð

snoopyzhao

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

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

¡ï ¡ï
ºØÒÇ(½ð±Ò+5): ¶àлָ½Ì£¡ 2011-08-23 14:53:23
ÓàÔó³É(½ð±Ò+2): лл²ÎÓëÓ¦Öú£¡ 2011-08-24 16:20:44
ÒýÓûØÌû:
1Â¥: Originally posted by ºØÒÇ at 2011-08-23 09:57:29:
´íÎóÌáʾ
POSITION.F(31): error #6358: Constants and expressions are invalid in read-only I/O lists.   [MSP]
      read(15,*) (msp(i),i=1,10)
Ç°ÃæmspÓж¨Òå
character msp(10)
  read(15,*) (msp(i ...

character ÀàÐÍÊý¾Ý²»ÄÜÓà read(15, *) ÕâÖÖ·½Ê½À´¶ÁÈ¡°É¡£ * Ö»ÊʺÏÓÚÊýÖµ£¬¶øÇÒ¶Á³öÀ´µÄ¸ñʽֻÄÜÊǸ¡µãÊý¡­¡­

Èç¹ûÄãµÄ msp ÊÇÓÃÀ´´æµ¥¸ö×Ö·û£¨10¸ö£¿£©£¬ÄÇôÄã¿ÉÒÔÓÃÏÂÃæµÄ´úÂëÊÔÊÔ£¨Ã»ÓвâÊÔ£¬²»±£Ö¤ÕýÈ·¡­¡­£©
CODE:
read(15, '(a1)') (msp(i), i=1,10)

2Â¥2011-08-23 14:31:59
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

snoopyzhao

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

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

¡ï ¡ï
ºØÒÇ(½ð±Ò+5): ¶àлָ½Ì£¡ 2011-08-23 14:53:34
ÓàÔó³É(½ð±Ò+2): ÐÁ¿àÁË£¡ 2011-08-24 16:20:55
¿´À´ÎÒÇ°ÃæµÄ˵·¨Ã»ÓÐÇÐÈëÒªµã£¬ÄãµÄÎÊÌâÔÚÓÚ
      real x y z
      character tile cellp cella cellb cellz msp(10)  
ÉÏÃæÁ½Ðи÷¸ö±äÁ¿Ö®¼äÓ¦¸ÃÓöººÅÀ´·Ö¸ô¡­¡­

µ±È»£¬×Ö·ûµÄÈ·²»ÄÜÓà read(15,*) ÕâÖÖ·½Ê½À´¶Á¡­¡­
3Â¥2011-08-23 14:47:28
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ºØÒÇ

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

¡ï
ÓàÔó³É(½ð±Ò+1): лл²ÎÓëÓ¦Öú£¡ 2011-08-24 16:21:05
²»ÊǶÁÈ¡µÄÔ­Òò
ÔÚ¶¨ÒåµÄʱºòÓá°£º£º¡±ºó¾Í¿ÉÒÔÁË
à½À²À²à½À²À²à½
4Â¥2011-08-23 14:55:38
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] һ־Ըɽ´ó07»¯Ñ§ 332·Ö ËÄÁù¼¶Òѹý ±¾¿ÆÉ½¶«Ë«·Ç Çóµ÷¼Á£¡ +3 ²»ÏëÀíÄã 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[¿¼ÑÐ] 265Çóµ÷¼Á +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[¿¼ÑÐ] »¯Ñ§Çóµ÷¼Á +4 ÁÙÔó¾³llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[¿¼ÑÐ] 279·ÖÇóµ÷¼Á Ò»Ö¾Ô¸211 +11 chaojifeixia 2026-03-19 12/600 2026-03-21 01:49 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 330Çóµ÷¼Á +4 С²Ä»¯±¾¿Æ 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸Î人Àí¹¤²ÄÁϹ¤³Ìר˶µ÷¼Á +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖÐÄÏ»¯Ñ§£¨0703£©×Ü·Ö337Çóµ÷¼Á +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[¿¼ÑÐ] Çóµ÷¼ÁÒ»Ö¾Ô¸ÄϾ©º½¿Õº½Ìì´óѧ289·Ö +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤×¨Ë¶µ÷¼Á +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[¿¼ÑÐ] 085410È˹¤ÖÇÄÜר˶317Çóµ÷¼Á£¨0854¶¼¿ÉÒÔ£© +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by ²»168
[¿¼ÑÐ] 266Çóµ÷¼Á +5 ÑôÑôÍÛÈû 2026-03-14 10/500 2026-03-19 15:08 by ÑôÑôÍÛÈû
[¿¼ÑÐ] Ò»Ö¾Ô¸¸£´ó288Óлú»¯Ñ§£¬Çóµ÷¼Á +3 Сľ³æ200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[¿¼ÑÐ] ²ÄÁÏר˶306Ó¢Ò»Êý¶þ +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[¿¼ÑÐ] ²ÄÁÏ£¬·ÄÖ¯£¬ÉúÎ0856¡¢0710£©£¬»¯Ñ§ÕÐÉúÀ² +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[¿¼ÑÐ] 334Çóµ÷¼Á +3 Ö¾´æ¸ßÔ¶ÒâÔÚ»úÐ 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[¿¼ÑÐ] ²ÄÁÏר˶326Çóµ÷¼Á +6 Ä«ìÏæ¦Ý· 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[¿¼ÑÐ] 327Çóµ÷¼Á +6 ʰ¹âÈÎȾ 2026-03-15 11/550 2026-03-15 22:47 by ʰ¹âÈÎȾ
[¿¼ÑÐ] 26¿¼ÑÐÒ»Ö¾Ô¸ÖйúʯÓÍ´óѧ(»ª¶«)305·ÖÇóµ÷¼Á +3 ¼ÎÄêÐÂ³Ì 2026-03-15 3/150 2026-03-15 13:58 by ¹þ¹þ¹þ¹þºÙºÙºÙ
[¿¼ÑÐ] 289Çóµ÷¼Á +4 ÕâôÃû×ÖÕ¦Ñù 2026-03-14 6/300 2026-03-14 18:58 by userper
[¿¼ÑÐ] ¸´ÊÔµ÷¼Á +3 ºôºô£¿~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û