[MITgcm-support] Problem with EXF package for global 1deg run

Abhisek Chakraborty abhisek.sac at gmail.com
Thu Jun 14 11:08:16 EDT 2012


Respected Dr. Matthew, many thanks for your kind guidance.


However, when applying the EXF package for 1deg*1deg global ocean run,
I am getting a compiler error showing:

"ABHI at SAMUDRA/oceansat3/abhisek/MODEL/TEMP/exf]mpiexec -n 8 mitgcm_uv
forrtl: error (65): floating invalid
Image              PC                Routine            Line        Source
aa                 4000000000104382  Unknown               Unknown  Unknown
aa                 40000000000D0150  Unknown               Unknown  Unknown
aa                 400000000045F090  Unknown               Unknown  Unknown
aa                 40000000003F2020  Unknown               Unknown  Unknown
aa                 4000000000482690  Unknown               Unknown  Unknown
aa                 4000000000482980  Unknown               Unknown  Unknown
aa                 4000000000365600  Unknown               Unknown  Unknown
aa                 4000000000005350  Unknown               Unknown  Unknown
libc.so.6.1        200000000045FC20  Unknown               Unknown  Unknown
aa                 4000000000005180  Unknown               Unknown  Unknown

=====================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   EXIT CODE: 134
=   CLEANING UP REMAINING PROCESSES
=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
=====================================================================================
APPLICATION TERMINATED WITH THE EXIT STRING: Aborted (signal 6)
Execution time ->   3:56.83 real    966.424 user    26.700 sys
ABHI at SAMUDRA/oceansat3/abhisek/MODEL/TEMP/exf]"

The  "data" files (data, data.exf, data.cal, data.pkg,
data.diagnostics) as well as required header files (which I have
modified) are attached here.

Can anyone please guide me to overcome such errors?

Thanking in advance,
Abhisek


On 6/14/12, Matthew Mazloff <mmazloff at ucsd.edu> wrote:
> Hi Abhisek
>
> Its not necessary for EXF.
>
> Its necessary because you are compiling the profiles package.
>
> You could remove
> profiles
> from
> packages.conf in the code folder, and then neither file will be needed
>
> Or you could just delete the line
>   profilesfiles(1)= 'prof_exemple',
> in data.profiles and then the profile package will be unused and won't
> look for prof_example.nc
>
> -Matt
>
>
>
> On Jun 13, 2012, at 9:09 AM, Abhisek Chakraborty wrote:
>
>> Dear All,
>>
>> I am doing some experiments with EXF package supplied with MITgcm. I
>> found the example case "verification/global_with_exf" (MITgcm 63m). In
>> the "input" directory there are two files "data.profiles" and
>> "prof_exemple.nc". I didn't find anything regarding these files in the
>> manual. However I found that the example case runs only with the
>> presence of these files, otherwise (if these files are removed) the
>> model stops, showing an error with missing input profile data.
>>
>> Can somebody please explain the requirements of these files? Also if
>> the profile information is necessary for EXF package, please explain
>> how to prepare the "prof_exemple.nc" for other applications using EXF
>> package.
>>
>> Thanking in advance,
>>
>> Abhisek
>>
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
-------------- next part --------------
# ******************************
# Model parameters
# Continuous equation parameters
# ******************************
 &PARM01
 tRef= 25.0,24.5,24.0 , 23.5, 23.0 ,22.5, 22.0 , 21.5,21.0 ,20.5, 20.0 ,
       19.5, 19.0 , 18.5, 18.0 ,17.5, 17.0 , 16.0 , 15.0 ,
       14.0 , 13.0 , 12.0 , 11.0 , 10.0 ,
        9.0 ,  8.0 ,  7.0 ,  6.0,   5.0 ,
        4.0 ,  3.0 ,  2.0 ,
 sRef=34.58, 34.62, 34.68, 34.72, 34.73, 
       34.65, 34.75, 34.82, 34.87, 34.90,
       34.90, 34.86, 34.78, 34.69, 34.60,
       34.58, 34.62, 34.68, 34.72, 34.73,
       34.74, 34.73, 34.73, 34.72, 34.72,
       34.71, 34.70, 34.69,34.78, 34.69, 34.60,34.72,
 no_slip_sides=.false.,
 no_slip_bottom=.TRUE.,
 viscAr=1.E-3,
 viscAh=1.E4,
 diffKhT=1.E2,
 diffKrT=1.E-5,
 diffKhS=1.E2,
 diffKrS=1.E-5,
 rhonil=1035.,
 rhoConstFresh=1000.,
 beta=1.E-11,
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 gravity=9.81,
 gBaro=9.81,
 rigidLid=.FALSE.,
 implicitFreeSurface=.true.,
 eosType='JMD95Z',
 ivdc_kappa=100.,
 implicitDiffusion=.TRUE.,
 allowFreezing=.TRUE.,
 exactConserv=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
 useCDscheme=.TRUE.,
