[MITgcm-support] EXF Interpolation problem.
Martin Losch
Martin.Losch at awi.de
Wed Jan 20 05:15:53 EST 2016
Hi Josue,
I think you need to specify
uwind_lat_inc = 91*0.25,
see verification/global_with_exf/input/data.exf for an example.
M.
> On 19 Jan 2016, at 23:16, Josué Martínez Moreno <josue_pingu at ciencias.unam.mx> wrote:
>
> 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
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list