²é¿´: 388  |  »Ø¸´: 1

zml90h

Òø³æ (СÓÐÃûÆø)

[ÇóÖú] ÇëÎÊһϣ¬Ó¦¸Ã´ÓÃüÁî´°¿ÚÊäÈëʲô£¬¸ñʽÊÇʲôÑù×ӵġ£

ÏÂÃæÒ»¶ÎFortran´úÂë»áÔÚÃüÁî´°¿ÚÌáʾÓû§ÊäÈëתÖáºÍת½Ç£¬µ«ÊÇÎÒ²»Ã÷°×ÊäÈëµÄ¸ñʽÊÇʲô£¬ÒòΪÎÒ¸Õ¸Õ½Ó´¥µ½FortranдµÄ³ÌÐò£¬ËùÒÔ¿´²»¶®¡£Òò´ËÔÚÕâÀïÇëÇó¸÷λ³æÓÑÃǰïСµÜ¿´¿´¡£Ð»Ð»ÁË¡£Ö÷ÒªÊÇǰ¼¸ÐС£

c  define rotation of lattice axes
c
    2 write(*,'(a)',advance='NO') ' axis (x/y/z/q) and degree : '
      read(*,*,end=3) axis,degree
      degr = 3.1415926535898d0 * degree / 180.0d0
      if(axis(1:1).eq.'q' .or. axis(1:1).eq.'Q') goto 3
      if(axis(1:1).eq.'x' .or. axis(1:1).eq.'X') then
         rotation(1,1) = one
         rotation(2,1) = zero
         rotation(3,1) = zero
         rotation(1,2) = zero
         rotation(2,2) = dcos(degr)
         rotation(3,2) = - dsin(degr)
         rotation(1,3) = zero
         rotation(2,3) = dsin(degr)
         rotation(3,3) = dcos(degr)
         call matrix_product(tmp_rotation,rotation,all_rotation)
         tmp_rotation = all_rotation
      elseif(axis(1:1).eq.'y' .or. axis(1:1).eq.'Y') then
         rotation(1,1) = dcos(degr)
         rotation(2,1) = zero
         rotation(3,1) = dsin(degr)
         rotation(1,2) = zero
         rotation(2,2) = one
         rotation(3,2) = zero
         rotation(1,3) = - dsin(degr)
         rotation(2,3) = zero
         rotation(3,3) = dcos(degr)
         call matrix_product(tmp_rotation,rotation,all_rotation)
         tmp_rotation = all_rotation
      elseif(axis(1:1).eq.'z' .or. axis(1:1).eq.'Z') then
         rotation(1,1) = dcos(degr)
         rotation(2,1) = - dsin(degr)
         rotation(3,1) = zero
         rotation(1,2) = dsin(degr)
         rotation(2,2) = dcos(degr)
         rotation(3,2) = zero
         rotation(1,3) = zero
         rotation(2,3) = zero
         rotation(3,3) = one
         call matrix_product(tmp_rotation,rotation,all_rotation)
         tmp_rotation = all_rotation
      endif
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

×öÒ»¸öÓµÓÐÕÜѧͷÄԵļ¼ÊõÈË
ÒÑÔÄ   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zml90h

Òø³æ (СÓÐÃûÆø)

ÈçºÎ·ÅÆúÇóÖú°¡ ÎÊÌâÒѾ­½â¾öÁË
×öÒ»¸öÓµÓÐÕÜѧͷÄԵļ¼ÊõÈË
2Â¥2013-12-19 11:19:41
ÒÑÔÄ   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ zml90h µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÖпƼ¼´óѧ£¬080502£¬354·ÖÇóµ÷¼Á +4 ÊØºòϦÑôCF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[¿¼ÑÐ] 344Çóµ÷¼Á +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[¿¼ÑÐ] 295Çóµ÷¼Á +3 Ò»Ö¾Ô¸¾©Çø211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[¿¼ÑÐ] 302Çóµ÷¼Á +10 ºôºôºô¡£¡£¡£¡£ 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[¿¼ÑÐ] 280Çóµ÷¼Á +6 ¹¾ààÏþÏþ 2026-03-18 7/350 2026-03-18 11:25 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] ÉúÎïѧ071000 329·ÖÇóµ÷¼Á +3 ÎÒ°®ÉúÎïÉúÎﰮΠ2026-03-17 3/150 2026-03-18 10:12 by macy2011
[¿¼ÑÐ] 268Çóµ÷¼Á +6 ¼òµ¥µã0 2026-03-17 6/300 2026-03-18 09:04 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] Áº³ÉΰÀÏʦ¿ÎÌâ×é»¶Ó­ÄãµÄ¼ÓÈë +8 һѼѼӴ 2026-03-14 10/500 2026-03-17 15:07 by һѼѼӴ
[¿¼ÑÐ] Ò»Ö¾Ô¸ËÕÖÝ´óѧ²ÄÁϹ¤³Ì£¨085601£©×¨Ë¶ÓпÆÑо­ÀúÈýÏî¹ú½±Á½¸öʵÓÃÐÍרÀûÒ»ÏîÊ¡¼¶Á¢Ïî +6 ´ó»ðɽС»ðɽ 2026-03-16 8/400 2026-03-17 15:05 by ÎÞи¿É»÷111
[¿¼ÑÐ] ²ÄÁϹ¤³Ìר˶274Ò»Ö¾Ô¸211Çóµ÷¼Á +6 Ñ¦ÔÆÅô 2026-03-15 6/300 2026-03-17 11:05 by ѧԱh26Tkc
[¿¼ÑÐ] ²ÄÁÏ080500µ÷¼ÁÇóÊÕÁô +3 Ò»¿Åmeteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[¿¼ÑÐ] ѧ˶285Çóµ÷¼Á +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[¿¼ÑÐ] Çó²ÄÁϵ÷¼Á 085600Ó¢Ò»Êý¶þ×Ü·Ö302 ǰÈý¿Æ235 ¾«Í¨»úÆ÷ѧϰ Ò»Ö¾Ô¸¹þ¹¤´ó +4 ÁÖyaxin 2026-03-12 4/200 2026-03-13 22:04 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 0856²ÄÁÏÓ뻯¹¤301Çóµ÷¼Á +5 ÞÈÊø¹â 2026-03-13 5/250 2026-03-13 22:00 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] ¿¼Ñе÷¼Á +4 ·Ò´ï46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[¿¼ÑÐ] ¹¤¿Æ²ÄÁÏ085601 279Çóµ÷¼Á +8 À§ÓÚÐdz¿ 2026-03-12 10/500 2026-03-13 15:42 by ms629
[¿¼ÑÐ] 274Çóµ÷¼Á +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[¿¼ÑÐ] һ־Ըɽ´ó07»¯Ñ§ 332·Ö ËÄÁù¼¶Òѹý ±¾¿ÆÉ½¶«Ë«·Ç Çóµ÷¼Á£¡ +3 ²»ÏëÀíÄã 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[¿¼ÑÐ] 290Çóµ÷¼Á +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û