        SUBROUTINE LOADGRIDA
C
C  Use of the WINTERACTER module is compulsory
C
!!      USE WINTERACTER
      IMPLICIT NONE
C
      INCLUDE 'RESOURCE.INC'
!        CALL WGRIDLABELROW(IDF_GRID1,1,'C1')
!        CALL WGRIDLABELROW(IDF_GRID1,2,'C2')
!        CALL WGRIDLABELROW(IDF_GRID1,3,'C3')
!        CALL WGRIDLABELROW(IDF_GRID1,4,'C4')
!        CALL WGRIDLABELROW(IDF_GRID1,5,'C5')
!        CALL WGRIDLABELROW(IDF_GRID1,6,'C6')
!        CALL WGRIDLABELROW(IDF_GRID1,7,'C7')
!        CALL WGRIDLABELROW(IDF_GRID1,8,'C8')
!        CALL WGRIDLABELROW(IDF_GRID1,9,'C9')
!        CALL WGRIDLABELROW(IDF_GRID1,10,'C10')
!        CALL WGRIDLABELROW(IDF_GRID1,11,'C11')
!        CALL WGRIDLABELROW(IDF_GRID1,12,'C12')
!        CALL WGRIDLABELROW(IDF_GRID1,13,'C13')
!        CALL WGRIDLABELROW(IDF_GRID1,14,'C14')
!        CALL WGRIDLABELROW(IDF_GRID1,15,'C15')
!        CALL WGRIDLABELROW(IDF_GRID1,16,'C16')
!        CALL WGRIDLABELROW(IDF_GRID1,17,'C17')
!        CALL WGRIDLABELROW(IDF_GRID1,18,'C18')
!        CALL WGRIDLABELROW(IDF_GRID1,19,'C19')
!        CALL WGRIDLABELROW(IDF_GRID1,20,'C20')
!        CALL WGRIDLABELROW(IDF_GRID1,21,'C21')
!        CALL WGRIDLABELROW(IDF_GRID1,22,'C22')
!        CALL WGRIDLABELROW(IDF_GRID1,23,'C23')
!        CALL WGRIDLABELROW(IDF_GRID1,24,'C24')
!        CALL WGRIDLABELROW(IDF_GRID1,25,'C25')
!        CALL WGRIDLABELROW(IDF_GRID1,26,'C26')
!        CALL WGRIDLABELROW(IDF_GRID1,27,'C27')
!        CALL WGRIDLABELROW(IDF_GRID1,28,'C28')
!        CALL WGRIDLABELROW(IDF_GRID1,29,'C29')
!        CALL WGRIDLABELROW(IDF_GRID1,30,'C30')
!        CALL WGRIDLABELROW(IDF_GRID1,31,'C31')
!        CALL WGRIDLABELROW(IDF_GRID1,32,'C32')
!        CALL WGRIDLABELROW(IDF_GRID1,33,'C33')
!        CALL WGRIDLABELROW(IDF_GRID1,34,'C34')
!        CALL WGRIDLABELROW(IDF_GRID1,35,'C35')
!        CALL WGRIDLABELROW(IDF_GRID1,36,'C36')
!        CALL WGRIDLABELROW(IDF_GRID1,37,'C37')
!        CALL WGRIDLABELROW(IDF_GRID1,38,'C38')
!        CALL WGRIDLABELROW(IDF_GRID1,39,'C39')
!        CALL WGRIDLABELROW(IDF_GRID1,40,'C40')
!        CALL WGRIDLABELROW(IDF_GRID1,41,'C41')
!        CALL WGRIDLABELROW(IDF_GRID1,42,'C42')
!        CALL WGRIDLABELROW(IDF_GRID1,43,'C43')
!        CALL WGRIDLABELROW(IDF_GRID1,44,'C44')
!        CALL WGRIDLABELROW(IDF_GRID1,45,'C45')
!        CALL WGRIDLABELROW(IDF_GRID1,46,'C46')
!        CALL WGRIDLABELROW(IDF_GRID1,47,'C47')
!        CALL WGRIDLABELROW(IDF_GRID1,48,'C48')
!        CALL WGRIDLABELROW(IDF_GRID1,49,'C49')
!        CALL WGRIDLABELROW(IDF_GRID1,50,'C50')
!        CALL WGRIDLABELROW(IDF_GRID1,51,'C51')
!        CALL WGRIDLABELROW(IDF_GRID1,52,'C52')
!        CALL WGRIDLABELROW(IDF_GRID1,53,'C53')
!        CALL WGRIDLABELROW(IDF_GRID1,54,'C54')
!        CALL WGRIDLABELROW(IDF_GRID1,55,'C55')
!        CALL WGRIDLABELROW(IDF_GRID1,56,'C56')
!        CALL WGRIDLABELROW(IDF_GRID1,57,'C57')
!        CALL WGRIDLABELROW(IDF_GRID1,58,'C58')
!        CALL WGRIDLABELROW(IDF_GRID1,59,'C59')
!        CALL WGRIDLABELROW(IDF_GRID1,60,'C60')
!        CALL WGRIDLABELROW(IDF_GRID1,61,'C61')
!        CALL WGRIDLABELROW(IDF_GRID1,62,'C62')
!        CALL WGRIDLABELROW(IDF_GRID1,63,'C63')
!        CALL WGRIDLABELROW(IDF_GRID1,64,'C64')
!        CALL WGRIDLABELROW(IDF_GRID1,65,'C65')
!        CALL WGRIDLABELROW(IDF_GRID1,66,'C66')
!        CALL WGRIDLABELROW(IDF_GRID1,67,'C67')
!        CALL WGRIDLABELROW(IDF_GRID1,68,'C68')
!        CALL WGRIDLABELROW(IDF_GRID1,69,'C69')
!        CALL WGRIDLABELROW(IDF_GRID1,70,'C70')
!        CALL WGRIDLABELROW(IDF_GRID1,71,'C71')
!        CALL WGRIDLABELROW(IDF_GRID1,72,'C72')
!        CALL WGRIDLABELROW(IDF_GRID1,73,'C73')
!        CALL WGRIDLABELROW(IDF_GRID1,74,'C74')
!        CALL WGRIDLABELROW(IDF_GRID1,75,'C75')
!        CALL WGRIDLABELROW(IDF_GRID1,76,'C76')
!        CALL WGRIDLABELROW(IDF_GRID1,77,'C77')
!        CALL WGRIDLABELROW(IDF_GRID1,78,'C78')
!        CALL WGRIDLABELROW(IDF_GRID1,79,'C79')
!        CALL WGRIDLABELROW(IDF_GRID1,80,'C80')
!        CALL WGRIDLABELROW(IDF_GRID1,81,'C81')
!        CALL WGRIDLABELROW(IDF_GRID1,82,'C82')
!        CALL WGRIDLABELROW(IDF_GRID1,83,'C83')
!        CALL WGRIDLABELROW(IDF_GRID1,84,'C84')
!        CALL WGRIDLABELROW(IDF_GRID1,85,'C85')
!        CALL WGRIDLABELROW(IDF_GRID1,86,'C86')
!        CALL WGRIDLABELROW(IDF_GRID1,87,'C87')
!        CALL WGRIDLABELROW(IDF_GRID1,88,'C88')
!        CALL WGRIDLABELROW(IDF_GRID1,89,'C89')
!        CALL WGRIDLABELROW(IDF_GRID1,90,'C90')
!        CALL WGRIDLABELROW(IDF_GRID1,91,'C91')
!        CALL WGRIDLABELROW(IDF_GRID1,92,'C92')
!        CALL WGRIDLABELROW(IDF_GRID1,93,'C93')
!        CALL WGRIDLABELROW(IDF_GRID1,94,'C94')
!        CALL WGRIDLABELROW(IDF_GRID1,95,'C95')
!        CALL WGRIDLABELROW(IDF_GRID1,96,'C96')
                RETURN
                END
        SUBROUTINE LOADGRIDB1(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,8
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=49,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=9,48
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=9,48
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'R^0')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'R^1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'R^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'R^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'R^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'R^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'R^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'R^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'R^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'R^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'R^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'R^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'R^12')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,'R^13')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,'R^14')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,'R^15')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,'R^16')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'R^17')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'R^18')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,'R^19')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,'R^20')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,'R^21')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,'R^22')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,'R^23')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,'R^24')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,'R^25')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,'R^26')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,'R^27')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,'R^28')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,'R^29')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,'R^30')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,'R^31')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,'R^32')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,'R^33')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,'R^34')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,'R^35')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,'R^36')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,'R^37')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,'R^38')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,'R^39')
                RETURN
                END
        SUBROUTINE LOADGRIDB2(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=67,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,66
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,66
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'1.0')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'R*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'R*SIN(T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'(R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'(2R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'(R^2)*SIN(2T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'(R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'(3R^3-2R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'(3R^3-2R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'(R^3)*SIN(3T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'(R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'(4R^4-3R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'(6R^4-6R^2+1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'(4R^4-3R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'(R^4)*SIN(4T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'(R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'(5R^5-4R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'(10R^5-12R^3+3R)*COS(T)
!!     1  ')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'(10R^5-12R^3+3R)*SIN(T)
!!     1  ')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'(5R^5-4R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'(R^5)*SIN(5T)')
!
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,
!!     1  '(R^6)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,
!!     1  '(6R^6-5R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,
!!     1  '(15R^6-20R^4-6R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,
 !    1  '(20R^6-30R^4+12R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,
!!     1  '(15R^6-20R^4-6R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,
!!     1  '(6R^6-5R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,
!!     1  '(R^6)*SIN(6T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,
!!     1  '(R^7)*COS(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,
!!     1  '(7R^7-6R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,
!!     1  '(21R^7-30R^5+10R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,
!!     1  '(35R^7-60R^5+30R^3-4R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,
!!     1  '(35R^7-60R^5+30R^3-4R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,
!!     1  '(21R^7-30R^5+10R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,
!!     1  '(7R^7-6R^5)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,
!!     1  '(R^7)*SIN(7T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,
!!     1  '(R^8)*COS(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,
!!     1  '(8R^8-7R^6)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,
!!     1  '(28R^8-42R^6+15R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,
!!     1  '(56R^8-105R^6+60R^4-10R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,
!!     1  '(70R^8-140R^6+90R^4-20R^2+1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,
!!     1  '(56R^8-105R^6+60R^4-10R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,
!!     1  '(28R^8-42R^6+15R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,
!!     1  '(8R^8-7R^6)*SIN(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,
!!     1  '(R^8)*SIN(8T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,
!!     1  '(R^9)*COS(9T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,
!!     1  '(9R^9-8R^7)*COS(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,
!!     1  '(36R^9-56R^7+21R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,49,
!!     1  '(84R^9-168R^7+105R^5-20R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,50,
!!     1  '(126R^9-280R^7+2105R^5-60R^3+5R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,51,
!!     1  '(126R^9-280R^7+2105R^5-60R^3+5R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,52,
!!     1  '(84R^9-168R^7+105R^5-20R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,53,
!!     1  '(36R^9-56R^7+21R^5)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,54,
!!     1  '(9R^9-8R^7)*SIN(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,55,
!!     1  '(R^9)*SIN(9T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,56,
!!     1  '(R^10)*COS(10T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,57,
!!     1  '(10R^10-9R^8)*COS(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,58,
!!     1  '(45R^10-72R^8+28R^6)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,59,
!!     1  '(120R^10-252R^8+168R^6-35R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,60,
!!     1  '(210R^10-504R^8+420R^6-140R^4-15R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,61,
!!     1  '(252R^10-630R^8+560R^6-210R^4+30R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,62,
!!     1  '(210R^10-504R^8+420R^6-140R^4-15R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,63,
!!     1  '(120R^10-252R^8+168R^6-35R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,64,
!!     1  '(45R^10-72R^8+28R^6)*SIN(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,65,
!!     1  '(10R^10-9R^8)*SIN(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,66,
!!     1  '(R^10)*SIN(10T)')
                RETURN
                END
        SUBROUTINE LOADGRIDB3(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=38,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,37
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,37
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'1.0')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'R*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'R*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'(2R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'(R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'(R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'(3R^3-2R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'(3R^3-2R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'(6R^4-6R^2+1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'(R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'(R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'(4R^4-3R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'(4R^4-3R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'(10R^5-12R^3+3R)*COS(T)
!!     1  ')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'(10R^5-12R^3+3R)*SIN(T)
!!     1  ')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,
!!     1  '(20R^6-30R^4+12R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'(R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'(R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'(5R^5-4R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'(5R^5-4R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,
!!     1  '(15R^6-20R^4-6R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,
!!     1  '(15R^6-20R^4-6R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,
!!     1  '(35R^7-60R^5+30R^3-4R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,
!!     1  '(35R^7-60R^5+30R^3-4R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,
!!     1  '(70R^8-140R^6+90R^4-20R^2+1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'(R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'(R^5)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,
!!     1  '(6R^6-5R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,
!!     1  '(6R^6-5R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,
!!     1  '(21R^7-30R^5+10R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,
!!     1  '(21R^7-30R^5+10R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,
!!     1  '(56R^8-105R^6+60R^4-10R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,
!!     1  '(56R^8-105R^6+60R^4-10R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,
!!     1  '(126R^9-280R^7+2105R^5-60R^3+5R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,
!!     1  '(126R^9-280R^7+2105R^5-60R^3+5R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,
!!     1  '(252R^10-630R^8+560R^6-210R^4+30R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,
!!     1  '(924R^12-2772R^10+3150R^8-1680R^6+420R^4-42R^2+1)')
                RETURN
                END

        SUBROUTINE LOADGRIDB4(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=4,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,3
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,3
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'P-V SAG in Lens Units')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'X - Periodicity')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'Y - Periodicity')
                RETURN
                END
        SUBROUTINE LOADGRIDB5(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*25,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,96
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
        CALL ITOAA(I,AJ)
        CFNAME='User-defined C'//AJ
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,CFNAME)
                        END DO
                RETURN
                END
        SUBROUTINE LOADGRIDB6(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,8
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=49,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=9,48
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=9,48
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'R^0')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'R^1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'R^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'R^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'R^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'R^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'R^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'R^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'R^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'R^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'R^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'R^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'R^12')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,'R^13')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,'R^14')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,'R^15')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,'R^16')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'R^17')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'R^18')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,'R^19')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,'R^20')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,'R^21')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,'R^22')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,'R^23')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,'R^24')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,'R^25')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,'R^26')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,'R^27')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,'R^28')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,'R^29')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,'R^30')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,'R^31')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,'R^32')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,'R^33')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,'R^34')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,'R^35')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,'R^36')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,'R^37')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,'R^38')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,'R^39')
                RETURN
                END
        SUBROUTINE LOADGRIDB7(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=92,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,91
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,91
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'X')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'Y^')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'X^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'XY')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'X^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'X^2*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'X*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'X^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'X^3*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'X^2*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'X*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'X^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'X^4*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'X^3*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'X^2*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'X*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,'X^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,'X^5*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,'X^4*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,'X^3*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'X^2*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'X*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,'Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,'X^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,'X^6*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,'X^5*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,'X^4*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,'X^3*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,'X^2*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,'X*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,'Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,'X^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,'X^7*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,'X^6*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,'X^5*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,'X^4*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,'X^3*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,'X^2*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,'X*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,'Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,'X^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,'X^8*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,'X^7*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,49,'X^6*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,50,'X^5*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,51,'X^4*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,52,'X^3*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,53,'X^2*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,54,'X*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,55,'Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,56,'X^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,57,'X^9*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,58,'X^8*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,59,'X^7*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,60,'X^6*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,61,'X^5*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,62,'X^4*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,63,'X^3*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,64,'X^2*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,65,'X*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,66,'Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,67,'X^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,68,'X^10*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,69,'X^9*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,70,'X^8*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,71,'X^7*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,72,'X^6*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,73,'X^5*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,74,'X^4*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,75,'X^3*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,76,'X^2*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,77,'X*Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,78,'Y^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,79,'X^12')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,80,'X^11*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,81,'X^10*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,82,'X^9*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,83,'X^8*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,84,'X^7*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,85,'X^6*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,86,'X^5*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,87,'X^4*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,88,'X^3*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,89,'X^2*Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,90,'X*Y^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,91,'Y^12')
                RETURN
                END
        SUBROUTINE LOADGRIDB8(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=92,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,91
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,91
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'X')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'Y^')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'X^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'XY')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'X^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'X^2*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'X*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'X^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'X^3*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'X^2*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'X*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'X^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'X^4*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'X^3*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'X^2*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'X*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,'X^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,'X^5*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,'X^4*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,'X^3*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'X^2*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'X*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,'Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,'X^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,'X^6*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,'X^5*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,'X^4*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,'X^3*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,'X^2*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,'X*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,'Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,'X^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,'X^7*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,'X^6*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,'X^5*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,'X^4*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,'X^3*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,'X^2*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,'X*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,'Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,'X^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,'X^8*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,'X^7*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,49,'X^6*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,50,'X^5*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,51,'X^4*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,52,'X^3*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,53,'X^2*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,54,'X*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,55,'Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,56,'X^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,57,'X^9*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,58,'X^8*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,59,'X^7*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,60,'X^6*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,61,'X^5*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,62,'X^4*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,63,'X^3*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,64,'X^2*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,65,'X*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,66,'Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,67,'X^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,68,'X^10*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,69,'X^9*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,70,'X^8*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,71,'X^7*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,72,'X^6*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,73,'X^5*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,74,'X^4*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,75,'X^3*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,76,'X^2*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,77,'X*Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,78,'Y^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,79,'X^12')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,80,'X^11*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,81,'X^10*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,82,'X^9*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,83,'X^8*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,84,'X^7*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,85,'X^6*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,86,'X^5*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,87,'X^4*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,88,'X^3*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,89,'X^2*Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,90,'X*Y^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,91,'Y^12')
                RETURN
                END
        SUBROUTINE LOADGRIDB9(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=67,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,66
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,66
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'1.0')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'R*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'R*SIN(T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'(R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'(2R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'(R^2)*SIN(2T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'(R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'(3R^3-2R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'(3R^3-2R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'(R^3)*SIN(3T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'(R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'(4R^4-3R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'(6R^4-6R^2+1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'(4R^4-3R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'(R^4)*SIN(4T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'(R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'(5R^5-4R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'(10R^5-12R^3+3R)*COS(T)
!!     1  ')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'(10R^5-12R^3+3R)*SIN(T)
!!     1  ')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'(5R^5-4R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'(R^5)*SIN(5T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,
!!     1  '(R^6)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,
!!     1  '(6R^6-5R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,
!!     1  '(15R^6-20R^4-6R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,
!!     1  '(20R^6-30R^4+12R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,
!!     1  '(15R^6-20R^4-6R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,
!!     1  '(6R^6-5R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,
!!     1  '(R^6)*SIN(6T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,
!!     1  '(R^7)*COS(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,
!!     1  '(7R^7-6R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,
!!     1  '(21R^7-30R^5+10R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,
!!     1  '(35R^7-60R^5+30R^3-4R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,
!!     1  '(35R^7-60R^5+30R^3-4R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,
!!     1  '(21R^7-30R^5+10R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,
!!     1  '(7R^7-6R^5)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,
!!     1  '(R^7)*SIN(7T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,
!!     1  '(R^8)*COS(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,
!!     1  '(8R^8-7R^6)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,
!!     1  '(28R^8-42R^6+15R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,
!!     1  '(56R^8-105R^6+60R^4-10R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,
!!     1  '(70R^8-140R^6+90R^4-20R^2+1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,
!!     1  '(56R^8-105R^6+60R^4-10R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,
!!     1  '(28R^8-42R^6+15R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,
!!     1  '(8R^8-7R^6)*SIN(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,
!!     1  '(R^8)*SIN(8T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,
!!     1  '(R^9)*COS(9T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,
!!     1  '(9R^9-8R^7)*COS(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,
!!     1  '(36R^9-56R^7+21R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,49,
!!     1  '(84R^9-168R^7+105R^5-20R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,50,
!!     1  '(126R^9-280R^7+2105R^5-60R^3+5R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,51,
!!     1  '(126R^9-280R^7+2105R^5-60R^3+5R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,52,
!!     1  '(84R^9-168R^7+105R^5-20R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,53,
!!     1  '(36R^9-56R^7+21R^5)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,54,
!!     1  '(9R^9-8R^7)*SIN(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,55,
!!     1  '(R^9)*SIN(9T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,56,
!!     1  '(R^10)*COS(10T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,57,
!!     1  '(10R^10-9R^8)*COS(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,58,
!!     1  '(45R^10-72R^8+28R^6)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,59,
!!     1  '(120R^10-252R^8+168R^6-35R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,60,
!!     1  '(210R^10-504R^8+420R^6-140R^4-15R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,61,
!!     1  '(252R^10-630R^8+560R^6-210R^4+30R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,62,
!!     1  '(210R^10-504R^8+420R^6-140R^4-15R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,63,
!!     1  '(120R^10-252R^8+168R^6-35R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,64,
!!     1  '(45R^10-72R^8+28R^6)*SIN(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,65,
!!     1  '(10R^10-9R^8)*SIN(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,66,
!!     1  '(R^10)*SIN(10T)')
                RETURN
                END
        SUBROUTINE LOADGRIDB10(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=38,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,37
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,37
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'1.0')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'R*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'R*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'(2R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'(R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'(R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'(3R^3-2R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'(3R^3-2R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'(6R^4-6R^2+1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'(R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'(R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'(4R^4-3R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'(4R^4-3R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'(10R^5-12R^3+3R)*COS(T)
!!     1  ')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'(10R^5-12R^3+3R)*SIN(T)
!!     1  ')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,
!!     1  '(20R^6-30R^4+12R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'(R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'(R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'(5R^5-4R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'(5R^5-4R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,
!!     1  '(15R^6-20R^4-6R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,
!!     1  '(15R^6-20R^4-6R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,
!!     1  '(35R^7-60R^5+30R^3-4R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,
!!     1  '(35R^7-60R^5+30R^3-4R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,
!!     1  '(70R^8-140R^6+90R^4-20R^2+1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'(R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'(R^5)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,
!!     1  '(6R^6-5R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,
!!     1  '(6R^6-5R^4)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,
!!     1  '(21R^7-30R^5+10R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,
!!     1  '(21R^7-30R^5+10R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,
!!     1  '(56R^8-105R^6+60R^4-10R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,
!!     1  '(56R^8-105R^6+60R^4-10R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,
!!     1  '(126R^9-280R^7+2105R^5-60R^3+5R)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,
!!     1  '(126R^9-280R^7+2105R^5-60R^3+5R)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,
!!     1  '(252R^10-630R^8+560R^6-210R^4+30R^2-1)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,
!!     1  '(924R^12-2772R^10+3150R^8-1680R^6+420R^4-42R^2+1)')
                RETURN
                END
        SUBROUTINE LOADGRIDB11(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*25,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,96
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
        CALL ITOAA(I,AJ)
        CFNAME='User-defined C'//AJ
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,CFNAME)
                        END DO
                RETURN
                END
        SUBROUTINE LOADGRIDB12a(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=11,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,10
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
        FTFL01(11,J)=0.0D0
                        DO I=1,11
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Order Number')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!!     1  'Constuction Wavelength (microns)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!!     1  'X - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!!     1  'Y - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!!     1  'Z - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,
!!     1  '+/- 1 (real or imaginary source point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,
!!     1  'X - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,
!!     1  'Y - Coordinate - Refeference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,
!!     1  'Z - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10
!!     1  ,'+/- 1 (real or imaginary refeference point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Added Phase Code')
                RETURN
                END
        SUBROUTINE LOADGRIDB12b(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=11,11
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=12,21
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=22,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,10
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
        FTFL01(11,J)=1.0D0
                        DO I=1,21
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Order Number')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!!     1  'Constuction Wavelength (microns)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!!     1  'X - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!!     1  'Y - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!!     1  'Z - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,
!!     1  '+/- 1 (real or imaginary source point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,
!!     1  'X - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,
!!     1  'Y - Coordinate - Refeference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,
!!     1  'Z - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10
!!     1  ,'+/- 1 (real or imaginary refeference point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Added Phase Code')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,
!!     1  'Phase Term R^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,
!!     1  'Phase Term R^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,
!!     1  'Phase Term R^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,
!!     1  'Phase Term R^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,
!!     1  'Phase Term R^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,
!!     1  'Phase Term R^12')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,
!!     1  'Phase Term R^14')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,
!!     1  'Phase Term R^16')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,
!!     1  'Phase Term R^18')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,
!!     1  'Phase Term R^20')
                RETURN
                END
        SUBROUTINE LOADGRIDB12c(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=11,11
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=12,89
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=90,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,10
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
        FTFL01(11,J)=2.0D0
                        DO I=1,89
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Order Number')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!!     1  'Constuction Wavelength (microns)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!!     1  'X - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!!     1  'Y - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!!     1  'Z - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,
!!     1  '+/- 1 (real or imaginary source point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,
!!     1  'X - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,
!!     1  'Y - Coordinate - Refeference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,
!!     1  'Z - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10
!!     1  ,'+/- 1 (real or imaginary refeference point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Added Phase Code')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'X')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'Y^')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'X^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'XY')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'X^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'X^2*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'X*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,'X^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,'X^3*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,'X^2*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,'X*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'X^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,'X^4*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,'X^3*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,'X^2*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,'X*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,'Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,'X^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,'X^5*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,'X^4*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,'X^3*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,'X^2*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,'X*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,'Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,'X^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,'X^6*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,'X^5*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,'X^4*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,'X^3*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,'X^2*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,'X*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,'Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,'X^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,49,'X^7*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,50,'X^6*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,51,'X^5*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,52,'X^4*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,53,'X^3*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,54,'X^2*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,55,'X*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,56,'Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,57,'X^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,58,'X^8*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,59,'X^7*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,60,'X^6*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,61,'X^5*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,62,'X^4*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,63,'X^3*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,64,'X^2*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,65,'X*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,66,'Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,67,'X^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,68,'X^9*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,69,'X^8*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,70,'X^7*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,71,'X^6*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,72,'X^5*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,73,'X^4*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,74,'X^3*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,75,'X^2*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,76,'X*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,77,'Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,78,'X^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,79,'X^10*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,80,'X^9*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,81,'X^8*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,82,'X^7*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,83,'X^6*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,84,'X^5*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,85,'X^4*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,86,'X^3*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,87,'X^2*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,88,'X*Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,89,'Y^11')
                RETURN
                END
        SUBROUTINE LOADGRIDB12d(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=11,11
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=12,89
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=90,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,10
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
        FTFL01(11,J)=3.0D0
                        DO I=1,89
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Order Number')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!!     1  'Constuction Wavelength (microns)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!!     1  'X - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!!     1  'Y - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!!     1  'Z - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,
!!     1  '+/- 1 (real or imaginary source point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,
!!     1  'X - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,
!!     1  'Y - Coordinate - Refeference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,
!!     1  'Z - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10
!!     1  ,'+/- 1 (real or imaginary refeference point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Added Phase Code')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'|X|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'Y^')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'|X|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'|X|Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'|X|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'|X|^2*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'|X|*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,'|X|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,'|X|^3*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,'|X|^2*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,'|X|*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'|X|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,'|X|^4*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,'|X|^3*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,'|X|^2*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,'|X|*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,'Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,'|X|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,'|X|^5*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,'|X|^4*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,'|X|^3*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,'|X|^2*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,'|X|*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,'Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,'|X|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,'|X|^6*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,'|X|^5*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,'|X|^4*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,'|X|^3*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,'|X|^2*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,'|X|*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,'Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,'|X|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,49,'|X|^7*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,50,'|X|^6*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,51,'|X|^5*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,52,'|X|^4*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,53,'|X|^3*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,54,'|X|^2*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,55,'|X|*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,56,'Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,57,'|X|^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,58,'|X|^8*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,59,'|X|^7*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,60,'|X|^6*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,61,'|X|^5*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,62,'|X|^4*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,63,'|X|^3*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,64,'|X|^2*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,65,'|X|*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,66,'Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,67,'|X|^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,68,'|X|^9*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,69,'|X|^8*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,70,'|X|^7*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,71,'|X|^6*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,72,'|X|^5*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,73,'|X|^4*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,74,'|X|^3*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,75,'|X|^2*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,76,'|X|*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,77,'Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,78,'|X|^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,79,'|X|^10*Y')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,80,'|X|^9*Y^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,81,'|X|^8*Y^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,82,'|X|^7*Y^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,83,'|X|^6*Y^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,84,'|X|^5*Y^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,85,'|X|^4*Y^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,86,'|X|^3*Y^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,87,'|X|^2*Y^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,88,'|X|*Y^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,89,'Y^11')
                RETURN
                END
        SUBROUTINE LOADGRIDB12e(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=11,11
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=12,89
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=90,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,10
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
        FTFL01(11,J)=4.0D0
                        DO I=1,89
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Order Number')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!!     1  'Constuction Wavelength (microns)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!!     1  'X - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!!     1  'Y - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!!     1  'Z - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,
!!     1  '+/- 1 (real or imaginary source point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,
!!     1  'X - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,
!!     1  'Y - Coordinate - Refeference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,
!!     1  'Z - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10
!!     1  ,'+/- 1 (real or imaginary refeference point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Added Phase Code')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'X')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'|Y|^')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'X^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'X|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'X^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'X^2*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'X*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,'X^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,'X^3*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,'X^2*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,'X*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'X^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,'X^4*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,'X^3*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,'X^2*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,'X*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,'|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,'X^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,'X^5*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,'X^4*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,'X^3*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,'X^2*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,'X*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,'|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,'X^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,'X^6*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,'X^5*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,'X^4*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,'X^3*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,'X^2*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,'X*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,'|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,'X^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,49,'X^7*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,50,'X^6*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,51,'X^5*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,52,'X^4*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,53,'X^3*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,54,'X^2*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,55,'X*|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,56,'|Y|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,57,'X^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,58,'X^8*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,59,'X^7*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,60,'X^6*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,61,'X^5*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,62,'X^4*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,63,'X^3*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,64,'X^2*|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,65,'X*|Y|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,66,'|Y|^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,67,'X^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,68,'X^9*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,69,'X^8*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,70,'X^7*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,71,'X^6*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,72,'X^5*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,73,'X^4*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,74,'X^3*|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,75,'X^2*|Y|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,76,'X*|Y|^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,77,'|Y|^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,78,'X^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,79,'X^10*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,80,'X^9*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,81,'X^8*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,82,'X^7*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,83,'X^6*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,84,'X^5*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,85,'X^4*|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,86,'X^3*|Y|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,87,'X^2*|Y|^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,88,'X*|Y|^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,89,'|Y|^11')
                RETURN
                END
        SUBROUTINE LOADGRIDB12f(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=11,11
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=12,89
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=90,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,10
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
        FTFL01(11,J)=5.0D0
                        DO I=1,89
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Order Number')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!!     1  'Constuction Wavelength (microns)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!!     1  'X - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!!     1  'Y - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!!     1  'Z - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,
!!     1  '+/- 1 (real or imaginary source point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,
!!     1  'X - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,
!!     1  'Y - Coordinate - Refeference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,
!!     1  'Z - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10
!!     1  ,'+/- 1 (real or imaginary refeference point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Added Phase Code')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'|X|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'|Y|^')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'|X|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'|X||Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'|X|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'|X|^2*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'|X|*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,'|X|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,'|X|^3*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,'|X|^2*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,'|X|*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,'|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,'|X|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,'|X|^4*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,'|X|^3*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,'|X|^2*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,'|X|*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,'|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,'|X|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,'|X|^5*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,'|X|^4*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,'|X|^3*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,'|X|^2*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,'|X|*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,'|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,'|X|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,'|X|^6*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,'|X|^5*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,'|X|^4*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,'|X|^3*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,'|X|^2*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,'|X|*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,'|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,'|X|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,49,'|X|^7*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,50,'|X|^6*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,51,'|X|^5*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,52,'|X|^4*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,53,'|X|^3*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,54,'|X|^2*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,55,'|X|*|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,56,'|Y|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,57,'|X|^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,58,'|X|^8*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,59,'|X|^7*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,60,'|X|^6*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,61,'|X|^5*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,62,'|X|^4*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,63,'|X|^3*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,64,'|X|^2*|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,65,'|X|*|Y|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,66,'|Y|^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,67,'|X|^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,68,'|X|^9*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,69,'|X|^8*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,70,'|X|^7*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,71,'|X|^6*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,72,'|X|^5*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,73,'|X|^4*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,74,'|X|^3*|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,75,'|X|^2*|Y|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,76,'|X|*|Y|^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,77,'|Y|^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,78,'|X|^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,79,'|X|^10*|Y|')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,80,'|X|^9*|Y|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,81,'|X|^8*|Y|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,82,'|X|^7*|Y|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,83,'|X|^6*|Y|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,84,'|X|^5*|Y|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,85,'|X|^4*|Y|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,86,'|X|^3*|Y|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,87,'|X|^2*|Y|^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,88,'|X|*|Y|^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,89,'|Y|^11')
                RETURN
                END
        SUBROUTINE LOADGRIDB12g(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=11,11
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=12,32
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=33,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,10
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
        FTFL01(11,J)=6.0D0
                        DO I=1,32
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Order Number')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!!     1  'Constuction Wavelength (microns)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!!     1  'X - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!!     1  'Y - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!!     1  'Z - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,
!!     1  '+/- 1 (real or imaginary source point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,
!!     1  'X - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,
!!     1  'Y - Coordinate - Refeference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,
!!     1  'Z - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10
!!     1  ,'+/- 1 (real or imaginary refeference point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Added Phase Code')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,
!!     1  'Phase Term |R|^0')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,
!!     1  'Phase Term |R|^1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,
!!     1  'Phase Term |R|^2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,
!!     1  'Phase Term |R|^3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,
!!     1  'Phase Term |R|^4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,
!!     1  'Phase Term |R|^5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,
!!     1  'Phase Term |R|^6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,
!!     1  'Phase Term |R|^7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,
!!     1  'Phase Term |R|^8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,
!!     1  'Phase Term |R|^9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,
!!     1  'Phase Term |R|^10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,
!!     1  'Phase Term |R|^11')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,
!!     1  'Phase Term |R|^12')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,
!!     1  'Phase Term |R|^13')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,
!!     1  'Phase Term |R|^14')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,
!!     1  'Phase Term |R|^15')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,
!!     1  'Phase Term |R|^16')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,
!!     1  'Phase Term |R|^17')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,
!!     1  'Phase Term |R|^18')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,
!!     1  'Phase Term |R|^19')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,
!!     1  'Phase Term |R|^20')
                RETURN
                END
        SUBROUTINE LOADGRIDB12h(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*31,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=11,11
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=12,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,10
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
        FTFL01(11,J)=7.0D0
                        DO I=1,96
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Order Number')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!!     1  'Constuction Wavelength (microns)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!!     1  'X - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!!     1  'Y - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!!     1  'Z - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,
!!     1  '+/- 1 (real or imaginary source point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,
!!     1  'X - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,
!!     1  'Y - Coordinate - Refeference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,
!!     1  'Z - Coordinate - Reference Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10
!!     1  ,'+/- 1 (real or imaginary refeference point)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Added Phase Code')
                        DO I=12,96
        CALL ITOAA(I,AJ)
        CFNAME='User-defined Phase C'//AJ
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,CFNAME)
                        END DO
                RETURN
                END
        SUBROUTINE LOADGRIDB13(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=15,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,14
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,14
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Order Number')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!!     1  'Constuction Wavelength(microns)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!!     1  'X - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!!     1  'Y - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!!     1  'Z - Coordinate - Source Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'Trace Direction, Source')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'Source Real/Imaginary')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,
!     1  'X - Coordinate - Ref. Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,
!     1  'Y - Coordinate - Ref. Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,
!     1  'Z - Coordinate - Ref. Point')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Trace Direction Ref.')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'Ref. Real/Imaginary')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'Source Configuration #')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'Ref. Configuration #')
                RETURN
                END
        SUBROUTINE LOADGRIDB14(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=67,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,66
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,66
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'1.0')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'R*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'R*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'(R^2)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'(R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'(R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'(R^3)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'(R^3)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'(R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'(R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'(R^4)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'(R^4)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'(R^4)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'(R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'(R^4)*SIN(4T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'(R^5)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'(R^5)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'(R^5)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'(R^5)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'(R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'(R^5)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,
!     1  '(R^6)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,
!     1  '(R^8)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,
!     1  '(R^10)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,
!     1  '(R^6)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,
!     1  '(R^6)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,
!     1  '(R^6)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,
!     1  '(R^6)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,
!     1  '(R^6)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,
!     1  '(R^6)*SIN(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,
!     1  '(R^7)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,
!     1  '(R^7)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,
!     1  '(R^7)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,
!     1  '(R^7)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,
!     1  '(R^7)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,
!     1  '(R^7)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,
!     1  '(R^7)*COS(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,
!     1  '(R^7)*SIN(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,
!     1  '(R^8)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,
!     1  '(R^8)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,
!     1  '(R^8)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,
!     1  '(R^8)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,
!     1  '(R^8)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,
!     1  '(R^8)*SIN(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,
!     1  '(R^8)*COS(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,
!     1  '(R^8)*SIN(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,
!     1  '(R^12)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,
!     1  '(R^14)')
                RETURN
                END
        SUBROUTINE LOADGRIDB15(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=67,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,66
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,66
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'1.0')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'R*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'R*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'(R^2)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'(R^2)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'(R^2)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'(R^3)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'(R^3)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'(R^3)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'(R^3)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'(R^4)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'(R^4)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'(R^4)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'(R^4)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'(R^4)*SIN(4T)')

!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'(R^5)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'(R^5)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'(R^5)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,19,'(R^5)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,20,'(R^5)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,21,'(R^5)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,22,
!     1  '(R^6)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,23,
!     1  '(R^8)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,24,
!     1  '(R^10)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,25,
!     1  '(R^6)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,26,
!     1  '(R^6)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,27,
!     1  '(R^6)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,28,
!     1  '(R^6)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,29,
!     1  '(R^6)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,30,
!     1  '(R^6)*SIN(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,31,
!     1  '(R^7)*COS(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,32,
!     1  '(R^7)*SIN(T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,33,
!     1  '(R^7)*COS(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,34,
!     1  '(R^7)*SIN(3T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,35,
!     1  '(R^7)*COS(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,36,
!     1  '(R^7)*SIN(5T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,37,
!     1  '(R^7)*COS(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,38,
!     1  '(R^7)*SIN(7T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,39,
!     1  '(R^8)*COS(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,40,
!     1  '(R^8)*SIN(2T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,41,
!     1  '(R^8)*COS(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,42,
!     1  '(R^8)*SIN(4T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,43,
!     1  '(R^8)*COS(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,44,
!     1  '(R^8)*SIN(6T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,45,
!     1  '(R^8)*COS(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,46,
!     1  '(R^8)*SIN(8T)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,47,
!     1  '(R^12)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,48,
!     1  '(R^14)')
                RETURN
                END
        SUBROUTINE LOADGRIDB16(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=12,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,11
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,11
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
                        END DO
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Fresnel -1, C1')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'Fresnel -1, C2')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'Fresnel -1, C3')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'Fresnel -1, C4')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'Fresnel -1, C5')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'Fresnel -1, C6')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'Fresnel -1, C7')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'Fresnel -1, C8')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'Fresnel -1, C9')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'Fresnel -1, C10')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Fresnel -1, C11')
                RETURN
                END
        SUBROUTINE LOADGRIDB17(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*25,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,96
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
        CALL ITOAA(I,AJ)
        CFNAME='User-defined C'//AJ
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,CFNAME)
                        END DO
                RETURN
                END
        SUBROUTINE LOADGRIDB18(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=19,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,2)
                        END DO
                        DO I=1,18
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,18
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
                        END DO
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,'Grazing Coefficient C1')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,'Grazing Coefficient C2')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,'Grazing Coefficient C3')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,'Grazing Coefficient C4')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,'Grazing Coefficient C5')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,6,'Grazing Coefficient C6')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,7,'Grazing Coefficient C7')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,8,'Grazing Coefficient C8')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,9,'Grazing Coefficient C9')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,10,'Grazing Coefficient C10')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,11,'Grazing Coefficient C11')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,12,'Grazing Coefficient C12')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,13,'Grazing Coefficient C13')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,14,'Grazing Coefficient C14')!
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,15,'Grazing Coefficient C15')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,16,'Grazing Coefficient C16')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,17,'Grazing Coefficient C17')
!      CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,18,'Grazing Coefficient C18')
                RETURN
                END
        SUBROUTINE LOADGRIDB19(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*25,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,4
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=5,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,0)
                        END DO
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,1,FTFL01(1,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,
!     1  'File ID Number (1 to 50)')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,2,FTFL01(2,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!     1  'Dimension of Apodization Grid')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,3,FTFL01(3,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!     1  'Clear Aperture Semi-Extent')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,4,FTFL01(4,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!     1  'Appodization Scale Factor (default = 1.0)')
                RETURN
                END
        SUBROUTINE LOADGRIDB20(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*25,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,5
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=6,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,0)
                        END DO
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,1,FTFL01(1,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,
!     1  'File ID Number (1 to 50)')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,2,FTFL01(2,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!     1  'Dimension of Phase Grid')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,3,FTFL01(3,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!     1  'Clear Aperture Semi-Extent')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,4,FTFL01(4,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!     1  'UNITS: 0=waves, 1=lens units, 2=microns')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,5,FTFL01(5,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,5,
!     1  'Phase Scale Factor (default = 1.0)')
                RETURN
                END
        SUBROUTINE LOADGRIDB21(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*25,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,96
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
        CALL ITOAA(I,AJ)
        CFNAME='User-defined C'//AJ
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,CFNAME)
                        END DO
                RETURN
                END
        SUBROUTINE LOADGRIDB22(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*25,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,4
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=5,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,0)
                        END DO
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,1,FTFL01(1,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,
!     1  'File ID Number (1 to 50)')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,2,FTFL01(2,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,
!     1  'Dimension of Sag Grid')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,3,FTFL01(3,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,3,
!     1  'Clear Aperture Semi-Extent')
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,4,FTFL01(4,J)
!     1  ,'(G23.15)')
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,4,
!     1  'Sag. Scale Factor (default = 1.0)')
                RETURN
                END
        SUBROUTINE LOADGRIDB23(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*28,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,1,FTFL01(1,J)
!     1  ,'(G23.15)')
        CFNAME='Number of Sag Points'
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,1,CFNAME)
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,2,FTFL01(2,J)
!     1  ,'(G23.15)')
        CFNAME='Reference Radius Value'
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,2,CFNAME)
                        DO I=3,96
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
        CALL ITOAA(I-2,AJ)
        CFNAME='Cubic Spline Sag Point # '//AJ
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,CFNAME)
                        END DO
                RETURN
                END
        SUBROUTINE LOADGRIDB24(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
        INTEGER I,J
        CHARACTER CFNAME*25,AJ*3
C
      INCLUDE 'RESOURCE.INC'
      INCLUDE 'DATLEN.INC'
                        DO I=1,96
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,' ')
!        CALL WGRIDCLEARCELL(IDF_GRID1,2,I)
                        END DO
                        DO I=1,96
