Forcing file not read

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
Manel_Grifoll
Posts: 5
Joined: Fri Mar 23, 2007 12:35 pm
Location: Laboratori d'Engineyria Maritima/Universitat Politècnica de Catalunya (LIM/UPC)

Forcing file not read

#1 Post by Manel_Grifoll » Fri Mar 26, 2010 12:30 pm

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

Post Reply