source: LMDZ.3.3/trunk/libf/dyn3d/testharm.F @ 2

Last change on this file since 2 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 985 bytes
Line 
1      PROGRAM testharm
2      implicit none
3
4#include "dimensions.h"
5#include "paramet.h"
6
7      integer i,j,ii,jj
8      real pi
9      real psi(jjp1,jjp1,iim) , fphi(iim,iim)
10      character file*20, str2*2, str2b*2
11
12#include "serre.h"
13#include "fxyprim.h"
14
15      pi=2.*asin(1.)
16
17      call harmspher( fphi,psi )
18      print*,'FIN de harmspher'
19
20      do i=1,iim
21         write(str2,'(i2.2)') i
22         file='Phi'//str2
23         open(10,file=file,form='formatted')
24         print*,'i,file',i,file
25         do ii=1,iim
26            write(10,*) fx(float(ii))*180./pi,fphi(ii,i)
27         enddo
28         close(10)
29      enddo
30
31      do i=1,iim
32         do j=1,jjp1
33            write(str2,'(i2.2)') i
34            write(str2b,'(i2.2)') j
35            file='Psi.'//str2//'.'//str2b
36            open(10,file=file,form='formatted')
37            do jj=1,jjp1
38               write(10,*) fy(float(jj))*180./pi,psi(jj,j,i)
39            enddo
40            close(10)
41         enddo
42      enddo
43
44      end
Note: See TracBrowser for help on using the repository browser.