source: trunk/libf/cosp/ini_histhfCOSP.h @ 1

Last change on this file since 1 was 1, checked in by emillour, 14 years ago

Import initial LMDZ5
File size: 8.3 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        CALL histvert(nid_hf_cosp,"height","height","m",Nlevout,vgrid%z,nvert)
21        print*,'Ok height Nlevout, height =',Nlevout,vgrid%z
22        CALL histvert(nid_hf_cosp,"height_mlev","height_mlev","m",Nlevlmdz,vgrid%mz,nvertm)
23        print*,'Ok height_mlev Nlevout height_mlev =',Nlevout,vgrid%mz
24!        CALL histvert(nid_hf_cosp,"presnivs","Vertical levels","mb",Nlevout,presnivs,nvert)
25
26        CALL histvert(nid_hf_cosp,"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp)
27
28        CALL histvert(nid_hf_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp)
29
30        CALL histvert(nid_hf_cosp,"column","column","count",Ncolumns,column_ax,nvertcol)
31
32! Sorties LIDAR
33       if (cfg%Llidar_sim) then
34         if (cfg%Lcllcalipso) then
35         CALL histdef(nid_hf_cosp, "cllcalipso", &
36                     "Lidar Low-level Cloud Fraction", "1", &
37                     iim, jj_nb,nhori,1,1,1,-99,32, &
38                     "ave(X)", zout,zstohf)
39         endif
40         if (cfg%Lclhcalipso) then
41         CALL histdef(nid_hf_cosp, "clhcalipso", &
42                     "Lidar High-level Cloud Fraction", "1", &
43                     iim, jj_nb,nhori,1,1,1,-99,32, &
44                     "ave(X)", zout,zstohf)
45         endif
46         if (cfg%Lclmcalipso) then
47         CALL histdef(nid_hf_cosp, "clmcalipso", &
48                     "Lidar Mid-level Cloud Fraction", "1", &
49                     iim, jj_nb,nhori,1,1,1,-99,32, &
50                     "ave(X)", zout,zstohf)
51         endif
52         if (cfg%Lcltcalipso) then
53         CALL histdef(nid_hf_cosp, "cltcalipso", &
54                     "Lidar Total Cloud Fraction", "1", &
55                     iim, jj_nb,nhori,1,1,1,-99,32, &
56                     "ave(X)", zout,zstohf)
57         endif
58         if (cfg%Lclcalipso) then
59         CALL histdef(nid_hf_cosp, "clcalipso", &
60                     "Lidar Cloud Fraction (532 nm)", "1", &
61                     iim,jj_nb,nhori, Nlevout,1,Nlevout,nvert, 32, &
62                     "ave(X)", zout,zstohf)
63         endif
64           if (cfg%Lcfad_lidarsr532) then
65              do ii=1,SR_BINS
66               CALL histdef(nid_hf_cosp, "cfad_lidarsr532_"//chcol(ii), &
67                           "Lidar Scattering Ratio CFAD (532 nm)","1", &
68                           iim,jj_nb,nhori, Nlevout,1,Nlevout,nvert, 32, &
69                           "ave(X)", zout,zstohf)   
70              enddo
71           endif
72           if (cfg%Lparasol_refl) then
73            CALL histdef(nid_hf_cosp, "parasol_refl", &
74                        "PARASOL-like mono-directional reflectance","1", &
75                        iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, &
76                        "ave(X)", zout,zstohf)   
77           endif
78           if (cfg%Latb532) then
79            do ii=1,Ncolumns
80             CALL histdef(nid_hf_cosp, "atb532_"//chcol(ii), &
81                         "Lidar Attenuated Total Backscatter (532 nm)","1", &
82                         iim,jj_nb,nhori, Nlevlmdz,1,Nlevlmdz,nvertm, 32, &
83                         "ave(X)", zout,zstohf)
84            enddo
85           endif
86           if (cfg%Lbeta_mol532) then
87            CALL histdef(nid_hf_cosp, "beta_mol532", &
88                        "Lidar Molecular Backscatter (532 nm)","m-1 sr-1", &
89                        iim,jj_nb,nhori, Nlevlmdz,1,Nlevlmdz,nvertm, 32, &
90                         "ave(X)", zout,zstohf)
91           endif
92        endif ! Lidar
93
94! Sorties RADAR
95!Attention A FAIRE
96!        if (cfg%Lradar_sim) then
97!         print*,'Ecriture sorties Radar'
98!          if (cfg%Lcfad_dbze94) then
99!              print*,'Ecriture de cfad_dbze94.nc '
100!              A revoir l axe vertical Nlvgrid
101!               do ii=1,DBZE_BINS
102!                   dbze_ax(ii) = CFAD_ZE_MIN + CFAD_ZE_WIDTH*(ii - 0.5)
103!               enddo
104!               call write_netcdf4d('cfad_dbze94.nc',use_vgrid,nlon,nlat,Nlevout,DBZE_BINS, &
105!                                   x,y,out_levs,dbze_ax,i,ndays,time,stradar%cfad_ze)
106!          endif
107!          if (cfg%Lclcalipso2) then
108!               call write_netcdf3d('clcalipso2.nc',use_vgrid,'clcalipso2', &
109!                              nlon,nlat,Nlevout,x,y,out_levs,i,ndays,time,stradar%lidar_only_freq_cloud)
110!          endif
111!          if (cfg%Ldbze94) then
112!             do ii=1,Ncolumns
113!                xcol(ii)=float(i)
114!             enddo
115!             call write_netcdf4d('dbze94.nc',use_vgrid,nlon,nlat,Nlevout,Ncolumns, &
116!                                 x,y,out_levs,xcol,i,ndays,time,sgradar%Ze_tot)
117!          endif
118!          if (cfg%Lcltlidarradar) then
119!             call write_netcdf2d('cltlidarradar.nc','cltlidarradar', &
120!                                 nlon,nlat,x,y,i,ndays,time,stradar%radar_lidar_tcc)
121!          endif
122!        endif  ! Radar
123
124! Sorties MISR
125!Attention A FAIRE
126!        if (cfg%Lmisr_sim) then
127!         print*,'Ecriture sorties Misr'
128!            call write_netcdf4d('clMISR.nc',use_vgrid,nlon,nlat,MISR_N_CTH,7, &
129!                                x,y,MISR_CTH,ISCCP_TAU,i,ndays,time,misr%fq_MISR)
130!        endif
131
132! Sorties ISCCP
133        if (cfg%Lisccp_sim) then
134         if (cfg%Lclisccp2) then
135            do ii=1,7
136             CALL histdef(nid_hf_cosp, "clisccp2_"//chcol(ii), &
137                         "Cloud Fraction as Calculated by the ISCCP Simulator","1", &
138                         iim,jj_nb,nhori,7,1,7,nvertisccp, 32, &
139                         "ave(X)", zout,zstohf)
140            enddo
141          endif
142          if (cfg%Lboxtauisccp) then
143            CALL histdef(nid_hf_cosp, "boxtauisccp", &
144                         "Optical Depth in Each Column as Calculated by the ISCCP Simulator","1", &
145                         iim,jj_nb,nhori,Ncolumns,1,Ncolumns,nvertcol, 32, &
146                         "ave(X)", zout,zstohf)
147          endif
148          if (cfg%Lboxptopisccp) then
149            CALL histdef(nid_hf_cosp, "boxptopisccp", &
150                         "Cloud Top Pressure in Each Column as Calculated by the ISCCP Simulator","Pa", &
151                         iim,jj_nb,nhori,Ncolumns,1,Ncolumns,nvertcol, 32, &
152                         "ave(X)", zout,zstohf)
153          endif
154          if (cfg%Ltclisccp) then
155           CALL histdef(nid_hf_cosp, "tclisccp", &
156                     "Total Cloud Fraction as Calculated by the ISCCP Simulator", "1", &
157                     iim, jj_nb,nhori,1,1,1,-99,32, &
158                     "ave(X)", zout,zstohf)
159          endif
160          if (cfg%Lctpisccp) then
161            CALL histdef(nid_hf_cosp, "ctpisccp", &
162                     "Mean Cloud Top Pressure as Calculated by the ISCCP Simulator", "Pa", &
163                     iim, jj_nb,nhori,1,1,1,-99,32, &
164                     "ave(X)", zout,zstohf)
165          endif
166          if (cfg%Ltauisccp) then
167           CALL histdef(nid_hf_cosp, "tauisccp", &
168                     "Optical Depth as Calculated by the ISCCP Simulator", "1", &
169                     iim, jj_nb,nhori,1,1,1,-99,32, &
170                     "ave(X)", zout,zstohf)
171          endif
172          if (cfg%Lalbisccp) then
173           CALL histdef(nid_hf_cosp, "albisccp", &
174                     "Mean Cloud Albedo as Calculated by the ISCCP Simulator", "1", &
175                     iim, jj_nb,nhori,1,1,1,-99,32, &
176                     "ave(X)", zout,zstohf) 
177          endif
178          if (cfg%Lmeantbisccp) then
179            CALL histdef(nid_hf_cosp, "meantbisccp", &
180             " Mean all-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator","K", &
181             iim, jj_nb,nhori,1,1,1,-99,32, &
182             "ave(X)", zout,zstohf)
183          endif
184          if (cfg%Lmeantbclrisccp) then
185           CALL histdef(nid_hf_cosp, "meantbclrisccp", &
186            "Mean clear-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator","K", &
187             iim, jj_nb,nhori,1,1,1,-99,32, &
188             "ave(X)", zout,zstohf) 
189          endif
190        endif ! Isccp
191
192
193        CALL histend(nid_hf_cosp)
194!$OMP END MASTER
195!$OMP BARRIER
Note: See TracBrowser for help on using the repository browser.