diff --git a/Inst_eddies/input_ssh.f90 b/Inst_eddies/input_ssh.f90 index bf4a75fc875cf9f0962e588155b71a01ad2cad75..a98d8950791af3f971b1b896d8a1efaf497dfeed 100644 --- a/Inst_eddies/input_ssh.f90 +++ b/Inst_eddies/input_ssh.f90 @@ -48,7 +48,7 @@ contains real:: step_deg(2) = huge(0.) ! longitude and latitude steps, in degrees real lon_max, lat_max ! longitude and latitude, in degrees integer copy - namelist /grid_nml/ corner_deg, step_deg, nlon, nlat + namelist /grid_nml/ corner_deg, step_deg, nlon, nlat, periodic !---------------------------------------------------------------------- diff --git a/Overlap/read_grid.F90 b/Overlap/read_grid.F90 index 6b06c7a6c0282ede1d88f53e086a3ef9027a90f1..d3c454290c25028f3e7d77badddc3cc1da6a156e 100644 --- a/Overlap/read_grid.F90 +++ b/Overlap/read_grid.F90 @@ -6,7 +6,7 @@ module read_grid_m ! size of ssh array in input NetCDF file, assuming no repeated point ! if the grid is global - logical, protected:: periodic ! grid is periodic in longitude + logical, protected:: periodic = .false. ! grid is periodic in longitude real, protected:: corner(2) ! longitude and latitude of the corner of the whole grid in input @@ -38,7 +38,7 @@ contains real:: step_deg(2) = [huge(0.), huge(0.)] ! longitude and latitude steps, in degrees - namelist /grid_nml/ corner_deg, step_deg, nlon, nlat + namelist /grid_nml/ corner_deg, step_deg, nlon, nlat, periodic !----------------------------------------------------------------- @@ -48,12 +48,6 @@ contains action = "read", position = "rewind") read(unit, nml = grid_nml) close(unit) - - ! As we are requiring the grid spacing to be uniform, the value of - ! "periodic" may be deduced from the values of step_deg(1) and nlon: - periodic = nint(360. / step_deg(1)) == nlon - print *, "periodic = ", periodic - corner = corner_deg * deg_to_rad step = step_deg * deg_to_rad end if