Opened 7 years ago

Closed 7 years ago

#1102 closed Bug (invalid)

mask send to IOM is wrong (ln_mskland=TRUE)

Reported by: marti Owned by: nemo
Priority: low Milestone:
Component: OCE Version: release-3.5
Severity: Keywords:
Cc:

Description

iom.F90, line 1094 and following :

  IF ( ln_mskland ) THEN
   ! mask land points, keep values on coast line -> specific mask for U, V and W points
     SELECT CASE ( cdgrd )
        CASE('T')   ;   zmask(:,:,:)       = tmask(:,:,:)
....
   CALL iom_set_domain_attr( "grid_"//cdgrd       , mask = zmask(:,:,1) /= 0. )
   CALL iom_set_grid_attr  ( "grid_"//cdgrd//"_3D", mask = zmask(:,:,:) /= 0. )

and zmask has dimension (jpi,jpj,jpk).

This should be done like lon/lat :

CALL iom_set_domain_attr("grid_"//cdgrd, lonvalue =    RESHAPE(plon(nldi:nlei, nldj:nlej),(/ ni*nj /)),   & 
   &                                     latvalue = RESHAPE(plat(nldi:nlei, nldj:nlej),(/ ni*nj /)))

Commit History (0)

(No commits)

Change History (1)

comment:1 Changed 7 years ago by smasson

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

same as #1103

Note: See TracTickets for help on using tickets.