气象ts评分_给大家分享一个格点插值到站点然后TS评分的程序

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看，没有帐号？立即注册

x

!****************************************************************************

!

!  PROGRAM: T639格点资料插值成站点资料

!

!****************************************************************************

module pra_mod

implicit none

integer,parameter :: Nlon=90                               ! T639 E0-180  2

integer,parameter :: Nlat=45                               ! T639 N0-90   2

integer,parameter :: Nsta_max=5000                         ! 最多5000站

real,parameter :: yuezhi=0.1                               ! 有无降水的阈值  mm

character(len = *), parameter :: dir_T639='C:\Users\Tianle\Desktop\插值检验\t639-05\'        ! input original_data

character(len = *), parameter :: file1_t639_name='16112102.027'       ! input original_data 体现那一天

character(len = *), parameter :: file2_t639_name='16112202.027'       ! input original_data

character(len = *), parameter :: dir_station='C:\Users\Tianle\Desktop\插值检验\'             ! input original_data

character(len = *), parameter :: filename_station='STATIONS.DAT'      ! input original_data

character(len = *), parameter :: dir_obs='C:\Users\Tianle\Desktop\插值检验\sur-05-24\'       ! input original_data

character(len = *), parameter :: file1_obs_name='16112205.000'        ! input original_data 文件日期

character(len = *), parameter :: file2_obs_name='16112305.000'        ! input original_data

endmodule pra_mod

module data_mod

use pra_mod

implicit none

real lon_t639(0:Nlon)

real lat_t639(0:Nlat)

real Rain_t639(0:Nlon,0:N