/[lmdze]/trunk/bibio/Writefield/write_field1d.f
ViewVC logotype

Annotation of /trunk/bibio/Writefield/write_field1d.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 108 - (hide annotations)
Tue Sep 16 14:00:41 2014 UTC (9 years, 9 months ago) by guez
File size: 1145 byte(s)
Imported writefield from LMDZ. Close at the end of gcm the files which
were created by writefiled (not done in LMDZ).

Removed procedures for the output of Grads files. Removed calls to
dump2d. In guide, replaced calls to wrgrads by calls to writefield.

In vlspltqs, removed redundant programming of saturation
pressure. Call foeew from module FCTTRE instead.

Bug fix in interpre: size of w exceeding size of correponding actual
argument wg in advtrac.

In leapfrog, call guide until the end of the run, instead of six hours
before the end.

Bug fix in readsulfate_preind: type of arguments.

1 guez 108 module write_field1d_m
2    
3     implicit none
4    
5     contains
6    
7     subroutine write_field1D(name,Field)
8    
9     use int2str_m, only: int2str
10    
11     integer, parameter :: MaxDim=1
12     character(len=*) :: name
13     real, dimension(:) :: Field
14     integer, dimension(MaxDim) :: Dim
15     integer :: i,nb
16     integer, parameter :: id=10
17     integer, parameter :: NbCol=4
18     integer :: ColumnSize
19     integer :: pos
20     character(len=255) :: form
21     character(len=255) :: MaxLen
22    
23    
24     open(unit=id,file=name//'.field',form='formatted',status='replace')
25     write (id,'("----- Field '//name//'",//)')
26     Dim=shape(Field)
27     MaxLen=int2str(len(trim(int2str(Dim(1)))))
28     ColumnSize=20+6+3+len(trim(int2str(Dim(1))))
29     Nb=0
30     Pos=2
31     do i=1,Dim(1)
32     nb=nb+1
33    
34     if (MOD(nb,NbCol)==0) then
35     form='(t'//trim(int2str(pos))// ',i'//trim(MaxLen) //'," ---> ",g22.16,/)'
36     Pos=2
37     else
38     form='(t'//trim(int2str(pos))// ',i'//trim(MaxLen) //'," ---> ",g22.16," | ",)'
39     Pos=Pos+ColumnSize
40     endif
41     write (id,form,advance='no') i,Field(i)
42     enddo
43    
44     close(id)
45    
46     end subroutine write_field1D
47    
48     end module write_field1d_m

  ViewVC Help
Powered by ViewVC 1.1.21