Re: Write columns

On 04/13/2012 08:59 AM, Dennis Shea wrote:
> Your setup for write_matrix is write_matrix (p, fmtx, False)
>
> ftxt = "params.soi"
> system("/bin/rm -f "+ftxt)
>
> opt = True
> opt@fout = ftxt
> opt@title = "Sand Silt Clay Porosity FC WP bexp AEP SHC"
> write_matrix (p, fmtx, opt)
>
> =========================================
>
>
>
> On 04/12/2012 05:45 PM, Guilherme Martins wrote:
>> Hi users,
>>
>> I'm using the asciiwrite to save a file, but when I use this command my
>> file is only a vector. I want the result like the below, that is,
>> columns. My script is below.
>>
>> Thanks,
>>
>> ; 0 1 2 3 4 5 6 7 8
>> ; Sand Silt Clay Porosity FC WP bexp AEP SHC
>> ; 0.92 0.05 0.03 0.437 0.091 0.033 1.7 0.07 5.83300E-05
>> ; 0.81 0.12 0.07 0.437 0.125 0.055 2.1 0.09 1.69720E-05
>> ; 0.65 0.25 0.10 0.453 0.207 0.095 3.1 0.15 7.19440E-06
>> ; 0.42 0.40 0.18 0.463 0.270 0.117 4.5 0.11 3.66670E-06
>> ; 0.20 0.65 0.15 0.501 0.330 0.133 4.7 0.21 1.88890E-06
>>
>> begin
>>
>> arquivo = asciiread("params.soi-EXER.04.QUEST.02.txt",(/11,9/),"float")
>>
>> do j = 1,4
>>
>> if (j .eq. 1) then
>> fator = 0.8
>> end if
>>
>> if (j .eq. 2) then
>> fator = 1.2
>> end if
>>
>> if (j .eq. 3) then
>> fator = 0.5
>> end if
>>
>> if (j .eq. 4) then
>> fator = 1.5
>> end if
>>
>> linhas = 11
>> colunas = 9
>>
>> sand = arquivo(:,0) * fator
>> silt = arquivo(:,1) * fator
>> clay = arquivo(:,2) * fator
>> porosity = arquivo(:,3) * fator
>> fc = arquivo(:,4) * fator
>> wp = arquivo(:,5) * fator
>> bexp = arquivo(:,6) * fator
>> aep = arquivo(:,7) * fator
>> shc = arquivo(:,8) * fator
>>
>> p = new( (/linhas,colunas/), typeof(sand) )
>>
>> p(:,0) = sand
>> p(:,1) = silt
>> p(:,2) = clay
>> p(:,3) = porosity
>> p(:,4) = fc
>> p(:,5) = wp
>> p(:,6) = bexp
>> p(:,7) = aep
>> p(:,8) = shc
>>
>> fmtx = "
>> f4.2,2x,f4.2,2x,f4.2,2x,f5.3,2x,f5.3,2x,f5.3,2x,f4.1,2x,f4.2,2x,e11.5"
>>
>> write_matrix (p, fmtx, False)
>> asciiwrite("params.soi."+fator+".txt",p)
>>
>> end do
>>
>> end
>>
>> --
>> *Guilherme Martins*
>> http://guilhermemartins.hd1.com.br/
>>
>>
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Apr 13 09:21:55 2012