Opened 6 years ago

Closed 6 years ago

#1355 closed Bug (fixed)

Out of range in PISCES/p4zopt

Reported by: omamce Owned by: cetlod
Priority: low Milestone:
Component: TOP Version: release-3.6
Severity: Keywords:
Cc:

Description

I have a crash in p4zopt in config ORCA1/LIM2/PISCES

forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
opa.xx             000000000146E569  Unknown               Unknown  Unknown
opa.xx             000000000146CE3E  Unknown               Unknown  Unknown
opa.xx             0000000001425452  Unknown               Unknown  Unknown
opa.xx             00000000013B9693  Unknown               Unknown  Unknown
opa.xx             00000000013BFEAB  Unknown               Unknown  Unknown
libpthread.so.0    00002AFA81390710  Unknown               Unknown  Unknown
opa.xx             0000000000B9710B  p4zopt_mp_p4z_opt         179  p4zopt.f90
opa.xx             0000000000EE43BB  p4zbio_mp_p4z_bio         153  p4zbio.f90
opa.xx             0000000000C16E92  p4zsms_mp_p4z_sms         107  p4zsms.f90
opa.xx             0000000000DBCA74  trcsms_mp_trc_sms          50  trcsms.f90
opa.xx             00000000006A5C1B  trcstp_mp_trc_stp         121  trcstp.f90
opa.xx             000000000055080B  step_mp_stp_              230  step.f90
opa.xx             0000000000511FE5  nemogcm_mp_nemo_g         153  nemogcm.f90
opa.xx             0000000000511F2A  MAIN__                     18  nemo.f90
opa.xx             0000000000511ED6  Unknown               Unknown  Unknown
libc.so.6          00002AFA815BCD1D  Unknown               Unknown  Unknown
opa.xx             0000000000511DC9  Unknown               Unknown  Unknown

The irgb index can have values LE 0 or GT 61. I have done a small patch, but I'm not sure it's convenient :

48    INTEGER, PARAMETER :: jp_rgb = 61

and (sea line 177)

174             DO ji = 1, jpi
175                zchl = ( trn(ji,jj,jk,jpnch) + trn(ji,jj,jk,jpdch) + rtrn ) * 1.e6
176                zchl = MIN(  10. , MAX( 0.05, zchl )  )
177                irgb = MAX ( 1, MIN ( jp_rgb, NINT( 41 + 20.* LOG10( zchl ) + rtrn )))
178                !
179                zekb(ji,jj,jk) = xkrgb(1,irgb) * e3t_0(ji,jj,jk)
180                zekg(ji,jj,jk) = xkrgb(2,irgb) * e3t_0(ji,jj,jk)
181                zekr(ji,jj,jk) = xkrgb(3,irgb) * e3t_0(ji,jj,jk)
182             END DO
183          END DO
184       END DO

Olivier

Commit History (0)

(No commits)

Change History (2)

comment:1 Changed 6 years ago by clevy

  • Owner changed from NEMO team to cetlod

comment:2 Changed 6 years ago by cetlod

  • Resolution set to fixed
  • Status changed from new to closed

OF course irgb can have values LE 0 or GT 61 when the values of trn(ji,jj,jk,jpnch)+trn(ji,jj,jk,jpdch) become too high.
This happens mostly when the model exploded. The proposed patch will simply hide the problem

Note: See TracTickets for help on using tickets.