source: LMDZ4/trunk/libf/phylmd/ini_coord_REGDYN.h @ 5434

Last change on this file since 5434 was 644, checked in by Laurent Fairhead, 20 years ago

Synchronisation avec tous les diagnostiques de Ionela IM
Inclusion du slab ocean IM
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.9 KB
RevLine 
[644]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.