source: trunk/LMDZ.COMMON/libf/dyn3d/interp_horiz.F @ 3666

Last change on this file since 3666 was 3630, checked in by afalco, 5 months ago

Pluto: larger size for array for High Resolution.
AF

File size: 4.3 KB
RevLine 
[1]1c
2c $Id: interp_horiz.F 1403 2010-07-01 09:02:53Z fairhead $
3c
4      subroutine interp_horiz (varo,varn,imo,jmo,imn,jmn,lm,
[3630]5     &  rlonuo,rlatvo,rlonun,rlatvn)
[1]6
7c===========================================================
8c  Interpolation Horizontales des variables d'une grille LMDZ
9c (des points SCALAIRES au point SCALAIRES)
10c  dans une autre grille LMDZ en conservant la quantite
11c  totale pour les variables intensives (/m2) : ex : Pression au sol
12c
13c Francois Forget (01/1995)
14c===========================================================
15
[3630]16      IMPLICIT NONE
[1]17
18c   Declarations:
19c ==============
20c
21c  ARGUMENTS
22c  """""""""
[3630]23
[1]24       integer imo, jmo ! dimensions ancienne grille (input)
25       integer imn,jmn  ! dimensions nouvelle grille (input)
26
27       real rlonuo(imo+1)     !  Latitude et
28       real rlatvo(jmo)       !  longitude des
[3630]29       real rlonun(imn+1)     !  bord des
[1]30       real rlatvn(jmn)     !  cases "scalaires" (input)
31
32       integer lm ! dimension verticale (input)
33       real varo (imo+1, jmo+1,lm) ! var dans l'ancienne grille (input)
34       real varn (imn+1,jmn+1,lm) ! var dans la nouvelle grille (output)
35
36c Autres variables
37c """"""""""""""""
38       real airetest(imn+1,jmn+1)
39       integer ii,jj,l
40
41       real airen (imn+1,jmn+1) ! aire dans la nouvelle grille
42c    Info sur les ktotal intersection entre les cases new/old grille
43       integer kllm, k, ktotal
[3630]44       parameter (kllm = 1400*200*10)
[1]45       integer iik(kllm), jjk(kllm),jk(kllm),ik(kllm)
46       real intersec(kllm)
47       real R
48       real totn, tots
49
50       logical firstcall, firsttest, aire_ok
51       save firsttest
52       data firsttest /.true./
53       data aire_ok /.true./
54
55
56
57
[3630]58
[1]59c initialisation
60c --------------
61c Si c'est le premier appel, on prepare l'interpolation
62c en calculant pour chaque case autour d'un point scalaire de la
63c nouvelle grille, la surface  de intersection avec chaque
64c    case de l'ancienne grille.
65
66
67        call iniinterp_horiz (imo,jmo,imn,jmn ,kllm,
68     &       rlonuo,rlatvo,rlonun,rlatvn,
69     &          ktotal,iik,jjk,jk,ik,intersec,airen)
70
71      do l=1,lm
72       do jj =1 , jmn+1
73        do ii=1, imn+1
74          varn(ii,jj,l) =0.
75        end do
76       end do
[3630]77      end do
78
[1]79c Interpolation horizontale
80c -------------------------
81c boucle sur toute les ktotal intersections entre les cases
82c de l'ancienne et la  nouvelle grille
83c
84      PRINT *, 'ktotal 1 = ', ktotal
[3630]85
[1]86      do k=1,ktotal
87        do l=1,lm
[3630]88         varn(iik(k),jjk(k),l) = varn(iik(k),jjk(k),l)
[1]89     &        + varo(ik(k), jk(k),l)*intersec(k)/airen(iik(k),jjk(k))
90        end do
91      end do
92
93c Une seule valeur au pole pour les variables ! :
94c -----------------------------------------------
95       do l=1, lm
96         totn =0.
97         tots =0.
98           do ii =1, imn+1
99             totn = totn + varn(ii,1,l)
100             tots = tots + varn (ii,jmn+1,l)
[3630]101           end do
[1]102           do ii =1, imn+1
103             varn(ii,1,l) = totn/REAL(imn+1)
104             varn(ii,jmn+1,l) = tots/REAL(imn+1)
[3630]105           end do
[1]106       end do
107
[3630]108
[1]109c---------------------------------------------------------------
[3630]110c  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST
[1]111!!       if (.not.(firsttest)) goto 99
112!!       firsttest = .false.
113!! !     write (*,*) 'INTERP. HORIZ. : TEST SUR LES AIRES:'
114!!       do jj =1 , jmn+1
115!!         do ii=1, imn+1
116!!           airetest(ii,jj) =0.
117!!         end do
[3630]118!!       end do
[1]119!!       PRINT *, 'ktotal = ', ktotal
120!!       PRINT *, 'jmn+1 =', jmn+1, 'imn+1', imn+1
[3630]121!!
[1]122!!       do k=1,ktotal
[3630]123!!          airetest(iik(k),jjk(k))= airetest(iik(k),jjk(k)) +intersec(k)
[1]124!!       end DO
[3630]125!!
126!!
[1]127!!       PRINT *, 'fin boucle'
128!!       do jj =1 , jmn+1
129!!        do ii=1, imn+1
130!!          r = airen(ii,jj)/airetest(ii,jj)
131!!          if ((r.gt.1.001).or.(r.lt.0.999)) then
132!! !             write (*,*) '********** PROBLEME D'' AIRES !!!',
133!! !     &                   ' DANS L''INTERPOLATION HORIZONTALE'
134!! !             write(*,*)'ii,jj,airen,airetest',
135!! !     &          ii,jj,airen(ii,jj),airetest(ii,jj)
136!!              aire_ok = .false.
137!!          end if
138!!        end do
139!!       end do
140!! !      if (aire_ok) write(*,*) 'INTERP. HORIZ. : AIRES OK'
141!!  99   continue
142
143c FIN TEST  FIN TEST  FIN TEST  FIN TEST  FIN TEST  FIN TEST  FIN TEST
144c---------------------------------------------------------------
145
146
147
148
149
150
151
152
153        return
154        end
Note: See TracBrowser for help on using the repository browser.