Opened 2 years ago

Last modified 15 months ago

#564 new defect

Problem with DGVM on unstructured grid in lps_establish

Reported by: jgipsl Owned by: somebody
Priority: major Milestone: Not scheduled yet
Component: Biogeochemical processes Version: trunc
Keywords: Cc:

Description

Following section in lpj_establish will not work on an unstructured grid (using DYNMICO):

   ! from how many sides is the grid box invaded (separate x and y directions
   ! because resolution may be strongly anisotropic)
   ! For the moment we only look into 4 direction but that can be expanded (JP) 
   nfrontx = 0
   IF ( neighbours(i,3) .GT. 0 ) THEN
      IF ( everywhere(neighbours(i,3),j) .GT. 1.-min_stomate ) nfrontx = nfrontx+1
   ENDIF
   IF ( neighbours(i,7) .GT. 0 ) THEN
      IF ( everywhere(neighbours(i,7),j) .GT. 1.-min_stomate ) nfrontx = nfrontx+1
   ENDIF

   nfronty = 0
   IF ( neighbours(i,1) .GT. 0 ) THEN
     IF ( everywhere(neighbours(i,1),j) .GT. 1.-min_stomate ) nfronty = nfronty+1
   ENDIF
   IF ( neighbours(i,5) .GT. 0 ) THEN
     IF ( everywhere(neighbours(i,5),j) .GT. 1.-min_stomate ) nfronty = nfronty+1
   ENDIF

   everywhere(i,j) = &
       everywhere(i,j) + migrate(j) * dt/one_year * &
      ( nfrontx / resolution(i,1) + nfronty / resolution(i,2) )

Usage of variables resolution and neighbours must be looked trough carefully.

Change History (1)

comment:1 Changed 15 months ago by luyssaert

  • Milestone set to Not scheduled yet
Note: See TracTickets for help on using tickets.