Opened 12 years ago

Closed 12 years ago

Last modified 13 months ago

#247 closed Bug (fixed)

Pb with the use of indices ji+1,jj+1in limsbc_2.F90 subroutine

Reported by: arnaud.caubel@… Owned by: ctlod
Priority: high Milestone:
Component: OCE Version: release-2
Severity: Keywords:
Cc:

Description

In the case of the use of key_vectopt_loop, we have in limsbc_2.F90 :

Do jj=2,jpjm1

Do ji=1,jpj


zutau = 0.5 * ( ztio_u(ji+1,jj) + ztio_u(ji+1,jj+1)

It seems there is a pb because the dimensions of ztio_u are (jpi,jpj).
I have noticed that because the model crashes (Invalid operation) on SX8 with -P stack and init stack=nan as compilation option.

Commit History (0)

(No commits)

Change History (4)

comment:1 Changed 12 years ago by ctlod

  • Owner changed from NEMO team to ctlod
  • Status changed from new to assigned

comment:2 in reply to: ↑ description ; follow-up: Changed 12 years ago by ctlod

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

Replying to arnaud.caubel@lsce.ipsl.fr:

In the case of the use of key_vectopt_loop, we have in limsbc_2.F90 :

Do jj=2,jpjm1

Do ji=1,jpj


zutau = 0.5 * ( ztio_u(ji+1,jj) + ztio_u(ji+1,jj+1)

It seems there is a pb because the dimensions of ztio_u are (jpi,jpj).
I have noticed that because the model crashes (Invalid operation) on SX8 with -P stack and init stack=nan as compilation option.

One way to solve this problem is to set ji loop indices from 2 to jpjm1 instead fs_2 to fs_jpim1 as currently done.

comment:3 in reply to: ↑ 2 Changed 12 years ago by ctlod

Replying to ctlod:

Replying to arnaud.caubel@lsce.ipsl.fr:

In the case of the use of key_vectopt_loop, we have in limsbc_2.F90 :

Do jj=2,jpjm1

Do ji=1,jpj


zutau = 0.5 * ( ztio_u(ji+1,jj) + ztio_u(ji+1,jj+1)

It seems there is a pb because the dimensions of ztio_u are (jpi,jpj).
I have noticed that because the model crashes (Invalid operation) on SX8 with -P stack and init stack=nan as compilation option.

One way to solve this problem is to set ji loop indices from 2 to jpjm1 instead fs_2 to fs_jpim1 as currently done.

correction, you should read: "from 2 to jpim1" instead "from 2 to jpjm1"


comment:4 Changed 13 months ago by dford

In 11251:

Branch to add Kd490 observation operator, and improve name checking when reading feedback files. See Met Office utils ticket 247.

Note: See TracTickets for help on using tickets.