[MITgcm-support] [EXTERNAL] How to select interpolation method in EXF_INTERP_UV
Menemenlis, Dimitris (US 329B)
dimitris.menemenlis at jpl.nasa.gov
Mon May 1 20:38:55 EDT 2023
Why don’t you try one of the standard, verification_other forcing fields as a starting point,
e.g., from here: https://github.com/MITgcm/verification_other/tree/master/global_oce_cs32
to make sure that your set-up is OK otherwise.
D.
On May 1, 2023, at 11:45 AM, Sabri MUTLU (MAM-ÇTÜE) <sabri.mutlu at tubitak.gov.tr> wrote:
Dear Dimitris,
Let me give you some details about my run. My domain area is a very small bay compared to the oceans. It is between the following coordinates, 29.2612 E --> 29.9476 E & 40.6560 N --> 40.8200 N.
I try to force this area with ERA5 data (res: 0.25 deg x 0.25 deg). Grid points are between 29 E --> 30 E and 40.5 N --> 41 N (5 x 3 points).
In theory, my forcing data is fully covering the domain area, so the forcing data can be interpolated into my domain without any problem.
Have no info in STDERR files. I have some info in the STDOUT file as follows:
(PID.TID 0000.0001) // End OBCS MONITOR field statistics
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) DEBUG_MSG: LEAVING S/R OBCS_MONITOR
(PID.TID 0000.0001) DEBUG_MSG: LEAVING S/R INITIALISE_VARIA
(PID.TID 0000.0001) DEBUG_MSG: ENTERED S/R MAIN_DO_LOOP
(PID.TID 0000.0001) DEBUG_MSG: CALLING S/R FORWARD_STEP
(PID.TID 0000.0001) DEBUG_MSG: ENTERED S/R FORWARD_STEP
(PID.TID 0000.0001) DEBUG_MSG: CALLING S/R LOAD_FIELDS_DRIVER
(PID.TID 0000.0001) DEBUG_MSG: ENTERED S/R LOAD_FIELDS_DRIVER
(PID.TID 0000.0001) DEBUG_MSG: CALLING S/R EXF_GETFORCING
(PID.TID 0000.0001) EXF_SET_UV: processing fields "uwind" & "vwind"
(PID.TID 0000.0001) EXF_SET_UV: files: u10_c1.bin & v10_c1.bin
(PID.TID 0000.0001) EXF_SET_UV: myIter, count0, count1: 0 1 2
(PID.TID 0000.0001) EXF_SET_UV: first, changed, fac: T T 0.100000000E+01
(PID.TID 0000.0001) EXF_SET_UV: fields "uwind" & "vwind", it= 0
(PID.TID 0000.0001) EXF_SET_UV: loading rec= 1 from file: "u10_c1.bin"
(PID.TID 0000.0001) EXF_SET_UV: loading rec= 1 from file: "v10_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: u10_c1.bin
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: v10_c1.bin
(PID.TID 0000.0001) EXF_INTERP_UV: fileU="u10_c1.bin", rec= 1 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_UV: fields "uwind" & "vwind", it= 0
(PID.TID 0000.0001) EXF_SET_UV: loading rec= 2 from file: "u10_c1.bin"
(PID.TID 0000.0001) EXF_SET_UV: loading rec= 2 from file: "v10_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: u10_c1.bin
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: v10_c1.bin
(PID.TID 0000.0001) EXF_INTERP_UV: fileU="u10_c1.bin", rec= 2 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: processing field "atemp", file: t2m_c1.bin
(PID.TID 0000.0001) EXF_SET_FLD: myIter, count0, count1: 0 1 2
(PID.TID 0000.0001) EXF_SET_FLD: first, changed, fac: T T 0.100000000E+01
(PID.TID 0000.0001) EXF_SET_FLD: field "atemp", it= 0, loading rec= 1
(PID.TID 0000.0001) EXF_SET_FLD: from file: "t2m_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: t2m_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="t2m_c1.bin", rec= 1 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: field "atemp", it= 0, loading rec= 2
(PID.TID 0000.0001) EXF_SET_FLD: from file: "t2m_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: t2m_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="t2m_c1.bin", rec= 2 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: processing field "aqh", file: aqh_c1.bin
(PID.TID 0000.0001) EXF_SET_FLD: myIter, count0, count1: 0 1 2
(PID.TID 0000.0001) EXF_SET_FLD: first, changed, fac: T T 0.100000000E+01
(PID.TID 0000.0001) EXF_SET_FLD: field "aqh", it= 0, loading rec= 1
(PID.TID 0000.0001) EXF_SET_FLD: from file: "aqh_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: aqh_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="aqh_c1.bin", rec= 1 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: field "aqh", it= 0, loading rec= 2
(PID.TID 0000.0001) EXF_SET_FLD: from file: "aqh_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: aqh_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="aqh_c1.bin", rec= 2 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: processing field "precip", file: tpr_c1.bin
(PID.TID 0000.0001) EXF_SET_FLD: myIter, count0, count1: 0 1 2
(PID.TID 0000.0001) EXF_SET_FLD: first, changed, fac: T T 0.100000000E+01
(PID.TID 0000.0001) EXF_SET_FLD: field "precip", it= 0, loading rec= 1
(PID.TID 0000.0001) EXF_SET_FLD: from file: "tpr_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: tpr_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="tpr_c1.bin", rec= 1 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: field "precip", it= 0, loading rec= 2
(PID.TID 0000.0001) EXF_SET_FLD: from file: "tpr_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: tpr_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="tpr_c1.bin", rec= 2 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: processing field "swdown", file: dws_c1.bin
(PID.TID 0000.0001) EXF_SET_FLD: myIter, count0, count1: 0 1 2
(PID.TID 0000.0001) EXF_SET_FLD: first, changed, fac: T T 0.100000000E+01
(PID.TID 0000.0001) EXF_SET_FLD: field "swdown", it= 0, loading rec= 1
(PID.TID 0000.0001) EXF_SET_FLD: from file: "dws_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: dws_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="dws_c1.bin", rec= 1 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: field "swdown", it= 0, loading rec= 2
(PID.TID 0000.0001) EXF_SET_FLD: from file: "dws_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: dws_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="dws_c1.bin", rec= 2 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: processing field "lwdown", file: dwl_c1.bin
(PID.TID 0000.0001) EXF_SET_FLD: myIter, count0, count1: 0 1 2
(PID.TID 0000.0001) EXF_SET_FLD: first, changed, fac: T T 0.100000000E+01
(PID.TID 0000.0001) EXF_SET_FLD: field "lwdown", it= 0, loading rec= 1
(PID.TID 0000.0001) EXF_SET_FLD: from file: "dwl_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: dwl_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="dwl_c1.bin", rec= 1 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
(PID.TID 0000.0001) EXF_SET_FLD: field "lwdown", it= 0, loading rec= 2
(PID.TID 0000.0001) EXF_SET_FLD: from file: "dwl_c1.bin"
(PID.TID 0000.0001) EXF_INTERP_READ: opening file: dwl_c1.bin
(PID.TID 0000.0001) EXF_INTERP: file="dwl_c1.bin", rec= 2 , x-Per,P.Sym= F F
(PID.TID 0000.0001) S.edge (j=-1,0,1) : proc= 0.0 0.0 0.0, yIn= 40.000000 40.250000 40.500000
(PID.TID 0000.0001) N.edge (j=+0,+1,+2) proc= 0.0 0.0 0.0, yIn= 40.751000 40.752000 40.753000
Best.
S.
________________________________
Kimden: "Menemenlis, Dimitris (US 329B)" <dimitris.menemenlis at jpl.nasa.gov>
Kime: "MITgcm Support" <mitgcm-support at mitgcm.org>
Gönderilenler: 1 Mayıs Pazartesi 2023 21:03:27
Konu: Re: [MITgcm-support] [EXTERNAL] How to select interpolation method in EXF_INTERP_UV
Any other info in STDOUT or STDERR?
If not try running with "debugMode=.TRUE.,” in eedata and “debuglevel=5,” in data
D.
On May 1, 2023, at 10:58 AM, Sabri MUTLU (MAM-ÇTÜE) <sabri.mutlu at tubitak.gov.tr> wrote:
Thank you Dimitris,
I am getting following error in the beginning of my run: STOP ABNORMAL END: S/R EXF_INTERP_UV
Unfortunately, changing the interpolation method doesn't solve the problem.
What could be the cause of my problem then?
Best.
S.
_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org
http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org
https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!O3ZFbreP1wHO7Es_3QowefxRta8bPN6kqWJO-WbwSYwhHFypK9hvkACQBxU-LvBh8btg5AMCbGQX8gfPit-gm-J6oZ9CgVw25d-X6sQ$
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230502/d83407dd/attachment-0001.html>
More information about the MITgcm-support
mailing list