How to write nx which bigger than 1000?

[yaqiang]

From ARL data format document we know the "I3" was used as string format.
(3I3) 1) Number x points, 2) Number y points, 3) Number levels

Now I have a data with the number points in x and y directions as 1280 and 640. So how to write x number to ARL data format? How to read out 1280 if only write 280 in the data file?

[yaqiang]

I found the answer from "pakset.f".

Code: Select all

!==>adjust grid parameters to deal with large grid sizes

  IF(GV(NG)%NXG.GE.1000.OR.GV(NG)%NYG.GE.1000)THEN
!    header record does not support grids of more than 999 
!    therefore in those situations the grid number is
!    converted to character to represent the 1000s digit
!    e.g. @(64)=<1000, A(65)=1000, B(66)=2000, etc

     KNX=GV(NG)%NXG/1000
     KNY=GV(NG)%NYG/1000
     WRITE(GV(NG)%IGC,'(A2)') CHAR(KNX+64)//CHAR(KNY+64)
     GV(NG)%XGPT=.TRUE.
  ELSE 
     WRITE(GV(NG)%IGC,'(I2)') GV(NG)%IG
     GV(NG)%XGPT=.FALSE. 
  END IF