source: LMDZ4/branches/LMDZ4_AR5/libf/cosp/ini_histhfCOSP.h @ 5373

Last change on this file since 5373 was 1928, checked in by idelkadi, 11 years ago
  • Correction du bug dans l'interface avec la physique : On ne veut pas que la distinction entre les nuages convectifs et stratiformes soit prise en compte dans les calculs Cosp. Dans l'interface avec la physique, sont passees en entree pour Cosp, les quantites totales (stratiforme + convective) du contenus en eau et d'autres variables. La fraction nuageuse convective calculee dans la physique est passee en entree pour Cosp dans la version buggee. Elle est remise a 0 dans cette version corrigee.
  • Mise a jour pour ISCCP : la fraction d'ensoleillement calculee par LMDZ est passee en argument pour Isccp
  • Rajouts d'autre variables diagnostiques : sunlit (=0 si nuit et =1 si jour) parasol_crefl (reflectance integree)
File size: 8.9 KB
Line 
1! Abderrahmane Idelkadi Septebmre 2009
2! Sorties journalieres de COSP
3! Pour l'instant sorties Lidar et ISCCP
4!
5! sorties par jour
6!
7!$OMP MASTER
8        zstohf = ecrit_hf
9        zout = freq_COSP
10!
11!       PRINT*, 'La frequence de sortie hf3d est de ', ecrit_hf
12!
13
14        idayref = day_ref
15        CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
16
17        CALL histbeg_phy("histhfCOSP",itau_phy,zjulian,dtime,nhori,nid_hf_cosp) 
18
19! Definition de l'axe vertical
20       if (use_vgrid) then
21        CALL histvert(nid_hf_cosp,"height","height","m",Nlevout,vgrid%z,nvert)
22       else
23        CALL histvert(nid_hf_cosp,"presnivs","Vertical levels","mb",Nlevout,presnivs,nvert,"down")
24       endif
25
26        CALL histvert(nid_hf_cosp,"height_mlev","height_mlev","m",Nlevlmdz,vgrid%mz,nvertm)
27
28        CALL histvert(nid_hf_cosp,"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp)
29
30        CALL histvert(nid_hf_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp,"down")
31
32        CALL histvert(nid_hf_cosp,"column","column","count",Ncolumns,column_ax(1:Ncolumns),nvertcol)
33
34! Sorties LIDAR
35       if (cfg%Llidar_sim) then
36         if (cfg%Lcllcalipso) then
37         CALL histdef(nid_hf_cosp, "cllcalipso", &
38                     "Lidar Low-level Cloud Fraction", "1", &
39                     iim, jj_nb,nhori,1,1,1,-99,32, &
40                     "ave(X)", zout,zstohf)
41         endif
42         if (cfg%Lclhcalipso) then
43         CALL histdef(nid_hf_cosp, "clhcalipso", &
44                     "Lidar High-level Cloud Fraction", "1", &
45                     iim, jj_nb,nhori,1,1,1,-99,32, &
46                     "ave(X)", zout,zstohf)
47         endif
48         if (cfg%Lclmcalipso) then
49         CALL histdef(nid_hf_cosp, "clmcalipso", &
50                     "Lidar Mid-level Cloud Fraction", "1", &
51                     iim, jj_nb,nhori,1,1,1,-99,32, &
52                     "ave(X)", zout,zstohf)
53         endif
54         if (cfg%Lcltcalipso) then
55         CALL histdef(nid_hf_cosp, "cltcalipso", &
56                     "Lidar Total Cloud Fraction", "1", &
57                     iim, jj_nb,nhori,1,1,1,-99,32, &
58                     "ave(X)", zout,zstohf)
59         endif
60         if (cfg%Lclcalipso) then
61         CALL histdef(nid_hf_cosp, "clcalipso", &
62                     "Lidar Cloud Fraction (532 nm)", "1", &
63                     iim,jj_nb,nhori, Nlevout,1,Nlevout,nvert, 32, &
64                     "ave(X)", zout,zstohf)
65         endif
66           if (cfg%Lcfad_lidarsr532) then
67              do ii=1,SR_BINS
68               CALL histdef(nid_hf_cosp, "cfad_lidarsr532_"//chcol(ii), &
69                           "Lidar Scattering Ratio CFAD (532 nm)","1", &
70                           iim,jj_nb,nhori, Nlevout,1,Nlevout,nvert, 32, &
71                           "ave(X)", zout,zstohf)   
72              enddo
73           endif
74           if (cfg%Lparasol_refl) then
75            CALL histdef(nid_hf_cosp, "parasol_refl", &
76                        "PARASOL-like mono-directional reflectance","1", &
77                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
78                        "ave(X)", zout,zstohf)   
79            CALL histdef(nid_hf_cosp, "parasol_crefl", &
80                        "PARASOL-like mono-directional reflectance","1", &
81                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
82                        "ave(X)", zout,zstohf)
83            CALL histdef(nid_hf_cosp, "Ncrefl", &
84                        "1 ou 0 si cldt > o <= 0 ","1", &
85                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
86                        "ave(X)", zout,zstohf)
87           endif
88           if (cfg%Latb532) then
89            do ii=1,Ncolumns
90             CALL histdef(nid_hf_cosp, "atb532_"//chcol(ii), &
91                         "Lidar Attenuated Total Backscatter (532 nm)","1", &
92                         iim,jj_nb,nhori, Nlevlmdz,1,Nlevlmdz,nvertm, 32, &
93                         "ave(X)", zout,zstohf)
94            enddo
95           endif
96           if (cfg%Lbeta_mol532) then
97            CALL histdef(nid_hf_cosp, "beta_mol532", &
98                        "Lidar Molecular Backscatter (532 nm)","m-1 sr-1", &
99                        iim,jj_nb,nhori, Nlevlmdz,1,Nlevlmdz,nvertm, 32, &
100                         "ave(X)", zout,zstohf)
101           endif
102        endif ! Lidar
103
104! Sorties RADAR
105!Attention A FAIRE
106!        if (cfg%Lradar_sim) then
107!         print*,'Ecriture sorties Radar'
108!          if (cfg%Lcfad_dbze94) then
109!              print*,'Ecriture de cfad_dbze94.nc '
110!              A revoir l axe vertical Nlvgrid
111!               do ii=1,DBZE_BINS
112!                   dbze_ax(ii) = CFAD_ZE_MIN + CFAD_ZE_WIDTH*(ii - 0.5)
113!               enddo
114!               call write_netcdf4d('cfad_dbze94.nc',use_vgrid,nlon,nlat,Nlevout,DBZE_BINS, &
115!                                   x,y,out_levs,dbze_ax,i,ndays,time,stradar%cfad_ze)
116!          endif
117!          if (cfg%Lclcalipso2) then
118!               call write_netcdf3d('clcalipso2.nc',use_vgrid,'clcalipso2', &
119!                              nlon,nlat,Nlevout,x,y,out_levs,i,ndays,time,stradar%lidar_only_freq_cloud)
120!          endif
121!          if (cfg%Ldbze94) then
122!             do ii=1,Ncolumns
123!                xcol(ii)=float(i)
124!             enddo
125!             call write_netcdf4d('dbze94.nc',use_vgrid,nlon,nlat,Nlevout,Ncolumns, &
126!                                 x,y,out_levs,xcol,i,ndays,time,sgradar%Ze_tot)
127!          endif
128!          if (cfg%Lcltlidarradar) then
129!             call write_netcdf2d('cltlidarradar.nc','cltlidarradar', &
130!                                 nlon,nlat,x,y,i,ndays,time,stradar%radar_lidar_tcc)
131!          endif
132!        endif  ! Radar
133
134! Sorties MISR
135!Attention A FAIRE
136!        if (cfg%Lmisr_sim) then
137!         print*,'Ecriture sorties Misr'
138!            call write_netcdf4d('clMISR.nc',use_vgrid,nlon,nlat,MISR_N_CTH,7, &
139!                                x,y,MISR_CTH,ISCCP_TAU,i,ndays,time,misr%fq_MISR)
140!        endif
141
142! Sorties ISCCP
143        if (cfg%Lisccp_sim) then
144
145          CALL histdef(nid_hf_cosp, "sunlit", &
146                     "0 si nuit 1 si jour", "1", &
147                     iim, jj_nb,nhori,1,1,1,-99,32, &
148                     "ave(X)", zout,zstomth)
149         if (cfg%Lclisccp2) then
150            do ii=1,7
151             CALL histdef(nid_hf_cosp, "clisccp2_"//chcol(ii), &
152                         "Cloud Fraction as Calculated by the ISCCP Simulator","1", &
153                         iim,jj_nb,nhori,7,1,7,nvertisccp, 32, &
154                         "ave(X)", zout,zstohf)
155            enddo
156          endif
157          if (cfg%Lboxtauisccp) then
158            CALL histdef(nid_hf_cosp, "boxtauisccp", &
159                         "Optical Depth in Each Column as Calculated by the ISCCP Simulator","1", &
160                         iim,jj_nb,nhori,Ncolumns,1,Ncolumns,nvertcol, 32, &
161                         "ave(X)", zout,zstohf)
162          endif
163          if (cfg%Lboxptopisccp) then
164            CALL histdef(nid_hf_cosp, "boxptopisccp", &
165                         "Cloud Top Pressure in Each Column as Calculated by the ISCCP Simulator","Pa", &
166                         iim,jj_nb,nhori,Ncolumns,1,Ncolumns,nvertcol, 32, &
167                         "ave(X)", zout,zstohf)
168          endif
169          if (cfg%Ltclisccp) then
170           CALL histdef(nid_hf_cosp, "tclisccp", &
171                     "Total Cloud Fraction as Calculated by the ISCCP Simulator", "1", &
172                     iim, jj_nb,nhori,1,1,1,-99,32, &
173                     "ave(X)", zout,zstohf)
174          endif
175          if (cfg%Lctpisccp) then
176            CALL histdef(nid_hf_cosp, "ctpisccp", &
177                     "Mean Cloud Top Pressure as Calculated by the ISCCP Simulator", "Pa", &
178                     iim, jj_nb,nhori,1,1,1,-99,32, &
179                     "ave(X)", zout,zstohf)
180          endif
181          if (cfg%Ltauisccp) then
182           CALL histdef(nid_hf_cosp, "tauisccp", &
183                     "Optical Depth as Calculated by the ISCCP Simulator", "1", &
184                     iim, jj_nb,nhori,1,1,1,-99,32, &
185                     "ave(X)", zout,zstohf)
186          endif
187          if (cfg%Lalbisccp) then
188           CALL histdef(nid_hf_cosp, "albisccp", &
189                     "Mean Cloud Albedo as Calculated by the ISCCP Simulator", "1", &
190                     iim, jj_nb,nhori,1,1,1,-99,32, &
191                     "ave(X)", zout,zstohf) 
192          endif
193          if (cfg%Lmeantbisccp) then
194            CALL histdef(nid_hf_cosp, "meantbisccp", &
195             " Mean all-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator","K", &
196             iim, jj_nb,nhori,1,1,1,-99,32, &
197             "ave(X)", zout,zstohf)
198          endif
199          if (cfg%Lmeantbclrisccp) then
200           CALL histdef(nid_hf_cosp, "meantbclrisccp", &
201            "Mean clear-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator","K", &
202             iim, jj_nb,nhori,1,1,1,-99,32, &
203             "ave(X)", zout,zstohf) 
204          endif
205        endif ! Isccp
206
207
208        CALL histend(nid_hf_cosp)
209!$OMP END MASTER
210!$OMP BARRIER
Note: See TracBrowser for help on using the repository browser.