Skip to content

Latest commit

 

History

History
42 lines (23 loc) · 1.57 KB

README.md

File metadata and controls

42 lines (23 loc) · 1.57 KB

fortran Harmonic Analysis 3D

This library contains Fortran scripts to do an harmonic analysis in 3D

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

  1. Open param.f90 and modify line 5,8,9,10,11,12 with your netcdf file (file_in), time (t_NAME), lon (x_NAME), lat (y_NAME), depth (z_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 :: z_NAME="depth"

character(len=*),parameter :: temp_NAME="variable"

also modify line 14 and 18 with your lon grid number (nx), lat grid number (ny), depth grid number (nz), 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

  1. Open jobcomp and modify NETCDFLIB and NETCDFINC path with your netcdf path

  2. Run ./jobcomp or in slurm: sbatch openmp_harmonic.sh