source: LMDZ5/branches/testing/libf/phylmd/ini_coord_REGDYN.h @ 5446

Last change on this file since 5446 was 1910, checked in by Laurent Fairhead, 11 years ago

Merged trunk changes r1860:1909 into testing branch

  • 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.9 KB
Line 
1c
2c $Header$
3c
4       nsrf=3
5       DO nreg=1, nbregdyn
6       DO i=1, klon
7
8c       IF (debut) THEN
9         IF(rlon(i).LT.0.) THEN
10           rlonPOS(i)=rlon(i)+360.
11         ELSE
12           rlonPOS(i)=rlon(i) 
13         ENDIF
14c       ENDIF
15
16        pct_ocean(i,nreg)=0
17
18c test si c'est 1 point d'ocean
19        IF(pctsrf(i,nsrf).EQ.1.) THEN
20
21         IF(nreg.EQ.1) THEN
22
23c TROP
24          IF(rlat(i).GE.-30.AND.rlat(i).LE.30.) THEN
25           pct_ocean(i,nreg)=1
26          ENDIF
27
28c PACIFIQUE NORD
29          ELSEIF(nreg.EQ.2) THEN
30           IF(rlat(i).GE.40.AND.rlat(i).LE.60.) THEN
31            IF(rlonPOS(i).GE.160..AND.rlonPOS(i).LE.235.) THEN
32             pct_ocean(i,nreg)=1
33            ENDIF
34           ENDIF
35c CALIFORNIE ST-CU
36         ELSEIF(nreg.EQ.3) THEN
37          IF(rlonPOS(i).GE.220..AND.rlonPOS(i).LE.250.) THEN
38           IF(rlat(i).GE.15.AND.rlat(i).LE.35.) THEN
39            pct_ocean(i,nreg)=1
40           ENDIF
41          ENDIF
42c HAWAI
43        ELSEIF(nreg.EQ.4) THEN
44         IF(rlonPOS(i).GE.180..AND.rlonPOS(i).LE.220.) THEN
45          IF(rlat(i).GE.15.AND.rlat(i).LE.35.) THEN
46           pct_ocean(i,nreg)=1
47          ENDIF
48         ENDIF
49c WARM POOL
50        ELSEIF(nreg.EQ.5) THEN
51         IF(rlonPOS(i).GE.70..AND.rlonPOS(i).LE.150.) THEN
52          IF(rlat(i).GE.-5.AND.rlat(i).LE.20.) THEN
53           pct_ocean(i,nreg)=1
54          ENDIF
55         ENDIF
56        ENDIF !nbregdyn
57c TROP
58c        IF(rlat(i).GE.-30.AND.rlat(i).LE.30.) THEN
59c         pct_ocean(i)=.TRUE.
60c         WRITE(*,*) 'pct_ocean =',i, rlon(i), rlat(i)
61c          ENDIF !lon
62c         ENDIF !lat
63
64        ENDIF !pctsrf
65       ENDDO !klon
66       ENDDO !nbregdyn
67cIM 190504      ENDIF !ok_regdyn
68 
69cIM somme de toutes les nhistoW BEG
70      IF (debut) THEN
71      DO nreg = 1, nbregdyn
72       DO k = 1, kmaxm1
73        DO l = 1, lmaxm1
74         DO iw = 1, iwmax
75          nhistoWt(k,l,iw,nreg)=0.
76         ENDDO !iw
77        ENDDO !l
78       ENDDO !k
79      ENDDO !nreg
80      ENDIF !(debut) THEN
81cIM 190504 BEG
Note: See TracBrowser for help on using the repository browser.