[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