# useNHMTerms=.TRUE.,
 readBinaryPrec=32,
 writeBinaryPrec=32,
 hFacMin=0.05,
 hFacMinDr = 50,
# globalfiles=.TRUE.,
# usesinglecpuio=.TRUE.,
&

# **************************
# Elliptic solver parameters
# **************************
 &PARM02
 cg2dMaxIters=500,
 cg2dTargetResidual=1.E-10,
 &

# ************************
 Time stepping parameters
# ************************
 &PARM03
 nIter0=      0,
 nTimeSteps = 90,
# 100 years of integration will yield a reasonable flow field
# startTime =          0.,
# endTime   = 3110400000.,
 deltaTmom = 10.,
 tauCD =     321428.,
 deltaTtracer= 86400.,
 deltaTClock = 86400.,
 deltaTfreesurf= 86400.,
 abEps = 0.1,
 pChkptFreq= 1728000.,
 dumpFreq=   2592000.,
 taveFreq=   2592000.,
 monitorFreq=1.,
 &

# Gridding parameters
 &PARM04
 usingSphericalPolarGrid=.TRUE.,
 delR= 10.,10.,10.,20., 25.,25.,25.,25.,50.,50.,
       50.,100.,100.,100.,100.,100.,100.,100.,100.,100.,100.,
       100.,100.,250.,250.,500.,500.,500.,500.,500.,500.,500.,  
ygOrigin=-79.5,
dySpacing=1.,
 dxSpacing=1.,
&

# Input datasets
 &PARM05
 bathyFile=      'bathymetry.bin',
 hydrogThetaFile='lev_t.bin',
 hydrogSaltFile= 'lev_s.bin',
 the_run_name=   'global_oce_latlon',
&
-------------- next part --------------
#
# *******************
# Calendar Parameters
# *******************
 &CAL_NML
 TheCalendar='model',
 startDate_1=20110101,
 startDate_2=000000,
 &
-------------- next part --------------
# Diagnostic Package Choices
#--------------------
#  dumpAtLast (logical): always write output at the end of simulation (default=F)
#  diag_mnc   (logical): write to NetCDF files (default=useMNC)
#--for each output-stream:
#  fileName(n) : prefix of the output file name (max 80c long) for outp.stream n
#  frequency(n):< 0 : write snap-shot output every |frequency| seconds
#               > 0 : write time-average output every frequency seconds
#  timePhase(n)     : write at time = timePhase + multiple of |frequency|
#    averagingFreq  : frequency (in s) for periodic averaging interval
#    averagingPhase : phase     (in s) for periodic averaging interval
#    repeatCycle    : number of averaging intervals in 1 cycle
#  levels(:,n) : list of levels to write to file (Notes: declared as REAL)
#                when this entry is missing, select all common levels of this list
#  fields(:,n) : list of selected diagnostics fields (8.c) in outp.stream n
#                (see "available_diagnostics.log" file for the full list of diags)
#  missing_value(n) : missing value for real-type fields in output file "n"
#  fileFlags(n)     : specific code (8c string) for output file "n"
#--------------------
# This example dumps EXF diagnostics as snapshot after 30 days
 &DIAGNOSTICS_LIST
  fields(1:4,1) = 'EXFtaux ','EXFtauy ','EXFqnet ','EXFempmr',
   fileName(1) = 'diagsEXF',
  frequency(1) = -864000.,
  fields(1:2,2)  = 'THETA   ','SALT    ',
   levels(1,2) = 1.,
   fileName(2) = 'diagSurf',
  frequency(2) = -864000.,
  fields(1:2,3)  = 'UVEL    ', 'VVEL    ',
   levels(1,3) = 2.,
   fileName(3) = 'diag_vel_k2',
  frequency(3) = -864000.,
 &

