Skip to content
Snippets Groups Projects
Commit a4a01f50 authored by GUEZ Lionel's avatar GUEZ Lionel
Browse files

Read `input_ssh_nml` before

Read `input_ssh_nml` before reading the NetCDF files. Preparing to
read from the namelist whether the grid is uniform.
parent 086d52f5
No related branches found
No related tags found
No related merge requests found
......@@ -57,6 +57,9 @@ contains
!----------------------------------------------------------------------
periodic = .false. ! default value
print *, "input_ssh: Enter input_ssh_nml:"
read(unit = *, nml = input_ssh_nml)
call nf95_open(ssh_fname, nf95_nowrite, ncid)
call nf95_find_coord(ncid, dimid = dimid, varid = varid, &
std_name = "longitude")
......@@ -77,9 +80,6 @@ contains
write(unit, nml = grid_nml)
close(unit)
step = step_deg * deg_to_rad
periodic = .false. ! default value
print *, "input_ssh: Enter input_ssh_nml:"
read(unit = *, nml = input_ssh_nml)
if (periodic) call assert(4 * max_radius(1) < nlon, "input_ssh: max_radius")
copy = merge(max_radius(1), 0, periodic)
allocate(ssh(1 - copy:nlon + copy, nlat), u(1 - copy:nlon + copy, nlat), &
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment