source: LMDZ4/trunk/libf/phy_IPCC_AR4/ini_histISCCP.h @ 1068

Last change on this file since 1068 was 868, checked in by Laurent Fairhead, 17 years ago

Preparation du remplacement de la physique utilisee pour l'exercice IPCC_AR4
par la version de la physique avec thermique. On garde le repertoire phylmd
pour un petit moment pour que les utilisateurs ne soient pas trop perdus ...
phy_IPCC_AR4 = phylmd
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 9.4 KB
Line 
1!
2! $Header$
3!
4      IF (ok_isccp) THEN
5c
6c$OMP MASTER
7      ndex2d = 0
8      ndex3d = 0
9c
10c pour les champs instantannes, il faut mettre la meme valeur pour
11c zout et zsto.
12c dtime est passe par ailleurs a histbeg
13c zstophy = frequence de stockage des champs tous les pdt physiques
14c zout = frequence d'ecriture des champs
15cIM 300505     zstophy = dtime
16c appel du simulateur toutes les 3heures
17!IM on lit la frequence d'appel dans physiq.def
18!         zcals(1) = dtime *6.  !toutes les 3h (en s)
19          zcals(1) = freq_ISCCP !toutes les freq_ISCCP secondes
20        DO n=1, napisccp
21          zcalh(n) = zcals(n)/3600. !stoutes les Xh (en heures)
22        ENDDO !n
23c
24c ecriture 8 fois par jour
25c       zout = dtime * REAL(NINT(86400./dtime*ecrit_isccp))
26c ecriture toutes les 2h (12 fois par jour)
27c       zout = dtime * 4.
28c ecriture toutes les 1/2 h (48 fois par jour)
29c       zout = dtime
30c
31c       IF(freqout_isccp.EQ.1.) THEN
32c ecriture jounaliere
33!IM on ecrit les resultats du simulateur ISCCP toutes les
34! ecrit_ISCCP secondes      zout_isccp(1) = ecrit_day !(en s)
35          zout_isccp(1) = ecrit_ISCCP !(en s)
36c ecriture mensuelle
37c         zout = dtime * ecrit_mth !(en s)
38        DO n=1, napisccp
39          zoutj(n)=zout_isccp(n)/86400. !(en jours)
40c
41c le nombre de sous-colonnes ncol : ncol=(100.*zcalh)/zoutd
42          ncol(n)=NINT((100.*zcalh(n))/zoutj(n))
43          IF(ncol(n).GT.ncolmx) THEN
44           PRINT*,'Warning: Augmenter le nombre colonnes du simulateur'
45           PRINT*,'         ISCCP ncol=', ncol,' ncolmx=',ncolmx
46c          PRINT*,'n ncol',n,ncol(n)
47           CALL abort
48          ENDIF
49c
50        DO l=1, ncol(n)
51          vertlev(l,n)=float(l)
52        ENDDO !ncol
53c
54        ENDDO !n
55
56c       PRINT*, 'La frequence de sortie ISCCP est de ', ecrit_isccp
57c
58        idayref = day_ref
59        CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
60c       write(*,*)'ISCCP ', itau_phy, zjulian
61c
62c
63c definition coordonnees lon,lat en globale
64c
65cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
66cym        DO i = 1, iim
67cym          zx_lon(i,1) = rlon(i+1)
68cym          zx_lon(i,jjmp1) = rlon(i+1)
69cym        ENDDO
70
71cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
72c
73cIM BEG region
74cym Desole dans un premier temps le mode region ne marchera pas
75cym Il faudra voir dans un second temps pour l'implementer
76cym Mais cela posera des problemes au niveau de la reconstruction
77
78          imin_ins=1
79          imax_ins=iim
80          jmin_ins=1
81          jmax_ins=jjmp1
82cym          do i=1,iim-1
83cym             if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i
84cym             if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1
85cym          enddo
86cym          do j=1,jjmp1
87cym             if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j
88cym             if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j
89cym          enddo
90c
91          print*,'On stoke le fichier histISCCP sur ',
92     s   imin_ins,imax_ins,jmin_ins,jmax_ins
93cym          print*,'On stoke le fichier histISCCP instantanne sur ',
94cym     s   zx_lon(imin_ins,1),zx_lon(imax_ins,1),
95cym     s   zx_lat(1,jmin_ins),zx_lat(1,jmax_ins)
96cIM END region
97c
98        IF(1.EQ.0) THEN
99cym         CALL histbeg("histISCCP.nc", iim,zx_lon(:,1),jjmp1,zx_lat(1,:),
100cym     .                 1, iim, 1, jjmp1,
101cym     .                 itau_phy, zjulian, dtime,
102cym     .                 nhori, nid_isccp)
103         CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime,
104     .                 nhori, nid_isccp)
105        ENDIF !(1.EQ.0) THEN
106c
107cym         CALL histbeg("histISCCP.nc", iim,zx_lon(:,1),
108cym     .                 jjmp1,zx_lat(1,:),
109cym     .                 imin_ins,imax_ins-imin_ins+1,
110cym     .                 jmin_ins,jmax_ins-jmin_ins+1,
111cym     .                 itau_phy, zjulian, dtime,
112cym     .                 nhori, nid_isccp)
113
114         CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime,
115     .                 nhori, nid_isccp)
116c
117        IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
118         CALL histvert(nid_isccp, "cldtopres","Cloud Top Pressure","mb",
119     .                 lmaxm1, cldtopres, nvert,'down')
120        ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
121         CALL histvert(nid_isccp,"cldtopres3","Cloud Top Pressure","mb",
122     .                 lmax3, cldtopres3, nvert3,'down')
123        ENDIF
124        DO n=1, napisccp
125         CALL histvert(nid_isccp, "Nbcol"//verticaxe(n),
126     .        "Nb of Column"//verticaxe(n),"1",
127     .        ncol(n), vertlev(:,n), nvlev(n),'up')
128        ENDDO
129c
130        IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
131c
132c variables a ecrire
133c
134         DO n=1, napisccp
135c
136         DO k=1, kmaxm1
137          CALL histdef(nid_isccp, "cldISCCP_"//taulev(k)//verticaxe(n),
138     .                "LMDZ ISCCP cld", "%",
139     .                iim, jj_nb,nhori,lmaxm1,1,lmaxm1,nvert,32,
140     .                "ave(X)", zcals(n),zout_isccp(n))
141         ENDDO
142c
143         CALL histdef(nid_isccp, "nsunlit"//verticaxe(n),
144     .                "Nb of calls with sunlit ", "%",
145     .                iim, jj_nb,nhori,1,1,1,-99,32,
146     .                "ave(X)", zcals(n),zout_isccp(n))
147c
148         ENDDO
149c
150        ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
151c
152         DO n=1, napisccp
153c
154c         print*,'n=',n,' avant histdef(..LMDZ ISCCP cld'
155c
156          DO k=1, kmaxm1
157           DO l=1, lmaxm1
158c
159           CALL histdef(nid_isccp, pclev(l)//taulev(k)//verticaxe(n),
160     .                "LMDZ ISCCP cld "//cnameisccp(l,k), "%",
161     .                iim, jj_nb,nhori,1,1,1,-99,32,
162     .                "ave(X)", zcals(n),zout_isccp(n))
163c
164           ENDDO
165          ENDDO
166c
167c         print*,'n=',n,' avant histdef(..Nb of calls sunlit'
168          CALL histdef(nid_isccp, "nsunlit"//verticaxe(n),
169     .                "Nb of calls with sunlit ", "%",
170     .                iim, jj_nb,nhori,1,1,1,-99,32,
171     .                "ave(X)", zcals(n),zout_isccp(n))
172c
173c 9types de nuages ISCCP-D2
174          CALL histdef(nid_isccp, "cirr",
175     .                "Cirrus lk ISCCP-D2", "%",
176     .                iim, jj_nb,nhori,1,1,1,-99,32,
177     .                "ave(X)", zcals(n),zout_isccp(n))
178c
179          CALL histdef(nid_isccp, "cist",
180     .                "CiSt lk ISCCP-D2", "%",
181     .                iim, jj_nb,nhori,1,1,1,-99,32,
182     .                "ave(X)", zcals(n),zout_isccp(n))
183c
184          CALL histdef(nid_isccp, "deep",
185     .                "Deep lk ISCCP-D2", "%",
186     .                iim, jj_nb,nhori,1,1,1,-99,32,
187     .                "ave(X)", zcals(n),zout_isccp(n))
188c
189          CALL histdef(nid_isccp, "alcu",
190     .                "AlCu lk ISCCP-D2", "%",
191     .                iim, jj_nb,nhori,1,1,1,-99,32,
192     .                "ave(X)", zcals(n),zout_isccp(n))
193c
194          CALL histdef(nid_isccp, "alst",
195     .                "AlSt lk ISCCP-D2", "%",
196     .                iim, jj_nb,nhori,1,1,1,-99,32,
197     .                "ave(X)", zcals(n),zout_isccp(n))
198c
199          CALL histdef(nid_isccp, "nist",
200     .                "NiSt lk ISCCP-D2", "%",
201     .                iim, jj_nb,nhori,1,1,1,-99,32,
202     .                "ave(X)", zcals(n),zout_isccp(n))
203c
204          CALL histdef(nid_isccp, "cumu",
205     .                "Cumu lk ISCCP-D2", "%",
206     .                iim, jj_nb,nhori,1,1,1,-99,32,
207     .                "ave(X)", zcals(n),zout_isccp(n))
208c
209          CALL histdef(nid_isccp, "stcu",
210     .                "StCu lk ISCCP-D2", "%",
211     .                iim, jj_nb,nhori,1,1,1,-99,32,
212     .                "ave(X)", zcals(n),zout_isccp(n))
213c
214          CALL histdef(nid_isccp, "stra",
215     .                "Stra lk ISCCP-D2", "%",
216     .                iim, jj_nb,nhori,1,1,1,-99,32,
217     .                "ave(X)", zcals(n),zout_isccp(n))
218c
219c 3_epaisseurs_optiques x3_pressions_au_sommet_des_nuages  types de nuages
220          CALL histdef(nid_isccp, "thin",
221     .                "Opt. thin ISCCP-D2 like clouds", "%",
222     .                iim, jj_nb,nhori,lmax3,1,lmax3,nvert3,32,
223     .                "ave(X)", zcals(n),zout_isccp(n))
224c
225          CALL histdef(nid_isccp, "mid",
226     .                "Opt. intermediate ISCCP-D2 like clouds", "%",
227     .                iim, jj_nb,nhori,lmax3,1,lmax3,nvert3,32,
228     .                "ave(X)", zcals(n),zout_isccp(n))
229c
230          CALL histdef(nid_isccp, "thick",
231     .                "Opt. thick ISCCP-D2 like clouds", "%",
232     .                iim, jj_nb,nhori,lmax3,1,lmax3,nvert3,32,
233     .                "ave(X)", zcals(n),zout_isccp(n))
234c
235c        IF(1.EQ.0) THEN
236c        IF(n.EQ.3) THEN
237c        IF(n.EQ.1) THEN
238c
239cIM 070905 BEG
240         IF(1.EQ.0) THEN
241          print*,'n=',n,' avant histdef(..boxptop axe'
242cIM verif boxptop
243          CALL histdef(nid_isccp,"boxptop"//verticaxe(n),
244     .                "Boxptop axe"//verticaxe(n), "mb",
245     .                iim, jj_nb,nhori,
246     .                ncol(n),1,ncol(n),nvlev(n),32,
247cIM  .                ncolmx,1,ncolmx,nvlev,32,
248cIM  .                "inst(X)",dtime,dtime)
249     .                "ave(X)",zcals(n),zout_isccp(n))
250         ENDIF !(1.EQ.0) THEN
251cIM 070905 END
252c        ENDIF !(n.EQ.3) THEN
253c       ENDIF !(1.EQ.0) THEN
254c
255c         print*,'n=',n,' avant histdef(..seed axe'
256          CALL histdef(nid_isccp, "seed"//verticaxe(n),
257     .                "seed axe"//verticaxe(n), "-",
258     .                iim, jj_nb,nhori,1,1,1,-99,32,
259cIM  .                "inst(X)", dtime,dtime)
260     .                "ave(X)", zcals(n),zout_isccp(n))
261c
262         ENDDO !n
263        ENDIF
264        CALL histend(nid_isccp)
265c
266c$OMP END MASTER
267      ENDIF ! ok_isccp
Note: See TracBrowser for help on using the repository browser.