[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