关于CRTM的一些疑问

李璐含

carrol 发表于 2016-5-26 17:57
我是换了里面的数据，运行结果应该会变的，还要看你模拟的是什么波长，廓线对可见光波段的影响不是很大

好的，谢谢

李璐含

carrol 发表于 2016-5-26 17:57
我是换了里面的数据，运行结果应该会变的，还要看你模拟的是什么波长，廓线对可见光波段的影响不是很大

我有几个问题想请教：
1、在安装的时候里面有一步是link to the library，你在这一步的时候是改了环境变量吗？因为如果是超算的话我们没有进入home的权限，你用的是自己的电脑吗？是把环境变量里添加了自己的路径，还是直接把include和lib的两个文件夹复制到自己的home下面？
2、运行的时候是一定要所有的example都要一起运行吗？还是进入到某一个自己想要运行的下面（如：运行forward中的Example1_Simple就要进到这个文件夹里然后运行就可以了）
3、运行的语句是make还是make test，还是别的，因为我没看明白手册里真正运行的时候用的是什么，只是make test时都成功了而已，等到make的时候就失败了，所以想问一下运行的具体语句
4、改变数据是改test里Load_Sfc_Data.inc和Load_Atm_Data.inc这两个里面的数据吗？
5、运行的时候要改变每一个example文件夹里的Example1_Simple.f90吗？

问题有点多，不知道是不是打扰您了。谢谢您

carrol

李璐含 发表于 2016-6-2 13:23
我有几个问题想请教：
1、在安装的时候里面有一步是link to the library，你在这一步的时候是改了环境变 ...

1.我们是在服务器上运算的，环境变量加过路径。
2.直接进入到文件夹里，不用全部运行
3.进入到文件夹里，make一下，然后运行，手册里面make test是把test里面的都编译了。你需要在REL-2.1.3目录下. configure/ifort.setup一下，然后进入文件目录里面make 。
4.一般是改变那两个文件
5.你想运行哪个就改变那个.f90文件就行

李璐含

carrol 发表于 2016-6-2 14:07
1.我们是在服务器上运算的，环境变量加过路径。
2.直接进入到文件夹里，不用全部运行
3.进入到文件夹里 ...

不好意思，又有些问题想要向你请教，问题有点麻烦。
     1、当我改变Load_Atm_Data，不改变gfortran中sensorID时运行之后output并未发生改变，具体情况如下：
fortran程序中这一部分不变
! Get sensor id from user
  ! -----------------------
  WRITE( *,'(/5x,"Enter sensor id [hirs4_n18, amsua_metop-a, or mhs_n18]: ")',ADVANCE='NO' )
  READ( *,'(a)' ) Sensor_Id
  Sensor_Id = ADJUSTL(Sensor_Id)
  WRITE( *,'(//5x,"Running CRTM for ",a," sensor...")' ) TRIM(Sensor_Id)
运行之后，没有任何变化。
后来我将sensorID=‘crisB1_npp’，再次运行的时候，例子的文件夹里比之前多生成了一个bin文件（crisB1_npp.RTSolution.bin），这个我看results这个文件夹里也有，我觉得应该是生成一个bin文件才算是运行成功了，但是之前为什么不能生成bin文件呢，而且你有没有添加过什么步骤，可以同时选择多个sensorID的。
      2、在以上的基础上，我分别再一个单独的example的文件夹里make了一次，但是只运行的界面如下：
==============================================
Building forward Example5_ClearSky using:
  Compiler: ifort
  Switches:  -c  -O2  -convert big_endian  -free  -assume byterecl  
==============================================
我回到文件夹里发现，除了.o文件和另外一个和fortran同一名字的文件更新了以外，并没有生成新的output文件以及bin文件，即输出数据没发生变化。但是我跳回test文件夹所在的主目录之后 make test，在回到我想要的文件夹里就发现生成了新的output文件和bin文件。所以我想问一下，你在单独的例子文件夹里是否只用了单纯的make语句，make后面有没有添加文件夹的名字之类的。
        3、在正演模拟亮温的时候，你们选用的是再分析数据吗？还是什么其他的数据，是根据数据本身分层来调节了Load_Atm_Data的分层吗？
   谢谢你！

carrol

李璐含 发表于 2016-6-2 19:06
不好意思，又有些问题想要向你请教，问题有点麻烦。
     1、当我改变Load_Atm_Data，不改变gfortran中s ...

