[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