DMP_TOOL modification

[jenniewaters]

Context

The current DMP_TOOL for calculating the restoration coefficients for use with the tra_dmp module in NEMO calculates the distance to coast (using the surface only) for every vertical level. Subsequently the tool is very inefficient for large model domains.

Proposal

• Include functionality to allow the distance to coast to be read in rather than calculated
• Only perform the distance to coast calculation once.

[jenniewaters]

My proposed code changes can be seen here:
​http://fcm3/projects/NEMO.xm/browser/branches/UKMO/dev_r5518_DMP_TOOLS

[jenniewaters]

In 10778:

Error: Failed to load processor CommitTicketReference

No macro or processor named 'CommitTicketReference' found

[jenniewaters]

Replying to jenniewaters:

Context

The current DMP_TOOL for calculating the restoration coefficients for use with the tra_dmp module in NEMO calculates the distance to coast (using the surface only) for every vertical level. Subsequently the tool is very inefficient for large model domains.

Proposal

• Include functionality to allow the distance to coast to be read in rather than calculated
• Only perform the distance to coast calculation once.

Future developments

Get the DMP_TOOL to calculate the distance to coast correctly for each vertical level. This could be achieved by reading in the land mask for the current vertical level in cofdis in coast_dist.F90.

[jenniewaters]

I ran both the trunk@5518 and my version of the DMP_TOOLS with namelist_025 (attached to the ticket). The resulting resto.nc files are identical, but the trunk version took 2hours 45 minutes to run, while the new version took less the 3 minutes.

[jenniewaters]

The namelist_025_readdistcoast is used to generate the new resto.nc file used in the PS43 testing. The file differs slightly from the versions produced using namelist_025 as the distance to coast is depth dependent. It took less than 30 seconds to produce this file.
resto_plots.png shows the level 53 resto field generated using namelist_025_readdistcoast (right plot) and the level 53 difference between resto fields generated using namelist_025_readdistcoast and namelist_025. The only differences are around sea mounts, islands and coastlines, as we would expect.