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?
How to write nx which bigger than 1000?
Re: How to write nx which bigger than 1000?
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