[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