[MITgcm-devel] problem in pkg cal ?

[heimbach at mit.edu]

Hi Jean-Michel,

I'll have a look at it in the next few days.

As a rule of thumb, this package is quite
fragile (or messy?).
Also too many computations based on Integer
(has partly been corrected).

C'Novim Godom from Moscow
-Patrick



Quoting Jean-Michel Campin <jmc at ocean.mit.edu>:

> Hi,
> 
> I started to run all the verification exp(s) with ifc -C (= check-bound
> option) on my laptop. I get few fails that I can fix easily, but
> one in pkg cal seems to be a real problem, and affects exp. 
> global_with_exf & lab_sea :
> 
> --> beginning of output
> (PID.TID 0000.0001) // Calendar configuration  >>> END <<<
> (PID.TID 0000.0001) //
> =======================================================
> (PID.TID 0000.0001)
> (PID.TID 0000.0001) Processing namelist file data.exf ...
> 
> Run-Time Error 406: Array bounds exceeded
>    In Procedure: cal_timepassed
> Diagnostics Entered From Subroutine cal_timepassed  Line  725
> Entered From Subroutine exf_readparms  Line  1776
> Entered From Subroutine packages_readparms  Line  1475
> Entered From Subroutine initialise_fixed  Line  1468
> Entered From Subroutine the_model_main  Line  1961
> Entered From MAIN PROGRAM   Line  3616
> End of diagnostics
> 
> <-- end of output
> 
> And when I add this write:
>     150               write(0,'(A,2I6)') 'mi,li=',mi,li
>     151               ndays = (ndaysnoleap - 1) + cal_IsLeap( yi, mythid ) -
>     152      &                cdi - ndaymonth(mi,li)
> to the file cal_timepassed.F
> with ndaymonth being declared: integer ndaymonth(12,2)
> I get (in global_with_exf):
> 
> mi,li=    10     1
> mi,li=    10     1
> mi,li=    10     1
> mi,li=    10     1
> mi,li=    10     1
> mi,li=     0     1
> Run-Time Error 406: Array bounds exceeded
>    In Procedure: cal_timepassed
> Diagnostics Entered From Subroutine cal_timepassed  Line  730
> Entered From Subroutine exf_readparms  Line  1889
> Entered From Subroutine packages_readparms  Line  1467
> Entered From Subroutine initialise_fixed  Line  1466
> Entered From Subroutine the_model_main  Line  1966
> Entered From MAIN PROGRAM   Line  3607
> 
> and a silimar problem in lab_sea:
> mi,li=    10     1
> mi,li=    10     1
> mi,li=    12     1
> mi,li=    10     1
> mi,li=    12     1
> mi,li=    10     1
> mi,li=    12     1
> mi,li=    10     1
> mi,li=    12     1
> mi,li=    10     1
> mi,li=    12     1
> mi,li=    10     1
> mi,li=    12     1
> mi,li=    10     1
> mi,li=    12     1
> mi,li=    10     1
> mi,li=    12     1
> mi,li=     0     1
> Run-Time Error 406: Array bounds exceeded
>    In Procedure: cal_timepassed
> Diagnostics Entered From Subroutine cal_timepassed  Line  726
> Entered From Subroutine exf_readparms  Line  1776
> Entered From Subroutine packages_readparms  Line  1475
> Entered From Subroutine initialise_fixed  Line  1468
> Entered From Subroutine the_model_main  Line  1961
> Entered From MAIN PROGRAM   Line  3616
> End of diagnostics
> 
> I don't have time now to look in details to this package cal, 
> but it might be interesting to figure out what's going on 
> in this S/R.
> 
> See you,
> 
> Jean-Michel
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-devel
>