source: LMDZ5/trunk/libf/misc/formcoord.F @ 5446

Last change on this file since 5446 was 1907, checked in by lguez, 11 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.1 KB
Line 
1!
2! $Header$
3!
4      subroutine formcoord(unit,n,x,a,rev,text)
5      implicit none
6      integer n,unit,ndec
7      logical rev
8      real x(n),a
9      character*4 text
10
11      integer i,id,i1,i2,in
12      real dx,dxmin
13
14      if(rev) then
15         id=-1
16         i1=n
17         i2=n-1
18         in=1
19         write(unit,3000) text(1:1)
20      else
21         id=1
22         i1=1
23         i2=2
24         in=n
25      endif
26
27      if (n.lt.2) then
28         ndec=1
29         write(unit,1000) text,n,x(1)*a
30      else
31         dxmin=abs(x(2)-x(1))
32         do i=2,n-1
33            dx=abs(x(i+1)-x(i))
34            if (dx.lt.dxmin) dxmin=dx
35         enddo
36
37         ndec=-log10(dxmin)+2
38         if(mod(n,6).eq.1) then
39            write(unit,1000) text,n,x(i1)*a
40            write(unit,2000) (x(i)*a,i=i2,in,id)
41         else
42            write(unit,1000) text,n
43            write(unit,2000) (x(i)*a,i=i1,in,id)
44         endif
45      endif
46
471000  format(a4,2x,i4,' LEVELS',43x,f12.2)
482000  format(6f12.2)
49c1000  format(a4,2x,i4,' LEVELS',43x,f12.<ndec>)
50c2000  format(6f12.<ndec>)
513000  format('FORMAT ',a1,'REV')
52      return
53
54      end
Note: See TracBrowser for help on using the repository browser.