Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Sat Jul 20, 2019 3:59 pm




Post new topic Reply to topic  [ 1 post ] 

All times are UTC

Author Message
 Post subject: Forcing file not read
PostPosted: Fri Mar 26, 2010 12:30 pm 
Offline

Joined: Fri Mar 23, 2007 12:35 pm
Posts: 5
Location: Laboratori d'Engineyria Maritima/Universitat Polit├Ęcnica de Catalunya (LIM/UPC)
Dear/s:

I am runnig an implementation with an unique forcing file, as it is shown below:

netcdf roms_frc {
dimensions:
xi_u = 176 ;
eta_u = 130 ;
xi_v = 177 ;
eta_v = 129 ;
xi_rho = 177 ;
eta_rho = 130 ;
xi_psi = 176 ;
eta_psi = 129 ;
sms_time = 12 ;
shf_time = 12 ;
swf_time = 12 ;
sst_time = 12 ;
srf_time = 12 ;
sss_time = 12 ;
tide_period = 10 ;
variables:
double sms_time(sms_time) ;
sms_time:long_name = "surface momentum stress time" ;
sms_time:units = "days" ;
sms_time:cycle_length = 360. ;
double shf_time(shf_time) ;
shf_time:long_name = "surface heat flux time" ;
shf_time:units = "days" ;
shf_time:cycle_length = 360. ;
double swf_time(swf_time) ;
swf_time:long_name = "surface freshwater flux time" ;
swf_time:units = "days" ;
swf_time:cycle_length = 360. ;
double sst_time(sst_time) ;
sst_time:long_name = "sea surface temperature time" ;
sst_time:units = "days" ;
sst_time:cycle_length = 360. ;
double sss_time(sss_time) ;
sss_time:long_name = "sea surface salinity time" ;
sss_time:units = "days" ;
sss_time:cycle_length = 360. ;
double srf_time(srf_time) ;
srf_time:long_name = "solar shortwave radiation time" ;
srf_time:units = "days" ;
srf_time:cycle_length = 360. ;
double sustr(sms_time, eta_u, xi_u) ;
sustr:long_name = "surface u-momentum stress" ;
sustr:units = "Newton meter-2" ;
double svstr(sms_time, eta_v, xi_v) ;
svstr:long_name = "surface v-momentum stress" ;
svstr:units = "Newton meter-2" ;
double shflux(shf_time, eta_rho, xi_rho) ;
shflux:long_name = "surface net heat flux" ;
shflux:units = "Watts meter-2" ;
double swflux(swf_time, eta_rho, xi_rho) ;
swflux:long_name = "surface freshwater flux (E-P)" ;
swflux:units = "centimeter day-1" ;
swflux:positive = "net evaporation" ;
swflux:negative = "net precipitation" ;
double SST(sst_time, eta_rho, xi_rho) ;
SST:long_name = "sea surface temperature" ;
SST:units = "Celsius" ;
double SSS(sss_time, eta_rho, xi_rho) ;
SSS:long_name = "sea surface salinity" ;
SSS:units = "PSU" ;
double dQdSST(sst_time, eta_rho, xi_rho) ;
dQdSST:long_name = "surface net heat flux sensitivity to SST" ;
dQdSST:units = "Watts meter-2 Celsius-1" ;
double swrad(srf_time, eta_rho, xi_rho) ;
swrad:long_name = "solar shortwave radiation" ;
swrad:units = "Watts meter-2" ;
swrad:positive = "downward flux, heating" ;
swrad:negative = "upward flux, cooling" ;
double tide_period(tide_period) ;
tide_period:long_name = "Tide angular period" ;
tide_period:units = "Hours" ;
double tide_Ephase(tide_period, eta_rho, xi_rho) ;
tide_Ephase:long_name = "Tidal elevation phase angle" ;
tide_Ephase:units = "Degrees" ;
double tide_Eamp(tide_period, eta_rho, xi_rho) ;
tide_Eamp:long_name = "Tidal elevation amplitude" ;
tide_Eamp:units = "Meter" ;
double tide_Cmin(tide_period, eta_rho, xi_rho) ;
tide_Cmin:long_name = "Tidal current ellipse semi-minor axis" ;
tide_Cmin:units = "Meter second-1" ;
double tide_Cmax(tide_period, eta_rho, xi_rho) ;
tide_Cmax:long_name = "Tidal current, ellipse semi-major axis" ;
tide_Cmax:units = "Meter second-1" ;
double tide_Cangle(tide_period, eta_rho, xi_rho) ;
tide_Cangle:long_name = "Tidal current inclination angle" ;
tide_Cangle:units = "Degrees between semi-major axis and East" ;
double tide_Cphase(tide_period, eta_rho, xi_rho) ;
tide_Cphase:long_name = "Tidal current phase angle" ;
tide_Cphase:units = "Degrees" ;
(........).

The tidal variables are reading correctly, but the others ones it seems that an incorrect values are read (at least does not coincide with .nc file).

*****************************************
GET_2DFLD - surface u-momentum stress, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.54176873E-05 Max = 5.72925972E-05)
GET_2DFLD - surface v-momentum stress, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -9.24406932E-05 Max = -3.40081785E-05)
GET_2DFLD - solar shortwave radiation flux, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.63061913E-05 Max = 1.75408070E-05)
GET_2DFLD - surface net heat flux, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.88596504E-05 Max = -2.89604253E-05)
GET_2DFLD - surface net freswater flux, (E-P), t = 15 00:00:00
(Rec=0001, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.56206398E-08 Max = 2.83432340E-08)
*********************************************
For example, the maximum of sustr (from .nc file)) for index=1 is 0.0587 and the minimum is -0.0104, the somthing is wrong. :(

The options file are:

#define UV_ADV
#define UV_COR
#define UV_QDRAG

#define DJ_GRADPS
#undef TS_A4HADVECTION
#undef TS_A4VADVECTION

#define NONLIN_EOS
#define SALINITY
#define MASKING
#define CURVGRID

#define SOLVE3D
#define TS_DIF2
#define MIX_GEO_TS
#define DIFF_GRID

#define SPHERICAL
#define SPLINES
#define AVERAGES
#define AVERAGES_AKV
#define AVERAGES_AKT
#define AVERAGES_AKS

#define VISC_GRID
#define SPONGE

#define LMD_MIXING
#ifdef LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_SKPP
# define LMD_BKPP
# define LMD_NONLOCAL
#endif

/* pels bottoms bhflux,bwflux,..etc*/
#define ANA_BSFLUX
#define ANA_BTFLUX

#define ZCLIMATOLOGY
#define ZCLM_NUDGING
#define ANA_M2OBC
#define ANA_FSOBC

#define M2CLIMATOLOGY
#define M3CLIMATOLOGY
#define TCLIMATOLOGY
#define M2CLM_NUDGING
#define M3CLM_NUDGING
#define TCLM_NUDGING

#define NORTHERN_WALL
#define WESTHERN_WALL

#define EST_FSCHAPMAN
#define EST_M2FLATHER
#define EST_M3GRADIENT
#define EST_TGRADIENT

#define SOUTH_FSCHAPMAN
#define SOUTH_M2FLATHER
#define SOUTH_M3GRADIENT
#define SOUTH_TGRADIENT

#define SSH_TIDES
#define UV_TIDES

#ifdef SSH_TIDES
# define ADD_FSOBC
#endif

#ifdef UV_TIDES
# undef FSOBC_REDUCED
# define ADD_M2OBC
#endif
**********************************************************

Do you thing somthing is incorredted formulated?

Thanks

_________________
Manel Grifoll
manel.grifoll@upc.edu
Laboratori d'Enginyeria Maritima
Universitat Politecnica de Catalunya


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 4 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group