[MITgcm-support] Query regarding "exf package".
himansu pradhan
oceancalling at gmail.com
Thu Dec 29 08:06:40 EST 2011
Hi MITGcm users,
I have been trying to use "exf package ". To be in detail : I was using
"winds stress files" in "data". Now i want to use this "wind stress files"
in in "data.exf" . My "data.exf" looks like this:
#........................data.exf........................
# External Forcing Data
# *********************
&EXF_NML_01
exf_iprec = 64,
readstressOnCgrid=.TRUE.,
&
&EXF_NML_02
ustressstartdate1=20110101,
ustressstartdate2=000000,
ustressperiod=21600,
#
vstressstartdate1=20110101,
vstressstartdate2=000000,
vstressperiod=21600,
#
ustressfile = 'pUstress.bin',
vstressfile = 'pVstress.bin',
#
&
....................................
and "data" file's "PARM03" and "PARM05" looks like :
....data..........
&PARM03
startTime=0,
endTime=21600,
deltaT=10.0,
abEps=0.1,
pChkptFreq=21600,
dumpFreq=3600,
# for time dependent open boundary conditions, activate the following 3
lines:
# periodicExternalForcing=.TRUE.,
# externForcingPeriod=10800.,
# externForcingCycle =43200.,
&
# Input datasets
&PARM05
bathyFile='top_MOD.bin',
hydrogThetaFile='15lvlT_may.bin',
hydrogSaltFile='15lvlS_may.bin',
# zonalWindFile='pUstress.bin',
# meridWindFile='pVstress.bin',
&
........................
I get the following error : ???????????????????/
[SVBabu at SVBabu input]$ ./mitgcmuv > 29dec1
** WARNINGS ** INI_THETA: found 32 wet grid-pts with theta=0
identically.
** WARNINGS ** INI_SALT: found 3999 wet grid-pts with salt=0
identically.
At line 2158 of file mdsio_read_field.f (unit = 9, file = 'pUstress.bin')
Fortran runtime error: Non-existing record number
??????????????
to note further : i use the same u/v -stress files in the "data.exf". And
my stress files are of single time slab. The model runs fine when these
stress files are used in "data". my data files are of precession : real 8 .
here is the matlab prog :
% creating BInary file 4m TExt file
ieee='b';accuracy='real*8';
u = importdata('Wu2D.txt');
W1 = -1* ((1.225)*(0.0012)*(u.*u));
fid=fopen('pUstress.bin','wb',ieee);fwrite(fid,W1',accuracy);fclose(fid);
Any suggestations/help
Thanking you in advance,
Himansu kesari pradhan
India.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20111229/b9f9d848/attachment.htm>
More information about the MITgcm-support
mailing list