def run():
print("Started model2roms")
confM2R = configM2R.Model2romsConfig()
if confM2R.createatmosforcing or confM2R.createoceanforcing:
if confM2R.createoceanforcing:
model2roms.convertMODEL2ROMS(confM2R)
clim2bry.writebry(confM2R)
# if confM2R.createAtmosForcing:
# atmosForcing.createAtmosFileUV(confM2R)
if confM2R.decimategridfile:
decimateGrid.createGrid(
confM2R.grdROMS,
"/Users/trondkr/Projects/KINO/GRID/kino_1600m_18072015.nc",
"/Users/trondkr/Projects/KINO/GRID/kino_1600m_18072015v2.nc", 2)
if confM2R.extractstations:
print(
"Running in station mode and extracting pre-defined station locations"
)
IOstation.getStationData(confM2R)
print('Finished ' + time.ctime(time.time()))
def run():
confM2R = configM2R.Model2romsConfig()
if confM2R.createatmosforcing or confM2R.createoceanforcing:
if confM2R.createoceanforcing:
model2roms.convertMODEL2ROMS(confM2R)
clim2bry.writebry(confM2R)
# if confM2R.createAtmosForcing:
# atmosForcing.createAtmosFileUV(confM2R)
print('Finished ' + time.ctime(time.time()))
def run():
print("Started model2roms")
confM2R = configM2R.Model2romsConfig()
if confM2R.createatmosforcing or confM2R.createoceanforcing:
if confM2R.createoceanforcing:
model2roms.convertMODEL2ROMS(confM2R)
clim2bry.writebry(confM2R)
# if confM2R.createAtmosForcing:
# atmosForcing.createAtmosFileUV(confM2R)
if confM2R.decimategridfile:
decimateGrid.createGrid(confM2R.grdROMS, "/Users/trondkr/Projects/KINO/GRID/kino_1600m_18072015.nc",
"/Users/trondkr/Projects/KINO/GRID/kino_1600m_18072015v2.nc", 2)
if confM2R.extractstations:
print("Running in station mode and extracting pre-defined station locations")
IOstation.getStationData(confM2R)
print('Finished ' + time.ctime(time.time()))
def main():
print '\n--------------------------\n'
print 'Started ' + time.ctime(time.time())
# EDIT ===================================================================
# Set show_progress to "False" if you do not want to see the progress
# indicator for horizontal interpolation.
show_progress = True
# Set compileAll to True if you want automatic re-compilation of all the
# fortran files necessary to run model2roms. Options are "gfortran" or "ifort". Edit
# compile.py to add other Fortran compilers.
compileAll = False
if compileAll is True:
import compile
compile.compileAll("ifort")
# Extract time-series of data for given longitude/latitude
extractStations = False
# Define a set of longitude/latitude positions with names to extract into
# station files (using extractStations)
if (extractStations):
stationNames = [
'NorthSea', 'Iceland', 'EastandWestGreenland', 'Lofoten',
'Georges Bank'
]
lonlist = [2.4301, -22.6001, -47.0801, 13.3801, -67.2001]
latlist = [54.5601, 63.7010, 60.4201, 67.5001, 41.6423]
# Create the bry, init, and clim files for a given grid and input data
createOceanForcing = True
# Create atmospheric forcing for the given grid
createAtmosForcing = False
# Create a smaller resolution grid based on your original. Decimates every second for
# each time run
decimateGridfile = False
# Write ice values to file (for Arctic regions)
writeIce = True
# Use ESMF for the interpolation. This requires that you have ESMF and ESMPy installed (import ESMF)
useESMF = True
# Apply filter to smooth the 2D fields after interpolation (time consuming)
useFilter = True
# Format to write the ouput to: 'NETCDF4', 'NETCDF4_CLASSIC', 'NETCDF3_64BIT', or 'NETCDF3_CLASSIC'
# Using NETCDF4 automatically turns on compression of files (ZLIB)
myformat = 'NETCDF4'
# Frequency of the input data: usually monthly
timeFrequencyOfInputData = "month" #, "month", "hour"
# Subset input data. If you have global data you may want to seubset these to speed up reading. Make
# sure that your input data are cartesian (0-360 or -180:180, -90:90)
subsetIndata = False
if subsetIndata:
subset = defineSubsetForIndata()
# Set the input data MODEL indatatype
indatatype = 'SODA'
indatatype = 'SODAMONTHLY'
indatatype = 'WOAMONTHLY'
indatatype = 'NORESM'
indatatype = 'GLORYS'
#indatatype = 'NS8KM'
#indatatype = 'NS8KMZ'
# GRIDTYPES ------------------------------------------------------------------------------
# Define what grid type you wnat to interpolate to
outgrid = "NS8KM"
#outgrid = "REGSCEN"
#outgrid = "KINO"
# outgrid= "AA10KM"
# Define what grid type you wnat to interpolate from: Can be Z for SIGMA for ROMS
# vertical coordinate system or ZLEVEL
ingridtype = "SIGMA"
ingridtype = "ZLEVEL"
outgridtype = "ROMS"
# PATH TO DATA ----------------------------------------------------------------------------
# Define the path to the input data
if indatatype == 'SODA':
modelpath = "/Volumes/MacintoshHD2/Datasets/SODA/"
if indatatype == 'SODAMONTHLY':
modelpath = "/Volumes/MacintoshHD2/Datasets/SODAMonthly/"
if indatatype == 'GLORYS':
modelpath = "/Volumes/MacintoshHD2/Datasets/GLOBAL_REANALYSIS_PHYS_001_009/"
modelpath = "/Users/trondkr/Projects/is4dvar/GLORYS2V3/"
modelpath = "/work/shared/imr/NS8KM/FORCING/GLORYS2V3/ftp.myocean.mercator-ocean.fr/Core/GLOBAL_REANALYSIS_PHYS_001_009/"
if indatatype == 'NORESM':
modelpath = "/Users/trondkr/Projects/RegScen/NRCP45AERCN_f19_g16_CLE_01/"
#modelpath = "/work/users/trondk/REGSCEN/NRCP45AERCN_f19_g16_CLE_01/"
#modelpath = "/work/users/trondk/REGSCEN/N20TRAERCN/"
if createAtmosForcing:
atmospath = "/Users/trondkr/Projects/RegScen/model2roms/TESTFILES/"
if indatatype == 'NS8KM':
modelpath = "/Users/trondkr/Projects/is4dvar/grid2lonlat/RESULTS/"
if indatatype == 'NS8KMZ':
modelpath = "/Users/trondkr/Projects/NOWMAPS/delivery-COPERNICUS/"
#modelpath = "/work/shared/imr/NS8KM/Z-LEVEL-ASSIMILATION2010-2013/"
if indatatype == 'WOAMONTHLY':
modelpath = "/Users/trondkr/Projects/is4dvar/createSSS/"
# PATH TO GRID -----------------------------------------------------------------------------
# Define the path to the grid file
if outgrid == "NS8KM":
romsgridpath = "/Users/trondkr/Projects/is4dvar/Grid/nordsjoen_8km_smoothed02022015.nc"
romsgridpath = "/work/shared/imr/NS8KM/FORCING/GRID/nordsjoen_8km_grid_hmax20m_v4.nc"
if outgrid == "KINO":
romsgridpath = "/work/users/trondk/KINO/GRID/kino_1600m_19022016_vf20.nc"
#romsgridpath = "/Users/trondkr/Projects/KINO/GRID/kino_1600m_07082015_vf20.nc"
if outgrid == "REGSCEN":
romsgridpath = "/Users/trondkr/Projects/RegScen/Grid/AA_10km_grid_noest.nc"
#romsgridpath = "/Users/trondkr/Projects/is4dvar/Grid/nordsjoen_8km_grid_hmax20m_v3.nc"
#romsgridpath = "/work/users/trondk/REGSCEN/GRID/AA_10km_grid_noest.nc"
if outgrid == "GREENLAND":
romsgridpath = "/Users/trondkr/Projects/RegScen/Grid/Sermilik_grid_4000m.nc"
romsgridpath = "/Users/trondkr/Projects/RegScen/model2roms/polarlowr_grid.nc"
if outgrid == "AA10KM":
romsgridpath = "/work-common/shared/imr/AA10KM/GRID/AA_10km_grid.nc"
if indatatype == 'WOAMONTHLY': isClimatology = True
else: isClimatology = False
# DETAILS -----------------------------------------------------------------------------------
# Define the period to create forcing for
start_year = 2013
end_year = 2014
start_month = 11
end_month = 12
start_day = 15
end_day = 15
if (int(calendar.monthrange(start_year, start_month)[1]) < start_day):
start_day = int(calendar.monthrange(start_year, start_month)[1])
if (int(calendar.monthrange(end_year, end_month)[1]) < end_day):
end_day = int(calendar.monthrange(end_year, end_month)[1])
startdate = datetime(start_year, start_month, start_day)
enddate = datetime(end_year, end_month, end_day)
years = [start_year + year for year in xrange(end_year + 1 - start_year)]
# Define what and name of variables to include in the forcing files
# -> myvars is the name model2roms uses to identify variables
# -> varNames is the name of the variable found in the NetCDF input files
if indatatype == 'SODA':
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
varNames = ['TEMP', 'SALT', 'SSH', 'U', 'V']
if indatatype == 'SODAMONTHLY':
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
varNames = ['temp', 'salt', 'ssh', 'u', 'v']
if indatatype == 'NS8KM':
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
varNames = ['votemper', 'vosaline', 'zeta', 'vozocrtx', 'vomecrty']
if indatatype == 'NS8KMZ':
fileNameIn, readFromOneFile = model2roms.getNS8KMZfilename(
startdate.year, startdate.month, startdate.day, "S", modelpath)
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
if (readFromOneFile):
varNames = ['temp', 'salt', 'zeta', 'u_eastward', 'v_northward']
else:
varNames = ['votemper', 'vosaline', 'zeta', 'vozocrtx', 'vomecrty']
if indatatype == 'GLORYS':
if (writeIce):
myvars = [
'temperature', 'salinity', 'ssh', 'uvel', 'vvel', 'uice',
'vice', 'aice', 'hice'
]
varNames = [
'votemper', 'vosaline', 'sossheig', 'vozocrtx', 'vomecrty',
'iicevelu', 'iicevelv', 'ileadfra', 'iicethic'
]
else:
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
varNames = [
'votemper', 'vosaline', 'sossheig', 'vozocrtx', 'vomecrty'
]
if indatatype == 'WOAMONTHLY':
myvars = ['temperature', 'salinity']
varNames = ['t_an', 's_an']
if indatatype == 'NORESM':
myvars = [
'temperature', 'salinity', 'ssh', 'uvel', 'vvel', 'ageice', 'uice',
'vice', 'aice', 'hice', 'snow_thick'
]
varNames = [
'templvl', 'salnlvl', 'sealv', 'uvellvl', 'vvellvl', 'iage',
'uvel', 'vvel', 'aice', 'hi', 'hs'
]
# NO EDIT BELOW ====================================================================================================
abbreviation = defineAbbreviation(outgrid)
climName, initName, bryName = defineOutputFilenames(
abbreviation, start_year, end_year, start_month, end_month, start_day,
end_day, indatatype)
if isClimatology is True:
climName = abbreviation + '_' + str(indatatype) + '_climatology.nc'
# Create the grid object for the output grid
grdROMS = grd.grdClass(romsgridpath, "ROMS", outgridtype, useESMF, 'ocean',
outgrid)
grdROMS.vars = myvars
if (useESMF):
import ESMF
manager = ESMF.Manager(logkind=ESMF.LogKind.MULTI, debug=True)
if createOceanForcing:
showInfo(myvars, romsgridpath, climName, initName, bryName, start_year,
end_year, isClimatology, useESMF, myformat)
model2roms.convertMODEL2ROMS(years,
startdate,
enddate,
timeFrequencyOfInputData,
climName,
initName,
modelpath,
romsgridpath,
myvars,
varNames,
show_progress,
indatatype,
outgridtype,
isClimatology,
writeIce,
useESMF,
useFilter,
myformat,
subsetIndata,
outgrid,
subset=None)
clim2bry.writeBry(grdROMS, start_year, bryName, climName, writeIce,
indatatype, myformat)
if createAtmosForcing:
atmosForcing.createAtmosFileUV(grdROMS, modelpath, atmospath,
startdate, enddate, useESMF, myformat,
abbreviation, indatatype, gridtype,
show_progress)
#if decimateGridfile:
#decimateGrid.createGrid(grdROMS, "/Users/trondkr/Projects/KINO/GRID/kino_1600m_18072015.nc", "/Users/trondkr/Projects/KINO/GRID/kino_1600m_18072015v2.nc", 2)
if extractStations:
print "Running in station mode and extracting pre-defined station locations"
IOstation.getStationData(years, IDS, modelpath, latlist, lonlist,
stationNames)
print 'Finished ' + time.ctime(time.time())
def main():
print '\n--------------------------\n'
print 'Started ' + time.ctime(time.time())
# EDIT ===================================================================
# Set show_progress to "False" if you do not want to see the progress
# indicator for horizontal interpolation.
show_progress = True
# Set compileAll to True if you want automatic re-compilation of all the
# fortran files necessary to run model2roms. Options are "gfortran" or "ifort". Edit
# compile.py to add other Fortran compilers.
compileAll = False
if compileAll is True:
import compile; compile.compileAll("ifort")
# Extract time-series of data for given longitude/latitude
extractStations = False
# Define a set of longitude/latitude positions with names to extract into
# station files (using extractStations)
if (extractStations):
stationNames = ['NorthSea', 'Iceland', 'EastandWestGreenland', 'Lofoten', 'Georges Bank']
lonlist = [2.4301, -22.6001, -47.0801, 13.3801, -67.2001]
latlist = [54.5601, 63.7010, 60.4201, 67.5001, 41.6423]
# Create the bry, init, and clim files for a given grid and input data
createOceanForcing = True
# Create atmospheric forcing for the given grid
createAtmosForcing = False
# Create a smaller resolution grid based on your original. Decimates every second for
# each time run
decimateGridfile = False
# Write ice values to file (for Arctic regions)
writeIce = False
# Use ESMF for the interpolation. This requires that you have ESMF and ESMPy installed (import ESMF)
useESMF = True
# Apply filter to smooth the 2D fields after interpolation (time consuming)
useFilter = True
# Format to write the ouput to: 'NETCDF4', 'NETCDF4_CLASSIC', 'NETCDF3_64BIT', or 'NETCDF3_CLASSIC'
# Using NETCDF4 automatically turns on compression of files (ZLIB)
myformat='NETCDF4'
# Frequency of the input data: usually monthly
timeFrequencyOfInputData = "day" #, "month", "hour"
# Subset input data. If you have global data you may want to seubset these to speed up reading. Make
# sure that your input data are cartesian (0-360 or -180:180, -90:90)
subsetIndata = False
if subsetIndata:
subset = defineSubsetForIndata()
# Set the input data MODEL indatatype
indatatype = 'SODA'
indatatype = 'SODAMONTHLY'
indatatype = 'WOAMONTHLY'
indatatype = 'NORESM'
indatatype = 'GLORYS'
#indatatype = 'NS8KM'
#indatatype = 'NS8KMZ'
# GRIDTYPES ------------------------------------------------------------------------------
# Define what grid type you wnat to interpolate to
outgrid = "NS8KM"
#outgrid = "REGSCEN"
#outgrid = "KINO"
# Define what grid type you wnat to interpolate from: Can be Z for SIGMA for ROMS
# vertical coordinate system or ZLEVEL
ingridtype = "SIGMA"
ingridtype = "ZLEVEL"
outgridtype="ROMS"
# PATH TO DATA ----------------------------------------------------------------------------
# Define the path to the input data
if indatatype == 'SODA':
modelpath = "/Volumes/MacintoshHD2/Datasets/SODA/"
if indatatype == 'SODAMONTHLY':
modelpath = "/Volumes/MacintoshHD2/Datasets/SODAMonthly/"
if indatatype == 'GLORYS':
modelpath = "/Volumes/MacintoshHD2/Datasets/GLOBAL_REANALYSIS_PHYS_001_009/"
modelpath = "/Users/trondkr/Projects/is4dvar/GLORYS2V3/"
modelpath = "/work/shared/imr/NS8KM/FORCING/GLORYS2V3/ftp.myocean.mercator-ocean.fr/Core/GLOBAL_REANALYSIS_PHYS_001_009/"
if indatatype == 'NORESM':
modelpath = "/Users/trondkr/Projects/RegScen/NRCP45AERCN_f19_g16_CLE_01/"
#modelpath = "/work/users/trondk/REGSCEN/NRCP45AERCN_f19_g16_CLE_01/"
#modelpath = "/work/users/trondk/REGSCEN/N20TRAERCN/"
if createAtmosForcing:
atmospath = "/Users/trondkr/Projects/RegScen/model2roms/TESTFILES/"
if indatatype == 'NS8KM':
modelpath = "/Users/trondkr/Projects/is4dvar/grid2lonlat/RESULTS/"
if indatatype == 'NS8KMZ':
modelpath = "/Users/trondkr/Dropbox/deliveryFOFINAL/"
modelpath = "/work/shared/imr/NS8KM/deliveryFO/FINAL/"
if indatatype == 'WOAMONTHLY':
modelpath = "/Users/trondkr/Projects/is4dvar/createSSS/"
# PATH TO GRID -----------------------------------------------------------------------------
# Define the path to the grid file
if outgrid == "NS8KM":
romsgridpath = "/Users/trondkr/Projects/is4dvar/Grid/nordsjoen_8km_smoothed02022015.nc"
romsgridpath = "/work/users/trondk/NS8km/FORCING/GRID/nordsjoen_8km_grid_hmax20m_v3.nc"
if outgrid == "KINO":
romsgridpath = "/work/users/trondk/KINO/GRID/kino_1600m_07082015_vf20.nc"
#romsgridpath = "/Users/trondkr/Projects/KINO/GRID/kino_1600m_07082015_vf20.nc"
if outgrid == "REGSCEN":
romsgridpath = "/Users/trondkr/Projects/RegScen/Grid/AA_10km_grid_noest.nc"
#romsgridpath = "/Users/trondkr/Projects/is4dvar/Grid/nordsjoen_8km_grid_hmax20m_v3.nc"
#romsgridpath = "/work/users/trondk/REGSCEN/GRID/AA_10km_grid_noest.nc"
if outgrid == "GREENLAND":
romsgridpath="/Users/trondkr/Projects/RegScen/Grid/Sermilik_grid_4000m.nc"
romsgridpath="/Users/trondkr/Projects/RegScen/model2roms/polarlowr_grid.nc"
if indatatype == 'WOAMONTHLY': isClimatology = True
else: isClimatology = False
# DETAILS -----------------------------------------------------------------------------------
# Define the period to create forcing for
start_year = 2012
end_year = 2013
start_month = 11
end_month = 12
start_day = 15
end_day = 15
if (int(calendar.monthrange(start_year, start_month)[1]) < start_day):
start_day = int(calendar.monthrange(start_year, start_month)[1])
if (int(calendar.monthrange(end_year, end_month)[1]) < end_day):
end_day = int(calendar.monthrange(end_year, end_month)[1])
startdate = datetime(start_year, start_month, start_day)
enddate = datetime(end_year, end_month, end_day)
years = [start_year+year for year in xrange(end_year+1-start_year)]
# Define what and name of variables to include in the forcing files
# -> myvars is the name model2roms uses to identify variables
# -> varNames is the name of the variable found in the NetCDF input files
if indatatype == 'SODA':
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
varNames = ['TEMP', 'SALT', 'SSH', 'U', 'V']
if indatatype == 'SODAMONTHLY':
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
varNames = ['temp', 'salt', 'ssh', 'u', 'v']
if indatatype == 'NS8KM':
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
varNames = ['votemper', 'vosaline', 'zeta', 'vozocrtx', 'vomecrty']
if indatatype == 'NS8KMZ':
fileNameIn, readFromOneFile = model2roms.getNS8KMZfilename(startdate.year, startdate.month, startdate.day, "S", modelpath)
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
if (readFromOneFile):
varNames = ['temp', 'salt', 'zeta', 'u_eastward', 'v_northward']
else:
varNames = ['votemper', 'vosaline', 'zeta', 'vozocrtx', 'vomecrty']
if indatatype == 'GLORYS':
if (writeIce):
myvars = ['temperature','salinity', 'ssh', 'uvel', 'vvel','uice','vice','aice','hice']
varNames = ['votemper', 'vosaline', 'sossheig', 'vozocrtx', 'vomecrty','iicevelu', 'iicevelv', 'ileadfra', 'iicethic']
else:
myvars = ['temperature', 'salinity', 'ssh', 'uvel', 'vvel']
varNames = ['votemper', 'vosaline', 'sossheig', 'vozocrtx', 'vomecrty']
if indatatype == 'WOAMONTHLY':
myvars = ['temperature','salinity']
varNames = ['t_an', 's_an']
if indatatype == 'NORESM':
myvars = ['temperature','salinity', 'ssh', 'uvel', 'vvel','ageice','uice','vice','aice','hice','snow_thick']
varNames = ['templvl','salnlvl','sealv', 'uvellvl', 'vvellvl','iage', 'uvel', 'vvel', 'aice', 'hi', 'hs']
# NO EDIT BELOW ====================================================================================================
abbreviation = defineAbbreviation(outgrid)
climName,initName,bryName = defineOutputFilenames(abbreviation,start_year,end_year,start_month,end_month,start_day,end_day,indatatype)
if isClimatology is True:
climName=abbreviation + '_' + str(indatatype) + '_climatology.nc'
# Create the grid object for the output grid
grdROMS = grd.grdClass(romsgridpath, "ROMS", outgridtype, useESMF,'ocean', outgrid)
grdROMS.vars=myvars
if (useESMF):
import ESMF
manager = ESMF.Manager(logkind = ESMF.LogKind.MULTI, debug = True)
if createOceanForcing:
showInfo(myvars, romsgridpath, climName, initName, bryName, start_year, end_year, isClimatology, useESMF, myformat)
model2roms.convertMODEL2ROMS(years, startdate, enddate, timeFrequencyOfInputData, climName, initName, modelpath, romsgridpath, myvars, varNames, show_progress,
indatatype, outgridtype, isClimatology, writeIce, useESMF, useFilter, myformat, subsetIndata, outgrid, subset=None)
clim2bry.writeBry(grdROMS, start_year, bryName, climName, writeIce, indatatype, myformat)
if createAtmosForcing:
atmosForcing.createAtmosFileUV(grdROMS,modelpath,atmospath,startdate,enddate,useESMF,
myformat,abbreviation,indatatype,gridtype,show_progress)
#if decimateGridfile:
#decimateGrid.createGrid(grdROMS, "/Users/trondkr/Projects/KINO/GRID/kino_1600m_18072015.nc", "/Users/trondkr/Projects/KINO/GRID/kino_1600m_18072015v2.nc", 2)
if extractStations:
print "Running in station mode and extracting pre-defined station locations"
IOstation.getStationData(years, IDS, modelpath, latlist, lonlist, stationNames)
print 'Finished ' + time.ctime(time.time())