Hi ROMS users
I'm running a simulation wherein I provide the surface forcing i.e. stresses, E-P and net heat flux directly (i.e. without using Bulk Formulae). I am providing them after interpolating from data sets to ROMS application grid. I was wondering if it is possible to provide these forcing parameters on a lat-lon grid instead of providing them on ROMS grid in direct application and let the model do interpolation on ROMS grid. I know it is possible to do so if you use BULK_FLUX. If it is possible, do the lat-lon grid need to have specific structure? Where can I find the ctl file for such a file? Thanks in advance.
interpolation to roms grid
-
- Posts: 27
- Joined: Mon Jan 27, 2014 9:50 pm
- Location: Indian Institute of Science
Re: interpolation to roms grid
Yes, you can do what you are asking. The interpolation code is the same, whether for winds or wind stress. Here's a typical "ncdump -h" of a forcing file that ROMS can interpolate from:
Code: Select all
netcdf MERRA_swrad_3hours_2016 {
dimensions:
lon = 576 ;
lat = 361 ;
srf_time = UNLIMITED ; // (2192 currently)
variables:
double lon(lon) ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
double lat(lat) ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
double srf_time(srf_time) ;
srf_time:units = "days since 1900-01-01 00:00:00" ;
srf_time:calendar = "LEAP" ;
float swrad(srf_time, lat, lon) ;
swrad:_FillValue = 1.e+15f ;
swrad:missing_value = 1.e+15f ;
swrad:long_name = "surface_incoming_shortwave_flux" ;
swrad:units = "W m-2" ;
swrad:coordinates = "lon lat" ;
// global attributes:
:Author = "<module>" ;
:Created = "2016-11-14 15:15:02" ;
:title = "MERRA-2 dataset. Modern Era Retrospective-analysis" ;
:history = "Mon Nov 14 16:31:47 2016: ncrcat MERRA_swrad_3hours_2016_01.nc MERRA_swrad_3hours_2016_02.nc MERRA_swrad_3hours_2016_03.nc MERRA_swrad_3hours_2016_04.nc MERRA_swrad_3hours_2016_05.nc MERRA_swrad_3hours_2016_06.nc MERRA_swrad_3hours_2016_07.nc MERRA_swrad_3hours_2016_08.nc MERRA_swrad_3hours_2016_09.nc MERRA_swrad_3hours_2016.nc" ;
:nco_openmp_thread_number = 1 ;
}