[MITgcm-support] [EXTERNAL] How to select interpolation method in EXF_INTERP_UV

Sabri MUTLU (MAM-ÇTÜE) sabri.mutlu at tubitak.gov.tr
Mon May 1 14:45:31 EDT 2023


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 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230501/9dd877c3/attachment-0001.html>


More information about the MITgcm-support mailing list