[MITgcm-support] EXF Interpolation problem.

Josué Martínez Moreno josue_pingu at ciencias.unam.mx
Tue Jan 19 17:16:18 EST 2016


Hi everyone,

I'm trying to build on Ubuntu 14.04, gfortran 4.8 and Intel i7 (8-core)
(Compiled on MPI), an oceanic experiment of the Gulf of México, with a
external forcing package (Atmospheric Forcing).

But I get the following error:

(PID.TID 0001.0001) *** ERROR *** EXF_INTERP_UV:
fileU="u_wind_10m_NCEP_2014_2015.bin", rec=     1
(PID.TID 0001.0001) *** ERROR *** EXF_INTERP_UV: input grid must encompass
output grid.
(PID.TID 0001.0001) *** ERROR *** i,j,bi,bj=       1       1     1     1 ,
xG,yG=  2.640750E+02  1.812500E+01
(PID.TID 0001.0001) *** ERROR *** nxIn=      192 , x_in(0,nxIn+1)=
 2.398750E+02  2.881250E+02
(PID.TID 0001.0001) *** ERROR *** nyIn=       92 , y_in(0,nyIn+1)=
 9.875000E+00  1.492500E+01

I send to you a piece of my Matlab code to generate the binary files (u_wind
& v_wind ). The domain that I take each 6 hours from the NCEP-NARR data was:

Lon=[240:0.25:287.75];
Lat=[10:0.25:32.75]';

size(Lon) =  192
size(lat2) =  92

griddata(lat_ncep,lon_ncep,u_wind,Lat,Lon,'linear');

And finally the configuration of data.exf:
#
&EXF_NML_02
#
 uwindstartdate1 = 20140101,
 uwindstartdate2 = 000000,
 uwindperiod = 21600.0,
#
&EXF_NML_04
#
 uwind_lon0    = 240.125,
 uwind_lon_inc = 0.25,
 uwind_lat0    = 10.125,
 uwind_lat_inc = 0.25,
 uwind_nlon    = 192,
 uwind_nlat    = 92,
#

Something to notice here is that the nxIn values are seemingly correct
(239.875 to
288.125 degrees) but the y_in values are completely wrong (9.875 to 14.925
degrees).

Could you help me?

Have a great day and thank you!
Universidad Nacional Autónoma de México
Muchas gracias
Josué Martínez
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160119/e2533bec/attachment.htm>


More information about the MITgcm-support mailing list