[MITgcm-support] Query regarding "exf package".

Menemenlis, Dimitris (3248) Dimitris.Menemenlis at jpl.nasa.gov
Thu Dec 29 10:29:40 EST 2011


Himansu, based on startdate in data.cal,
make sure that there are forcing records available for your
first integration time step.

Specifically, based on data.exf below, your startdate in
data.cal has to be on or after 
 startdate1=20110101,
 startdate2=000000,

Also make sure that your initial temperature and salinity
files are consistent with your model bathymetry.

Good luck

Dimitris Menemenlis

On Dec 29, 2011, at 5:06 AM, himansu pradhan wrote:

> 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.
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list