[MITgcm-support] monthly forcing data

[Martin Losch]

Just a wild guess: 

in EXF_INTERP_SIZE.h exf_max_nlon/nlat are set to 520/260, maybe that’s too small for your input file (1440/721). Try increasing that or use EXF_INTERP_USE_DYNALLOC (if that’s not already defined, by default).

OR if EXF_INTERP_USE_DYNALLOC is already defined, there may be an issue with the dynamic allocation of the arrays used for reading the input, in that case you could undefined this flag and then really increase exf_max_nlon/nlat. 


Martin

> On 3. Apr 2022, at 07:52, ESTANISLAO GAVILAN PASCUAL-AHU <e.gavilan at hhu.edu.cn> wrote:
> 
> Hi Martin,
> 
> This is the specific error
> ===================================================================================
> =   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
> =   PID 122635 RUNNING AT b1u23n3
> =   EXIT CODE: 11
> =   CLEANING UP REMAINING PROCESSES
> =   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
> ===================================================================================
>    Intel(R) MPI Library troubleshooting guide:
>       https://software.intel.com/node/561764
> ===================================================================================
> 
> After reading the code it seems the input file dimesions are wrong. However, I have checked the dimesions and they should be fine. 
> These are the dimensions setting for the 0.25 resolution
>  runoff_lon0       = 0.0D0,
>  runoff_lon_inc    = 0.25D0,
>  runoff_lat0       = -90.D0,
>  runoff_lat_inc    = 720*0.25D0,
>  runoff_nlon       = 1440,
>  runoff_nlat       = 721,
> This is the one for 1 degree resolution
> runoff_lon0       = 0.5D0,
> runoff_lon_inc    = 1.D0,
> runoff_lat0       = -89.5D0,
> runoff_lat_inc    = 179*1.D0,
> runoff_nlon       = 360,
> runoff_nlat       = 180,
> 
> I have also tried using zero matrix and using the flag -12. Every time that I used the 0.25 resolution the model cannot open the file in the routine exf_interp_read.F.
> I have also used the debug model true,  and this is the last line I get.
> 
> (PID.TID 0000.0001) EXF_SET_FLD: processing field "runoff",  file: runoff.bin
> (PID.TID 0000.0001) EXF_SET_FLD:   myIter, count0, count1:         0     12      1
> (PID.TID 0000.0001) EXF_SET_FLD:   first, changed, fac:   T   F    0.49798387096774194
> (PID.TID 0000.0001) EXF_SET_FLD: field "runoff", it=         0, loading rec=    12
> (PID.TID 0000.0001) EXF_SET_FLD:   from file: "runoff.bin"
> 
> Maybe, I am missing something obvious, but I cannot see it. Thank again for all your help
> 
> Kind regards,
> 
> Estanislao
> 
> _______________________________________________
> 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/20220404/29c03dff/attachment.html>