포트란 고수님들 디버깅 도와주세요.. ㅠ
글쓴이: hjchoi0620 / 작성시간: 월, 2022/08/08 - 3:14오후
안녕하세요.
저는 포트란으로 구성된 모델을 돌리는 대학원생입니다.
포트란을 혼자서 공부하면서 모델을 돌리고 있는데요.
디버깅 과정에서 어떤식으로 접근해서 봐야하는지 감을 아예 못 잡고 있어서,
포트란 잘하시는 분들이 접근 방법이라도 알려주시면 좋을 것 같아서 글을 올립니다.
모델이 잘 돌아가다가 중간에 멈췄는데, 아래의 메세지를 뱉었습니다.
Error termination. Backtrace:
At line 132 of file /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/p4zopt.f90
Fortran runtime error: Index '0' of dimension 2 of array 'rkrgb' below lower bound of 1
Error termination. Backtrace:
#0 0x14f0cb20ed21 in ???
#1 0x14f0cb20f869 in ???
#2 0x14f0cb20fee6 in ???
#3 0x555f56bbf3a2 in __p4zopt_MOD_p4z_opt
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/p4zopt.f90:132
#4 0x555f575f5bf2 in __p4zbio_MOD_p4z_bio
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/p4zbio.f90:121
#5 0x555f56ccdd8d in __p4zsms_MOD_p4z_sms
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/p4zsms.f90:163
#6 0x555f570a3c6b in __trcsms_pisces_MOD_trc_sms_pisces
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/trcsms_pisces.f90:81
#7 0x555f5709f11a in __trcsms_MOD_trc_sms
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/trcsms.f90:92
#8 0x555f55d83b27 in __trcstp_MOD_trc_stp
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/trcstp.f90:144
#9 0x555f55b959a9 in __step_MOD_stp
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/step.f90:237
#10 0x555f55af1b7c in __nemogcm_MOD_nemo_gcm
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/BLD/ppsrc/nemo/nemogcm.f90:192
#11 0x555f55aecf91 in nemo
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/WORK/nemo.f90:18
#12 0x555f55aecfca in main
at /media/cmlws/Data1/hjc/NEMO/r4.0.6/cfgs/ORCA_test4/WORK/nemo.f90:11제가 보기에는 p4zopt.f90의 rkrgb라는 변수의 문제인 듯 해서, 저 코드 위에 write들을 썼는데
DO jk = 1, jpkm1
DO jj = 1, jpj
DO ji = 1, jpi
zchl = ( zchl3d(ji,jj,jk) + rtrn ) * 1.e6
zchl = MIN( 10. , MAX( 0.05, zchl ) )
irgb = NINT( 41 + 20.* LOG10( zchl ) + rtrn )
!
write(*,*) ' irgb = ',irgb
write(*,*) ' rkrgb = ',rkrgb(1,irgb)
write(*,*) ' e3t_n = ',e3t_n(ji,jj,jk)
ekb(ji,jj,jk) = rkrgb(1,irgb) * e3t_n(ji,jj,jk) !!!! << line 132
ekg(ji,jj,jk) = rkrgb(2,irgb) * e3t_n(ji,jj,jk)
ekr(ji,jj,jk) = rkrgb(3,irgb) * e3t_n(ji,jj,jk)
END DO
END DO
END DO를 썼더니 처음에 썼던 에러는 안나오는데, 이번에는 그냥 모델이 돌지 않고 에러가 났던 같은 지점에서 가만히 멈춰 있더라구요..
모델을 컴파일할 때 주었던 옵션은 아래와 같고, gfortran을 썼습니다..
%FCFLAGS -fdefault-real-8 -funroll-all-loops -cpp -fcray-pointer -ffree-line-length-none -fbounds-check -g -O0 -fbacktrace
디버깅을 하는 좋은 방법이 있으시면 알려주셔도 좋습니다. 더 궁금하신 점이 있으시면 알려주세요. 도움 부탁드립니다!!
감사합니다.
Forums:
댓글 달기