[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