[MITgcm-support] EXF Interpolation problem
Josué Martínez Moreno
josue_pingu at ciencias.unam.mx
Fri Jan 22 12:23:55 EST 2016
Hi everyone,
That solved the problem. Thanks Martin.
Have a great day and thank you for all!
Universidad Nacional Autonoma de Mexico
Muchas gracias
Josue Martinez
On Wed, Jan 20, 2016 at 11:00 AM, <mitgcm-support-request at mitgcm.org> wrote:
> Send MITgcm-support mailing list submissions to
> mitgcm-support at mitgcm.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mitgcm.org/mailman/listinfo/mitgcm-support
> or, via email, send a message with subject or body 'help' to
> mitgcm-support-request at mitgcm.org
>
> You can reach the person managing the list at
> mitgcm-support-owner at mitgcm.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of MITgcm-support digest..."
>
>
> Today's Topics:
>
> 1. EXF Interpolation problem. (Josu? Mart?nez Moreno)
> 2. Re: EXF Interpolation problem. (Martin Losch)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 19 Jan 2016 16:16:18 -0600
> From: Josu? Mart?nez Moreno <josue_pingu at ciencias.unam.mx>
> To: mitgcm-support at mitgcm.org
> Subject: [MITgcm-support] EXF Interpolation problem.
> Message-ID:
> <
> CAFuuELQPDNYJFLxG9rUoyXy33NKrrfntZeurJA+nHUbZG69TNw at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> 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.html
> >
>
> ------------------------------
>
> Message: 2
> Date: Wed, 20 Jan 2016 11:15:53 +0100
> From: Martin Losch <Martin.Losch at awi.de>
> To: MITgcm Support <mitgcm-support at mitgcm.org>
> Subject: Re: [MITgcm-support] EXF Interpolation problem.
> Message-ID: <6396F058-D5E3-40FC-8219-6C6F5644AA11 at awi.de>
> Content-Type: text/plain; charset="utf-8"
>
> 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
>
>
>
>
> ------------------------------
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> End of MITgcm-support Digest, Vol 151, Issue 12
> ***********************************************
>
--
Muchas gracias
Josué
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160122/4a33c2e3/attachment.htm>
More information about the MITgcm-support
mailing list