不知道你修改了里面的程序没有，直接转里面的Example时，它的结果一般都被输出到了屏幕上
  ! 6. **** OUTPUT THE RESULTS TO SCREEN ****
  !
  ! User should read the user guide or the source code of the routine
  ! CRTM_RTSolution_Inspect in the file CRTM_RTSolution_Define.f90 to
  ! select the needed variables for outputs.  These variables are contained
  ! in the structure RTSolution.
  DO m = 1, N_PROFILES
    WRITE( *,'(//7x,"Profile ",i0," output for ",a )') m, SENSOR_ID
    DO l = 1, n_Channels
      WRITE( *, '(/5x,"Channel ",i0," results")') RTSolution(l,m)%Sensor_Channel
      CALL CRTM_RTSolution_Inspect(RTSolution(l,m))
    END DO
  END DO

你说的生成bin文件那部分是下面的程序在做的
  ! 8. **** COMPARE RTSolution RESULTS TO SAVED VALUES ****
  !
  !   Step 8 is not part of the example to show how to use CRTM.
  !   It is to check the user results against the results in the CRTM package.
  !
  WRITE( *, '( /5x, "Comparing calculated results with saved ones..." )' )

  ! 8a. Create the output file if it does not exist
  ! -----------------------------------------------
  ! ...Generate a filename
  rts_File = RESULTS_PATH//SENSOR_ID//'.RTSolution.bin'
  ! ...Check if the file exists
  IF ( .NOT. File_Exists(rts_File) ) THEN
    Message = 'RTSolution save file does not exist. Creating...'
    CALL Display_Message( PROGRAM_NAME, Message, INFORMATION )
    ! ...File not found, so write RTSolution structure to file
    Error_Status = CRTM_RTSolution_WriteFile( rts_File, RTSolution, Quiet=.TRUE. )
    IF ( Error_Status /= SUCCESS ) THEN
      Message = 'Error creating RTSolution save file'
      CALL Display_Message( PROGRAM_NAME, Message, FAILURE )
      STOP
    END IF
这部分是另一个工作，没有做过。
一般情况下，要模拟亮温，辐射什么的是不会做这部分工作的，直接上面6.output the results 就行了。

修改了Load_Atm_Data文件，好像要先make clean一下，再make才生效，具体原因我不太清楚。Example5_ClearSky 里面就是没有上面说的8....那部分的工作，所以不会输出新的bin文件。你在test所在主目录make test后，是把所有的Example都编译了一遍。
一般，你只要进入你需要的目录中进行编译运行就行，make 不用加什么的。

我现在是用的再分析数据放到里面在模拟亮温，根据再分析资料进行了调节，不过出来的结果有些问题。

李璐含

carrol 发表于 2016-6-2 21:30
不知道你修改了里面的程序没有，直接转里面的Example时，它的结果一般都被输出到了屏幕上
  ! 6. **** O ...

你用的是ncep吗？如果要是提供的数据分好的气压层和水汽的分层有不一样怎么办？比如说如果气压分层是28层，水汽分层是14层，你是怎么处理的，还是你用的数据分层都是一样的

carrol

李璐含 发表于 2016-6-3 08:58
你用的是ncep吗？如果要是提供的数据分好的气压层和水汽的分层有不一样怎么办？比如说如果气压分层是28层 ...

我用的是ERA-interim的资料，分层是一样的，necp的资料我不太懂所以没用过，要不你插值试试

李璐含

carrol 发表于 2016-6-3 10:23
我用的是ERA-interim的资料，分层是一样的，necp的资料我不太懂所以没用过，要不你插值试试

1、就是我也不是特别明白应该选哪个数据，我是想模拟亮温和实测的亮温对比，我都看人家模拟向日葵8号上的AHI用的是ERA-interim的资料，咱们国家的模拟风云的用的是NCEP。你考虑的时候有针对性的选的吗？
2、你是不是说你替换数据的时候出现问题了，能说说你出现什么问题了吗？我替换数据的时候一直都不成功，想看看有什么共同性吗？
3、如果有些是28层，有些是14层，你插值了吗？用的什么方法？
4、数据里不是分为level_pressure 和layer_pressure ,前者是0：k后者是1：k，可是很多数据里level和layer的分层都是0:27，前者并没有比后者多一个数据，这种情况你出现过吗？怎么处理的？

SHERO桃桃

carrol 发表于 2016-6-2 14:07
1.我们是在服务器上运算的，环境变量加过路径。
2.直接进入到文件夹里，不用全部运行
3.进入到文件夹里 ...

能不能请问一下 link to library 这一步 环境变量是怎么加的

carrol

SHERO桃桃 发表于 2016-10-12 14:11
能不能请问一下 link to library 这一步 环境变量是怎么加的

我当时没有要改环境变量