source: LMDZ4/branches/LMDZ4-dev-20091210/libf/cosp/ini_histdayCOSP.h @ 5064

Last change on this file since 5064 was 1262, checked in by idelkadi, 15 years ago

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