anchor.diff:  !
     Write(us,*)"${GMT} pscontour -P -X4 -Y8 -I -JE0/-90/12 -R0/360/-90/-52 -K ", trim(adjustl(data_file)), & 
     	     " -C"//trim(adjustl(table_name)), " > ", trim(adjustl(ps_file)) 
  -  Write(us,*)"${GMT} pscontour -G15 -W1/255 -A -JE -R -O -K ", trim(adjustl(data_file)), " -C"//trim(adjustl(table_name)), " >> ", & 
  +  Write(us,*)"${GMT} pscontour -G15 -W1/255 -A -JE -R -O -K ", trim(adjustl(data_file)), &
  +              " -C"//trim(adjustl(table_name)), " >> ", & 
                 trim(adjustl(ps_file))
     Write(us,*)"${GMT} pscoast -R -JE -Di -Ba30f30g30 -W2/255/0/0 -A0 -O -K >> ", & 
                 trim(adjustl(ps_file))
  
positive.code:   	      trim(adjustl(ps_file))
    Write(us,*)"${GMT} pscontour -I -JJ -R -O -K ", trim(adjustl(data_file)), " -C"//trim(adjustl(table_name)), " >> ", & 
                trim(adjustl(ps_file)) 
    Write(us,*)"${GMT} pscontour -G16 -W1/255 -A -JJ -R -O -K ", trim(adjustl(data_file)), " -C"//trim(adjustl(table_name)), " >> ", & 
                trim(adjustl(ps_file))
    Write(us,*)"${GMT} pscoast -R -JJ -Di -B -W2/0 -A1000 -O -K >> ", & 
                trim(adjustl(ps_file))
  
negative.code:   	      trim(adjustl(ps_file))
    Write(us,*)"${GMT} pscontour -I -JJ -R -O -K ", trim(adjustl(data_file)), " -C"//trim(adjustl(table_name)), " >> ", & 
                trim(adjustl(ps_file)) 
    Write(us,*)"${GMT} pscontour -G16 -W1/255 -A -JJ -R -O -K ", trim(adjustl(data_file)), &
                " -C"//trim(adjustl(table_name)), " >> ", & 
                trim(adjustl(ps_file))
    Write(us,*)"${GMT} pscoast -R -JJ -Di -B -W2/0 -A1000 -O -K >> ", & 
                trim(adjustl(ps_file))
  