²é¿´: 3532  |  »Ø¸´: 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:29:55
ÓàÔó³É(½ð±Ò+1): лл²ÎÓëÓ¦Öú£¡ 2011-08-24 16:22:51
ÖÁÉÙ
      integer mnum Nanum
      character tile cellp
»¹ÊDz»ÕýÈ·µÄ°É¡­¡­
ÁíÍ⣬\ ËÆºõÒ²ÓÐÎÊÌâ¡­¡­
7Â¥2011-08-23 15:10:46
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤£¨0856£©304Çó BÇø µ÷¼Á +3 Çñgl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[¿¼ÑÐ] ³õÊÔ 317 +3 °ëÀ­Ô±û 2026-03-20 3/150 2026-03-21 13:25 by zhukairuo
[¿¼ÑÐ] Çóµ÷¼Á +3 Ê®Èý¼ÓÓÍ 2026-03-21 3/150 2026-03-21 13:18 by zhukairuo
[¿¼ÑÐ] 336Çóµ÷¼Á +3 rmc8866 2026-03-21 3/150 2026-03-21 11:16 by Çóµ÷¼Ázz
[¿¼ÑÐ] ²ÄÁÏ 271Çóµ÷¼Á +3 Õ¹ÐÅÔÃ_ 2026-03-21 3/150 2026-03-21 10:17 by ÍõЦÓîÒ»¶¨ÊÇÑо
[¿¼ÑÐ] Äϲý´óѧ²ÄÁÏר˶311·ÖÇóµ÷¼Á +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸Ìì½ò´óѧ»¯Ñ§¹¤ÒÕרҵ£¨081702£©315·ÖÇóµ÷¼Á +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄϾ©Àí¹¤´óѧ085701×ÊÔ´Óë»·¾³302·ÖÇóµ÷¼Á +4 ¿ûè÷ÎÀ¶Ó 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤ 322Çóµ÷¼Á +4 È»11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[¿¼ÑÐ] 350Çóµ÷¼Á +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[¿¼ÑÐ] 329Çóµ÷¼Á +9 ÏëÉÏѧ߹߹ 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[¿¼ÑÐ] Çóµ÷¼ÁÒ»Ö¾Ô¸ÄϾ©º½¿Õº½Ìì´óѧ289·Ö +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[¿¼²©] ÕÐÊÕ²©Ê¿1-2ÈË +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by ßÉßÉßÉßɽÐ
[¿¼ÑÐ] 081700»¯¹¤Ñ§Ë¶µ÷¼Á +3 ¡¾1¡¿ 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[¿¼ÑÐ] 328Çóµ÷¼Á£¬Ó¢ÓïÁù¼¶551£¬ÓпÆÑо­Àú +4 ÉúÎ﹤³Ìµ÷¼Á 2026-03-16 12/600 2026-03-19 11:10 by ÉúÎ﹤³Ìµ÷¼Á
[¿¼ÑÐ] 311Çóµ÷¼Á +11 ¶¬Ê®Èý 2026-03-15 12/600 2026-03-18 14:36 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] Ò»Ö¾Ô¸ËÕÖÝ´óѧ²ÄÁϹ¤³Ì£¨085601£©×¨Ë¶ÓпÆÑо­ÀúÈýÏî¹ú½±Á½¸öʵÓÃÐÍרÀûÒ»ÏîÊ¡¼¶Á¢Ïî +6 ´ó»ðɽС»ðɽ 2026-03-16 8/400 2026-03-17 15:05 by ÎÞи¿É»÷111
[¿¼ÑÐ] 275Çóµ÷¼Á +4 Ì«Ñô»¨ÌìÌ쿪ÐÄ 2026-03-16 4/200 2026-03-17 10:53 by ¹¦·ò·è¿ñ
[¿¼ÑÐ] 321Çóµ÷¼Á +5 ´óÃ×·¹£¡ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[¿¼ÑÐ] 0856Çóµ÷¼Á +3 ÁõÃÎ΢ 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û