HarmonicanalysisTideLoad.f90

!  HarmonicanalysisTideLoad.f90 
!
!  FUNCTIONS:
!  HarmonicanalysisTideLoad - Entry point of console application.
!
!****************************************************************************
       program HarmonicanalysisTideLoad
       implicit none
	 character*800::dtmfl,tidefl
	 integer knd
	 real*8::dk,itd
!---------------------------------------------------------------------
      !输入潮汐负荷类型knd
      !knd=0 surface atmosphere tide, =1 ocean tide
      knd=1!knd=0大气压潮，=1海洋潮汐
      !输入陆海地形球坐标格网文件名，用于分离陆海区域knd=0 or =1
      !The spatial resolution of the land-sea terrain grid should not be lower than that of the surface load grid.
      write(dtmfl,*)'sphetopo30m.dat'
      !输入系列分潮调和常数球坐标格网文件名文件,球坐标格网头文件第七个数为长整型ETideLoad系统格式时间
      !The seventh numerical value of spherical coordinate load ewh grid file header is the long integer time agreed by ETideLoad
      write(tidefl,*)'sphgridtidefls.txt'
      !输入迭代终止条件：残差标准差为原格值标准的dk倍
      !input the iteration condition parameters
      !Iteration termination condition: The standard deviation of the residual grid value is less than dk of the standard deviation
      !of the original grid value, or the difference of the residual standard deviation of the previous step iteration relative to 
      !the current step iteration is less than itd of the standard deviation of the original grid values.
      dk=1.d-3; itd=1.d-4
      write(*, *)"    Begin compulation......"
      write(*, *)"    name doodson  C00(+) relative error(+)%  C00(-) relative error(-)%"
      call tideharmanalysis(dtmfl,tidefl,knd,dk,itd)
      pause
      end