#--------------------
# Parameter for Diagnostics of per level statistics:
#--------------------
#  diagSt_mnc (logical): write stat-diags to NetCDF files (default=diag_mnc)
#  diagSt_regMaskFile : file containing the region-mask to read-in
#  nSetRegMskFile   : number of region-mask sets within the region-mask file
#  set_regMask(i)   : region-mask set-index that identifies the region "i"
#  val_regMask(i)   : region "i" identifier value in the region mask
#--for each output-stream:
#  stat_fName(n) : prefix of the output file name (max 80c long) for outp.stream n
#  stat_freq(n):< 0 : write snap-shot output every |stat_freq| seconds
#               > 0 : write time-average output every stat_freq seconds
#  stat_phase(n)    : write at time = stat_phase + multiple of |stat_freq|
#  stat_region(:,n) : list of "regions" (default: 1 region only=global)
#  stat_fields(:,n) : list of selected diagnostics fields (8.c) in outp.stream n
#                (see "available_diagnostics.log" file for the full list of diags)
#--------------------
 &DIAG_STATIS_PARMS
#stat_fields(1:5,1)  = 'ETAN    ','UVEL    ','VVEL    ','WVEL    ', 'THETA   ',
# stat_fName(1) = 'dynStDiag',
#  stat_freq(1) = -864000.,
# stat_phase(1) = 0.,
 &
-------------- next part --------------
#
# *********************
# External Forcing Data
# *********************
 &EXF_NML_01
 exf_iprec         = 32,
exf_verbose       = .TRUE.,
repeatPeriod = 31104000.0
useExfCheckRange=.FALSE.,
&

# *********************
 &EXF_NML_02
 lwdownstartdate1=20110101,
 lwdownstartdate2=000000,
 lwdownperiod=2592000.0,
#
 swdownstartdate1=20110101,
 swdownstartdate2=000000,
 swdownperiod=2592000.0,
#
 runoffstartdate1=20110101,
 runoffstartdate2=000000,
 runoffperiod=2592000.0,
#
 precipstartdate1=20110101,
 precipstartdate2=000000,
 precipperiod=2592000.0,
#
 evapstartdate1=20110101,
 evapstartdate2=000000,
 evapperiod=2592000.0,	
#
 aqhstartdate1=20110101,
 aqhstartdate2=000000,
 aqhperiod=2592000.0,
#
 atempstartdate1=20110101,
 atempstartdate2=000000,
 atempperiod=2592000.0,
#
 uwindstartdate1=20110101,
 uwindstartdate2=000000,
 uwindperiod=2592000.0,
#
 vwindstartdate1=20110101,
 vwindstartdate2=000000,
 vwindperiod=2592000.0,
#
 climsststartdate1  = 20110101,
 climsststartdate2  = 000000,
 climsstperiod      = 2592000.0,
# 2 months restoring timescale for temperature
 climsstTauRelax         = 5184000.0,
#
 climsssstartdate1  = 20110101,
 climsssstartdate2  = 000000,
 climsssperiod      = 2592000.0,
# 6 months restoring timescale for salinity
 climsssTauRelax         = 15552000.0,
#
 lwdownfile    = 'dlwrf.bin',
 swdownfile    = 'dswrf.bin',
 runofffile  = 'runoff.bin',
 precipfile  = 'prate.bin',
 evapfile = 'erate.bin'
 aqhfile = 'Q2.bin'
 atempfile = 'T2.bin'
 uwindfile = 'u10.bin'
 vwindfile = 'v10.bin'
 climsstfile  = 'lev_sst.bin',
 climsssfile  = 'lev_sss.bin',
#
 &

# *********************
 &EXF_NML_03
 &
# *********************
 &EXF_NML_04
 lwdown_lon0    = 0.9375
 lwdown_lon_inc = 1.875
 lwdown_lat0    = -87.6
 lwdown_lat_inc = 93*1.884
 lwdown_nlon    = 192
 lwdown_nlat    = 94
