计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2999|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[Fortran] fortran建立文件夹并按序建立文件的程序。

[复制链接 Copy URL]

132

帖子

0

威望

1025

eV
积分
1157

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
最近测试CCD(coupled cluster doubles approximation)计算效率时候发现双电子积分矩阵太大时(四个空轨道的双电子积分,在硬盘上大概有60个G),使用了‘direct’在硬盘上存储双电子积分,读取速度会变的很慢。也不知道是什么原因。。。
然后自己想尝试下在存储双电子积分时候直接多开一点文件。 在网上搜了搜教程之后自己凑了个代码,能实现创建文件夹并且在文件夹下按序列创建程序。但是能力有限,感觉写的特别丑。求大神指导怎样改进能写的漂亮点!
附上代码
==================
     program creat
      implicit none
      integer i
      character(len=50) folder, opt
      character(len=13) filename
      character(len=5) tmpc
      character(len=100) pathname
      folder= 'storage_'
      opt(1:9)='mkdir -p '
      opt(10:17)=folder
      filename(1:8)='test____'
      call system(adjustl(trim(opt)))  ! to create folder
      pathname(1:2)='./'
      pathname(3:10)=folder
      pathname(11:11)='/'
      write(6,*)'pathname', trim(pathname)
      do i = 1,100
         write(tmpc,'(i5.5)') i
         filename(1:8)='test____'
         filename(9:13)=tmpc
         write(6,*)'index of file', tmpc
         pathname(12:24)=filename
         open(unit=i+100,file=trim(pathname), form='unformatted',
     &          access='direct', recl=10000)
         close(i+100)
      enddo

      opt(1:6)='rm -r '
      opt(7:15)=folder
      call system(adjustl(trim(opt)))  ! to rm -r folder
      end

403

帖子

4

威望

2874

eV
积分
3357

Level 5 (御坂)

2#
发表于 Post on 2018-10-16 19:30:50 | 只看该作者 Only view this author
要不试试HDF或者NetCDF库?

132

帖子

0

威望

1025

eV
积分
1157

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2018-10-17 07:59:44 | 只看该作者 Only view this author
万里云 发表于 2018-10-16 19:30
要不试试HDF或者NetCDF库?

谢谢!下载下来学习下。
----------------------------------
https://www.hdfgroup.org/

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-23 18:58 , Processed in 0.272907 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list