anchor.diff:            !*** Write the field information
             VALUE_INDEX=1
             if(FIRST_NODE)THEN
  -             write(10,'( '' #Fields=3'' )')
  +             write(10,'( '' #Fields=5'' )')
                write(10,'('' 1) coordinates, coordinate, rectangular cartesian, #Components=3'')')
                do nj=1,3
                   if(nj.eq.1) write(10,'(2X,''x.  '')',advance="no")
  
positive.code:              !write(10,'('' 4) pressure, field, rectangular cartesian, #Components=1'')')
               !write(10,'(2X,''1.  '')',advance="no")
               !write(10,'(''Value index='',I1,'', #Derivatives='',I1)',advance="yes") VALUE_INDEX,0
               VALUE_INDEX=VALUE_INDEX+1
               write(10,'('' 3) pressure, field, rectangular cartesian, #Components=1'')')
               write(10,'(2X,''1.  '')',advance="no")
               write(10,'(''Value index='',I1,'', #Derivatives='',I1)',advance="yes") VALUE_INDEX,0
               !Compliance
               !write(10,'('' 5) compliance, field, rectangular cartesian, #Components=1'')')
               !write(10,'(2X,''1.  '')',advance="no")
               !write(10,'(''Value index='',I1,'', #Derivatives='',I1)',advance="yes") VALUE_INDEX,0
               !VALUE_INDEX=VALUE_INDEX+1
               !Pleural pressure
               !write(10,'('' 6) pleural pressure, field, rectangular cartesian, #Components=1'')')
               !write(10,'(2X,''1.  '')',advance="no")
               !write(10,'(''Value index='',I1,'', #Derivatives='',I1)',advance="yes") VALUE_INDEX,0
               !VALUE_INDEX=VALUE_INDEX+1
               !Tidal volume
               !write(10,'('' 7) tidal volume, field, rectangular cartesian, #Components=1'')')
               !write(10,'(2X,''1.  '')',advance="no")
               !write(10,'(''Value index='',I1,'', #Derivatives='',I1)',advance="yes") VALUE_INDEX,0
            endif !FIRST_NODE
            !***      write the node
            write(10,'(1X,''Node: '',I12)') np
            do nj=1,3
               write(10,'(2X,4(1X,F12.6))') (node_xyz(nj,np))      !Coordinates
            enddo !njj2
            !write(10,'(2X,4(1X,F12.6))') (unit_field(nu_vent,NOLIST)) !Ventilation
            !write(10,'(2X,4(1X,F12.6))') (unit_field(nu_vol,nolist))   !Volume (end expiration)
            !write(10,'(2X,4(1X,F12.6))') (unit_field(nu_press,nolist)) !Pressure
            write(10,'(2X,4(1X,F12.6))') (unit_field(nu_comp,nolist))  !Compliance (end exp)
  
negative.code:              write(10,'(2X,''1.  '')',advance="no")
               write(10,'(''Value index='',I1,'', #Derivatives='',I1)',advance="yes") VALUE_INDEX,0
               !VALUE_INDEX=VALUE_INDEX+1
               !Volume
               !write(10,'('' 3) volume, field, rectangular cartesian, #Components=1'')')
               !write(10,'(2X,''1.  '')',advance="no")
               !write(10,'(''Value index='',I1,'', #Derivatives='',I1)',advance="yes") VALUE_INDEX,0
  