#
 swdown_lon0    = 0.9375
 swdown_lon_inc = 1.875
 swdown_lat0    = -87.6
 swdown_lat_inc = 93*1.884
 swdown_nlon    = 192
 swdown_nlat    = 94
#
 runoff_lon0    = 0.9375
 runoff_lon_inc = 1.875
 runoff_lat0    = -87.6
 runoff_lat_inc = 93*1.884
 runoff_nlon    = 192
 runoff_nlat    = 94
#
 precip_lon0    = 0.9375
 precip_lon_inc = 1.875
 precip_lat0    = -87.6
 precip_lat_inc = 93*1.884
 precip_nlon    = 192
 precip_nlat    = 94
#
 evap_lon0    = 0.9375
 evap_lon_inc = 1.875
 evap_lat0    = -87.6
 evap_lat_inc = 93*1.884
 evap_nlon    = 192
 evap_nlat    = 94
#
 aqh_lon0    = 0.9375
 aqh_lon_inc = 1.875
 aqh_lat0    = -87.6
 aqh_lat_inc = 93*1.884
 aqh_nlon    = 192
 aqh_nlat    = 94
#
 atemp_lon0    = 0.9375
 atemp_lon_inc = 1.875
 atemp_lat0    = -87.6
 atemp_lat_inc = 93*1.884
 atemp_nlon    = 192
 atemp_nlat    = 94
#
 uwind_lon0    = 0.9375
 uwind_lon_inc = 1.875
 uwind_lat0    = -87.6
 uwind_lat_inc = 93*1.884
 uwind_nlon    = 192
 uwind_nlat    = 94
#
 climsst_lon0    = 2.0
 climsst_lon_inc = 4.0
 climsst_lat0    = -78.0
 climsst_lat_inc = 39*4.0
 climsst_nlon    = 90
 climsst_nlat    = 40
#
 climsss_lon0    = 2.0
 climsss_lon_inc = 4.0
 climsss_lat0    = -78.0
 climsss_lat_inc = 39*4.0
 climsss_nlon    = 90
 climsss_nlat    = 40
#
 &
-------------- next part --------------
# GM+Redi package parameters:

#-from MOM :
# GM_background_K: 	G & Mc.W  diffusion coefficient
# GM_maxSlope    :	max slope of isopycnals
# GM_Scrit       :	transition for scaling diffusion coefficient
# GM_Sd          :	half width scaling for diffusion coefficient
# GM_taper_scheme:	slope clipping or one of the tapering schemes
# GM_Kmin_horiz  :	horizontal diffusion minimum value 

#-Option parameters (needs to "define" options in GMREDI_OPTIONS.h")
# GM_isopycK     :	isopycnal diffusion coefficient (default=GM_background_K)
# GM_AdvForm     :	turn on GM Advective form       (default=Skew flux form)

 &GM_PARM01 
  GM_background_K    = 1.e+3,
  GM_taper_scheme    = 'gkw91',
  GM_maxSlope        = 1.e-2,
  GM_Kmin_horiz      = 50.,
  GM_Scrit           = 4.e-3,
  GM_Sd              = 1.e-3,
# GM_Visbeck_alpha   = 0.,
# GM_Visbeck_length  = 2.e+5,
# GM_Visbeck_depth   = 1.e+3,
# GM_Visbeck_maxval_K= 2.5e+3,
 &

-------------- next part --------------
A non-text attachment was scrubbed...
Name: data.pkg
Type: application/octet-stream
Size: 108 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120614/a858e684/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: available_diagnostics.log
Type: application/octet-stream
Size: 23803 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120614/a858e684/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: DIAGNOSTICS_SIZE.h
Type: text/x-chdr
Size: 1449 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120614/a858e684/attachment-0003.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: EXF_OPTIONS.h
Type: text/x-chdr
Size: 6877 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120614/a858e684/attachment-0004.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: packages.conf
Type: application/octet-stream
Size: 202 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120614/a858e684/attachment-0005.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SIZE.h
Type: text/x-chdr
Size: 2189 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120614/a858e684/attachment-0005.bin>


More information about the MITgcm-support mailing list