R语言计算地球固体潮汐数据

• 利用R语言基于Fortan ETERNA 3.4计算地球固体潮汐

本文章算是一个笔记，给出利用jkennel在github上发行的earthtide库实现对于地球固体潮汐的计算，步骤稍显粗略。

1. 首先安装Anaconda环境，详细步骤可以参考【Anaconda教程01】怎么安装Anaconda3
   
2. 安装Rstudio环境【Anaconda教程01】怎么安装Anaconda3
   
3. 配置Rstudio的包镜像

Tools–>Global options–>Package–>CRAN mirror
修改至合肥

5. 安装earthtide库
   

6. 运行程序

library(earthtide)

# One month of hourly data
tide <-  Earthtide$
  new(utc = as.POSIXct("2015-01-01", tz = "UTC") + 0:(24*31) * 3600,
      latitude = 52.3868,
      longitude = 9.7144,
      wave_groups = data.frame(start = 0.0, end = 6.0))$
  predict(method = "gravity", astro_update = 1)$    # compute gravity
  lod_tide()$                                       # LOD tide column
  pole_tide()$                                      # pole tide column
  tide()                                            # return result


plot(tide[1:2],type = 'l',lwd=2)