[MITgcm-support] strange hflux value in MON exf_hflux
Chun-Yan Zhou
c.zhou at dundee.ac.uk
Mon Jan 14 14:52:44 EST 2013
Hi, there,
I use hflux in data.exf, with the max and min 679 and -247,
however, the exf_hflux_max in the output is 1200, which is much more than the biggest number in the file I provide. Besides, I did not provide the sflux file, the exf_sflux_max still had value.
MON exf_hflux_max = 1.2212290649558E+03
(PID.TID 0000.0001) %MON exf_hflux_min = 2.7547387234863E+00
(PID.TID 0000.0001) %MON exf_hflux_mean = 5.1822075200978E+02
(PID.TID 0000.0001) %MON exf_hflux_sd = 3.0564424766070E+02
(PID.TID 0000.0001) %MON exf_hflux_del2 = 5.4699216688885E-02
(PID.TID 0000.0001) %MON exf_sflux_max = 3.2260977333231E-07
(PID.TID 0000.0001) %MON exf_sflux_min = 1.0412672397707E-09
(PID.TID 0000.0001) %MON exf_sflux_mean = 1.3920299907157E-07
(PID.TID 0000.0001) %MON exf_sflux_sd = 7.6904666535711E-08
(PID.TID 0000.0001) %MON exf_sflux_del2 = 1.5054136398581E-11
Maybe there is another factor contributing the exf_hflux_ma?
Any help is much appreciated.
data.exf
#
# *********************
# External Forcing Data
# *********************
&EXF_NML_01
repeatPeriod= 31104000,
exf_iprec = 64,
#exf_verbose = .TRUE.,
&
# *********************
&EXF_NML_02
uwindstartdate1=20020101,
uwindstartdate2=000000,
uwindperiod=2592000.0,
#
vwindstartdate1=20020101,
vwindstartdate2=000000,
vwindperiod=2592000.0,
#
#
hfluxstartdate1=20020101,
hfluxstartdate2=000000,
hfluxperiod=2592000.0,
#
climsststartdate1 = 20020101,
climsststartdate2 = 000000,
climsstperiod = 2592000.0,
climsssstartdate1 = 20020101,
climsssstartdate2 = 000000,
climsssperiod = 2592000.0,
uwindfile = 'windxgf.bin',
vwindfile = 'windygf.bin',
# 0.5 months restoring timescale for temperature
climsstTauRelax = 129600,
hfluxfile = 'qnet.bin',
climsstfile = 'sst.bin',
climsssTauRelax = 129600,
climsssfile = 'sss.bin',
#
&
# *********************
&EXF_NML_03
&
# *********************
&EXF_NML_04
#
&
# *********************
&EXF_NML_OBCS
obcsSstartdate1 = 20020101,
obcsSstartdate2 = 000000,
obcsSperiod = 2592000,
obcsEstartdate1 = 20020101,
obcsEstartdate2 = 000000,
obcsEperiod = 2592000,
&
EXF_OPTION.h
C Bulk formulae related flags.
#define ALLOW_ATM_TEMP
#define ALLOW_ATM_WIND
#define ALLOW_DOWNWARD_RADIATION
#define ALLOW_RUNOFF
#if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND))
# define ALLOW_BULKFORMULAE
# undef ALLOW_BULK_LARGEYEAGER04
#endif
C Zenith Angle/Albedo related flags.
#ifdef ALLOW_DOWNWARD_RADIATION
# undef ALLOW_ZENITHANGLE
#endif
C Use ocean_emissivity*lwdwon in lwFlux. This flag should be define
C unless to reproduce old results (obtained with inconsistent old code)
#ifdef ALLOW_DOWNWARD_RADIATION
# define EXF_LWDOWN_WITH_EMISSIVITY
#endif
C Relaxation to monthly climatologies.
#define ALLOW_CLIMSST_RELAXATION
#define ALLOW_CLIMSSS_RELAXATION
C Use spatial interpolation to interpolate
C forcing files from input grid to model grid.
#undef USE_EXF_INTERPOLATION
C for interpolated vector fields, rotate towards model-grid axis
C using old rotation formulae (instead of grid-angles)
#define ALLOW_ADDWAVE
C input wave paremeters
#define ALLOW_SEDSURF
C input sediment surf input gtracer
#undef EXF_USE_OLD_VEC_ROTATION
C for interpolation around N & S pole, use the old formulation
C (no pole symmetry, single vector-comp interp, reset to 0 zonal-comp @ N.pole)
#undef EXF_USE_OLD_INTERP_POLE
#define EXF_INTERP_USE_DYNALLOC
#if ( defined (EXF_INTERP_USE_DYNALLOC) && defined (USING_THREADS) )
# define EXF_IREAD_USE_GLOBAL_POINTER
#endif
#endif /* ndef ALLOW_AUTODIFF_TAMC */
#endif /* ALLOW_EXF */
#endif /* EXF_OPTIONS_H */
Chunyan
The University of Dundee is a registered Scottish Charity, No: SC015096
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130114/619f0352/attachment.htm>
More information about the MITgcm-support
mailing list