Bulk surface treament for era5

[nataraj2]

This PR implements a bulk surface treatment for the surface fluxes for momentum, sensible, and latent heat fluxes. Smagorinsky 3d LES, Kessler microphysics are used. No PBL scheme. This is intended only as a testing routine, and not for long term use.

1. The animations show the velocity magnitude (km/hr) with and without the bulk surface treatment.
2. The images show the hurricane track and the intensification ie. max speed vs time (in knots), with and without the bulf surface treatment.

Observations:

• The sensible and latent heat flux source terms in the bulk surface treatment in the ocean are critical for intensification.
• Even with increased values of drag coefficient over land, the decay on land was not captured. The decay on land is caused by dissipation mechanisms. The PBL scheme that enable vertical mixing and conveys the effect of surface drag vertically, is missing.