Fortran代码逆向推导输入文件结构及逻辑问询

代码逐行解析 #

1. 文件打开

   OPEN(10, FILE=TRIM(10)//'.txt', STATUS='OLD')
   

   打开通道号为10的已存在文件，文件名由某个字符串变量（代码中TRIM(10)疑似笔误，应为具体变量名）拼接.txt后缀得到，STATUS='OLD'要求文件必须预先存在。

2. 跳过文件头

   READ(10, *)     ! header
   

   读取并丢弃文件第一行，该行作为文件头信息，不参与后续数据处理。

3. 温度数组（TEMP）处理

   • 读取温度数量标识NB_TEMP，并根据其绝对值分配数组存储空间：

     READ(10, *) NB_TEMP
     ALLOCATE(TEMP(abs(NB_TEMP)))
     

   • 若NB_TEMP < 0：
     将NB_TEMP转为绝对值，读取温度起始值TEMP_0和步长TEMP_D，通过循环生成NB_TEMP个连续温度值存入数组：

     NB_TEMP = ABS(NB_TEMP)
     READ(10, *) TEMP_0, TEMP_D
     TEMP(1) = TEMP_0
     DO I=2, NB_TEMP
       TEMP(I) = TEMP(I-1) + TEMP_D
     ENDDO
     

   • 若NB_TEMP > 0：
     直接读取NB_TEMP个温度值，直接存入数组：

     READ(10,*) TEMP(:)
     

4. 角度数组（PRS）处理
   逻辑与温度数组完全一致，仅增加了角度转弧度的步骤：

   • 读取角度数量标识NB_PRS，分配数组存储空间：

     READ(10, *) NB_PRS
     ALLOCATE(PRS(ABS(NB_PRS)))
     

   • 若NB_PRS < 0：
     转为绝对值后读取起始角度PRS_0和步长PRS_D，生成角度值后统一转为弧度（乘以PI/180）：

     NB_PRS = ABS(NB_PRS)
     READ(10, *) PRS_0, PRS_D
     PRS(1) = PRS_0*PI/180.
     DO I=2, NB_PRS
       PRS(I) = PRS(I-1) + PRS_D*PI/180.
     ENDDO
     

   • 若NB_PRS > 0：
     直接读取NB_PRS个角度值，再统一转为弧度：

     READ(10,*) PRS(:)
     DO I=1, NB_PRS
       PRS(I) = PRS(I)*PI/180.
     ENDDO
     

输入文件结构说明 #

输入文件为纯文本格式，结构分为固定的几个部分，NB_TEMP和NB_PRS的正负决定了后续行的格式，支持两种模式混合使用：

模式1：生成模式（NB_TEMP/NB_PRS为负） #

[任意文本的文件头行]
-3
25.0 5.0
-4
0.0 90.0

解释：

• 第一行：文件头（内容无限制，仅作标识）
• 第二行：-3表示要自动生成3个温度值
• 第三行：温度起始值25.0，步长5.0 → 最终数组为[25.0, 30.0, 35.0]
• 第四行：-4表示要自动生成4个角度值
• 第五行：角度起始值0.0，步长90.0 → 转为弧度后数组为[0.0, π/2, π, 3π/2]

模式2：直接输入模式（NB_TEMP/NB_PRS为正） #

[任意文本的文件头行]
3
25.0 30.0 35.0
4
0.0 90.0 180.0 270.0

解释：

• 第一行：文件头
• 第二行：3表示直接输入3个温度值
• 第三行：3个温度值（可空格分隔或每行一个）
• 第四行：4表示直接输入4个角度值
• 第五行：4个角度值（程序会自动转为弧度）

模式3：混合模式 #

[任意文本的文件头行]
-2
10.0 2.0
3
30.0 60.0 90.0

解释：

• 温度部分自动生成2个值：[10.0, 12.0]
• 角度部分直接输入3个值，转为弧度后存储

关键说明 #

• NB_PRS是从文件的新一行单独读取的，和NB_TEMP的值完全无关，不会重复读取同一个数值。
• 所有数值支持空格、换行等任意空白符分隔（Fortran自由格式读取会自动忽略空白）。

内容的提问来源于stack exchange，提问作者Jake Stravinsky