This library contains Fortran scripts to compute annual and semiannual harmonics using an harmonic fit
References
Wyrtki, K. (1965). The annual and semiannual variation of sea surface temperature in the North Pacific Ocean. Limnology and Oceanography, 10(3), 307-313.
Campos, F., Mosquera, K. (2020) Boletín científico El Niño, Instituto Geofísico del Perú, Vol. 7 Nº 04. pág. 12 -16 link
Getting Started
- Open param.f90 and modify line 5,8,9,10,11 with your netcdf file (file_in), time (t_NAME), lon (x_NAME), lat (y_NAME) and variable (temp_NAME):
character(len=*),parameter :: file_in="your_file.nc"
character(len=*),parameter :: t_NAME="time"
character(len=*),parameter :: y_NAME="latitude"
character(len=*),parameter :: x_NAME="longitude"
character(len=*),parameter :: temp_NAME="variable"
also modify line 13, 15 and 17 with your lon grid number (nx), lat grid number (ny), time steps (nt) missing_val , scale factor (sf_thetao) and add offset (af_thetao):
integer, parameter :: nx = 280, ny = 321, nt = 9647
integer i, j, k, ierr
real, parameter :: pi=3.1415927, missing_val=-32767, sf_thetao=1, af_thetao=0
-
Open jobcomp and modify NETCDFLIB and NETCDFINC path with your netcdf path
-
Run ./jobcomp or in slurm: sbatch openmp_harmonic.sh