!        CALL WGRIDSTATECELL(IDF_GRID1,2,I,1)
                        END DO
                        DO I=1,96
!        CALL WGRIDPUTCELLDOUBLE(IDF_GRID1,2,I,FTFL01(I,J)
!     1  ,'(G23.15)')
        CALL ITOAA(I,AJ)
        CFNAME='User-defined C'//AJ
!        CALL WGRIDPUTCELLSTRING(IDF_GRID1,1,I,CFNAME)
                        END DO
                RETURN
                END

        SUBROUTINE DOCOEFFS96(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,15,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C15,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,16,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C16,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,17,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C17,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,18,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C18,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,19,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C19,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,20,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C20,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,21,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C21,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,22,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C22,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,23,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C23,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,24,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C24,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,25,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C25,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,26,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C26,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,27,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C27,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,28,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C28,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,29,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C29,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,30,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C30,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,31,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C31,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,32,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C32,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,33,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C33,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,34,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C34,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,35,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C35,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,36,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C36,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,37,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C37,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,38,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C38,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,39,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C39,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,40,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C40,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,41,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C41,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,42,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C42,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,43,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C43,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,44,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C44,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,45,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C45,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,46,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C46,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,47,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C47,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,48,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C48,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,49,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C49,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,50,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C50,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,51,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C51,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,52,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C52,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,53,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C53,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,54,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C54,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,55,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C55,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,56,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C56,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,57,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C57,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,58,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C58,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,59,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C59,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,60,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C60,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,61,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C61,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,62,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C62,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,63,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C63,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,64,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C64,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,65,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C65,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,66,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C66,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,67,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C67,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,68,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C68,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,69,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C69,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,70,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C70,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,71,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C71,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,72,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C72,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,73,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C73,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,74,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C74,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,75,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C75,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,76,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C76,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,77,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C77,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,78,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C78,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,79,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C79,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,80,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C80,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,81,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C81,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,82,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C82,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,83,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C83,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,84,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C84,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,85,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C85,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,86,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C86,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,87,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C87,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,88,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C88,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,89,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C89,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,90,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C90,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,91,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C91,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,92,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C92,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,93,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C93,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,94,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C94,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,95,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C95,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,96,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C96,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS48(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,15,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C15,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,16,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C16,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,17,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C17,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,18,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C18,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,19,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C19,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,20,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C20,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,21,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C21,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,22,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C22,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,23,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C23,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,24,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C24,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,25,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C25,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,26,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C26,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,27,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C27,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,28,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C28,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,29,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C29,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,30,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C30,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,31,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C31,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,32,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C32,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,33,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C33,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,34,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C34,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,35,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C35,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,36,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C36,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,37,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C37,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,38,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C38,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,39,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C39,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,40,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C40,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,41,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C41,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,42,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C42,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,43,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C43,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,44,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C44,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,45,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C45,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,46,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C46,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,47,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C47,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,48,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C48,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS31(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,15,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C15,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,16,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C16,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,17,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C17,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,18,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C18,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,19,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C19,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,20,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C20,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,21,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C21,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,22,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C22,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,23,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C23,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,24,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C24,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,25,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C25,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,26,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C26,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,27,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C27,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,28,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C28,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,29,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C29,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,30,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C30,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,31,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C31,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS18(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,15,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C15,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,16,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C16,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,17,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C17,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,18,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C18,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS14(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS11(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS3(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS4(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS5(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS37(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,15,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C15,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,16,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C16,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,17,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C17,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,18,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C18,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,19,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C19,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,20,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C20,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,21,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C21,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,22,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C22,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,23,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C23,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,24,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C24,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,25,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C25,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,26,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C26,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,27,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C27,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,28,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C28,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,29,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C29,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,30,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C30,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,31,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C31,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,32,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C32,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,33,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C33,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,34,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C34,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,35,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C35,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,36,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C36,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,37,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C37,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS89(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,15,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C15,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,16,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C16,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,17,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C17,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,18,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C18,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,19,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C19,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,20,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C20,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,21,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C21,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,22,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C22,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,23,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C23,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,24,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C24,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,25,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C25,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,26,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C26,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,27,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C27,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,28,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C28,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,29,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C29,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,30,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C30,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,31,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C31,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,32,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C32,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,33,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C33,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,34,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C34,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,35,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C35,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,36,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C36,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,37,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C37,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,38,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C38,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,39,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C39,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,40,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C40,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,41,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C41,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,42,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C42,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,43,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C43,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,44,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C44,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,45,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C45,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,46,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C46,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,47,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C47,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,48,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C48,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,49,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C49,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,50,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C50,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,51,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C51,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,52,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C52,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,53,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C53,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,54,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C54,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,55,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C55,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,56,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C56,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,57,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C57,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,58,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C58,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,59,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C59,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,60,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C60,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,61,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C61,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,62,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C62,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,63,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C63,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,64,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C64,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,65,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C65,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,66,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C66,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,67,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C67,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,68,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C68,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,69,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C69,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,70,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C70,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,71,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C71,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,72,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C72,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,73,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C73,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,74,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C74,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,75,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C75,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,76,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C76,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,77,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C77,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,78,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C78,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,79,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C79,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,80,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C80,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,81,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C81,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,82,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C82,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,83,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C83,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,84,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C84,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,85,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C85,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,86,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C86,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,87,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C87,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,88,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C88,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,89,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C89,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS66(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,15,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C15,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,16,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C16,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,17,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C17,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,18,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C18,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,19,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C19,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,20,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C20,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,21,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C21,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,22,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C22,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,23,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C23,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,24,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C24,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,25,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C25,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,26,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C26,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,27,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C27,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,28,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C28,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,29,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C29,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,30,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C30,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,31,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C31,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,32,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C32,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,33,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C33,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,34,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C34,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,35,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C35,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,36,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C36,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,37,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C37,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,38,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C38,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,39,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C39,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,40,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C40,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,41,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C41,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,42,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C42,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,43,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C43,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,44,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C44,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,45,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C45,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,46,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C46,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,47,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C47,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,48,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C48,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,49,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C49,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,50,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C50,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,51,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C51,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,52,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C52,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,53,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C53,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,54,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C54,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,55,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C55,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,56,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C56,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,57,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C57,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,58,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C58,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,59,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C59,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,60,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C60,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,61,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C61,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,62,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C62,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,63,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C63,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,64,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C64,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,65,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C65,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,66,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C66,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS39(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,15,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C15,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,16,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C16,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,17,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C17,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,18,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C18,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,19,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C19,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,20,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C20,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,21,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C21,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,22,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C22,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,23,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C23,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,24,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C24,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,25,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C25,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,26,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C26,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,27,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C27,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,28,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C28,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,29,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C29,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,30,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C30,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,31,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C31,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,32,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C32,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,33,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C33,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,34,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C34,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,35,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C35,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,36,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C36,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,37,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C37,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,38,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C38,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,39,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C39,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
        SUBROUTINE DOCOEFFS91(J)
C
C  Use of the WINTERACTER module is compulsory
C
!      USE WINTERACTER
      IMPLICIT NONE
      REAL*8 DW1
      INTEGER J
      CHARACTER AW1*23,AJ*3
      INCLUDE 'DATMAI.INC'
      INCLUDE 'DATLEN.INC'
      INCLUDE 'RESOURCE.INC'
        CALL ITOAA(J,AJ)
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,1,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C1,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,2,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C2,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,3,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C3,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,4,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C4,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,5,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C5,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,6,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C6,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,7,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C7,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,8,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C8,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,9,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C9,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,10,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C10,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,11,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C11,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,12,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C12,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,13,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C13,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,14,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C14,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,15,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C15,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,16,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C16,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,17,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C17,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,18,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C18,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,19,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C19,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,20,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C20,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,21,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C21,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,22,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C22,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,23,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C23,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,24,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C24,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,25,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C25,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,26,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C26,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,27,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C27,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,28,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C28,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,29,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C29,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,30,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C30,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,31,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C31,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,32,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C32,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,33,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C33,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,34,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C34,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,35,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C35,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,36,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C36,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,37,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C37,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,38,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C38,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,39,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C39,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,40,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C40,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,41,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C41,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,42,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C42,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,43,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C43,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,44,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C44,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,45,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C45,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,46,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C46,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,47,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C47,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,48,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C48,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,49,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C49,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,50,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C50,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,51,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C51,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,52,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C52,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,53,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C53,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,54,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C54,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,55,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C55,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,56,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C56,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,57,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C57,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,58,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C58,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,59,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C59,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,60,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C60,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,61,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C61,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,62,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C62,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,63,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C63,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,64,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C64,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,65,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C65,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,66,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C66,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,67,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C67,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,68,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C68,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,69,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C69,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,70,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C70,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,71,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C71,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,72,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C72,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,73,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C73,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,74,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C74,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,75,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C75,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,76,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C76,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,77,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C77,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,78,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C78,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,79,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C79,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,80,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C80,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,81,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C81,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,82,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C82,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,83,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C83,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,84,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C84,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,85,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C85,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,86,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C86,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,87,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C87,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,88,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C88,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,89,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C89,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,90,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C90,'//AJ//','//AW1
        CALL PROCES
!        CALL WGRIDGETCELLDOUBLE(IDF_GRID1,2,91,DW1)
        CALL DTOA23(DW1,AW1)
        INPUT='C91,'//AJ//','//AW1
        CALL PROCES
                        RETURN
                        END
